dppdiv-par-sse: MbRandom.cpp:1742: unsigned int MbRandom::categoricalRv(const double*, unsigned int): Assertion `u < 1.0e-6' failed.

[Mark Derbyshire]

Dear dppdiv users and developers,

I'm currently trying to run dppdiv-par on a large concatenated amino acid alignment with 107 taxa and 37621 sites.

The computer cluster I am using runs the Cray Linux Environment with the SLURM system for job scheduling.

I have a rooted newick tree with no branch or inner node labels.

I am running these commands:

#SBATCH --nodes=1
#SBATCH --time=10:00:00

export PATH=/group/y95/mderbyshire/local/general/src/dppdiv-parallel-0.0.1/:$PATH
export OMP_NUM_THREADS=6

srun dppdiv-par -in combined_alignment.txt -out test1 -tre test.nwk -cal calibrations.txt

Unfortunately, this produces the following error after one MCMC sample (sometimes before one):

==> Finished with Model
==> Starting MCMC
     1 -- -1404533.164 -> -1404443.867
dppdiv-par: MbRandom.cpp:1742: unsigned int MbRandom::categoricalRv(const double*, unsigned int): Assertion `u < 1.0e-6' failed.

My first question is, are you able to determine what is causing this error?

My second question is, is there any way to tell how many threads the software is using? If I launch multiple MPI tasks the executable spawns multiple times, so I've specified just the number of nodes, hoping that dppdiv will use the maximum available number of cores (24).

Any help would be much appreciated.

Kind regards,

Mark