Hi,
there is another TOV solver available in the "reprimand" thorn. It only solves the NS equations, however, but does not set up any initial data. It provides ordinary C++ functions for EOS handling and TOV solving, but there is no cactus interface. To set up initial data, you would need to write a small thorn or replace the TOV solver part in the available TOVSolver thorn with calls to reprimand. The thorn is a based on the standalone C++ library also named reprimand (which also has a python interface that can be pip-installed). https://github.com/wokast/RePrimAnd It is documented here https://wokast.github.io/RePrimAnd/index.html
It might be well suited for your special use case since it is designed for dealing also with discontinuous EOS, as described in this article: https://arxiv.org/abs/2404.11346
Cheers, Wolfgang.
On Thu, 2024-10-31 at 20:58 +0000, CJ Osakwe wrote:
Hello,
I am trying to set up a stable hybrid star (with a quark matter core and hadronic crust, each with their own equation of state) in the Einstein Toolkit. My initial approach was to adapt the TOVSolver thorn to achieve this, but the code is removing the density discontinuity at the hadronic-quark matter interface when it interpolates from 1D to 3D. As well the discontinuity doesn't seem to be reflected in the pressure anyway. I'm wondering if there is a thorn that is better suited for setting up a stable hybrid star?
Cheers, CJ Osakwe PhD candidate, Department of Physics and Astronomy University of Calgary
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users