EOS with domain dependency. More...
#include <eos.h>
Public Member Functions | |
MEos (int ndom_i, const Eos **mono_eos_i) | |
Standard constructor. | |
MEos (const Eos &eos1, const Eos &eos2) | |
Constructor for 2 domains. | |
MEos (const Eos &eos1, const Eos &eos2, const Eos &eos3) | |
Constructor for 3 domains. | |
MEos (const Eos &eos1, const Eos &eos2, const Eos &eos3, const Eos &eos4) | |
Constructor for 4 domains. | |
MEos (const MEos &) | |
Copy constructor. | |
virtual | ~MEos () |
Destructor. | |
void | operator= (const MEos &) |
Assignment to another MEos . | |
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 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 from the log-enthalpy. | |
virtual double | der_ener_ent_p (double ent, const Param *par=0x0) const |
Computes the logarithmic derivative from the log-enthalpy. | |
virtual double | der_press_ent_p (double ent, const Param *par=0x0) const |
Computes the logarithmic derivative from the log-enthalpy. | |
const char * | get_name () const |
Returns the EOS name. | |
void | set_name (const char *name_i) |
Sets the EOS name. | |
virtual bool | operator== (const Eos &) const =0 |
Comparison operator (egality). | |
virtual bool | operator!= (const Eos &) const =0 |
Comparison operator (difference). | |
Cmp | nbar_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the baryon density field from the log-enthalpy field and extra parameters. | |
Scalar | nbar_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the baryon density field from the log-enthalpy field and extra parameters. | |
Cmp | ener_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the total energy density from the log-enthalpy and extra parameters. | |
Scalar | ener_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the total energy density from the log-enthalpy and extra parameters. | |
Cmp | press_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the pressure from the log-enthalpy and extra parameters. | |
Scalar | press_ent (const Scalar &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the pressure from the log-enthalpy and extra parameters. | |
Cmp | der_nbar_ent (const Cmp &ent, int nzet, int l_min=0, const Param *par=0x0) const |
Computes the logarithmic derivative 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 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 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 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 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 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 | |
MEos (FILE *) | |
Constructor from a binary file (created by the function sauve(FILE*) ). | |
MEos (ifstream &) | |
Constructor from a formatted file. | |
virtual ostream & | operator>> (ostream &) const |
Operator >>. | |
void | calcule (const Cmp &thermo, int nzet, int l_min, double(Eos::*fait)(double, const Param *) const, const Param *par, Cmp &resu) const |
General computational method for Cmp 's. | |
void | calcule (const Scalar &thermo, int nzet, int l_min, double(Eos::*fait)(double, const Param *) const, const Param *par, Scalar &resu) const |
General computational method for Scalar 's. | |
Protected Attributes | |
const Eos ** | mono_eos |
Array (upon the domains) containing the various EOS. | |
int | ndom |
Number of domains. | |
bool | constructed_from_file |
Indicates wether the EOS has been constructed from a file. | |
char | name [100] |
EOS name. | |
Friends | |
Eos * | Eos::eos_from_file (FILE *) |
Eos * | Eos::eos_from_file (ifstream &) |
ostream & | operator<< (ostream &, const Eos &) |
Display. |
EOS with domain dependency.
Definition at line 2269 of file eos.h.
MEos::MEos | ( | int | ndom_i, | |
const Eos ** | mono_eos_i | |||
) |
MEos::MEos | ( | const MEos & | meos | ) |
MEos::MEos | ( | 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 129 of file meos.C.
References Eos::eos_from_file(), fread_be(), mono_eos, and ndom.
MEos::MEos | ( | 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&)
. The construction functions from a file
Definition at line 141 of file meos.C.
References Eos::eos_from_file(), mono_eos, and ndom.
MEos::~MEos | ( | ) | [virtual] |
Destructor.
Definition at line 158 of file meos.C.
References constructed_from_file, mono_eos, and ndom.
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.
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.
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 from the log-enthalpy and extra parameters.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
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 from the log-enthalpy and extra parameters.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
Definition at line 426 of file eos.C.
References Eos::calcule(), Eos::der_ener_ent_p(), and Cmp::get_mp().
double MEos::der_ener_ent_p | ( | double | ent, | |
const Param * | par = 0x0 | |||
) | const [virtual] |
Computes the logarithmic derivative from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H defined by Eq. (4) | |
par | possible extra parameters of the EOS |
Implements Eos.
Definition at line 301 of file meos.C.
References Eos::der_ener_ent_p(), Param::get_int(), and mono_eos.
Scalar Eos::der_nbar_ent | ( | const Scalar & | ent, | |
int | nzet, | |||
int | l_min = 0 , |
|||
const Param * | par = 0x0 | |||
) | const [inherited] |
Computes the logarithmic derivative from the log-enthalpy and extra parameters.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
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 from the log-enthalpy and extra parameters.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
Definition at line 403 of file eos.C.
References Eos::calcule(), Eos::der_nbar_ent_p(), and Cmp::get_mp().
double MEos::der_nbar_ent_p | ( | double | ent, | |
const Param * | par = 0x0 | |||
) | const [virtual] |
Computes the logarithmic derivative from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H defined by Eq. (4) | |
par | possible extra parameters of the EOS |
Implements Eos.
Definition at line 291 of file meos.C.
References Eos::der_nbar_ent_p(), Param::get_int(), and mono_eos.
Scalar Eos::der_press_ent | ( | const Scalar & | ent, | |
int | nzet, | |||
int | l_min = 0 , |
|||
const Param * | par = 0x0 | |||
) | const [inherited] |
Computes the logarithmic derivative from the log-enthalpy and extra parameters.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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. |
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 from the log-enthalpy and extra parameters.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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. |
Definition at line 448 of file eos.C.
References Eos::calcule(), Eos::der_press_ent_p(), and Cmp::get_mp().
double MEos::der_press_ent_p | ( | double | ent, | |
const Param * | par = 0x0 | |||
) | const [virtual] |
Computes the logarithmic derivative from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H defined by Eq. (4) | |
par | possible extra parameters of the EOS |
Implements Eos.
Definition at line 311 of file meos.C.
References Eos::der_press_ent_p(), Param::get_int(), and mono_eos.
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.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
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.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
Definition at line 359 of file eos.C.
References Eos::calcule(), Eos::ener_ent_p(), and Cmp::get_mp().
double MEos::ener_ent_p | ( | double | ent, | |
const Param * | par = 0x0 | |||
) | const [virtual] |
Computes the total energy density from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H defined by Eq. (4) | |
par | possible extra parameters of the EOS |
Implements Eos.
Definition at line 271 of file meos.C.
References Eos::ener_ent_p(), Param::get_int(), and mono_eos.
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:
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 297 of file eos_from_file.C.
Eos * Eos::eos_from_file | ( | FILE * | fich | ) | [static, inherited] |
Construction of an EOS from a binary file.
The file must have been created by the function sauve(FILE*)
.
Definition at line 165 of file eos_from_file.C.
References fread_be().
const char * Eos::get_name | ( | ) | const [inherited] |
int MEos::identify | ( | ) | const [virtual] |
Returns a number to identify the sub-classe of Eos
the object belongs to.
Implements Eos.
Definition at line 151 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.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
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.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
Definition at line 334 of file eos.C.
References Eos::calcule(), Cmp::get_mp(), and Eos::nbar_ent_p().
double MEos::nbar_ent_p | ( | double | ent, | |
const Param * | par = 0x0 | |||
) | const [virtual] |
Computes the baryon density from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H defined by Eq. (4) | |
par | possible extra parameters of the EOS |
Implements Eos.
Definition at line 260 of file meos.C.
References Param::get_int(), mono_eos, and Eos::nbar_ent_p().
virtual bool Eos::operator!= | ( | const Eos & | ) | const [pure virtual, inherited] |
Comparison operator (difference).
bool MEos::operator!= | ( | const Eos & | eos_i | ) | const [virtual] |
void MEos::operator= | ( | const MEos & | ) |
virtual bool Eos::operator== | ( | const Eos & | ) | const [pure virtual, inherited] |
Comparison operator (egality).
bool MEos::operator== | ( | const Eos & | eos_i | ) | const [virtual] |
Comparison operator (egality).
Definition at line 218 of file meos.C.
References identify(), Eos::identify(), mono_eos, and ndom.
ostream & MEos::operator>> | ( | ostream & | ost | ) | const [protected, virtual] |
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.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
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.
ent | [input, unit: ] log-enthalpy H defined by , where e is the (total) energy density, p the pressure, n the baryon density, and 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 |
Definition at line 381 of file eos.C.
References Eos::calcule(), Cmp::get_mp(), and Eos::press_ent_p().
double MEos::press_ent_p | ( | double | ent, | |
const Param * | par = 0x0 | |||
) | const [virtual] |
Computes the pressure from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H defined by Eq. (4) | |
par | possible extra parameters of the EOS |
Implements Eos.
Definition at line 281 of file meos.C.
References Param::get_int(), mono_eos, and Eos::press_ent_p().
void MEos::sauve | ( | FILE * | fich | ) | const [virtual] |
Save in a file.
Reimplemented from Eos.
Definition at line 186 of file meos.C.
References fwrite_be(), mono_eos, ndom, and Eos::sauve().
void Eos::set_name | ( | const char * | name_i | ) | [inherited] |
ostream& operator<< | ( | ostream & | , | |
const Eos & | ||||
) | [friend, inherited] |
Display.
bool MEos::constructed_from_file [protected] |
const Eos** MEos::mono_eos [protected] |
int MEos::ndom [protected] |