had_gibbs_nseos Class Reference

#include <had_gibbs_nseos.h>


Detailed Description

Create a Gibbs phase transition between two hadronic EOSs to create neutron star matter.

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 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_eoshep1
 The hadronic EOS for Phase 1.
hadronic_eoshep2
 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

Member Function Documentation

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

int set_solver ( mroot< void *, mm_funct< void * > > &  mr  )  [inline]

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.


The documentation for this class was generated from the following file:

Documentation generated with Doxygen and provided under the GNU Free Documentation License. See License Information for details.

Project hosting provided by SourceForge.net Logo, O2scl Sourceforge Project Page