<html>#2497: IllinoisGRMHD is incompatible with setting TmunuBase::stress_energy_at_RHS = "no"
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Gabriele Bozzola</td></tr>
<tr><td style='text-align:right'>   Status:</td><td>open</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>major</td></tr>
<tr><td style='text-align:right'>Component:</td><td>EinsteinToolkit thorn</td></tr>
</table>

<p>Comment (by Gabriele Bozzola):</p>
<p>Upon further inspection, the PR does not solve the problem. The group <code>SetTmunu</code> is not added to the scheduler when <code>TmunuBase::stress_energy_at_RHS = "no"</code></p>
<p>I don’t know how to fix this in a clean and easy way. The scheduler looks like this</p>
<div class="codehilite"><pre><span></span><span>            GROUP HydroBase_Con2Prim: Convert from conservative to primitive variables</span>
<span>          GROUP SetTmunu: Group for calculating the stress-energy tensor</span>
<span>            TmunuBase::TmunuBase_ZeroTmunu: Initialise the stress-energy tensor to zero</span>
<span>            GROUP AddToTmunu: Add to the stress-energy tensor here</span>
<span>               IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder.</span>
<span>               IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz: Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero.            </span>
</pre></div>


<p><code>IllinoisGRMHD</code> computes Tmunu when it is doing primitive-recovery (in <code>IllinoisGRMHD_convser_to_prims</code>), so, it has to be scheduled after <code>TmunuBase_ZeroTmunu</code>. However, the group <code>SetTmunu</code> does not exist at all when `<code>TmunuBase::stress_energy_at_RHS = "no"</code></p>
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2497/illinoisgrmhd-is-incompatible-with-setting'>https://bitbucket.org/einsteintoolkit/tickets/issues/2497/illinoisgrmhd-is-incompatible-with-setting</a></p>
</html>