LORENE
|
Equation of state base class. More...
#include <eos.h>
Public Member Functions | |
virtual | ~Eos () |
Destructor. | |
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. | |
virtual void | sauve (FILE *) const |
Save in a file. | |
virtual double | nbar_ent_p (double ent, const Param *par=0x0) const =0 |
Computes the baryon density from the log-enthalpy and extra parameters (virtual function implemented in the derived classes). | |
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. | |
virtual double | ener_ent_p (double ent, const Param *par=0x0) const =0 |
Computes the total energy density from the log-enthalpy and extra parameters (virtual function implemented in the derived classes). | |
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. | |
virtual double | press_ent_p (double ent, const Param *par=0x0) const =0 |
Computes the pressure from the log-enthalpy and extra parameters (virtual function implemented in the derived classes). | |
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. | |
virtual double | der_nbar_ent_p (double ent, const Param *par=0x0) const =0 |
Computes the logarithmic derivative ![]() | |
Cmp | der_nbar_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the logarithmic derivative ![]() | |
Scalar | der_nbar_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the logarithmic derivative ![]() | |
virtual double | der_ener_ent_p (double ent, const Param *par=0x0) const =0 |
Computes the logarithmic derivative ![]() | |
Cmp | der_ener_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the logarithmic derivative ![]() | |
Scalar | der_ener_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the logarithmic derivative ![]() | |
virtual double | der_press_ent_p (double ent, const Param *par=0x0) const =0 |
Computes the logarithmic derivative ![]() | |
Cmp | der_press_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the logarithmic derivative ![]() | |
Scalar | der_press_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the logarithmic derivative ![]() | |
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 () | |
Standard constructor. | |
Eos (const char *name_i) | |
Standard constructor with name. | |
Eos (const Eos &) | |
Copy constructor | |
Eos (FILE *) | |
Constructor from a binary file (created by the function sauve(FILE*) ). | |
Eos (ifstream &) | |
Constructor from a formatted file. | |
virtual ostream & | operator>> (ostream &) const =0 |
Operator >> | |
void | calcule (const Cmp &thermo, int nzet, int l_min, double(Eos::*fait)(double, const Param *) const, const Param *par, Cmp &resu) const |
General computational method for Cmp 's. | |
void | calcule (const Scalar &thermo, int nzet, int l_min, double(Eos::*fait)(double, const Param *) const, const Param *par, Scalar &resu) const |
General computational method for Scalar 's. | |
Protected Attributes | |
char | name [100] |
EOS name. | |
Friends | |
ostream & | operator<< (ostream &, const Eos &) |
Display. | |
|
protected |
|
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.C.
References name.
|
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 140 of file eos.C.
References name.
|
protected |
General computational method for Cmp
's.
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. |
|
protected |
General computational method for Scalar
's.
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. |
Cmp Lorene::Eos::der_ener_ent | ( | const Cmp & | ent, |
int | nzet, | ||
int | l_min = 0 , |
||
const Param * | par = 0x0 |
||
) | const |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 430 of file eos.C.
References Lorene::Cmp::get_mp().
Scalar Lorene::Eos::der_ener_ent | ( | const Scalar & | ent, |
int | nzet, | ||
int | l_min = 0 , |
||
const Param * | par = 0x0 |
||
) | const |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 440 of file eos.C.
References Lorene::Tensor::get_mp().
|
pure virtual |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
par | possible extra parameters of the EOS |
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_fitting, Lorene::Eos_mag, Lorene::Eos_multi_poly, and Lorene::Eos_tabul.
Cmp Lorene::Eos::der_nbar_ent | ( | const Cmp & | ent, |
int | nzet, | ||
int | l_min = 0 , |
||
const Param * | par = 0x0 |
||
) | const |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 407 of file eos.C.
References Lorene::Cmp::get_mp().
Scalar Lorene::Eos::der_nbar_ent | ( | const Scalar & | ent, |
int | nzet, | ||
int | l_min = 0 , |
||
const Param * | par = 0x0 |
||
) | const |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 417 of file eos.C.
References Lorene::Tensor::get_mp().
|
pure virtual |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
par | possible extra parameters of the EOS |
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_fitting, Lorene::Eos_mag, Lorene::Eos_multi_poly, and Lorene::Eos_tabul.
Cmp Lorene::Eos::der_press_ent | ( | const Cmp & | ent, |
int | nzet, | ||
int | l_min = 0 , |
||
const Param * | par = 0x0 |
||
) | const |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
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. |
Definition at line 452 of file eos.C.
References Lorene::Cmp::get_mp().
Scalar Lorene::Eos::der_press_ent | ( | const Scalar & | ent, |
int | nzet, | ||
int | l_min = 0 , |
||
const Param * | par = 0x0 |
||
) | const |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
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. |
Definition at line 462 of file eos.C.
References Lorene::Tensor::get_mp().
|
pure virtual |
Computes the logarithmic derivative
ent | [input, unit: ![]() ![]() ![]() |
par | possible extra parameters of the EOS |
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_fitting, Lorene::Eos_mag, Lorene::Eos_multi_poly, and Lorene::Eos_tabul.
Cmp Lorene::Eos::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.
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 363 of file eos.C.
References Lorene::Cmp::get_mp().
Scalar Lorene::Eos::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.
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 373 of file eos.C.
References Lorene::Tensor::get_mp().
Computes the total energy density from the log-enthalpy and extra parameters (virtual function implemented in the derived classes).
ent | [input, unit: ![]() ![]() ![]() |
par | possible extra parameters of the EOS |
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_consistent, Lorene::Eos_fitting, Lorene::Eos_mag, Lorene::Eos_multi_poly, and Lorene::Eos_tabul.
Construction of an EOS from a binary file.
The file must have been created by the function sauve(FILE*)
.
Definition at line 177 of file eos_from_file.C.
References Lorene::fread_be().
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:
Eos_poly
).Eos_poly_newt
).Eos_incomp
).Eos_incomp_newt
).Eos_Fermi
)MEos
)Eos_multi_poly
)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 314 of file eos_from_file.C.
Returns a number to identify the sub-classe of Eos
the object belongs to.
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_CompOSE, Lorene::Eos_consistent, Lorene::Eos_fit_SLy4, Lorene::Eos_fit_FPS, Lorene::Eos_fit_AkmalPR, Lorene::Eos_mag, Lorene::Eos_multi_poly, Lorene::Eos_SLy4, Lorene::Eos_FPS, Lorene::Eos_BPAL12, Lorene::Eos_AkmalPR, Lorene::Eos_BBB2, Lorene::Eos_BalbN1H1, and Lorene::Eos_GlendNH3.
Cmp Lorene::Eos::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.
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 338 of file eos.C.
References Lorene::Cmp::get_mp().
Scalar Lorene::Eos::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.
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 348 of file eos.C.
References Lorene::Tensor::get_mp().
Computes the baryon density from the log-enthalpy and extra parameters (virtual function implemented in the derived classes).
ent | [input, unit: ![]() ![]() ![]() |
par | possible extra parameters of the EOS |
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_consistent, Lorene::Eos_fitting, Lorene::Eos_mag, Lorene::Eos_multi_poly, and Lorene::Eos_tabul.
Comparison operator (difference)
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_CompOSE, Lorene::Eos_consistent, Lorene::Eos_fit_SLy4, Lorene::Eos_fit_FPS, Lorene::Eos_fit_AkmalPR, Lorene::Eos_mag, Lorene::Eos_multi_poly, Lorene::Eos_SLy4, Lorene::Eos_FPS, Lorene::Eos_BPAL12, Lorene::Eos_AkmalPR, Lorene::Eos_BBB2, Lorene::Eos_BalbN1H1, and Lorene::Eos_GlendNH3.
Comparison operator (egality)
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_CompOSE, Lorene::Eos_consistent, Lorene::Eos_fit_SLy4, Lorene::Eos_fit_FPS, Lorene::Eos_fit_AkmalPR, Lorene::Eos_mag, Lorene::Eos_multi_poly, Lorene::Eos_SLy4, Lorene::Eos_FPS, Lorene::Eos_BPAL12, Lorene::Eos_AkmalPR, Lorene::Eos_BBB2, Lorene::Eos_BalbN1H1, and Lorene::Eos_GlendNH3.
Operator >>
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_CompOSE, Lorene::Eos_consistent, Lorene::Eos_fit_SLy4, Lorene::Eos_fit_FPS, Lorene::Eos_fit_AkmalPR, Lorene::Eos_mag, Lorene::Eos_multi_poly, Lorene::Eos_SLy4, Lorene::Eos_FPS, Lorene::Eos_BPAL12, Lorene::Eos_AkmalPR, Lorene::Eos_BBB2, Lorene::Eos_BalbN1H1, and Lorene::Eos_GlendNH3.
Cmp Lorene::Eos::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.
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 385 of file eos.C.
References Lorene::Cmp::get_mp().
Scalar Lorene::Eos::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.
ent | [input, unit: ![]() ![]() ![]() |
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 |
Definition at line 395 of file eos.C.
References Lorene::Tensor::get_mp().
Computes the pressure from the log-enthalpy and extra parameters (virtual function implemented in the derived classes).
ent | [input, unit: ![]() ![]() ![]() |
par | possible extra parameters of the EOS |
Implemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_consistent, Lorene::Eos_fitting, Lorene::Eos_mag, Lorene::Eos_multi_poly, and Lorene::Eos_tabul.
Save in a file.
Reimplemented in Lorene::Eos_poly, Lorene::Eos_poly_newt, Lorene::Eos_incomp, Lorene::Eos_incomp_newt, Lorene::Eos_strange, Lorene::Eos_strange_cr, Lorene::Eos_Fermi, Lorene::MEos, Lorene::Eos_fitting, Lorene::Eos_mag, Lorene::Eos_multi_poly, and Lorene::Eos_tabul.
Definition at line 179 of file eos.C.
References Lorene::fwrite_be(), identify(), and name.