26 #ifndef O2SCL_COLD_NSTAR_H
27 #define O2SCL_COLD_NSTAR_H
29 #include <o2scl/eos_had_base.h>
30 #include <o2scl/tov_solve.h>
31 #include <o2scl/tov_solve.h>
32 #include <o2scl/table.h>
33 #include <o2scl/fermion.h>
34 #include <o2scl/root_cern.h>
35 #include <o2scl/mroot_cern.h>
36 #include <o2scl/mroot_hybrids.h>
37 #include <o2scl/eos_tov.h>
39 #ifndef DOXYGEN_NO_O2NS
199 void fixed(
double target_mass);
354 #ifndef DOXYGEN_INTERNAL
401 #ifndef DOXYGEN_NO_O2NS
eos_tov_interp def_eos_tov
Default EOS object for the TOV solver.
eos_had_base * hep
A pointer to the equation of state.
Naive static cold neutron star.
int set_root(root<> &rf)
Set the equation solver for the EOS.
void fixed(double target_mass)
Calculate the profile for a fixed gravitational mass.
int verbose
Verbosity parameter (default 0)
tov_solve * tp
A pointer to the TOV object.
o2_shared_ptr< table_units<> >::type get_results()
Return the results data table.
double deny_urca
The smallest density where Urca becomes disallowed.
void calc_nstar()
Calculate the M vs. R curve.
fermion * pp
A pointer to the proton.
fermion def_p
The default proton.
double acausal_ed
The energy density at which the EOS becomes acausal.
root_cern def_root
The default equation solver for the EOS.
int set_tov(tov_solve &ts)
Specify the object for solving the TOV equations.
bool well_formed
If true, the energy density of the EOS is monotonically increasing and the pressure is always positiv...
fermion def_n
The default neutron.
double solver_tol
Solver tolerance (default )
bool eos_set
True if equation of state has been set.
double solve_fun(double x)
Solve to ensure zero charge in -equilibrium.
An EOS for the TOV solver using simple linear interpolation and an optional crust EOS...
fermion_zerot fzt
Zero-temperature fermion thermodynamics.
Hadronic equation of state [abstract base].
o2_shared_ptr< table_units<> >::type eost
Storage for the EOS table.
fermion * np
A pointer to the neutron.
tov_solve def_tov
The default TOV equation solver.
Class to solve the Tolman-Oppenheimer-Volkov equations.
void calc_eos(double np_0=0.0)
Calculate the given equation of state.
double nb_end
The final baryon density (default 2.0)
double calc_urca(double np_0=0.0)
Compute the density at which the direct Urca process is allowe.
o2_shared_ptr< table_units<> >::type get_eos_results()
Get the eos table (after having called calc_eos())
double barn
The baryon density.
void set_eos(eos_had_base &he)
Set the equation of state.
int set_n_and_p(fermion &n, fermion &p)
Set the neutron and proton.
bool solver_success
If true, the last call of calc_eos() succeeded.
double acausal_pr
The pressure at which the EOS becomes acausal.
double acausal
The density at which the EOS becomes acausal.
double dnb
The baryon density stepsize (default 0.01)
double nb_start
The starting baryon density (default 0.05)
double allow_urca
The smallest density where Urca becomes allowed.
bool err_nonconv
If true, throw an exception if the calculation fails (default true)
bool include_muons
If true, include muons (default false)
root * rp
A pointer to the solver.
o2_shared_ptr< table_units<> >::type get_tov_results()
Get the results from the TOV (after having called calc_nstar())
double pressure_flat
The smallest baryon density where the pressure starts to decrease.