There appears to be some inconsistency with make_center, for a tabulated equation of state (p,h) is returned whereas, for the polytropic equation of state, it looks like p and the natural log of h is returned! I haven't got the EoS relationships at hand so this needs to be checked later.
There appears to be some inconsistency with make_center, for a tabulated equation of state (p,h) is returned whereas, for the polytropic equation of state, it looks like p and the natural log of h is returned! I haven't got the EoS relationships at hand so this needs to be checked later.