template EOS on Real type for inversion

[zingale]

We should track how many EOS calls it takes for us to converge in Newton-Raphson and see whether we can first invert the EOS in single precision and then do a final double precision iteration. This could improve performance if we are taking several iterations.