Program for obtaining static Neutron Star solutions in regularized 4D Einstein-Gauss-Bonnet gravity¹, for tabulated EoS², using the KEH/CST numerical scheme³. A C++ version of the Savitzky-Golay filter⁴ is used, along with some modules from the BOOST C++ library⁵. Both are included and implemented as header-only libraries. In addition the C++ code is wrapped using SWIG⁶, into a python library which can be imported and used in any way, as shown in swig/test.py.

Compile, using the provided makefile. Type "make" in the main directory, where the makefile is also located, to create the executionable file.

EGB.exe takes 7 inputs in-turn. The parameters are specified using the following flags:

1. -f eos_name (The EoS file name).
2. -c coupling (The EGB coupling constant in km^2).
3. -e central_pressure (The central pressure in dyn/cm^2/10^35).
4. -t relative_error (The relative error for the iteration scheme).
5. -m maximum_iterations (The maximum number of iterations).
6. -l relaxation_factor
7. -p print_option:
   • 0: Prints gravitational mass M and radius R.
   • 1: Prints (0) along with the distance, metric, scalar, energy density and pressure profiles.

For example,

./EGB -f ppsly4.cold.rns1.1.txt -c 20.0 -e 10.0 -t 1e-05 -m 100 -l 0.2 -p 0

typical input parameter space:

• central pressure ~ 0.1-30.0 (10^35 dyn/cm^2)
• coupling ~ 0-70 (km^2)