[Users] hybrid star setup in the TOVSolver thorn

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