missing method
gsl is more tolerant, because of boundary condutions, spline
has zero zecond derivatives, not test function

metaprogramming cleare, wilke mixture pass with pure species objects