Comment (by Roland Haas):

The current code has all kinds of things that seems to have support for piecwiese polytropes (at the very least for a single piece) eg the function compute_P_cold__eps_cold__dPcold_drho__eps_th__h__gamma_cold. Are all those incorrectly coded and only work for a Gamma-law?

I had a look at the notebooks, none tell me how to use the code to set up piecewise polytropes. They do go into details on how the piecewise polytropes may work showing first python code and then C code that (presumably) does the same thing. I could not find any information how to use it though. Admittedly I did not read all the 10+ jupyter notebooks so it could be in one the ones I did not read.

‌

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2514/broken-support-for-piecewise-polytropes-in