<div dir="auto">Hi Roland and Zach,<div dir="auto"><br></div><div dir="auto">Just saying that my students Pedro Espino and Gabriele worked on this. We already have a working code that is interfaced with ILGRMHD, but the interface is generic enough and should be straightforward to include it in other codes as well. It is a summer project for us. There are some final outstanding issues that we need to iron out, and as soon as we are happy with it, our goal is to work with you and share the code with the ET.</div><div dir="auto"><br></div><div dir="auto">Cheers,</div><div dir="auto"><br></div><div dir="auto">Vasilis</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 18, 2020, 9:04 AM Roland Haas <<a href="mailto:trac-noreply@einsteintoolkit.org">trac-noreply@einsteintoolkit.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">#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 href="http://lists.einsteintoolkit.org/pipermail/users/2020-June/007479.html" rel="nofollow noreferrer" target="_blank">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" target="_blank" rel="noreferrer">https://bitbucket.org/einsteintoolkit/tickets/issues/2414/include-con2prim-framework-by-siegel</a></p>
_______________________________________________<br>
Trac mailing list<br>
<a href="mailto:Trac@einsteintoolkit.org" target="_blank" rel="noreferrer">Trac@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/trac" rel="noreferrer noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/trac</a><br>
</blockquote></div>