<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 Vikram Manikantan):</p>
<p>Hi Leonardo, </p>
<p>In your suggested modification above, do you mean to replace the current lines that have ‘SetTmunu’ or just append this IF statement to the end of the file? See below for existing code in schedule.ccl.</p>
<div class="codehilite"><pre><span></span><span class="n">schedule</span> <span class="n">IllinoisGRMHD_conserv_to_prims</span> <span class="k">in</span> <span class="n">SetTmunu</span> <span class="k">after</span> <span class="p">(</span><span class="n">compute_B_and_Bstagger_from_A</span><span class="p">,</span> <span class="n">TmunuBase_ZeroTmunu</span><span class="p">)</span>
<span>{</span>
<span class="n">LANG</span><span class="p">:</span> <span class="k">C</span>
<span>}</span> <span class="ss">"Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder."</span>
<span class="n">schedule</span> <span class="n">IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz</span> <span class="k">in</span> <span class="n">SetTmunu</span> <span class="k">after</span> <span class="n">IllinoisGRMHD_conserv_to_prims</span>
<span>{</span>
<span class="o">#</span> <span class="n">We</span> <span class="n">must</span> <span class="n">sync</span> <span>{</span><span class="n">P</span><span class="p">,</span><span class="n">rho_b</span><span class="p">,</span><span class="n">vx</span><span class="p">,</span><span class="n">vy</span><span class="p">,</span><span class="n">vz</span><span>}</span> <span class="n">here</span><span class="p">.</span>
<span class="n">SYNC</span><span class="p">:</span> <span class="n">grmhd_primitives_allbutBi</span>
<span class="n">LANG</span><span class="p">:</span> <span class="k">C</span>
<span>}</span> <span class="ss">"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></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>