Hi,

can we change the parameters TOV_Gamma and TOV_K in EinsteinInitial/TOVSolver to be simple scalars instead of arrays? At the moment this makes it possible to setup initial data of multiple stars with different EOSs. However, I don't know of a code which would then be able to evolve this and how that should look like. Is anyone using this ability to choose different EOSs for multiple stars?

Frank