Code for: "One- and two-argument equation of state parametrizations with continuous sound speed for neutron star simulations". Pearson-Potekhin fitting scheme

This project contains a series of python scripts that performs an analytical fit to one- or two-argument equations of state (EoS, i.e. pressure as a function of mass-energy density, or pressure as a function of mass-energy density and electronic charge fraction) from the CompOSE database. We use the analytical respresentation described by Haensel & Potekhin 2004, A&A, 428, 191-197 and Pearson et al. 2018, MNRAS, 481, 2994-3026 to perform a fit to the high-density core region, while a polytropic model is added for the crust. We use the Generalised Piecewise Polytrope scheme described in Boyle et al. 2020, Phys. Rev. D, 102, 083027 to ensure continuity in pressure, sound speed and energy-density across the crust-core boundary.