xflouris / libpll

Phylogenetic Likelihood Library
GNU Affero General Public License v3.0
27 stars 6 forks source link

Bug with states <> 4 and invariant sites #83

Closed xflouris closed 8 years ago

xflouris commented 8 years ago

PLL_ATTRIB_ARCH_AVX fails when number of states is different than 4, and invariant sites are used in a second evaluation of log-likelihood. PLL_ATTRIB_ARCH_CPU and PLL_ATTRIB_PATTERN_TIP with any combination of AVX or CPU work fine. The problem appears only if AVX is used without tip pattern lookups.

The error can be reproduced in the unrooted example by changing states (e.g. to 5) and extending the subst_params array (e.g. to 10 ones).

xflouris commented 8 years ago

fixed