Lorene::Eos_tabul Class Reference
[Equations of state]

Base class for tabulated equations of state. More...

#include <eos_tabul.h>

Inheritance diagram for Lorene::Eos_tabul:
Lorene::Eos Lorene::Eos_AkmalPR Lorene::Eos_BalbN1H1 Lorene::Eos_BBB2 Lorene::Eos_BPAL12 Lorene::Eos_CompOSE Lorene::Eos_FPS Lorene::Eos_GlendNH3 Lorene::Eos_SLy4

List of all members.

Public Member Functions

virtual ~Eos_tabul ()
 Destructor.
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).
virtual int identify () const =0
 Returns a number to identify the sub-classe of Eos the object belongs to.
Cmp nbar_ent (const Cmp &ent, int nzet, int l_min=0, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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_tabul (const char *name_i, const char *table, const char *path)
 Standard constructor.
 Eos_tabul (const char *name_i, const char *file_name)
 Standard constructor from the full filename.
 Eos_tabul (const Eos_tabul &)
 Copy constructor.
 Eos_tabul (const char *name_i)
 Standard constructor with a name.
 Eos_tabul (FILE *)
 Constructor from a binary file (created by the function sauve(FILE*) ).
 Eos_tabul (ifstream &ist, const char *table)
 Constructor from a formatted file.
 Eos_tabul (ifstream &ist)
 Constructor from a formatted file.
virtual void read_table ()
 Reads the file containing the table and initializes in the arrays logh , logp and dlpsdlh .
virtual ostream & operator>> (ostream &) const =0
 Operator >>.
void calcule (const Cmp &thermo, int nzet, int l_min, double(Eos::*fait)(double, const Param *) 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, Param *par, Scalar &resu) const
 General computational method for Scalar 's.

Protected Attributes

string tablename
 Name of the file containing the tabulated data.
string authors
 Authors - reference for the table.
double hmin
 Lower boundary of the enthalpy interval.
double hmax
 Upper boundary of the enthalpy interval.
Tbllogh
 Table of $\log H$.
Tbllogp
 Table of $\log p$.
Tbldlpsdlh
 Table of $d\log P/d\log H$.
Tbllognb
 Table of $\log n_b$.
Tbldlpsdlnb
 Table of $d\log P/d\log nb$.
double * press
double * nb
double * ro
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 tabulated equations of state.

()

EoS data are to be stored in a formatted file in the following format. The first five lines are comments, preceded by hashes. Then is given the number of data lines. After, three other lines of comments, data are given in four columns. The first one is just a number, not used by Lorene. Second is the baryon number density in ${\rm fm}^{-3}$, then is the total energy density (including rest mass) and pressure, both in cgs units. Here is an example from the file eos_fps.d:

#  Date: Tue, 21 Nov 2000 17:24:31 +0100
#  From: xxx
#  FPS e.o.s.:   BPS below n.drip, then FPS. Supplied by N.Stergioulas
#  (June 1998) crust bottom at n=0.0957 , rho=1.60E14 1.E50,1.60E14
#
129  <-- Number of lines
#
#            n_B [fm^{-3}]       rho [g/cm^3]        p [dyn/cm^2]
#  
   6     0.2722041559E-13    0.4518590000E+02    0.1701030000E+15
   7     0.1279071657E-12    0.2123260000E+03    0.5817150000E+16  

Note that units used in the table file are not those of Lorene. The interpolation through the tables is a cubic Hermite interpolation, which is thermodynamically consistent, i.e. preserves the Gibbs-Duhem relation. It is defined in [Nozawa, Stergioulas, Gourgoulhon & Eriguchi, Astron. Astrophys. Suppl. Ser. 132 , 431 (1998)], and derives from a general technique presented in [Swesty, J. Comp. Phys. 127 , 118 (1996)].

Definition at line 173 of file eos_tabul.h.


Constructor & Destructor Documentation

Lorene::Eos_tabul::Eos_tabul ( const char *  name_i,
const char *  table,
const char *  path 
) [protected]

Standard constructor.

Parameters:
name_i Name of the equation of state
table Name of the file containing the EOS table
path Path to the directory containing the EOS file

Definition at line 148 of file eos_tabul.C.

References read_table(), and tablename.

Lorene::Eos_tabul::Eos_tabul ( const char *  name_i,
const char *  file_name 
) [protected]

Standard constructor from the full filename.

Parameters:
name_i Name of the equation of state
table Full name of the file containing the EOS table (including the absolute path).

Definition at line 160 of file eos_tabul.C.

References read_table(), and tablename.

Lorene::Eos_tabul::Eos_tabul ( const Eos_tabul  )  [protected]

Copy constructor.

Lorene::Eos_tabul::Eos_tabul ( const char *  name_i  )  [protected]

Standard constructor with a name.

Definition at line 206 of file eos_tabul.C.

Lorene::Eos_tabul::Eos_tabul ( 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 172 of file eos_tabul.C.

References read_table(), and tablename.

Lorene::Eos_tabul::Eos_tabul ( ifstream &  ist,
const char *  table 
) [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& ) .

Parameters:
ist input file stream containing a name as first line and the path to the directory containing the EOS file as second line
table Name of the file containing the EOS table

Definition at line 186 of file eos_tabul.C.

References read_table(), and tablename.

Lorene::Eos_tabul::Eos_tabul ( ifstream &  ist  )  [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& ) .

Parameters:
ist input file stream containing a name as first line and the full filename (including the path) containing the EOS file as second line
table Name of the file containing the EOS table

Definition at line 197 of file eos_tabul.C.

References read_table(), and tablename.

Lorene::Eos_tabul::~Eos_tabul (  )  [virtual]

Destructor.

Definition at line 216 of file eos_tabul.C.

References dlpsdlh, dlpsdlnb, logh, lognb, and logp.


Member Function Documentation

void Lorene::Eos::calcule ( const Scalar thermo,
int  nzet,
int  l_min,
double(Eos::*)(double, const Param *) const   fait,
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 282 of file eos.C.

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

void Lorene::Eos::calcule ( const Cmp thermo,
int  nzet,
int  l_min,
double(Eos::*)(double, const Param *) const   fait,
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 210 of file eos.C.

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

Scalar Lorene::Eos::der_ener_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
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 461 of file eos.C.

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

Cmp Lorene::Eos::der_ener_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
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 451 of file eos.C.

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

double Lorene::Eos_tabul::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
Returns:
dln(e)/dln(H)

Implements Lorene::Eos.

Definition at line 485 of file eos_tabul.C.

References der_nbar_ent_p(), ener_ent_p(), hmax, hmin, and press_ent_p().

Scalar Lorene::Eos::der_nbar_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
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 438 of file eos.C.

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

Cmp Lorene::Eos::der_nbar_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
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 428 of file eos.C.

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

double Lorene::Eos_tabul::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
Returns:
dln(n)/dln(H)

Implements Lorene::Eos.

Definition at line 462 of file eos_tabul.C.

References der_press_ent_p(), der_press_nbar_p(), hmax, and hmin.

Scalar Lorene::Eos::der_press_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
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 483 of file eos.C.

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

Cmp Lorene::Eos::der_press_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
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 473 of file eos.C.

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

double Lorene::Eos_tabul::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
Returns:
dln(p)/dln(H)

Implements Lorene::Eos.

Definition at line 505 of file eos_tabul.C.

References dlpsdlh, Lorene::Tbl::get_taille(), hmax, hmin, Lorene::log10(), logh, and logp.

double Lorene::Eos_tabul::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
Returns:
dln(p)/dln(n)

Definition at line 535 of file eos_tabul.C.

References dlpsdlnb, Lorene::Tbl::get_taille(), hmax, hmin, Lorene::log10(), and logh.

Scalar Lorene::Eos::ener_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
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 394 of file eos.C.

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

Cmp Lorene::Eos::ener_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
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 384 of file eos.C.

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

double Lorene::Eos_tabul::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
Returns:
energy density e [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Implements Lorene::Eos.

Reimplemented in Lorene::Eos_consistent.

Definition at line 409 of file eos_tabul.C.

References dlpsdlh, Lorene::Tbl::get_taille(), hmax, hmin, Lorene::log10(), logh, logp, and Lorene::pow().

Eos * Lorene::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 317 of file eos_from_file.C.

Eos * Lorene::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 180 of file eos_from_file.C.

References Lorene::fread_be().

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

Returns the EOS name.

Definition at line 176 of file eos.C.

References Lorene::Eos::name.

virtual int Lorene::Eos::identify (  )  const [pure virtual, inherited]
Scalar Lorene::Eos::nbar_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
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 369 of file eos.C.

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

Cmp Lorene::Eos::nbar_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
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 359 of file eos.C.

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

double Lorene::Eos_tabul::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
Returns:
baryon density n [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$]

Implements Lorene::Eos.

Reimplemented in Lorene::Eos_consistent.

Definition at line 381 of file eos_tabul.C.

References dlpsdlh, Lorene::exp(), Lorene::Tbl::get_taille(), hmax, hmin, Lorene::log10(), logh, logp, and Lorene::pow().

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

Comparison operator (difference).

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

Comparison operator (egality).

virtual ostream& Lorene::Eos::operator>> ( ostream &   )  const [protected, pure virtual, inherited]
Scalar Lorene::Eos::press_ent ( const Scalar ent,
int  nzet,
int  l_min = 0,
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 416 of file eos.C.

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

Cmp Lorene::Eos::press_ent ( const Cmp ent,
int  nzet,
int  l_min = 0,
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 406 of file eos.C.

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

double Lorene::Eos_tabul::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
Returns:
pressure p [unit: $\rho_{\rm nuc} c^2$], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Implements Lorene::Eos.

Reimplemented in Lorene::Eos_consistent.

Definition at line 437 of file eos_tabul.C.

References dlpsdlh, Lorene::Tbl::get_taille(), hmax, hmin, Lorene::log10(), logh, logp, and Lorene::pow().

void Lorene::Eos_tabul::read_table (  )  [protected, virtual]

Reads the file containing the table and initializes in the arrays logh , logp and dlpsdlh .

Reimplemented in Lorene::Eos_consistent.

Definition at line 241 of file eos_tabul.C.

References authors, dlpsdlh, dlpsdlnb, hmax, hmin, Lorene::log(), Lorene::log10(), logh, lognb, logp, Lorene::pow(), Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), and tablename.

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

Save in a file.

Reimplemented from Lorene::Eos.

Reimplemented in Lorene::Eos_CompOSE.

Definition at line 228 of file eos_tabul.C.

References tablename.

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

Sets the EOS name.

Definition at line 170 of file eos.C.

References Lorene::Eos::name.


Friends And Related Function Documentation

Eos* Eos::eos_from_file ( FILE *   )  [friend]
ostream& operator<< ( ostream &  ,
const Eos  
) [friend, inherited]

Display.


Member Data Documentation

string Lorene::Eos_tabul::authors [protected]

Authors - reference for the table.

Definition at line 182 of file eos_tabul.h.

Table of $d\log P/d\log H$.

Definition at line 197 of file eos_tabul.h.

Table of $d\log P/d\log nb$.

Definition at line 203 of file eos_tabul.h.

double Lorene::Eos_tabul::hmax [protected]

Upper boundary of the enthalpy interval.

Definition at line 188 of file eos_tabul.h.

double Lorene::Eos_tabul::hmin [protected]

Lower boundary of the enthalpy interval.

Definition at line 185 of file eos_tabul.h.

Table of $\log H$.

Definition at line 191 of file eos_tabul.h.

Table of $\log n_b$.

Definition at line 200 of file eos_tabul.h.

Table of $\log p$.

Definition at line 194 of file eos_tabul.h.

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

EOS name.

Definition at line 199 of file eos.h.

string Lorene::Eos_tabul::tablename [protected]

Name of the file containing the tabulated data.

Definition at line 180 of file eos_tabul.h.


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

Generated on 7 Dec 2019 for LORENE by  doxygen 1.6.1