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 _______________________________________________ Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users
