<html>#2414: include con2prim framework by Siegel, Moesta et al. into ET
<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></td></tr>
<tr><td style='text-align:right'> Type:</td><td>proposal</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 Roland Haas):</p>
<p>Relevant bits from today’s ET call minutes (<a data-is-external-link="true" href="http://lists.einsteintoolkit.org/pipermail/users/2020-June/007479.html" rel="nofollow">http://lists.einsteintoolkit.org/pipermail/users/2020-June/007479.html</a>):</p>
<blockquote>
<p>Other groups are ading these methods, like Gabriele Bozzola to IllinoisGRMHD, and Wolfgang Kastaun is also working on these methods. Wolfgang's method has<br />
a proof of solution to roots...if it has a solution. ZE, often with AMR the conservatives can be ill-valued. His code has a lot of dependencies, Boost, Ninja, etc.</p>
<p>ZE, the ideal interface for any GRMHD code is not a con2prim solver that does everything at all points, but one where you call it and it does it at one point. <br />
Convenient if you can do it in one OpenMP loop, so the point-by-point is better. Suppose it fails to find a root, if it fails fall-back to another method. Make the interface a point-wise method. RH agrees. But point-wise is maybe the<br />
other extreme, you may want that it takes a 1D array of data and then flags individual SIMD-parallelization a Newton-Raphson method. ZE, in some cases it is<br />
faster to us a quartic solver in some cases.</p>
<p>RH, if this is a summer project for other groups with an end time, suggest waiting until the end of summer to see what their results are.</p>
</blockquote>
<p></p>
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2414/include-con2prim-framework-by-siegel'>https://bitbucket.org/einsteintoolkit/tickets/issues/2414/include-con2prim-framework-by-siegel</a></p>
</html>