Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

crash in fortran code with new clang #44

Open
richfitz opened this issue Dec 15, 2023 · 0 comments
Open

crash in fortran code with new clang #44

richfitz opened this issue Dec 15, 2023 · 0 comments

Comments

@richfitz
Copy link
Owner

richfitz commented Dec 15, 2023

docker pull ghcr.io/r-hub/containers/clang17
docker run -v ${PWD}:/src -it --rm ghcr.io/r-hub/containers/clang17
apt-get update && apt-get -y install gdb git libfftw3-dev libgsl-dev vim
R -e 'options(repos = "https://cloud.r-project.org"); install.packages(c("Rcpp", "ape", "deSolve", "subplex", "Rcpp"))'

# Fresh
git clone https://github.com/richfitz/diversitree
git -C diversitree checkout cran-fixes

# There are lots of concerning warnings in about the problematic part
# of the fortran code here.
R CMD INSTALL diversitree

cat <<EOF > err.R
library(diversitree)
RNGkind(sample.kind = "Rounding")
pars <- c(0.1, 0.2, 0.03, 0.03, 0.01, 0.01)
set.seed(4)
phy <- tree.bisse(pars, max.t=30, x0=0)
unresolved <- data.frame(tip.label=I(c("sp25", "sp30", "sp40", "sp56", "sp20")),
                         Nc =c(10, 9, 6, 5, 2),
                         n0=0, n1=0)
unresolved.bd <- structure(unresolved[["Nc"]], names=unresolved[["tip.label"]])
lik.bisse3 <- make.bisse(phy, phy[["tip.state"]], unresolved)
lik.bisse3(pars)
EOF
Rscript err.R

Confirmed with images 8c7713c7033e and 4a6a9279096b

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant