Version 1.5.6.

Fix bug for log-likelihood of multifurcating tree with multithreading (reported by Filipe Romero)
Cibiv · Dec 4, 2017 · 23e4d13 · 23e4d13
1 parent 0b47ae8
commit 23e4d13
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/phylokernelnew.h b/phylokernelnew.h
@@ -1191,7 +1191,7 @@ void PhyloTree::computePartialLikelihoodGenericSIMD(TraversalInfo &info, size_t
         /*--------------------- multifurcating node ------------------*/
 
         // now for-loop computing partial_lh over all site-patterns
-        VectorClass *partial_lh_all = (VectorClass*) &buffer_partial_lh_ptr[block*VectorClass::size()*2*thread_id];
+        VectorClass *partial_lh_all = (VectorClass*) &buffer_partial_lh_ptr[(2*block+nstates)*VectorClass::size()*thread_id];
         double *vec_tip = (double*)&partial_lh_all[block];
 
         for (ptn = ptn_lower; ptn < ptn_upper; ptn+=VectorClass::size()) {