Eos_poly_newt Class Reference
[Equations of state]

Polytropic equation of state (Newtonian case). More...

#include <eos.h>

Inheritance diagram for Eos_poly_newt:
Eos_poly Eos

List of all members.

Public Member Functions

 Eos_poly_newt (double gamma, double kappa)
 Standard constructor.
 Eos_poly_newt (const Eos_poly_newt &)
 Copy constructor.
virtual ~Eos_poly_newt ()
 Destructor.
void operator= (const Eos_poly_newt &)
 Assignment to another Eos_poly_newt.
virtual bool operator== (const Eos &) const
 Comparison operator (egality).
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.
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 specific enthalpy.
virtual double ener_ent_p (double ent, const Param *par=0x0) const
 Computes the total energy density from the specific enthalpy.
virtual double press_ent_p (double ent, const Param *par=0x0) const
 Computes the pressure from the specific 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 specific 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 specific 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 specific enthalpy.
virtual bool operator== (const Eos &) const =0
 Comparison operator (egality).
virtual bool operator!= (const Eos &) const =0
 Comparison operator (difference).
double get_gam () const
 Returns the adiabatic index $\gamma$ (cf. Eq. (3)).
double get_kap () const
 Returns the pressure coefficient $\kappa$ (cf.
double get_m_0 () const
 Return the individual particule mass $m_0$ (cf.
double get_mu_0 () const
 Return the relativistic chemical potential at zero pressure [unit: $m_B c^2$, with $m_B = 1.66\ 10^{-27} \ {\rm kg}$].
const char * get_name () const
 Returns the EOS name.
void set_name (const char *name_i)
 Sets the EOS name.
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 Eos * eos_from_file (FILE *)
 Construction of an EOS from a binary file.
static Eos * eos_from_file (ifstream &)
 Construction of an EOS from a formatted file.

Protected Member Functions

 Eos_poly_newt (FILE *)
 Constructor from a binary file (created by the function sauve(FILE*) ).
 Eos_poly_newt (ifstream &)
 Constructor from a formatted file.
virtual ostream & operator>> (ostream &) const
 Operator >>.
void set_auxiliary ()
 Computes the auxiliary quantities gam1 , unsgam1 , gam1sgamkap from the values of gam and kap.
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

double gam
 Adiabatic index $\gamma$ (cf. Eq. (3)).
double kap
 Pressure coefficient $\kappa$ (cf.
double m_0
 Individual particule mass $m_0$ (cf.
double mu_0
 Relativistic chemical potential at zero pressure [unit: $m_B c^2$, with $m_B = 1.66\ 10^{-27} \ {\rm kg}$].
double gam1
 $\gamma-1$
double unsgam1
 $1/(\gamma-1)$
double gam1sgamkap
 $(\gamma-1) / (\gamma \kappa) m_0$
double rel_mu_0
 $\mu_0/m_0$
double ent_0
 Enthalpy at zero pressure ($\ln (\mu_0/m_0)$).
char name [100]
 EOS name.

Friends

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

Detailed Description

Polytropic equation of state (Newtonian case).

()

This equation of state (EOS) corresponds to identical non relativistic particles of rest mass is $m_0$, whose internal energy density $\epsilon$ is related to their numerical density n by

\[ \epsilon(n) = {\kappa \over \gamma-1} n^\gamma \ . \qquad\qquad (1) \]

The (non-relativistic) chemical potential is then

\[ \mu(n) := {d\epsilon\over dn} = {\kappa \gamma \over \gamma-1} n^{\gamma-1} \ . \qquad\qquad (2) \]

The pressure is given by the (zero-temperature) First Law of Thermodynamics: $p = \mu n - \epsilon$, so that

\[ p(n) = \kappa n^\gamma \ . \qquad\qquad (3) \]

The (non-relativistic) specific enthalpy is :

\[ h(n) := {\epsilon + p \over m_0 n} \ . \qquad\qquad (4) \]

According to the (zero-temperature) First Law of Thermodynamics, the specific enthalpy is related to the chemical potential by

\[ h = {\mu \over m_0} \ . \qquad\qquad (5) \]

From this expression and relation (2), the expression of the particle density in term of the specific enthalpy is

\[ n(h) = \left[ {\gamma-1\over \gamma} {m_0 \over \kappa} h \right] ^{1/(\gamma-1)} \ . \qquad\qquad (6) \]

The energy density and pressure as functions of H can then be obtained by inserting this relation into Eq. (1) and (3).

Definition at line 1034 of file eos.h.

Constructor & Destructor Documentation

Eos_poly_newt::Eos_poly_newt ( double  gamma,
double  kappa 
)

Standard constructor.

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

Parameters:
gamma adiabatic index $\gamma$ (cf. Eq. (3))
kappa pressure coefficient $\kappa$ (cf. Eq. (3)) [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 90 of file eos_poly_newt.C.

References Eos::set_name().

Eos_poly_newt::Eos_poly_newt ( const Eos_poly_newt &  eosi  ) 

Copy constructor.

Definition at line 100 of file eos_poly_newt.C.

Eos_poly_newt::Eos_poly_newt ( 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 105 of file eos_poly_newt.C.

Eos_poly_newt::Eos_poly_newt ( 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 109 of file eos_poly_newt.C.

Eos_poly_newt::~Eos_poly_newt (  )  [virtual]

Destructor.

Definition at line 116 of file eos_poly_newt.C.

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_poly_newt::der_ener_ent_p ( double  ent,
const Param *  par = 0x0 
) const [virtual]

Computes the logarithmic derivative $d\ln e/d\ln h$ from the specific enthalpy.

Parameters:
ent [input, unit: $c^2$] specific enthalpy H defined by Eq. (4)
Returns:
dln(e)/dln(h)

Reimplemented from Eos_poly.

Definition at line 273 of file eos_poly_newt.C.

References exp(), Eos_poly::gam, Eos_poly::gam1, Eos_poly::gam1sgamkap, Eos_poly::kap, Eos_poly::m_0, pow(), and Eos_poly::unsgam1.

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_poly_newt::der_nbar_ent_p ( double  ent,
const Param *  par = 0x0 
) const [virtual]

Computes the logarithmic derivative $d\ln n/d\ln h$ from the specific enthalpy.

Parameters:
ent [input, unit: $c^2$] specific enthalpy H defined by Eq. (4)
Returns:
dln(n)/dln(h)

Reimplemented from Eos_poly.

Definition at line 264 of file eos_poly_newt.C.

References Eos_poly::gam1.

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_poly_newt::der_press_ent_p ( double  ent,
const Param *  par = 0x0 
) const [virtual]

Computes the logarithmic derivative $d\ln p/d\ln h$ from the specific enthalpy.

Parameters:
ent [input, unit: $c^2$] specific enthalpy H defined by Eq. (4)
Returns:
dln(p)/dln(h)

Reimplemented from Eos_poly.

Definition at line 297 of file eos_poly_newt.C.

References Eos_poly::gam, and Eos_poly::gam1.

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_poly_newt::ener_ent_p ( double  ent,
const Param *  par = 0x0 
) const [virtual]

Computes the total energy density from the specific enthalpy.

Parameters:
ent [input, unit: $c^2$] specific enthalpy H defined by Eq. (4)
Returns:
energy density e [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Reimplemented from Eos_poly.

Definition at line 229 of file eos_poly_newt.C.

References Eos_poly::gam, Eos_poly::gam1sgamkap, Eos_poly::kap, Eos_poly::m_0, pow(), and Eos_poly::unsgam1.

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().

double Eos_poly::get_gam (  )  const [inherited]

Returns the adiabatic index $\gamma$ (cf. Eq. (3)).

Definition at line 252 of file eos_poly.C.

References Eos_poly::gam.

double Eos_poly::get_kap (  )  const [inherited]

Returns the pressure coefficient $\kappa$ (cf.

Eq. (3)) [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 256 of file eos_poly.C.

References Eos_poly::kap.

double Eos_poly::get_m_0 (  )  const [inherited]

Return the individual particule mass $m_0$ (cf.

Eq. (1)) [unit: $m_B = 1.66\ 10^{-27} \ {\rm kg}$].

Definition at line 260 of file eos_poly.C.

References Eos_poly::m_0.

double Eos_poly::get_mu_0 (  )  const [inherited]

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

Definition at line 264 of file eos_poly.C.

References Eos_poly::mu_0.

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

Returns the EOS name.

Definition at line 165 of file eos.C.

References Eos::name.

int Eos_poly_newt::identify (  )  const [virtual]

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

Reimplemented from Eos_poly.

Definition at line 121 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_poly_newt::nbar_ent_p ( double  ent,
const Param *  par = 0x0 
) const [virtual]

Computes the baryon density from the specific enthalpy.

Parameters:
ent [input, unit: $c^2$] specific enthalpy H defined by Eq. (4)
Returns:
baryon density n [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$]

Reimplemented from Eos_poly.

Definition at line 215 of file eos_poly_newt.C.

References Eos_poly::gam1sgamkap, pow(), and Eos_poly::unsgam1.

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

Comparison operator (difference).

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

Comparison operator (difference).

Reimplemented from Eos_poly.

Definition at line 179 of file eos_poly_newt.C.

References operator==().

void Eos_poly_newt::operator= ( const Eos_poly_newt &  eosi  ) 

Assignment to another Eos_poly_newt.

Reimplemented from Eos_poly.

Definition at line 125 of file eos_poly_newt.C.

References Eos_poly::gam, Eos_poly::kap, Eos_poly::m_0, Eos::name, Eos_poly::set_auxiliary(), and Eos::set_name().

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

Comparison operator (egality).

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

Comparison operator (egality).

Reimplemented from Eos_poly.

Definition at line 140 of file eos_poly_newt.C.

References Eos_poly::gam, identify(), Eos::identify(), Eos_poly::kap, and Eos_poly::m_0.

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

Operator >>.

Reimplemented from Eos_poly.

Definition at line 196 of file eos_poly_newt.C.

References Eos_poly::gam, and Eos_poly::kap.

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_poly_newt::press_ent_p ( double  ent,
const Param *  par = 0x0 
) const [virtual]

Computes the pressure from the specific enthalpy.

Parameters:
ent [input, unit: $c^2$] specific enthalpy H defined by Eq. (4)
Returns:
pressure p [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Reimplemented from Eos_poly.

Definition at line 247 of file eos_poly_newt.C.

References Eos_poly::gam, Eos_poly::gam1sgamkap, Eos_poly::kap, pow(), and Eos_poly::unsgam1.

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

Save in a file.

Reimplemented from Eos_poly.

Definition at line 190 of file eos_poly_newt.C.

void Eos_poly::set_auxiliary (  )  [protected, inherited]

Computes the auxiliary quantities gam1 , unsgam1 , gam1sgamkap from the values of gam and kap.

Definition at line 238 of file eos_poly.C.

References Eos_poly::ent_0, Eos_poly::gam, Eos_poly::gam1, Eos_poly::gam1sgamkap, Eos_poly::kap, log(), Eos_poly::m_0, Eos_poly::mu_0, Eos_poly::rel_mu_0, and Eos_poly::unsgam1.

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.

Reimplemented from Eos_poly.

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

Display.

Member Data Documentation

double Eos_poly::ent_0 [protected, inherited]

Enthalpy at zero pressure ($\ln (\mu_0/m_0)$).

Definition at line 780 of file eos.h.

double Eos_poly::gam [protected, inherited]

Adiabatic index $\gamma$ (cf. Eq. (3)).

Definition at line 754 of file eos.h.

double Eos_poly::gam1 [protected, inherited]

$\gamma-1$

Definition at line 776 of file eos.h.

double Eos_poly::gam1sgamkap [protected, inherited]

$(\gamma-1) / (\gamma \kappa) m_0$

Definition at line 778 of file eos.h.

double Eos_poly::kap [protected, inherited]

Pressure coefficient $\kappa$ (cf.

Eq. (3)) [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 761 of file eos.h.

double Eos_poly::m_0 [protected, inherited]

Individual particule mass $m_0$ (cf.

Eq. (1)) [unit: $m_B = 1.66\ 10^{-27} \ {\rm kg}$].

Definition at line 766 of file eos.h.

double Eos_poly::mu_0 [protected, inherited]

Relativistic chemical potential at zero pressure [unit: $m_B c^2$, with $m_B = 1.66\ 10^{-27} \ {\rm kg}$].

(standard value: 1)

Definition at line 772 of file eos.h.

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

EOS name.

Definition at line 186 of file eos.h.

double Eos_poly::rel_mu_0 [protected, inherited]

$\mu_0/m_0$

Definition at line 779 of file eos.h.

double Eos_poly::unsgam1 [protected, inherited]

$1/(\gamma-1)$

Definition at line 777 of file eos.h.

The documentation for this class was generated from the following files:
Generated on 7 Oct 2014 for LORENE by  doxygen 1.6.1