LORENE
|
Class for a magnetized (tabulated) equation of state. More...
#include <eos_mag.h>
Public Member Functions | |
virtual | ~Eos_mag () |
Destructor. 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... | |
virtual void | sauve (FILE *) const |
Save in a file. More... | |
double | mag_ent_p (double ent, const Param *par=0x0) const |
Computes the magnetisation. 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 from the log-enthalpy. More... | |
virtual double | der_ener_ent_p (double ent, const Param *par=0x0) const |
Computes the logarithmic derivative from the log-enthalpy. More... | |
virtual double | der_press_ent_p (double ent, const Param *par=0x0) const |
Computes the logarithmic derivative from the log-enthalpy. More... | |
virtual double | der_press_nbar_p (double ent, const Param *par=0x0) const |
Computes the logarithmic derivative from the log-enthalpy. More... | |
virtual double | csound_square_ent_p (double, const Param *) const |
Computes the sound speed squared 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 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 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 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 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 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 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 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_mag (const char *name_i, const char *table, const char *path) | |
Standard constructor. More... | |
Eos_mag (const char *name_i, const char *file_name) | |
Standard constructor from the full filename. More... | |
Eos_mag (const Eos_mag &) | |
Copy constructor. More... | |
Eos_mag (FILE *) | |
Constructor from a binary file (created by the function sauve(FILE*) ). More... | |
Eos_mag (ifstream &ist, const char *table) | |
Constructor from a formatted file. More... | |
Eos_mag (ifstream &ist) | |
Constructor from a formatted file. More... | |
void | read_table () |
Reads the file containing the table and initializes in the arrays logh , logp and dlpsdlh . 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 | |
string | tablename |
Name of the file containing the tabulated data. More... | |
double | hmin |
Lower boundary of the log-enthalpy interval. More... | |
double | hmax |
Upper boundary of the log-enthalpy interval. More... | |
double | Bmax |
Upper boundary of the magnetic field interval. More... | |
Tbl * | logp |
Table of . More... | |
Tbl * | logh |
Table of . More... | |
Tbl * | Bfield |
Table of . More... | |
Tbl * | dlpsdlh |
Table of . More... | |
Tbl * | dlpsdB |
Table of . More... | |
Tbl * | d2lp |
Table of . 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 &) |
Class for a magnetized (tabulated) equation of state.
()
This EOS depends on two variables : density (or log-enthalpy) and the magnetic field amplitude. The interpolation through the tables is a cubic Hermite interpolation, which is thermodynamically consistent, i.e. preserves the Gibbs-Duhem relation.
|
protected |
Standard constructor.
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 101 of file eos_mag.C.
References read_table(), and tablename.
|
protected |
Standard constructor from the full filename.
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 113 of file eos_mag.C.
References read_table().
|
protected |
Copy constructor.
|
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 123 of file eos_mag.C.
References read_table(), and tablename.
|
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& ) .
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 138 of file eos_mag.C.
References read_table(), and tablename.
|
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& ) .
ist | input file stream containing a name as first line and the full filename (including the path) containing the EOS file as second line |
Definition at line 148 of file eos_mag.C.
References read_table(), and tablename.
|
virtual |
|
protectedinherited |
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 213 of file eos.C.
References Lorene::Cmp::get_etat().
|
protectedinherited |
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 285 of file eos.C.
References Lorene::Scalar::get_etat().
|
inherited |
Computes the sound speed squared from the enthalpy with extra parameters.
ent | [input, unit: c^2] enthalpy |
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 499 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::csound_square_ent_p(), and Lorene::Tensor::get_mp().
|
inlinevirtual |
Computes the sound speed squared from the enthapy with extra parameters (virtual function implemented in the derived classes).
ent | [input, unit: c^2] enthalpy |
par | possible extra parameters of the EOS |
Implements Lorene::Eos.
Definition at line 293 of file eos_mag.h.
References Lorene::c_est_pas_fait().
|
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 454 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::der_ener_ent_p(), and Lorene::Cmp::get_mp().
|
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 464 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::der_ener_ent_p(), and Lorene::Tensor::get_mp().
|
virtual |
Computes the logarithmic derivative from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Implements Lorene::Eos.
Definition at line 523 of file eos_mag.C.
References Lorene::c_est_pas_fait().
|
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 431 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::der_nbar_ent_p(), and Lorene::Cmp::get_mp().
|
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 441 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::der_nbar_ent_p(), and Lorene::Tensor::get_mp().
|
virtual |
Computes the logarithmic derivative from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Implements Lorene::Eos.
Definition at line 511 of file eos_mag.C.
References Lorene::c_est_pas_fait().
|
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 476 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::der_press_ent_p(), and Lorene::Cmp::get_mp().
|
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 486 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::der_press_ent_p(), and Lorene::Tensor::get_mp().
|
virtual |
Computes the logarithmic derivative from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Implements Lorene::Eos.
Definition at line 536 of file eos_mag.C.
References Lorene::c_est_pas_fait().
|
virtual |
Computes the logarithmic derivative from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Definition at line 548 of file eos_mag.C.
References Lorene::c_est_pas_fait().
|
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 387 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::ener_ent_p(), and Lorene::Cmp::get_mp().
|
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 397 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Eos::ener_ent_p(), and Lorene::Tensor::get_mp().
|
virtual |
Computes the total energy density from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Implements Lorene::Eos.
Definition at line 376 of file eos_mag.C.
References Bfield, Bmax, d2lp, dlpsdB, dlpsdlh, Lorene::exp(), Lorene::Param::get_double_mod(), Lorene::Param::get_n_double_mod(), hmax, hmin, logh, and logp.
|
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().
|
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:
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 345 of file eos_from_file.C.
|
inherited |
|
virtual |
Returns a number to identify the sub-classe of Eos
the object belongs to.
Implements Lorene::Eos.
Definition at line 169 of file eos_from_file.C.
double Lorene::Eos_mag::mag_ent_p | ( | double | ent, |
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 362 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Cmp::get_mp(), and Lorene::Eos::nbar_ent_p().
|
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 372 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Tensor::get_mp(), and Lorene::Eos::nbar_ent_p().
|
virtual |
Computes the baryon density from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Implements Lorene::Eos.
Definition at line 333 of file eos_mag.C.
References Bfield, Bmax, d2lp, dlpsdB, dlpsdlh, Lorene::exp(), Lorene::Param::get_double_mod(), Lorene::Param::get_n_double_mod(), hmax, hmin, logh, and logp.
|
virtual |
Comparison operator (difference)
Implements Lorene::Eos.
Definition at line 199 of file eos_mag.C.
References operator==().
|
virtual |
Comparison operator (egality)
Implements Lorene::Eos.
Definition at line 186 of file eos_mag.C.
References identify(), and Lorene::Eos::identify().
|
protectedvirtual |
|
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 409 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Cmp::get_mp(), and Lorene::Eos::press_ent_p().
|
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 419 of file eos.C.
References Lorene::Eos::calcule(), Lorene::Tensor::get_mp(), and Lorene::Eos::press_ent_p().
|
virtual |
Computes the pressure from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Implements Lorene::Eos.
Definition at line 422 of file eos_mag.C.
References Bfield, Bmax, d2lp, dlpsdB, dlpsdlh, Lorene::Param::get_double_mod(), Lorene::Param::get_n_double_mod(), hmax, hmin, logh, and logp.
|
protected |
Reads the file containing the table and initializes in the arrays logh
, logp
and dlpsdlh
.
Definition at line 227 of file eos_mag.C.
References Bfield, Bmax, d2lp, dlpsdB, dlpsdlh, hmax, hmin, Lorene::log(), logh, logp, Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), and tablename.
|
virtual |
Save in a file.
Reimplemented from Lorene::Eos.
Definition at line 174 of file eos_mag.C.
References Lorene::Eos::sauve(), and tablename.
|
inherited |
|
friend |
The construction functions from a file.
|
protected |
|
protected |
|
protected |
|
protected |