<html>#2514: broken support for piecewise polytropes in IllinoisGRMHD
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Roland Haas</td></tr>
<tr><td style='text-align:right'>   Status:</td><td>new</td></tr>
<tr><td style='text-align:right'>Milestone:</td><td></td></tr>
<tr><td style='text-align:right'>  Version:</td><td>development version</td></tr>
<tr><td style='text-align:right'>     Type:</td><td>bug</td></tr>
<tr><td style='text-align:right'> Priority:</td><td>minor</td></tr>
<tr><td style='text-align:right'>Component:</td><td></td></tr>
</table>

<p>Comment (by Roland Haas):</p>
<p>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 <code>compute_P_cold__eps_cold__dPcold_drho__eps_th__h__gamma_cold</code>. Are <em>all</em> those incorrectly coded and only work for a Gamma-law?</p>
<p>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.</p>
<p>‌</p>
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2514/broken-support-for-piecewise-polytropes-in'>https://bitbucket.org/einsteintoolkit/tickets/issues/2514/broken-support-for-piecewise-polytropes-in</a></p>
</html>