LORENE
|
Equation of state for fitting the 1-parameter EoSs from the CompOSE database.
More...
#include <eos_compose_fit.h>
Public Member Functions | |
Eos_compose_fit (const string ¶m_file) | |
Constructor from a parameter file. More... | |
virtual | ~Eos_compose_fit () |
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... | |
void | write_lorene_table (const string &, int nlines=200) const |
Save into a table in Lorene format. More... | |
virtual double | nbar_ent_p (double ent, const Param *par=0x0) const |
Higher bound in density. 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 *par=0x0) 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_compose_fit (FILE *) | |
Constructor from a binary file (created by the function sauve(FILE*) ). More... | |
Eos_compose_fit (ifstream &) | |
Constructor from a formatted file. More... | |
void | read_and_compute (ifstream &) |
Reads the Compose data and makes the fit. More... | |
void | read_compose_data (int &nbp, Tbl *&logh, Tbl *&logp, Tbl *&loge, Tbl *&lognb, Tbl *&gam1) |
Reads Compose data and stores the values of thermodynamic quantities. More... | |
void | adiabatic_index_fit (int i_min, int i_max, const Tbl &logh_read, const Tbl &logp_read, const Tbl &loge_read, const Tbl &lognb_read, const Tbl &gam1_read) |
From the read values, makes the fit on the adiabatic index and deduces the other quantities from thermodynamic relations and definitions. More... | |
double | integrate_equations (double kappa, double mean_gam, const Scalar &gamma1, const Scalar &log_ent, const Scalar &ent, Scalar &YYY) |
Integrates the differential system giving all thermo quantities from the adiabatic index. More... | |
virtual ostream & | operator>> (ostream &) const |
Operator >> More... | |
const string & | get_tablename () const |
Returns the name (given in the parameter file, see the introduction of the class). More... | |
double | get_nbmin () const |
double | get_nbmax () const |
Lower bound in density. 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... | |
int | n_coefs |
Number of coeficients for polynomial regression. More... | |
double | nb_min |
Lower bound in baryon density, below which the EoS is assumed to be a polytrope. More... | |
double | nb_mid |
Middle bound in baryon density, above which which the EoS is determined from the polynomial fit to the adiabatic index. More... | |
double | nb_max |
Higher bound on the density, above which the EoS is assumed to be a polytrope. More... | |
double | hmin |
Values of enthalpy corresponding to nb_min and nb_max. More... | |
double | hmax |
const Eos_poly * | p_eos_low |
Pointer on a polytropic EoS for the description of low densities (nb<nb_min) More... | |
const Eos_poly * | p_eos_high |
Pointer on a polytropic EoS for the description of high densities (nb>nb_max) More... | |
const Mg3d * | mg |
Multi-grid defining the number of domains and spectral coefficients. More... | |
const Map_af * | mp |
Mapping in . More... | |
Scalar * | log_p |
Table of . More... | |
Scalar * | log_e |
Table of . More... | |
Scalar * | log_nb |
Table of . More... | |
Scalar * | log_cs2 |
Table of . More... | |
char | name [100] |
EOS name. More... | |
Private Member Functions | |
Eos_compose_fit (const Eos_compose_fit &) | |
Copy constructor (private to make Eos_compose_fit a non-copiable class) More... | |
Friends | |
Eos * | Eos::eos_from_file (FILE *) |
The construction functions from a file. More... | |
Eos * | Eos::eos_from_file (ifstream &) |
Equation of state for fitting the 1-parameter EoSs from the CompOSE database.
()
A polynomial fit is done ont the adidabatic index for the highest densities. At the lowest ones, a simple polytrope is fitted. In between, a linear interpolation is done in terms of adiabatic index. When built with Eos::eos_from_file()
, the file must be composed of the following lines:
7 Type of the EOS EoS fitted from CompOSE data /full/path/to/the/eos/table/ 0.01 1.e-9 # limiting density values 6 # polynomial degree for the regression 33 # number of Chebyshev point sin each domain
On the second line the name is just for output. The path to the directory containing CompOSE EoS (third line), does not contain the name of the tables, it is assumed to be called eos.nb
and eos.thermo
(see CompOSE documentation). The fourth line gives limit values in density for the fit: above the first value, the polynomial regression is done. Below the second value, a polytrope is fitted.
Definition at line 90 of file eos_compose_fit.h.
Lorene::Eos_compose_fit::Eos_compose_fit | ( | const string & | param_file | ) |
Constructor from a parameter file.
files_path | Absolute name (including path) of the parameter file. |
Definition at line 71 of file eos_compose_fit.C.
References read_and_compute().
|
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 86 of file eos_compose_fit.C.
References Lorene::fread_be(), hmin, log_cs2, log_e, log_nb, log_p, mg, mp, n_coefs, nb_max, nb_mid, nb_min, p_eos_high, p_eos_low, 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& ) .
Definition at line 129 of file eos_compose_fit.C.
References read_and_compute().
|
private |
Copy constructor (private to make Eos_compose_fit
a non-copiable class)
|
virtual |
Destructor.
Definition at line 142 of file eos_compose_fit.C.
References log_cs2, log_e, log_nb, log_p, mg, mp, p_eos_high, and p_eos_low.
|
protected |
From the read values, makes the fit on the adiabatic index and deduces the other quantities from thermodynamic relations and definitions.
Definition at line 200 of file eos_compose_fit_build.C.
References Lorene::Valeur::c, Lorene::Valeur::c_cf, Lorene::Valeur::coef(), Lorene::Valeur::coef_i(), Lorene::Scalar::domain(), Lorene::exp(), Lorene::Tbl::get_dim(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nzone(), integrate_equations(), Lorene::log(), log_cs2, log_e, log_nb, log_p, mg, mp, n_coefs, p_eos_high, Lorene::poly_regression(), Lorene::pow(), Lorene::Map::r, Lorene::Tbl::set(), Lorene::Scalar::set_domain(), Lorene::Tbl::set_etat_qcq(), Lorene::Scalar::set_spectral_va(), Lorene::Scalar::std_spectral_base(), and Lorene::Scalar::val_grid_point().
|
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().
|
virtual |
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 481 of file eos_compose_fit.C.
References Lorene::Eos_poly::csound_square_ent_p(), Lorene::exp(), hmin, Lorene::log(), log_cs2, p_eos_high, p_eos_low, and Lorene::Scalar::val_point().
|
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 446 of file eos_compose_fit.C.
References csound_square_ent_p(), ener_ent_p(), and press_ent_p().
|
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 437 of file eos_compose_fit.C.
References csound_square_ent_p().
|
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 459 of file eos_compose_fit.C.
References ener_ent_p(), and press_ent_p().
|
virtual |
Computes the logarithmic derivative from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Definition at line 472 of file eos_compose_fit.C.
References der_nbar_ent_p(), and der_press_ent_p().
|
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 395 of file eos_compose_fit.C.
References Lorene::Eos_poly::ener_ent_p(), Lorene::exp(), hmin, Lorene::log(), log_e, p_eos_high, p_eos_low, and Lorene::Scalar::val_point().
|
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 |
|
inlineprotected |
|
inlineprotected |
Returns the name (given in the parameter file, see the introduction of the class).
Definition at line 240 of file eos_compose_fit.h.
References tablename.
|
virtual |
Returns a number to identify the sub-classe of Eos
the object belongs to.
Implements Lorene::Eos.
Definition at line 151 of file eos_from_file.C.
|
protected |
Integrates the differential system giving all thermo quantities from the adiabatic index.
Definition at line 364 of file eos_compose_fit_build.C.
References Lorene::log(), Lorene::Scalar::primr(), Lorene::Scalar::std_spectral_base(), and Lorene::Scalar::val_grid_point().
|
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 |
Higher bound in density.
Computes the baryon density from the log-enthalpy.
ent | [input, unit: ] log-enthalpy H |
Implements Lorene::Eos.
Definition at line 374 of file eos_compose_fit.C.
References Lorene::exp(), hmin, Lorene::log(), log_nb, Lorene::Eos_poly::nbar_ent_p(), p_eos_high, p_eos_low, and Lorene::Scalar::val_point().
|
virtual |
Comparison operator (difference)
Implements Lorene::Eos.
Definition at line 169 of file eos_compose_fit.C.
References operator==().
|
virtual |
Comparison operator (egality)
Implements Lorene::Eos.
Definition at line 158 of file eos_compose_fit.C.
References identify(), and Lorene::Eos::identify().
|
protectedvirtual |
Operator >>
Implements Lorene::Eos.
Definition at line 218 of file eos_compose_fit.C.
References Lorene::exp(), hmin, mg, mp, n_coefs, nb_max, nb_mid, nb_min, p_eos_high, p_eos_low, tablename, and Lorene::Map_af::val_r_jk().
|
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 416 of file eos_compose_fit.C.
References Lorene::exp(), hmin, Lorene::log(), log_p, p_eos_high, p_eos_low, Lorene::Eos_poly::press_ent_p(), and Lorene::Scalar::val_point().
|
protected |
Reads the Compose data and makes the fit.
Definition at line 239 of file eos_compose_fit.C.
References Lorene::exp(), hmin, mg, n_coefs, nb_mid, nb_min, read_compose_data(), Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), and tablename.
|
protected |
Reads Compose data and stores the values of thermodynamic quantities.
Definition at line 63 of file eos_compose_fit_build.C.
References Lorene::compute_derivative(), Lorene::log(), nb_max, Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), and tablename.
|
virtual |
Save in a file.
Reimplemented from Lorene::Eos.
Definition at line 178 of file eos_compose_fit.C.
References Lorene::fwrite_be(), Lorene::Eos_poly::get_gam(), Lorene::Eos_poly::get_kap(), Lorene::Eos_poly::get_m_0(), Lorene::Eos_poly::get_mu_0(), hmin, log_cs2, log_e, log_nb, log_p, mg, mp, n_coefs, nb_max, nb_mid, nb_min, p_eos_high, p_eos_low, Lorene::Eos::sauve(), Lorene::Mg3d::sauve(), Lorene::Scalar::sauve(), Lorene::Map_af::sauve(), and tablename.
|
inherited |
void Lorene::Eos_compose_fit::write_lorene_table | ( | const string & | name_file, |
int | nlines = 200 |
||
) | const |
Save into a table in Lorene format.
Definition at line 301 of file eos_compose_fit.C.
References ener_ent_p(), Lorene::exp(), Lorene::Mg3d::get_nzone(), hmin, Lorene::log(), log_e, log_nb, log_p, mg, nbar_ent_p(), press_ent_p(), and tablename.
|
friend |
The construction functions from a file.
|
protected |
Values of enthalpy corresponding to nb_min and nb_max.
Definition at line 117 of file eos_compose_fit.h.
|
protected |
Table of .
Definition at line 141 of file eos_compose_fit.h.
|
protected |
Table of .
Definition at line 135 of file eos_compose_fit.h.
|
protected |
Table of .
Definition at line 138 of file eos_compose_fit.h.
|
protected |
Table of .
Definition at line 132 of file eos_compose_fit.h.
|
protected |
Multi-grid defining the number of domains and spectral coefficients.
Definition at line 126 of file eos_compose_fit.h.
|
protected |
Mapping in .
Definition at line 129 of file eos_compose_fit.h.
|
protected |
Number of coeficients for polynomial regression.
Note: this is in general different from the number of spectral coefficients used for the representation of thermodynamical fields.
Definition at line 103 of file eos_compose_fit.h.
|
protected |
Higher bound on the density, above which the EoS is assumed to be a polytrope.
Definition at line 114 of file eos_compose_fit.h.
|
protected |
Middle bound in baryon density, above which which the EoS is determined from the polynomial fit to the adiabatic index.
Definition at line 111 of file eos_compose_fit.h.
|
protected |
Lower bound in baryon density, below which the EoS is assumed to be a polytrope.
Definition at line 106 of file eos_compose_fit.h.
|
protected |
Pointer on a polytropic EoS for the description of high densities (nb>nb_max)
Definition at line 123 of file eos_compose_fit.h.
|
protected |
Pointer on a polytropic EoS for the description of low densities (nb<nb_min)
Definition at line 120 of file eos_compose_fit.h.
|
protected |
Name of the file containing the tabulated data.
Definition at line 97 of file eos_compose_fit.h.