#include <had_gibbs_nseos.h>
This computes neutron star matter assuming a Gibbs phase transition between two hadronic_eos 's. When possible it generally uses hadronic_eos::calc_e() over hadronic_eos::calc_p().
This class is not very sophisticated and the functions calc_e() and fixed_chi() will fail if the initial guess is not good enough.
Definition at line 46 of file had_gibbs_nseos.h.
Public Member Functions | |
virtual int | calc_e (double nb, fermion &n1, fermion &p1, thermo &th1, fermion &n2, fermion &p2, thermo &th2, double &chi, int type) |
Compute the properties of neutron star matter at. | |
virtual int | fixed_chi (double chi, fermion &n1, fermion &p1, thermo &th1, fermion &n2, fermion &p2, thermo &th2, double &nb, int type) |
Compute the properties of neutron star matter in the. | |
int | set_hadronic_eos1 (hadronic_eos &he) |
Set the Phase 1 hadronic equation of state. | |
int | set_hadronic_eos2 (hadronic_eos &he) |
Set the Phase 2 hadronic equation of state. | |
int | set_solver (mroot< void *, mm_funct< void * > > &mr) |
Change the solver. | |
Data Fields | |
gsl_mroot_hybrids< void *, mm_funct< void * > > | def_solver |
The default solver. | |
int | phase_guess |
The guess for the phase of the next call to calc_e() (default: phase1 ). | |
Static Public Attributes | |
Possible values for phase_guess | |
static const int | phase1 = 1 |
static const int | phase2 = 2 |
static const int | mixed_phase = 3 |
Possible values for type | |
static const int | nuclear_matter = 1 |
static const int | neutron_matter = 2 |
static const int | nstar_matter = 3 |
Protected Types | |
typedef struct had_gibbs_nseos::solvepar_s | solvepar |
typedef struct had_gibbs_nseos::chipar_s | chipar |
Protected Member Functions | |
int | phase1_nstar (size_t nv, const ovector_view &x, ovector_view &y, void *&pa) |
Solve for Phase 1 matter. | |
int | phase2_nstar (size_t nv, const ovector_view &x, ovector_view &y, void *&pa) |
Solve for Phase 2 matter. | |
int | gibbs_nstar (size_t nv, const ovector_view &x, ovector_view &y, void *&pa) |
Solve for mixed phase matter. | |
int | fixed_chi_nstar (size_t nv, const ovector_view &x, ovector_view &y, void *&pa) |
Solve for mixed phase matter at fixed chi. | |
int | gibbs_nuc (size_t nv, const ovector_view &x, ovector_view &y, void *&pa) |
Solve for mixed phase matter. | |
int | fixed_chi_nuc (size_t nv, const ovector_view &x, ovector_view &y, void *&pa) |
Solve for mixed phase matter at fixed chi. | |
int | gibbs_neut (size_t nv, const ovector_view &x, ovector_view &y, void *&pa) |
Solve for mixed phase matter. | |
int | fixed_chi_neut (size_t nv, const ovector_view &x, ovector_view &y, void *&pa) |
Solve for mixed phase matter at fixed chi. | |
Protected Attributes | |
hadronic_eos * | hep1 |
The hadronic EOS for Phase 1. | |
hadronic_eos * | hep2 |
The hadronic EOS for Phase 2. | |
fermion | elec |
The electron. | |
fermion | mu |
The muon. | |
mroot< void *, mm_funct< void * > > * | solverp |
The solver. | |
Storage for the solution vectors | |
ovector | gx |
ovector | px |
Data Structures | |
struct | chipar_s |
Useful structure for parameters to had_gibbs_nseos::fixed_chi_eqs() [protected]. More... | |
struct | solvepar_s |
Useful structure for parameters to phase1_eqs(), phase2_eqs(), and gibbs_eqs() [protected]. More... |
virtual int calc_e | ( | double | nb, | |
fermion & | n1, | |||
fermion & | p1, | |||
thermo & | th1, | |||
fermion & | n2, | |||
fermion & | p2, | |||
thermo & | th2, | |||
double & | chi, | |||
int | type | |||
) | [virtual] |
Compute the properties of neutron star matter at.
fixed density
virtual int fixed_chi | ( | double | chi, | |
fermion & | n1, | |||
fermion & | p1, | |||
thermo & | th1, | |||
fermion & | n2, | |||
fermion & | p2, | |||
thermo & | th2, | |||
double & | nb, | |||
int | type | |||
) | [virtual] |
Compute the properties of neutron star matter in the.
mixed phase for fixed chi
Change the solver.
In order to work properly, the solver needs to gracefully recover from a function which returns a non-zero value. If a gsl_mroot_hybrids object is used, then gsl_mroot_hybrids::shrink_step should be set to true
. This is done automatically for the default solver.
Definition at line 94 of file had_gibbs_nseos.h.
Documentation generated with Doxygen and provided under the GNU Free Documentation License. See License Information for details.
Project hosting provided by
,
O2scl Sourceforge Project Page