Possibly they are entirely working, but they are hard to debug because doing the numerical derivatives correctly is very tricky, especially shifting the mole fractions while keeping molar density constant. A new method might be needed for testing use only.

* Add methods for calculation of pure fluid critical points

* Pass the template arguments along

* Expose the values and Jacobian

And JSON I/O