#1752: GRHydro con2prim is not a projection ------------------------------------+--------------------------------------- Reporter: rhaas | Owner: Type: defect | Status: new Priority: minor | Milestone: Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: GRHydro ------------------------------------+---------------------------------------
Comment (by rhaas):
I am thinking of the con2prim method that Wolfgang Kastaun implemented in the current WhsikyThermal which I am told works well and does not rely on an initial guess.
Keeping around the initial guess would also work. I think the memory costs are not that high since you don't need copies of all primitives, only the combinations that go into con2prim as an initial guess (which should be at most 3 scalars for the 3d Newton-Raphson inversion).