Big speedup for gradient calculations

See https://en.cppreference.com/w/cpp/types/result_of

Switch critical tracing to use one of the molar concentrations as the fixed parameter when polishing

Expose more critical methods in pybind11 interface

Also not clear if this is the best way to include it.  Probably need the std::variant approach here too.