Thread-6" java.lang.Error: randomChoicePDF falls through -- negative, infinite or NaN components in input distribution, or all zeroes?

[Bill]

Dear BEAST users,

My phylogeography analysis (discrete, GMRF skyride) stopped a third of the way into the run with the following error below. I am using BEAST v1.8.1 on a Win7 PC 64-bit OS using a TESLA K20c video card (CUDA v4.2.1) and 192 GB RAM. The JAVA version is JRE6.

Any suggestions?

Bill

Exception in thread "Thread-6" java.lang.Error: randomChoicePDF falls through -- negative, infinite or NaN components in input distribution, or all zeroes?

at dr.math.MathUtils.randomChoicePDF(Unknown Source)
at dr.app.beagle.evomodel.treelikelihood.AncestralStateBeagleTreeLikelihood.drawChoice(Unknown Source)
at dr.app.beagle.evomodel.treelikelihood.AncestralStateBeagleTreeLikelihood.traverseSample(Unknown Source)
at dr.app.beagle.evomodel.treelikelihood.AncestralStateBeagleTreeLikelihood.redrawAncestralStates(Unknown Source)
at dr.app.beagle.evomodel.treelikelihood.AncestralStateBeagleTreeLikelihood.getStatesForNode(Unknown Source)
at dr.app.beagle.evomodel.treelikelihood.AncestralStateBeagleTreeLikelihood$1.getTraitString(Unknown Source)
at dr.evomodelxml.tree.TreeLoggerParser$2.getTraitString(Unknown Source)
at dr.evolution.tree.Tree$Utils.writeTreeTraits(Unknown Source)
at dr.evolution.tree.Tree$Utils.newick(Unknown Source)
at dr.evolution.tree.Tree$Utils.newick(Unknown Source)
at d3150000 -8168.8391 -2362.3891 -5806.4499 2.82668E6 4.0215E-9 9.24122E-8 11.0000 3.81 hours/million states
0 0.0
1 0.0
2 0.0
3 0.0
r.evolution.tree.Tree$Utils.newick(Unknown Source)
at dr.evomodel.tree.TreeLogger.log(Unknown Source)
at dr.inference.mcmc.MCMC$1.currentState(Unknown Source)
at dr.inference.markovchain.MarkovChain.fireCurrentModel(Unknown Source)
at dr.inference.markovchain.MarkovChain.runChain(Unknown Source)
at dr.inference.mcmc.MCMC.chain(Unknown Source)
at dr.inference.mcmc.MCMC.run(Unknown Source)
at java.lang.Thread.run(Unknown Source)
12.059499166666665 hours

[Matthew]

I got the same error.

[Andrew Rambaut]

Dear Matthew,

Are you using double precision with your TESLA card? Check in the output that BEAGLE is reporting it is using double

precision. If this is not the issue try using '-beagle_scaling always'.

Andrew

--
You received this message because you are subscribed to the Google Groups "beast-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to beast-users...@googlegroups.com.
To post to this group, send email to beast...@googlegroups.com.
Visit this group at http://groups.google.com/group/beast-users.
For more options, visit https://groups.google.com/d/optout.

[Kehinde Oyeniran]

Hi, 

I am having this error now running Beast 1.10.4 on the HPC cluster, and I can't figure out where it is coming from. I have now attached photos from the log files. It is as if the program is skipping the MCMC chain cos the job status shows it completed successfully with code 0 on the cluster. Kindly help.

[Neo One]

have you found any solution?

[Kehinde Oyeniran]

This error went away after I reinstalled beast and updated my java. I think the cluster also upgraded the beast. 

Cheers 

Kehinde Adewole Oyeniran, Ph.D.

Arturo-Fallaschi Postdoctoral Fellow | Plant Systems Biology | International Centre for Genetic Engineering and Biotechnology (ICGEB), Cape Town 7925, South Africa | Email: Kehinde....@icgeb.org, kehinde....@uct.ac.za |

Department of Biological Sciences | Bamidele Olumilua University of Education, Science and Technology Ikere (BOUESTI), Nigeria | Email: oyeniran...@bouesti.edu.ng | Blog:

www.electoralreform-ng.blogspot.com | Tel: +27(0)619507372,  +234(0)8038294995 |

“As you rise, you must lift the others up"—Justice Skweyiya (RIP)

To view this discussion visit https://groups.google.com/d/msgid/beast-users/1133f9ba-9f70-49e8-9ac3-c0614d7107e9n%40googlegroups.com.