PMSF on morphological data crashes with the `+ASC` flag

[HS6986]

Dear IQ-TREE Developers,

I'm trying a new approach I devised to model morphological data in IQ-TREE v3.0.1—applying the PMSF method using randomly generated frequency vectors considering state space heterogeneity, which might account for heterogeneity of state space and state frequencies across characters well. The first test analysis seemed to finish successfully, but the analysis with the +ASC flag crashed with the following error messages:

ERROR: phylokernelnew.h:3255: double PhyloTree::computeLikelihoodBranchGenericSIMD(PhyloNeighbor *, PhyloNode *, bool) [VectorClass = Vec4d, SAFE_NUMERIC = true, FMA = true, SITE_MODEL = false]: Assertion `all_prob_const < 1.0 && all_prob_const >= 0.0' failed.
ERROR: 
ERROR: *** IQ-TREE CRASHES WITH SIGNAL ABORTED
ERROR: *** For bug report please send to developers:
ERROR: ***    Log file: ASC.log
ERROR: ***    Alignment files (if possible)

I understand that the +ASC flag in mixture models is not currently implemented (#12). However, when +ASC is called in mixture models, IQ-TREE should stop with an error Mixture model +ASC is not supported yet. Contact author if needed.. As such, I suspect that +ASC in the PMSF method should be possible, but that due to some bug it is obstructed.

I would greatly appreciate it if you could confirm that and fix it if it is due to a bug.

Here's the folder of the analyses.

Thank you very much for your time and support.

[HS6986]

Dear IQ-TREE Developers @bqminh @thomaskf,

Hi, just following up — this issue has been open for several months without any response.

I'd really appreciate it if you could take a look when you have a chance.

Thanks for your attention.