Hi Roland and Zach,

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.

Cheers,

Vasilis

On Thu, Jun 18, 2020, 9:04 AM Roland Haas <trac-noreply@einsteintoolkit.org> wrote:
#2414: include con2prim framework by Siegel, Moesta et al. into ET
Reporter:Roland Haas
Status:new
Milestone:
Version:
Type:proposal
Priority:major
Component:EinsteinToolkit thorn

Comment (by Roland Haas):

Relevant bits from today’s ET call minutes (http://lists.einsteintoolkit.org/pipermail/users/2020-June/007479.html):

Other groups are ading these methods, like Gabriele Bozzola to IllinoisGRMHD, and Wolfgang Kastaun is also working on these methods.  Wolfgang's method has
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.

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. 
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
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
faster to us a quartic solver in some cases.

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.

‌

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2414/include-con2prim-framework-by-siegel

_______________________________________________
Trac mailing list
Trac@einsteintoolkit.org
http://lists.einsteintoolkit.org/mailman/listinfo/trac