LORENE
Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Friends | List of all members
Lorene::Eos_strange Class Reference
Equations of state

Strange matter EOS (MIT Bag model). More...

#include <eos.h>

Inheritance diagram for Lorene::Eos_strange:
Lorene::Eos

Public Member Functions

 Eos_strange (double n0_b60_i, double b60_i, double ent0_i, double eps_fit_i, double rho0_b60_i)
 Standard constructor. More...
 
 Eos_strange (const Eos_strange &)
 Copy constructor. More...
 
virtual ~Eos_strange ()
 Destructor. More...
 
void operator= (const Eos_strange &)
 Assignment to another Eos_strange. More...
 
virtual bool operator== (const Eos &) const
 Comparison operator (egality) More...
 
virtual bool operator!= (const Eos &) const
 Comparison operator (difference) More...
 
virtual int identify () const
 Returns a number to identify the sub-classe of Eos the object belongs to. More...
 
double get_n0_b60 () const
 Returns the baryon density at zero pressure divided by $B_{60}^{3/4}$ [unit: ${\rm fm}^{-3}$]. More...
 
double get_b60 () const
 Returns the bag constant [unit: $60\ {\rm MeV\, fm}^{-3}$]. More...
 
double get_ent0 () const
 Returns the log-enthalpy threshold for setting the energy density to a non zero value (should be negative). More...
 
double get_eps_fit () const
 Returns the fitting parameter $\epsilon_{\rm fit}$ related to the square of sound velocity by $c_s^2 = 1/3(1+\epsilon_{\rm fit})$. More...
 
double get_rho0_b60 () const
 Returns the energy density at zero pressure divided by $B_{60}$. More...
 
virtual void sauve (FILE *) const
 Save in a file. More...
 
virtual double nbar_ent_p (double ent, const Param *par=0x0) const
 Computes the baryon density from the log-enthalpy. More...
 
virtual double ener_ent_p (double ent, const Param *par=0x0) const
 Computes the total energy density from the log-enthalpy. More...
 
virtual double press_ent_p (double ent, const Param *par=0x0) const
 Computes the pressure from the log-enthalpy. More...
 
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. More...
 
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. More...
 
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. More...
 
virtual double csound_square_ent_p (double, const Param *) const
 Computes the sound speed squared $ c_s^2 = c^2 \frac{dp}{de}$ from the enthapy with extra parameters (virtual function implemented in the derived classes). More...
 
const char * get_name () const
 Returns the EOS name. More...
 
void set_name (const char *name_i)
 Sets the EOS name. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
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. More...
 
Scalar csound_square_ent (const Scalar &ent, int nzet, int l_min=0, Param *par=0x0) const
 Computes the sound speed squared $ c_s^2 = c^2 \frac{dp}{de}$ from the enthalpy with extra parameters. More...
 

Static Public Member Functions

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

Protected Member Functions

 Eos_strange (FILE *)
 Constructor from a binary file (created by the function sauve(FILE*) ). More...
 
 Eos_strange (ifstream &)
 Constructor from a formatted file. More...
 
void set_auxiliary ()
 Computes the auxiliary quantities n0 , rh0 , b34 and fach from the values of the other parameters. More...
 
virtual ostream & operator>> (ostream &) const
 Operator >> More...
 
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. More...
 
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. More...
 

Protected Attributes

double n0_b60
 Baryon density at zero pressure divided by $B_{60}^{3/4}$. More...
 
double b60
 Bag constant [unit: $60\ {\rm MeV\, fm}^{-3}$]. More...
 
double ent0
 Log-enthalpy threshold for setting the energy density to a non zero value (should be negative). More...
 
double eps_fit
 Fitting parameter $\epsilon_{\rm fit}$ related to the square of sound velocity by $c_s^2 = 1/3(1+\epsilon_{\rm fit})$. More...
 
double rho0_b60
 Energy density at zero pressure divided by $B_{60}$. More...
 
double n0
 Baryon density at zero pressure. More...
 
double rho0
 Energy density at zero pressure. More...
 
double b34
 $B_{60}^{3/4}$ More...
 
double fach
 Factor $(4+\epsilon_{\rm fit})/(1+\epsilon_{\rm fit})$. More...
 
char name [100]
 EOS name. More...
 

Friends

Eos * Eos::eos_from_file (FILE *)
 The construction functions from a file. More...
 
Eos * Eos::eos_from_file (ifstream &)
 

Detailed Description

Strange matter EOS (MIT Bag model).

This equation of state (EOS) corresponds to u,d,s degenerate symetric matter in the MIT bag model, according to approximate formula given in Zdunik, Astron. Astrophys. 359 , 311 (2000). ()

Definition at line 1999 of file eos.h.

Constructor & Destructor Documentation

◆ Eos_strange() [1/4]

Lorene::Eos_strange::Eos_strange ( double  n0_b60_i,
double  b60_i,
double  ent0_i,
double  eps_fit_i,
double  rho0_b60_i 
)

Standard constructor.

Parameters
n0_b60_iBaryon density at zero pressure divided by $B_{60}^{3/4}$ [unit: ${\rm fm}^{-3}$]
b60_iBag constant [unit: $60\ {\rm MeV\, fm}^{-3}$]
ent0_iLog-enthalpy threshold for setting the energy density to a non zero value (should be negative)
eps_fit_iFitting parameter $\epsilon_{\rm fit}$ related to the square of sound velocity by $c_s^2 = 1/3(1+\epsilon_{\rm fit})$ [cf. Zdunik, Astron. Astrophys. 359 , 311 (2000)]
rho0_b60_iEnergy density at zero pressure divided by $B_{60}$ [unit: ${\rm MeV\, fm^{-3}}$]

Definition at line 101 of file eos_strange.C.

References set_auxiliary().

◆ Eos_strange() [2/4]

Lorene::Eos_strange::Eos_strange ( const Eos_strange &  eos_i)

Copy constructor.

Definition at line 117 of file eos_strange.C.

References set_auxiliary().

◆ Eos_strange() [3/4]

Lorene::Eos_strange::Eos_strange ( 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 132 of file eos_strange.C.

References b60, ent0, eps_fit, Lorene::fread_be(), n0_b60, rho0_b60, and set_auxiliary().

◆ Eos_strange() [4/4]

Lorene::Eos_strange::Eos_strange ( 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 147 of file eos_strange.C.

References b60, ent0, eps_fit, n0_b60, rho0_b60, and set_auxiliary().

◆ ~Eos_strange()

Lorene::Eos_strange::~Eos_strange ( )
virtual

Destructor.

Definition at line 165 of file eos_strange.C.

Member Function Documentation

◆ calcule() [1/2]

void Lorene::Eos::calcule ( const Cmp &  thermo,
int  nzet,
int  l_min,
double(Eos::*)(double, const Param *) const  fait,
Param *  par,
Cmp &  resu 
) const
protectedinherited

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.
parpossible extra parameters of the EOS
resu[output] result of the computation.

Definition at line 213 of file eos.C.

References Lorene::Cmp::get_etat().

◆ calcule() [2/2]

void Lorene::Eos::calcule ( const Scalar &  thermo,
int  nzet,
int  l_min,
double(Eos::*)(double, const Param *) const  fait,
Param *  par,
Scalar &  resu 
) const
protectedinherited

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.
parpossible extra parameters of the EOS
resu[output] result of the computation.

Definition at line 285 of file eos.C.

References Lorene::Scalar::get_etat().

◆ csound_square_ent()

Scalar Lorene::Eos::csound_square_ent ( const Scalar &  ent,
int  nzet,
int  l_min = 0,
Param *  par = 0x0 
) const
inherited

Computes the sound speed squared $ c_s^2 = c^2 \frac{dp}{de}$ from the enthalpy with extra parameters.

Parameters
ent[input, unit: c^2] enthalpy
nzetnumber of domains where the derivative dln(e)/dln(H) is to be computed.
l_minindex 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.
parpossible extra parameters of the EOS
Returns
$c_s^2 $ [unit: c^2]

Definition at line 499 of file eos.C.

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

◆ csound_square_ent_p()

virtual double Lorene::Eos_strange::csound_square_ent_p ( double  ,
const Param *   
) const
inlinevirtual

Computes the sound speed squared $ c_s^2 = c^2 \frac{dp}{de}$ from the enthapy with extra parameters (virtual function implemented in the derived classes).

Parameters
ent[input, unit: c^2] enthalpy
parpossible extra parameters of the EOS
Returns
$c_s^2 $ [unit: c^2]

Implements Lorene::Eos.

Definition at line 2229 of file eos.h.

References Lorene::c_est_pas_fait().

◆ der_ener_ent() [1/2]

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
nzetnumber of domains where the derivative dln(e)/dln(H) is to be computed.
l_minindex 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.
parpossible extra parameters of the EOS
Returns
dln(e)/dln(H)

Definition at line 454 of file eos.C.

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

◆ der_ener_ent() [2/2]

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
nzetnumber of domains where the derivative dln(e)/dln(H) is to be computed.
l_minindex 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.
parpossible extra parameters of the EOS
Returns
dln(e)/dln(H)

Definition at line 464 of file eos.C.

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

◆ der_ener_ent_p()

double Lorene::Eos_strange::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 382 of file eos_strange.C.

References ent0, eps_fit, Lorene::exp(), and fach.

◆ der_nbar_ent() [1/2]

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
nzetnumber of domains where the derivative dln(n)/dln(H) is to be computed.
l_minindex 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.
parpossible extra parameters of the EOS
Returns
dln(n)/dln(H)

Definition at line 431 of file eos.C.

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

◆ der_nbar_ent() [2/2]

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
nzetnumber of domains where the derivative dln(n)/dln(H) is to be computed.
l_minindex 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.
parpossible extra parameters of the EOS
Returns
dln(n)/dln(H)

Definition at line 441 of file eos.C.

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

◆ der_nbar_ent_p()

double Lorene::Eos_strange::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 367 of file eos_strange.C.

References ent0, and eps_fit.

◆ der_press_ent() [1/2]

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
nzetnumber of domains where the derivative dln(p)/dln(H) is to be computed.
parpossible extra parameters of the EOS
l_minindex 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 476 of file eos.C.

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

◆ der_press_ent() [2/2]

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
nzetnumber of domains where the derivative dln(p)/dln(H) is to be computed.
parpossible extra parameters of the EOS
l_minindex 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 486 of file eos.C.

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

◆ der_press_ent_p()

double Lorene::Eos_strange::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 400 of file eos_strange.C.

References ent0, Lorene::exp(), and fach.

◆ ener_ent() [1/2]

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
nzetnumber of domains where the energy density is to be computed.
l_minindex 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.
parpossible 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 387 of file eos.C.

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

◆ ener_ent() [2/2]

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
nzetnumber of domains where the energy density is to be computed.
l_minindex 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.
parpossible 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 397 of file eos.C.

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

◆ ener_ent_p()

double Lorene::Eos_strange::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.

Definition at line 332 of file eos_strange.C.

References ent0, eps_fit, Lorene::exp(), fach, and rho0.

◆ eos_from_file() [1/2]

Eos * Lorene::Eos::eos_from_file ( FILE *  fich)
staticinherited

Construction of an EOS from a binary file.

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

Definition at line 193 of file eos_from_file.C.

References Lorene::fread_be().

◆ eos_from_file() [2/2]

Eos * Lorene::Eos::eos_from_file ( ifstream &  fich)
staticinherited

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 = Tabulated EOS in CompOSE format
  • 18 = magnetized (tabulated) equation of state
  • 19 = relativistic ideal Fermi gas at zero temperature (class Eos_Fermi)
  • 20 = Tabulated EOS in CompOSE format corrected for thermo. consistency
  • 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 345 of file eos_from_file.C.

◆ get_b60()

double Lorene::Eos_strange::get_b60 ( ) const
inline

Returns the bag constant [unit: $60\ {\rm MeV\, fm}^{-3}$].

Definition at line 2125 of file eos.h.

References b60.

◆ get_ent0()

double Lorene::Eos_strange::get_ent0 ( ) const
inline

Returns the log-enthalpy threshold for setting the energy density to a non zero value (should be negative).

Definition at line 2130 of file eos.h.

References ent0.

◆ get_eps_fit()

double Lorene::Eos_strange::get_eps_fit ( ) const
inline

Returns the fitting parameter $\epsilon_{\rm fit}$ related to the square of sound velocity by $c_s^2 = 1/3(1+\epsilon_{\rm fit})$.

[cf. Zdunik, Astron. Astrophys. 359 , 311 (2000)]

Definition at line 2136 of file eos.h.

References eps_fit.

◆ get_n0_b60()

double Lorene::Eos_strange::get_n0_b60 ( ) const
inline

Returns the baryon density at zero pressure divided by $B_{60}^{3/4}$ [unit: ${\rm fm}^{-3}$].

Definition at line 2122 of file eos.h.

References n0_b60.

◆ get_name()

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

Returns the EOS name.

Definition at line 179 of file eos.C.

References Lorene::Eos::name.

◆ get_rho0_b60()

double Lorene::Eos_strange::get_rho0_b60 ( ) const
inline

Returns the energy density at zero pressure divided by $B_{60}$.

[unit: $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$]

Definition at line 2141 of file eos.h.

References rho0_b60.

◆ identify()

int Lorene::Eos_strange::identify ( ) const
virtual

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

Implements Lorene::Eos.

Definition at line 147 of file eos_from_file.C.

◆ nbar_ent() [1/2]

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
nzetnumber of domains where the baryon density is to be computed.
l_minindex 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.
parpossible extra parameters of the EOS
Returns
baryon density [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$]

Definition at line 362 of file eos.C.

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

◆ nbar_ent() [2/2]

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
nzetnumber of domains where the baryon density is to be computed.
l_minindex 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.
parpossible extra parameters of the EOS
Returns
baryon density [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$]

Definition at line 372 of file eos.C.

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

◆ nbar_ent_p()

double Lorene::Eos_strange::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.

Definition at line 317 of file eos_strange.C.

References ent0, eps_fit, Lorene::exp(), and n0.

◆ operator!=()

bool Lorene::Eos_strange::operator!= ( const Eos &  eos_i) const
virtual

Comparison operator (difference)

Implements Lorene::Eos.

Definition at line 268 of file eos_strange.C.

References operator==().

◆ operator=()

void Lorene::Eos_strange::operator= ( const Eos_strange &  eosi)

Assignment to another Eos_strange.

Definition at line 175 of file eos_strange.C.

References b60, ent0, eps_fit, n0_b60, Lorene::Eos::name, rho0_b60, set_auxiliary(), and Lorene::Eos::set_name().

◆ operator==()

bool Lorene::Eos_strange::operator== ( const Eos &  eos_i) const
virtual

Comparison operator (egality)

Implements Lorene::Eos.

Definition at line 214 of file eos_strange.C.

References b60, ent0, eps_fit, Lorene::Eos::identify(), identify(), n0_b60, and rho0_b60.

◆ operator>>()

ostream & Lorene::Eos_strange::operator>> ( ostream &  ost) const
protectedvirtual

Operator >>

Implements Lorene::Eos.

Definition at line 290 of file eos_strange.C.

References b60, ent0, eps_fit, n0_b60, and rho0_b60.

◆ press_ent() [1/2]

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
nzetnumber of domains where the pressure is to be computed.
l_minindex 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.
parpossible 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 409 of file eos.C.

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

◆ press_ent() [2/2]

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
nzetnumber of domains where the pressure is to be computed.
l_minindex 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.
parpossible 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 419 of file eos.C.

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

◆ press_ent_p()

double Lorene::Eos_strange::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.

Definition at line 350 of file eos_strange.C.

References ent0, Lorene::exp(), fach, and rho0.

◆ sauve()

void Lorene::Eos_strange::sauve ( FILE *  fich) const
virtual

Save in a file.

Reimplemented from Lorene::Eos.

Definition at line 278 of file eos_strange.C.

References b60, ent0, eps_fit, Lorene::fwrite_be(), n0_b60, rho0_b60, and Lorene::Eos::sauve().

◆ set_auxiliary()

void Lorene::Eos_strange::set_auxiliary ( )
protected

Computes the auxiliary quantities n0 , rh0 , b34 and fach from the values of the other parameters.

Definition at line 194 of file eos_strange.C.

References b34, b60, eps_fit, fach, n0, n0_b60, Lorene::pow(), rho0, and rho0_b60.

◆ set_name()

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

Sets the EOS name.

Definition at line 173 of file eos.C.

References Lorene::Eos::name.

Friends And Related Function Documentation

◆ Eos::eos_from_file

Eos* Eos::eos_from_file ( FILE *  )
friend

The construction functions from a file.

Member Data Documentation

◆ b34

double Lorene::Eos_strange::b34
protected

$B_{60}^{3/4}$

Definition at line 2043 of file eos.h.

◆ b60

double Lorene::Eos_strange::b60
protected

Bag constant [unit: $60\ {\rm MeV\, fm}^{-3}$].

Definition at line 2011 of file eos.h.

◆ ent0

double Lorene::Eos_strange::ent0
protected

Log-enthalpy threshold for setting the energy density to a non zero value (should be negative).

Definition at line 2016 of file eos.h.

◆ eps_fit

double Lorene::Eos_strange::eps_fit
protected

Fitting parameter $\epsilon_{\rm fit}$ related to the square of sound velocity by $c_s^2 = 1/3(1+\epsilon_{\rm fit})$.

[cf. Zdunik, Astron. Astrophys. 359 , 311 (2000)]

Definition at line 2022 of file eos.h.

◆ fach

double Lorene::Eos_strange::fach
protected

Factor $(4+\epsilon_{\rm fit})/(1+\epsilon_{\rm fit})$.

Definition at line 2048 of file eos.h.

◆ n0

double Lorene::Eos_strange::n0
protected

Baryon density at zero pressure.

[unit: $0.1{\rm \ fm}^{-3}$ (Lorene's unit)]

Definition at line 2032 of file eos.h.

◆ n0_b60

double Lorene::Eos_strange::n0_b60
protected

Baryon density at zero pressure divided by $B_{60}^{3/4}$.

[unit: ${\rm fm}^{-3}$]

Definition at line 2008 of file eos.h.

◆ name

char Lorene::Eos::name[100]
protectedinherited

EOS name.

Definition at line 215 of file eos.h.

◆ rho0

double Lorene::Eos_strange::rho0
protected

Energy density at zero pressure.

[unit: $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$ (Lorene's unit)]

Definition at line 2038 of file eos.h.

◆ rho0_b60

double Lorene::Eos_strange::rho0_b60
protected

Energy density at zero pressure divided by $B_{60}$.

[unit: ${\rm MeV\, fm^{-3}}$]

Definition at line 2027 of file eos.h.

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