[Users] GRHydro_c2p_failure_problem

Roland Haas rhaas at illinois.edu
Thu Dec 16 09:41:04 CST 2021

Hello Jimmy,

sorry for the late response. Unfortunately there is not usually a
simple fix to con2prim failures since it involves a non-linear root
finding with constraints. 

You may find that sometimes it helps to enable MHD simulation by

HydroBase::Bvec_evolution_method = "GRHdyro"

even though you are evolving a magnetic field of strength zero. Setting
this parameter selects a different con2prim routine that can be more

It can also help to increase resolution in particular since the gallery
example is quite poorly resolved.

Typically you would like to have a resolution of at least 0.1Msun (so
150m or so) for production simulations.

You can find some parameter files and discussion for production
simulations on the website of the University of Parma:



> Hi,
>   I've successfully run the gallery BNS simulation and as a second step, I try to replace the single polytropic EoS with the 7-segment SLyPP one mentioned in the 2021_ETK_Workshop. Mainly I made two substitutions in the original par file, I changed the Thorn EOS_Omni parameters into the hybrid form and I changed the type of the GRHydro_eos_table into "Hybrid". However, the simulation aborted due to the GRHydro_c2p_failure, as can be seen in the output .err file:
>  ----------------------------------------------------------------------
>   while executing schedule bin CCTK_POSTSTEP, routine GRHydro::check_GRHydro_C2P_failed
>   in thorn GRHydro, file GRHydro_Con2Prim.F90:1787:
>   ->                   xyz location:     14.06250       -5.625000        1.125000  
>  ----------------------------------------------------------------------
> I noticed that the failure point is near the outer layer of the neutron star. For clearance, I uploaded my par file and the err file. I wonder how to address this issue and I'm looking forward to your help.
> Yours,
> Jimmy

My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://pgp.mit.edu .
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 833 bytes
Desc: OpenPGP digital signature
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20211216/e8d86e93/attachment.bin 

More information about the Users mailing list