Eos_multi_poly Class Reference
[Equations of state]

Base class for a multiple polytropic equation of state. More...

#include <eos_multi_poly.h>

Inheritance diagram for Eos_multi_poly:
Eos

List of all members.

Public Member Functions

 Eos_multi_poly (int npoly, double *gamma_i, double kappa0_i, double logP1_i, double *logRho_i, double *decInc_i)
 Standard constructor (sets m0 to 1).
 Eos_multi_poly (const Eos_multi_poly &)
 Copy constructor.
virtual ~Eos_multi_poly ()
 Destructor.
void operator= (const Eos_multi_poly &)
 Assignment to another Eos_multi_poly.
virtual bool operator== (const Eos &) const
 Read/write kappa.
virtual bool operator!= (const Eos &) const
 Comparison operator (difference).
virtual int identify () const
 Returns a number to identify the sub-classe of Eos the object belongs to.
const int & get_npeos () const
 Returns the number of polytropes npeos.
const double & get_gamma (int n) const
 Returns the adiabatic index $\gamma$.
const double & get_kappa0 () const
 Returns the pressure coefficient for the crust.
const double & get_logP1 () const
 Returns the exponent of the pressure at the fiducial density.
const double & get_logRho (int n) const
 Returns the exponent of fiducial densities.
const double & get_kappa (int n) const
 Returns the pressure coefficient $\kappa$ [unit: $\rho_{\rm nuc} c^2 / n_{\rm nuc}^\gamma$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$ and $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$.
const double & get_nbCrit (int n) const
 Returns the critical number density.
const double & get_entCrit (int n) const
 Returns the critical enthalpy.
virtual void sauve (FILE *) const
 Save in a file.
virtual double nbar_ent_p (double ent, const Param *par=0x0) const
 Computes the baryon density from the log-enthalpy.
virtual double ener_ent_p (double ent, const Param *par=0x0) const
 Computes the total energy density from the log-enthalpy.
virtual double press_ent_p (double ent, const Param *par=0x0) const
 Computes the pressure from the log-enthalpy.
virtual double der_nbar_ent_p (double ent, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy.
virtual double der_ener_ent_p (double ent, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy.
virtual double der_press_ent_p (double ent, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy.
virtual double der_press_nbar_p (double ent, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln p/d\ln n$ from the log-enthalpy.
const char * get_name () const
 Returns the EOS name.
void set_name (const char *name_i)
 Sets the EOS name.
virtual bool operator== (const Eos &) const =0
 Comparison operator (egality).
virtual bool operator!= (const Eos &) const =0
 Comparison operator (difference).
Cmp nbar_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the baryon density field from the log-enthalpy field and extra parameters.
Scalar nbar_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the baryon density field from the log-enthalpy field and extra parameters.
Cmp ener_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the total energy density from the log-enthalpy and extra parameters.
Scalar ener_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the total energy density from the log-enthalpy and extra parameters.
Cmp press_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the pressure from the log-enthalpy and extra parameters.
Scalar press_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the pressure from the log-enthalpy and extra parameters.
Cmp der_nbar_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy and extra parameters.
Scalar der_nbar_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy and extra parameters.
Cmp der_ener_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy and extra parameters.
Scalar der_ener_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy and extra parameters.
Cmp der_press_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy and extra parameters.
Scalar der_press_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const
 Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy and extra parameters.

Static Public Member Functions

static Eoseos_from_file (FILE *)
 Construction of an EOS from a binary file.
static Eoseos_from_file (ifstream &)
 Construction of an EOS from a formatted file.

Protected Member Functions

 Eos_multi_poly (FILE *)
 Constructor from a binary file (created by the function sauve(FILE*) ).
 Eos_multi_poly (ifstream &)
 Constructor from a formatted file.
void set_auxiliary ()
 Computes the auxiliary quantities.
virtual ostream & operator>> (ostream &) const
 Operator >>.
void calcule (const Cmp &thermo, int nzet, int l_min, double(Eos::*fait)(double, const Param *) const, const Param *par, Cmp &resu) const
 General computational method for Cmp 's.
void calcule (const Scalar &thermo, int nzet, int l_min, double(Eos::*fait)(double, const Param *) const, const Param *par, Scalar &resu) const
 General computational method for Scalar 's.

Protected Attributes

int npeos
 Number of polytropic equations of state.
double * gamma
 Array (size: npeos) of adiabatic index $\gamma$.
double kappa0
 Pressure coefficient for the crust [unit: $({\rm g/cm^3})^{1-\gamma_0}$].
double logP1
 Exponent of the pressure at the fiducial density $\rho_1$.
double * logRho
 Array (size: npeos - 1) of the exponent of fiducial densities.
double * kappa
 Array (size: npeos) of pressure coefficient $\kappa$ [unit: $\rho_{\rm nuc} c^2 / n_{\rm nuc}^\gamma$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$ and $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$.
double * nbCrit
 Array (size npeos - 1) of the number density at which the polytropic EOS changes its index and constant.
double * entCrit
 Array (size npeos - 1) of the critical enthalpy at which the polytropic EOS changes its index and constant.
double * decInc
 Array (size npeos - 1) of the percentage which detemines the terminating enthalpy for lower density and the starting enthalpy for higher density.
double m0
 Individual particule mass $m0$ [unit: $m_B = 1.66\ 10^{-27} \ {\rm kg}$].
double * mu0
 Array (size: npeos) of the relativistic chemical potential at zero pressure [unit: $m_B c^2$, with $m_B = 1.66\ 10^{-27} \ {\rm kg}$].
char name [100]
 EOS name.

Friends

EosEos::eos_from_file (FILE *)
 The construction functions from a file.
EosEos::eos_from_file (ifstream &)
ostream & operator<< (ostream &, const Eos &)
 Display.

Detailed Description

Base class for a multiple polytropic equation of state.

This equation of state mimics some realistic, tabulated EOSs. ()

Definition at line 80 of file eos_multi_poly.h.


Constructor & Destructor Documentation

Eos_multi_poly::Eos_multi_poly ( int  npoly,
double *  gamma_i,
double  kappa0_i,
double  logP1_i,
double *  logRho_i,
double *  decInc_i 
)

Standard constructor (sets m0 to 1).

The individual particle mass $m0$ is set to the mean baryon mass $m_B = 1.66\ 10^{-27} \ {\rm kg}$.

Parameters:
npoly number of polytropes
gamma_i array of adiabatic index $\gamma$
kappa0_i pressure coefficient for the crust
logP1_i exponent of the pressure at the fiducial density
logRho_i array of the exponent of fiducial densities
decInc_i array of percentage

Definition at line 87 of file eos_multi_poly.C.

References decInc, gamma, logRho, npeos, and set_auxiliary().

Eos_multi_poly::Eos_multi_poly ( const Eos_multi_poly eosmp  ) 

Copy constructor.

Definition at line 119 of file eos_multi_poly.C.

References decInc, entCrit, gamma, kappa, logRho, mu0, nbCrit, and npeos.

Eos_multi_poly::Eos_multi_poly ( FILE *  fich  )  [protected]

Constructor from a binary file (created by the function sauve(FILE*) ).

This constructor is protected because any EOS construction from a binary file must be done via the function Eos::eos_from_file(FILE*) .

Definition at line 169 of file eos_multi_poly.C.

References decInc, fread_be(), gamma, kappa0, logP1, logRho, m0, npeos, and set_auxiliary().

Eos_multi_poly::Eos_multi_poly ( ifstream &  fich  )  [protected]

Constructor from a formatted file.

This constructor is protected because any EOS construction from a formatted file must be done via the function Eos::eos_from_file(ifstream&) .

Definition at line 201 of file eos_multi_poly.C.

References decInc, gamma, kappa0, logP1, logRho, m0, npeos, and set_auxiliary().

Eos_multi_poly::~Eos_multi_poly (  )  [virtual]

Destructor.

Definition at line 235 of file eos_multi_poly.C.

References decInc, entCrit, gamma, kappa, logRho, mu0, and nbCrit.


Member Function Documentation

void Eos::calcule ( const Scalar thermo,
int  nzet,
int  l_min,
double(Eos::*)(double, const Param *) const   fait,
const Param par,
Scalar resu 
) const [protected, inherited]

General computational method for Scalar 's.

Parameters:
thermo [input] thermodynamical quantity (for instance the enthalpy field)from which the thermodynamical quantity resu is to be computed.
nzet [input] number of domains where resu is to be computed.
l_min [input] index of the innermost domain is which resu is to be computed [default value: 0]; resu is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
fait [input] pointer on the member function of class Eos which performs the pointwise calculation.
par possible extra parameters of the EOS
resu [output] result of the computation.

Definition at line 264 of file eos.C.

References Scalar::annule(), Valeur::c, Valeur::coef_i(), Tbl::get_etat(), Scalar::get_etat(), Tensor::get_mp(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Tbl::get_taille(), Valeur::set_etat_c_qcq(), Tbl::set_etat_qcq(), Mtbl::set_etat_qcq(), Scalar::set_etat_qcq(), Tbl::set_etat_zero(), Scalar::set_etat_zero(), Scalar::set_spectral_va(), Tbl::t, and Mtbl::t.

void Eos::calcule ( const Cmp thermo,
int  nzet,
int  l_min,
double(Eos::*)(double, const Param *) const   fait,
const Param par,
Cmp resu 
) const [protected, inherited]

General computational method for Cmp 's.

Parameters:
thermo [input] thermodynamical quantity (for instance the enthalpy field)from which the thermodynamical quantity resu is to be computed.
nzet [input] number of domains where resu is to be computed.
l_min [input] index of the innermost domain is which resu is to be computed [default value: 0]; resu is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
fait [input] pointer on the member function of class Eos which performs the pointwise calculation.
par possible extra parameters of the EOS
resu [output] result of the computation.

Definition at line 199 of file eos.C.

References Cmp::annule(), Valeur::c, Valeur::coef_i(), Tbl::get_etat(), Cmp::get_etat(), Cmp::get_mp(), Mg3d::get_nzone(), Tbl::get_taille(), Valeur::set_etat_c_qcq(), Tbl::set_etat_qcq(), Mtbl::set_etat_qcq(), Cmp::set_etat_qcq(), Tbl::set_etat_zero(), Cmp::set_etat_zero(), Tbl::t, Mtbl::t, and Cmp::va.

Scalar Eos::der_ener_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the derivative dln(e)/dln(H) is to be computed.
l_min index of the innermost domain is which the coefficient dln(n)/dln(H) is to be computed [default value: 0]; the derivative dln(e)/dln(H) is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
dln(e)/dln(H)

Definition at line 436 of file eos.C.

References Eos::calcule(), Eos::der_ener_ent_p(), and Tensor::get_mp().

Cmp Eos::der_ener_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the derivative dln(e)/dln(H) is to be computed.
l_min index of the innermost domain is which the coefficient dln(n)/dln(H) is to be computed [default value: 0]; the derivative dln(e)/dln(H) is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
dln(e)/dln(H)

Definition at line 426 of file eos.C.

References Eos::calcule(), Eos::der_ener_ent_p(), and Cmp::get_mp().

double Eos_multi_poly::der_ener_ent_p ( double  ent,
const Param par = 0x0 
) const [virtual]

Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H
par possible extra parameters of the EOS
Returns:
dln(e)/dln(H)

Implements Eos.

Definition at line 850 of file eos_multi_poly.C.

References decInc, der_press_ent_p(), der_press_nbar_p(), entCrit, exp(), gamma, kappa, log10(), m0, mu0, npeos, and pow().

Scalar Eos::der_nbar_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the derivative dln(n)/dln(H) is to be computed.
l_min index of the innermost domain is which the coefficient dln(n)/dln(H) is to be computed [default value: 0]; the derivative dln(n)/dln(H) is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
dln(n)/dln(H)

Definition at line 413 of file eos.C.

References Eos::calcule(), Eos::der_nbar_ent_p(), and Tensor::get_mp().

Cmp Eos::der_nbar_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the derivative dln(n)/dln(H) is to be computed.
l_min index of the innermost domain is which the coefficient dln(n)/dln(H) is to be computed [default value: 0]; the derivative dln(n)/dln(H) is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
dln(n)/dln(H)

Definition at line 403 of file eos.C.

References Eos::calcule(), Eos::der_nbar_ent_p(), and Cmp::get_mp().

double Eos_multi_poly::der_nbar_ent_p ( double  ent,
const Param par = 0x0 
) const [virtual]

Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H
par possible extra parameters of the EOS
Returns:
dln(n)/dln(H)

Implements Eos.

Definition at line 791 of file eos_multi_poly.C.

References decInc, der_press_ent_p(), der_press_nbar_p(), entCrit, exp(), gamma, m0, mu0, and npeos.

Scalar Eos::der_press_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the derivative dln(p)/dln(H) is to be computed.
par possible extra parameters of the EOS
l_min index of the innermost domain is which the coefficient dln(n)/dln(H) is to be computed [default value: 0]; the derivative dln(p)/dln(H) is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
Returns:
dln(p)/dln(H)

Definition at line 458 of file eos.C.

References Eos::calcule(), Eos::der_press_ent_p(), and Tensor::get_mp().

Cmp Eos::der_press_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the derivative dln(p)/dln(H) is to be computed.
par possible extra parameters of the EOS
l_min index of the innermost domain is which the coefficient dln(n)/dln(H) is to be computed [default value: 0]; the derivative dln(p)/dln(H) is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
Returns:
dln(p)/dln(H)

Definition at line 448 of file eos.C.

References Eos::calcule(), Eos::der_press_ent_p(), and Cmp::get_mp().

double Eos_multi_poly::der_press_ent_p ( double  ent,
const Param par = 0x0 
) const [virtual]

Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H
par possible extra parameters of the EOS
Returns:
dln(p)/dln(H)

Implements Eos.

Definition at line 983 of file eos_multi_poly.C.

References decInc, entCrit, exp(), gamma, kappa, log10(), m0, mu0, npeos, and pow().

double Eos_multi_poly::der_press_nbar_p ( double  ent,
const Param par = 0x0 
) const [virtual]

Computes the logarithmic derivative $d\ln p/d\ln n$ from the log-enthalpy.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H
par possible extra parameters of the EOS
Returns:
dln(p)/dln(n)

Definition at line 1081 of file eos_multi_poly.C.

References decInc, entCrit, gamma, log10(), and npeos.

Scalar Eos::ener_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the total energy density from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the energy density is to be computed.
l_min index of the innermost domain is which the energy density is to be computed [default value: 0]; the energy density is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
energy density [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Definition at line 369 of file eos.C.

References Eos::calcule(), Eos::ener_ent_p(), and Tensor::get_mp().

Cmp Eos::ener_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the total energy density from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the energy density is to be computed.
l_min index of the innermost domain is which the energy density is to be computed [default value: 0]; the energy density is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
energy density [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Definition at line 359 of file eos.C.

References Eos::calcule(), Eos::ener_ent_p(), and Cmp::get_mp().

double Eos_multi_poly::ener_ent_p ( double  ent,
const Param par = 0x0 
) const [virtual]

Computes the total energy density from the log-enthalpy.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H
par possible extra parameters of the EOS
Returns:
energy density e [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Implements Eos.

Definition at line 584 of file eos_multi_poly.C.

References decInc, entCrit, exp(), gamma, kappa, log10(), m0, mu0, npeos, and pow().

Eos * Eos::eos_from_file ( ifstream &  fich  )  [static, inherited]

Construction of an EOS from a formatted file.

The fist line of the file must start by the EOS number, according to the following conventions:

  • 1 = relativistic polytropic EOS (class Eos_poly ).
  • 2 = Newtonian polytropic EOS (class Eos_poly_newt ).
  • 3 = Relativistic incompressible EOS (class Eos_incomp ).
  • 4 = Newtonian incompressible EOS (class Eos_incomp_newt ).
  • 5 = Strange matter (MIT Bag model)
  • 6 = Strange matter (MIT Bag model) with crust
  • 10 = SLy4 (Douchin & Haensel 2001)
  • 11 = FPS (Friedman-Pandharipande + Skyrme)
  • 12 = BPAL12 (Bombaci et al. 1995)
  • 13 = AkmalPR (Akmal, Pandharipande & Ravenhall 1998)
  • 14 = BBB2 (Baldo, Bombaci & Burgio 1997)
  • 15 = BalbN1H1 (Balberg 2000)
  • 16 = GlendNH3 (Glendenning 1985, case 3)
  • 17 = Compstar (Tabulated EOS for 2010 CompStar school)
  • 18 = magnetized (tabulated) equation of state
  • 19 = relativistic ideal Fermi gas at zero temperature (class Eos_Fermi)
  • 100 = Multi-domain EOS (class MEos )
  • 110 = Multi-polytropic EOS (class Eos_multi_poly )
  • 120 = Fitted SLy4 (Shibata 2004)
  • 121 = Fitted FPS (Shibata 2004)
  • 122 = Fitted AkmalPR (Taniguchi 2005)

The second line in the file should contain a name given by the user to the EOS. The following lines should contain the EOS parameters (one parameter per line), in the same order than in the class declaration.

Definition at line 297 of file eos_from_file.C.

Eos * Eos::eos_from_file ( FILE *  fich  )  [static, inherited]

Construction of an EOS from a binary file.

The file must have been created by the function sauve(FILE*) .

Definition at line 165 of file eos_from_file.C.

References fread_be().

const double& Eos_multi_poly::get_entCrit ( int  n  )  const [inline]

Returns the critical enthalpy.

Definition at line 248 of file eos_multi_poly.h.

References entCrit, and npeos.

const double& Eos_multi_poly::get_gamma ( int  n  )  const [inline]

Returns the adiabatic index $\gamma$.

Definition at line 213 of file eos_multi_poly.h.

References gamma, and npeos.

const double& Eos_multi_poly::get_kappa ( int  n  )  const [inline]

Returns the pressure coefficient $\kappa$ [unit: $\rho_{\rm nuc} c^2 / n_{\rm nuc}^\gamma$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$ and $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$.

Definition at line 236 of file eos_multi_poly.h.

References kappa, and npeos.

const double& Eos_multi_poly::get_kappa0 (  )  const [inline]

Returns the pressure coefficient for the crust.

Definition at line 219 of file eos_multi_poly.h.

References kappa0.

const double& Eos_multi_poly::get_logP1 (  )  const [inline]

Returns the exponent of the pressure at the fiducial density.

Definition at line 222 of file eos_multi_poly.h.

References logP1.

const double& Eos_multi_poly::get_logRho ( int  n  )  const [inline]

Returns the exponent of fiducial densities.

Definition at line 225 of file eos_multi_poly.h.

References logRho, and npeos.

const char * Eos::get_name (  )  const [inherited]

Returns the EOS name.

Definition at line 165 of file eos.C.

References Eos::name.

const double& Eos_multi_poly::get_nbCrit ( int  n  )  const [inline]

Returns the critical number density.

Definition at line 242 of file eos_multi_poly.h.

References nbCrit, and npeos.

const int& Eos_multi_poly::get_npeos (  )  const [inline]

Returns the number of polytropes npeos.

Definition at line 210 of file eos_multi_poly.h.

References npeos.

int Eos_multi_poly::identify (  )  const [virtual]

Returns a number to identify the sub-classe of Eos the object belongs to.

Implements Eos.

Definition at line 153 of file eos_from_file.C.

Scalar Eos::nbar_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the baryon density field from the log-enthalpy field and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the baryon density is to be computed.
l_min index of the innermost domain is which the baryon density is to be computed [default value: 0]; the baryon density is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
baryon density [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$]

Definition at line 344 of file eos.C.

References Eos::calcule(), Tensor::get_mp(), and Eos::nbar_ent_p().

Cmp Eos::nbar_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the baryon density field from the log-enthalpy field and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the baryon density is to be computed.
l_min index of the innermost domain is which the baryon density is to be computed [default value: 0]; the baryon density is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
baryon density [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$]

Definition at line 334 of file eos.C.

References Eos::calcule(), Cmp::get_mp(), and Eos::nbar_ent_p().

double Eos_multi_poly::nbar_ent_p ( double  ent,
const Param par = 0x0 
) const [virtual]

Computes the baryon density from the log-enthalpy.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H
par possible extra parameters of the EOS
Returns:
baryon density n [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$]

Implements Eos.

Definition at line 486 of file eos_multi_poly.C.

References decInc, entCrit, exp(), gamma, kappa, log10(), m0, mu0, npeos, and pow().

virtual bool Eos::operator!= ( const Eos  )  const [pure virtual, inherited]

Comparison operator (difference).

bool Eos_multi_poly::operator!= ( const Eos eos_i  )  const [virtual]

Comparison operator (difference).

Definition at line 378 of file eos_multi_poly.C.

References operator==().

void Eos_multi_poly::operator= ( const Eos_multi_poly  ) 

Assignment to another Eos_multi_poly.

Definition at line 251 of file eos_multi_poly.C.

virtual bool Eos::operator== ( const Eos  )  const [pure virtual, inherited]

Comparison operator (egality).

bool Eos_multi_poly::operator== ( const Eos eos_i  )  const [virtual]

Read/write kappa.

Comparison operator (egality)

Definition at line 334 of file eos_multi_poly.C.

References gamma, get_gamma(), get_kappa(), get_npeos(), identify(), Eos::identify(), kappa, and npeos.

ostream & Eos_multi_poly::operator>> ( ostream &  ost  )  const [protected, virtual]

Operator >>.

Implements Eos.

Definition at line 412 of file eos_multi_poly.C.

References entCrit, gamma, kappa, logP1, logRho, mu0, nbCrit, npeos, and pow().

Scalar Eos::press_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the pressure from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the pressure is to be computed.
l_min index of the innermost domain is which the pressure is to be computed [default value: 0]; the pressure is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
pressure [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Definition at line 391 of file eos.C.

References Eos::calcule(), Tensor::get_mp(), and Eos::press_ent_p().

Cmp Eos::press_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
const Param par = 0x0 
) const [inherited]

Computes the pressure from the log-enthalpy and extra parameters.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right) $, where e is the (total) energy density, p the pressure, n the baryon density, and $m_B$ the baryon mass
nzet number of domains where the pressure is to be computed.
l_min index of the innermost domain is which the pressure is to be computed [default value: 0]; the pressure is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
par possible extra parameters of the EOS
Returns:
pressure [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Definition at line 381 of file eos.C.

References Eos::calcule(), Cmp::get_mp(), and Eos::press_ent_p().

double Eos_multi_poly::press_ent_p ( double  ent,
const Param par = 0x0 
) const [virtual]

Computes the pressure from the log-enthalpy.

Parameters:
ent [input, unit: $c^2$] log-enthalpy H
par possible extra parameters of the EOS
Returns:
pressure p [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Implements Eos.

Definition at line 693 of file eos_multi_poly.C.

References decInc, entCrit, exp(), gamma, kappa, log10(), m0, mu0, npeos, and pow().

void Eos_multi_poly::sauve ( FILE *  fich  )  const [virtual]

Save in a file.

Reimplemented from Eos.

Definition at line 388 of file eos_multi_poly.C.

References decInc, fwrite_be(), gamma, kappa0, logP1, logRho, and npeos.

void Eos_multi_poly::set_auxiliary (  )  [protected]

Computes the auxiliary quantities.

Definition at line 263 of file eos_multi_poly.C.

References entCrit, gamma, kappa, kappa0, log(), logP1, logRho, m0, mu0, nbCrit, npeos, and pow().

void Eos::set_name ( const char *  name_i  )  [inherited]

Sets the EOS name.

Definition at line 159 of file eos.C.

References Eos::name.


Friends And Related Function Documentation

Eos* Eos::eos_from_file ( FILE *   )  [friend]

The construction functions from a file.

ostream& operator<< ( ostream &  ,
const Eos  
) [friend, inherited]

Display.


Member Data Documentation

double* Eos_multi_poly::decInc [protected]

Array (size npeos - 1) of the percentage which detemines the terminating enthalpy for lower density and the starting enthalpy for higher density.

Definition at line 125 of file eos_multi_poly.h.

double* Eos_multi_poly::entCrit [protected]

Array (size npeos - 1) of the critical enthalpy at which the polytropic EOS changes its index and constant.

Definition at line 119 of file eos_multi_poly.h.

double* Eos_multi_poly::gamma [protected]

Array (size: npeos) of adiabatic index $\gamma$.

Definition at line 90 of file eos_multi_poly.h.

double* Eos_multi_poly::kappa [protected]

Array (size: npeos) of pressure coefficient $\kappa$ [unit: $\rho_{\rm nuc} c^2 / n_{\rm nuc}^\gamma$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$ and $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$.

Definition at line 109 of file eos_multi_poly.h.

double Eos_multi_poly::kappa0 [protected]

Pressure coefficient for the crust [unit: $({\rm g/cm^3})^{1-\gamma_0}$].

Definition at line 95 of file eos_multi_poly.h.

double Eos_multi_poly::logP1 [protected]

Exponent of the pressure at the fiducial density $\rho_1$.

Definition at line 98 of file eos_multi_poly.h.

double* Eos_multi_poly::logRho [protected]

Array (size: npeos - 1) of the exponent of fiducial densities.

Definition at line 101 of file eos_multi_poly.h.

double Eos_multi_poly::m0 [protected]

Individual particule mass $m0$ [unit: $m_B = 1.66\ 10^{-27} \ {\rm kg}$].

Definition at line 130 of file eos_multi_poly.h.

double* Eos_multi_poly::mu0 [protected]

Array (size: npeos) of the relativistic chemical potential at zero pressure [unit: $m_B c^2$, with $m_B = 1.66\ 10^{-27} \ {\rm kg}$].

(The value for the EOS which covers the lowest density: 1)

Definition at line 138 of file eos_multi_poly.h.

char Eos::name[100] [protected, inherited]

EOS name.

Definition at line 186 of file eos.h.

double* Eos_multi_poly::nbCrit [protected]

Array (size npeos - 1) of the number density at which the polytropic EOS changes its index and constant.

Definition at line 114 of file eos_multi_poly.h.

int Eos_multi_poly::npeos [protected]

Number of polytropic equations of state.

Definition at line 87 of file eos_multi_poly.h.


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

Generated on 7 Oct 2014 for LORENE by  doxygen 1.6.1