Check PhreeqcMatrix.subset() or removeNaN
Subsetting a PhreeqcMatrix gives strange results and it allows to go "out of bounds" on the .values without throwing a segfault and without returning NaNs.
MWE.tgz contains a minimal example to reproduce the issue.