Equation of state SLy4 (Douchin & Haensel 2001). More...

#include <eos_tabul.h>

Inheritance diagram for Lorene::Eos_SLy4:

Public Member Functions
	Eos_SLy4 (const char *path)
	Standard constructor. More...

virtual	~Eos_SLy4 ()
	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...

const string &	get_tablename () const

const Tbl &	get_logh () const

double	get_hmin () const
	Returns minimal value of log-enthalpy in the table. More...

double	get_hmax () const
	Returns maximal value of log-enthalpy in the table. More...

virtual void	sauve (FILE *) const
	Save in a file. More...

virtual double	nbar_ent_p (double ent, const Param *par=0x0) const
	Computes the baryon density from the log-enthalpy. More...

virtual double	ener_ent_p (double ent, const Param *par=0x0) const
	Computes the total energy density from the log-enthalpy. More...

virtual double	press_ent_p (double ent, const Param *par=0x0) const
	Computes the pressure from the log-enthalpy. More...

virtual double	der_nbar_ent_p (double ent, const Param *par=0x0) const
	Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy. More...

virtual double	der_ener_ent_p (double ent, const Param *par=0x0) const
	Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy. More...

virtual double	der_press_ent_p (double ent, const Param *par=0x0) const
	Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy. More...

virtual double	der_press_nbar_p (double ent, const Param *par=0x0) const
	Computes the logarithmic derivative $d\ln p/d\ln n$ from the log-enthalpy. More...

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

const char *	get_name () const
	Returns the EOS name. More...

void	set_name (const char *name_i)
	Sets the EOS name. More...

Cmp	nbar_ent (const Cmp &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the baryon density field from the log-enthalpy field and extra parameters. More...

Scalar	nbar_ent (const Scalar &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the baryon density field from the log-enthalpy field and extra parameters. More...

Cmp	ener_ent (const Cmp &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the total energy density from the log-enthalpy and extra parameters. More...

Scalar	ener_ent (const Scalar &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the total energy density from the log-enthalpy and extra parameters. More...

Cmp	press_ent (const Cmp &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the pressure from the log-enthalpy and extra parameters. More...

Scalar	press_ent (const Scalar &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the pressure from the log-enthalpy and extra parameters. More...

Cmp	der_nbar_ent (const Cmp &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy and extra parameters. More...

Scalar	der_nbar_ent (const Scalar &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy and extra parameters. More...

Cmp	der_ener_ent (const Cmp &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy and extra parameters. More...

Scalar	der_ener_ent (const Scalar &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy and extra parameters. More...

Cmp	der_press_ent (const Cmp &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy and extra parameters. More...

Scalar	der_press_ent (const Scalar &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy and extra parameters. More...

Scalar	csound_square_ent (const Scalar &ent, int nzet, int l_min=0, Param *par=0x0) const
	Computes the sound speed squared $c_s^2 = c^2 \frac{dp}{de}$ from the enthalpy with extra parameters. More...

Static Public Member Functions
static Eos *	eos_from_file (FILE *)
	Construction of an EOS from a binary file. More...

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

Protected Member Functions
	Eos_SLy4 (FILE *)
	Constructor from a binary file (created by the function `sauve(FILE*)` ). More...

	Eos_SLy4 (ifstream &)
	Constructor from a formatted file. More...

virtual ostream &	operator>> (ostream &) const
	Operator >> More...

virtual void	read_table ()
	Reads the file containing the table and initializes in the arrays `logh` , `logp` and `dlpsdlh` . 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...

string	authors
	Authors - reference for the table. More...

double	hmin
	Lower boundary of the enthalpy interval. More...

double	hmax
	Upper boundary of the enthalpy interval. More...

Tbl *	logh
	Table of $\log H$ . More...

Tbl *	logp
	Table of $\log p$ . More...

Tbl *	dlpsdlh
	Table of $d\log P/d\log H$ . More...

Tbl *	lognb
	Table of $\log n_b$ . More...

Tbl *	dlpsdlnb
	Table of $d\log P/d\log nb$ . More...

Tbl *	log_cs2
	Table of $\log c_s^2 = \log \left( c^2 \frac{d p}{d e} \right)$ . More...

double *	press

double *	nb

double *	ro

char	name [100]
	EOS name. More...

Private Member Functions
	Eos_SLy4 (const Eos_SLy4 &)
	Copy constructor (private to make `Eos_SLy4` 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 &)

Detailed Description

Equation of state SLy4 (Douchin & Haensel 2001).

()

Interior: neutrons, protons, electrons and muons described by the Skyrme Lyon 4 potential.

Inner crust: Douchin & Haensel 2001

Outer crust: Haensel & Pichon, Astron. Astrophys. 283 ,
313 (1994)

Definition at line 423 of file eos_tabul.h.

Constructor & Destructor Documentation

◆ Eos_SLy4() [1/4]

Lorene::Eos_SLy4::Eos_SLy4 ( const char * path )

Standard constructor.

Parameters

path	Path to the directory containing the EOS file

Definition at line 67 of file eos_sly4.C.

◆ Eos_SLy4() [2/4]

Lorene::Eos_SLy4::Eos_SLy4 ( 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 75 of file eos_sly4.C.

◆ Eos_SLy4() [3/4]

Lorene::Eos_SLy4::Eos_SLy4 ( ifstream & fich )

protected

Constructor from a formatted file.

This constructor is protected because any EOS construction from a formatted file must be done via the function Eos::eos_from_file(ifstream& ) .

Definition at line 81 of file eos_sly4.C.

◆ Eos_SLy4() [4/4]

Lorene::Eos_SLy4::Eos_SLy4 ( const Eos_SLy4 & )

private

Copy constructor (private to make Eos_SLy4 a non-copiable class)

◆ ~Eos_SLy4()

Lorene::Eos_SLy4::~Eos_SLy4 ( )

virtual

Destructor.

Definition at line 89 of file eos_sly4.C.

Member Function Documentation

◆ calcule() [1/2]

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

protectedinherited

General computational method for Cmp 's.

Parameters

thermo	[input] thermodynamical quantity (for instance the enthalpy field)from which the thermodynamical quantity `resu` is to be computed.
nzet	[input] number of domains where `resu` is to be computed.
l_min	[input] index of the innermost domain is which `resu` is to be computed [default value: 0]; `resu` is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
fait	[input] pointer on the member function of class `Eos` which performs the pointwise calculation.
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().

◆ calcule() [2/2]

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

protectedinherited

General computational method for Scalar 's.

Parameters

thermo	[input] thermodynamical quantity (for instance the enthalpy field)from which the thermodynamical quantity `resu` is to be computed.
nzet	[input] number of domains where `resu` is to be computed.
l_min	[input] index of the innermost domain is which `resu` is to be computed [default value: 0]; `resu` is computed only in domains whose indices are in [l_min,l_min+nzet-1] . In the other domains, it is set to zero.
fait	[input] pointer on the member function of class `Eos` which performs the pointwise calculation.
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().

◆ csound_square_ent()

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

inherited

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

Parameters

ent	[input, unit: c^2] enthalpy
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

Returns: [unit: c^2]

Definition at line 499 of file eos.C.

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

◆ csound_square_ent_p()

double Lorene::Eos_tabul::csound_square_ent_p	(	double	ent,
		const Param *
	)		const

virtualinherited

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

Parameters

ent	[input, unit: c^2] enthalpy
par	possible extra parameters of the EOS

Returns: [unit: c^2]

Implements Lorene::Eos.

Reimplemented in Lorene::Eos_consistent.

Definition at line 602 of file eos_tabul.C.

References Lorene::Tbl::get_taille(), Lorene::Eos_tabul::hmax, Lorene::Eos_tabul::hmin, Lorene::log10(), Lorene::Eos_tabul::log_cs2, Lorene::Eos_tabul::logh, Lorene::Eos_tabul::lognb, and Lorene::pow().

◆ der_ener_ent() [1/2]

Cmp Lorene::Eos::der_ener_ent	(	const Cmp &	ent,
		int	nzet,
		int	l_min = `0`,
		Param *	par = `0x0`
	)		const

inherited

Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: dln(e)/dln(H)

Definition at line 454 of file eos.C.

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

◆ der_ener_ent() [2/2]

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

inherited

Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: dln(e)/dln(H)

Definition at line 464 of file eos.C.

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

◆ der_ener_ent_p()

double Lorene::Eos_tabul::der_ener_ent_p	(	double	ent,
		const Param *	par = `0x0`
	)		const

virtualinherited

Computes the logarithmic derivative $d\ln e/d\ln H$ from the log-enthalpy.

Parameters

ent	[input, unit: ] log-enthalpy H

Returns: dln(e)/dln(H)

Implements Lorene::Eos.

Definition at line 526 of file eos_tabul.C.

References Lorene::Eos_tabul::der_nbar_ent_p(), Lorene::Eos_tabul::ener_ent_p(), Lorene::Eos_tabul::hmax, Lorene::Eos_tabul::hmin, and Lorene::Eos_tabul::press_ent_p().

◆ der_nbar_ent() [1/2]

Cmp Lorene::Eos::der_nbar_ent	(	const Cmp &	ent,
		int	nzet,
		int	l_min = `0`,
		Param *	par = `0x0`
	)		const

inherited

Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: dln(n)/dln(H)

Definition at line 431 of file eos.C.

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

◆ der_nbar_ent() [2/2]

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

inherited

Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: dln(n)/dln(H)

Definition at line 441 of file eos.C.

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

◆ der_nbar_ent_p()

double Lorene::Eos_tabul::der_nbar_ent_p	(	double	ent,
		const Param *	par = `0x0`
	)		const

virtualinherited

Computes the logarithmic derivative $d\ln n/d\ln H$ from the log-enthalpy.

Parameters

ent	[input, unit: ] log-enthalpy H

Returns: dln(n)/dln(H)

Implements Lorene::Eos.

Definition at line 503 of file eos_tabul.C.

References Lorene::Eos_tabul::der_press_ent_p(), Lorene::Eos_tabul::der_press_nbar_p(), Lorene::Eos_tabul::hmax, and Lorene::Eos_tabul::hmin.

◆ der_press_ent() [1/2]

Cmp Lorene::Eos::der_press_ent	(	const Cmp &	ent,
		int	nzet,
		int	l_min = `0`,
		Param *	par = `0x0`
	)		const

inherited

Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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.

Returns: dln(p)/dln(H)

Definition at line 476 of file eos.C.

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

◆ der_press_ent() [2/2]

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

inherited

Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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.

Returns: dln(p)/dln(H)

Definition at line 486 of file eos.C.

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

◆ der_press_ent_p()

double Lorene::Eos_tabul::der_press_ent_p	(	double	ent,
		const Param *	par = `0x0`
	)		const

virtualinherited

Computes the logarithmic derivative $d\ln p/d\ln H$ from the log-enthalpy.

Parameters

ent	[input, unit: ] log-enthalpy H

Returns: dln(p)/dln(H)

Implements Lorene::Eos.

Definition at line 546 of file eos_tabul.C.

References Lorene::Eos_tabul::dlpsdlh, Lorene::Tbl::get_taille(), Lorene::Eos_tabul::hmax, Lorene::Eos_tabul::hmin, Lorene::log10(), Lorene::Eos_tabul::logh, and Lorene::Eos_tabul::logp.

◆ der_press_nbar_p()

double Lorene::Eos_tabul::der_press_nbar_p	(	double	ent,
		const Param *	par = `0x0`
	)		const

virtualinherited

Computes the logarithmic derivative $d\ln p/d\ln n$ from the log-enthalpy.

Parameters

ent	[input, unit: ] log-enthalpy H

Returns: dln(p)/dln(n)

Definition at line 576 of file eos_tabul.C.

References Lorene::Eos_tabul::dlpsdlnb, Lorene::Tbl::get_taille(), Lorene::Eos_tabul::hmax, Lorene::Eos_tabul::hmin, Lorene::log10(), and Lorene::Eos_tabul::logh.

◆ ener_ent() [1/2]

Cmp Lorene::Eos::ener_ent	(	const Cmp &	ent,
		int	nzet,
		int	l_min = `0`,
		Param *	par = `0x0`
	)		const

inherited

Computes the total energy density from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: energy density [unit: $\rho_{\rm nuc} c^2$ ], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Definition at line 387 of file eos.C.

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

◆ ener_ent() [2/2]

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

inherited

Computes the total energy density from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: energy density [unit: $\rho_{\rm nuc} c^2$ ], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Definition at line 397 of file eos.C.

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

◆ ener_ent_p()

double Lorene::Eos_tabul::ener_ent_p	(	double	ent,
		const Param *	par = `0x0`
	)		const

virtualinherited

Computes the total energy density from the log-enthalpy.

Parameters

ent	[input, unit: ] log-enthalpy H

Returns: energy density e [unit: $\rho_{\rm nuc} c^2$ ], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Implements Lorene::Eos.

Reimplemented in Lorene::Eos_consistent.

Definition at line 450 of file eos_tabul.C.

References Lorene::Eos_tabul::dlpsdlh, Lorene::Tbl::get_taille(), Lorene::Eos_tabul::hmax, Lorene::Eos_tabul::hmin, Lorene::log10(), Lorene::Eos_tabul::logh, Lorene::Eos_tabul::logp, and Lorene::pow().

◆ eos_from_file() [1/2]

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

staticinherited

Construction of an EOS from a binary file.

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

Definition at line 193 of file eos_from_file.C.

References Lorene::fread_be().

◆ eos_from_file() [2/2]

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

staticinherited

Construction of an EOS from a formatted file.

The fist line of the file must start by the EOS number, according to the following conventions:

1 = relativistic polytropic EOS (class Eos_poly ).
2 = Newtonian polytropic EOS (class Eos_poly_newt ).
3 = Relativistic incompressible EOS (class Eos_incomp ).
4 = Newtonian incompressible EOS (class Eos_incomp_newt ).
5 = Strange matter (MIT Bag model)
6 = Strange matter (MIT Bag model) with crust
10 = SLy4 (Douchin & Haensel 2001)
- 11 = FPS (Friedman-Pandharipande + Skyrme)
12 = BPAL12 (Bombaci et al. 1995)
13 = AkmalPR (Akmal, Pandharipande & Ravenhall 1998)
14 = BBB2 (Baldo, Bombaci & Burgio 1997)
15 = BalbN1H1 (Balberg 2000)
- 16 = GlendNH3 (Glendenning 1985, case 3)
- 17 = Tabulated EOS in CompOSE format
18 = magnetized (tabulated) equation of state
19 = relativistic ideal Fermi gas at zero temperature (class Eos_Fermi)
20 = Tabulated EOS in CompOSE format corrected for thermo. consistency
100 = Multi-domain EOS (class MEos )
110 = Multi-polytropic EOS (class Eos_multi_poly )
120 = Fitted SLy4 (Shibata 2004)
121 = Fitted FPS (Shibata 2004)
122 = Fitted AkmalPR (Taniguchi 2005)

The second line in the file should contain a name given by the user to the EOS. The following lines should contain the EOS parameters (one parameter per line), in the same order than in the class declaration.

Definition at line 345 of file eos_from_file.C.

◆ get_hmax()

double Lorene::Eos_tabul::get_hmax ( ) const

inlineinherited

Returns maximal value of log-enthalpy in the table.

Definition at line 302 of file eos_tabul.h.

References Lorene::Eos_tabul::hmax.

◆ get_hmin()

double Lorene::Eos_tabul::get_hmin ( ) const

inlineinherited

Returns minimal value of log-enthalpy in the table.

Definition at line 299 of file eos_tabul.h.

References Lorene::Eos_tabul::hmin.

◆ get_name()

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

inherited

Returns the EOS name.

Definition at line 179 of file eos.C.

References Lorene::Eos::name.

◆ identify()

int Lorene::Eos_SLy4::identify ( ) const

virtual

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

Implements Lorene::Eos.

Definition at line 153 of file eos_from_file.C.

◆ nbar_ent() [1/2]

Cmp Lorene::Eos::nbar_ent	(	const Cmp &	ent,
		int	nzet,
		int	l_min = `0`,
		Param *	par = `0x0`
	)		const

inherited

Computes the baryon density field from the log-enthalpy field and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: baryon density [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$ ]

Definition at line 362 of file eos.C.

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

◆ nbar_ent() [2/2]

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

inherited

Computes the baryon density field from the log-enthalpy field and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: baryon density [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$ ]

Definition at line 372 of file eos.C.

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

◆ nbar_ent_p()

double Lorene::Eos_tabul::nbar_ent_p	(	double	ent,
		const Param *	par = `0x0`
	)		const

virtualinherited

Computes the baryon density from the log-enthalpy.

Parameters

ent	[input, unit: ] log-enthalpy H

Returns: baryon density n [unit: $n_{\rm nuc} := 0.1 \ {\rm fm}^{-3}$ ]

Implements Lorene::Eos.

Reimplemented in Lorene::Eos_consistent.

Definition at line 422 of file eos_tabul.C.

References Lorene::Eos_tabul::dlpsdlh, Lorene::exp(), Lorene::Tbl::get_taille(), Lorene::Eos_tabul::hmax, Lorene::Eos_tabul::hmin, Lorene::log10(), Lorene::Eos_tabul::logh, Lorene::Eos_tabul::logp, and Lorene::pow().

◆ operator!=()

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

virtual

Comparison operator (difference)

Implements Lorene::Eos.

Definition at line 114 of file eos_sly4.C.

References operator==().

◆ operator==()

bool Lorene::Eos_SLy4::operator== ( const Eos & eos_i ) const

virtual

Comparison operator (egality)

Implements Lorene::Eos.

Definition at line 101 of file eos_sly4.C.

References Lorene::Eos::identify(), and identify().

◆ operator>>()

ostream & Lorene::Eos_SLy4::operator>> ( ostream & ost ) const

protectedvirtual

Operator >>

Implements Lorene::Eos.

Definition at line 125 of file eos_sly4.C.

◆ press_ent() [1/2]

Cmp Lorene::Eos::press_ent	(	const Cmp &	ent,
		int	nzet,
		int	l_min = `0`,
		Param *	par = `0x0`
	)		const

inherited

Computes the pressure from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: pressure [unit: $\rho_{\rm nuc} c^2$ ], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Definition at line 409 of file eos.C.

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

◆ press_ent() [2/2]

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

inherited

Computes the pressure from the log-enthalpy and extra parameters.

Parameters

ent	[input, unit: ] log-enthalpy H defined by $H = c^2 \ln\left( {e+p \over m_B c^2 n} \right)$ , where e is the (total) energy density, p the pressure, n the baryon density, and 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

Returns: pressure [unit: $\rho_{\rm nuc} c^2$ ], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Definition at line 419 of file eos.C.

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

◆ press_ent_p()

double Lorene::Eos_tabul::press_ent_p	(	double	ent,
		const Param *	par = `0x0`
	)		const

virtualinherited

Computes the pressure from the log-enthalpy.

Parameters

ent	[input, unit: ] log-enthalpy H

Returns: pressure p [unit: $\rho_{\rm nuc} c^2$ ], where $\rho_{\rm nuc} := 1.66\ 10^{17} \ {\rm kg/m}^3$

Implements Lorene::Eos.

Reimplemented in Lorene::Eos_consistent.

Definition at line 478 of file eos_tabul.C.

References Lorene::Eos_tabul::dlpsdlh, Lorene::Tbl::get_taille(), Lorene::Eos_tabul::hmax, Lorene::Eos_tabul::hmin, Lorene::log10(), Lorene::Eos_tabul::logh, Lorene::Eos_tabul::logp, and Lorene::pow().

◆ read_table()

void Lorene::Eos_tabul::read_table ( )

protectedvirtualinherited

Reads the file containing the table and initializes in the arrays logh , logp and dlpsdlh .

Reimplemented in Lorene::Eos_consistent.

Definition at line 261 of file eos_tabul.C.

References Lorene::Eos_tabul::authors, Lorene::compute_derivative(), Lorene::Eos_tabul::dlpsdlh, Lorene::Eos_tabul::dlpsdlnb, Lorene::Eos_tabul::hmax, Lorene::Eos_tabul::hmin, Lorene::log(), Lorene::log10(), Lorene::Eos_tabul::log_cs2, Lorene::Eos_tabul::logh, Lorene::Eos_tabul::lognb, Lorene::Eos_tabul::logp, Lorene::pow(), Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), and Lorene::Eos_tabul::tablename.

◆ sauve()

void Lorene::Eos_tabul::sauve ( FILE * fich ) const

virtualinherited

Save in a file.

Reimplemented from Lorene::Eos.

Reimplemented in Lorene::Eos_CompOSE.

Definition at line 248 of file eos_tabul.C.

References Lorene::Eos::sauve(), and Lorene::Eos_tabul::tablename.

◆ set_name()

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

inherited

Sets the EOS name.

Definition at line 173 of file eos.C.

References Lorene::Eos::name.

Friends And Related Function Documentation

◆ Eos::eos_from_file

Eos* Eos::eos_from_file ( FILE * )

friend

The construction functions from a file.

Member Data Documentation

◆ authors

string Lorene::Eos_tabul::authors

protectedinherited

Authors - reference for the table.

Definition at line 194 of file eos_tabul.h.

◆ dlpsdlh

Tbl* Lorene::Eos_tabul::dlpsdlh

protectedinherited

Table of $d\log P/d\log H$ .

Definition at line 209 of file eos_tabul.h.

◆ dlpsdlnb

Tbl* Lorene::Eos_tabul::dlpsdlnb

protectedinherited

Table of $d\log P/d\log nb$ .

Definition at line 215 of file eos_tabul.h.

◆ hmax

double Lorene::Eos_tabul::hmax

protectedinherited

Upper boundary of the enthalpy interval.

Definition at line 200 of file eos_tabul.h.

◆ hmin

double Lorene::Eos_tabul::hmin

protectedinherited

Lower boundary of the enthalpy interval.

Definition at line 197 of file eos_tabul.h.

◆ log_cs2

Tbl* Lorene::Eos_tabul::log_cs2

protectedinherited

Table of $\log c_s^2 = \log \left( c^2 \frac{d p}{d e} \right)$ .

Definition at line 218 of file eos_tabul.h.

◆ logh

Tbl* Lorene::Eos_tabul::logh

protectedinherited

Table of $\log H$ .

Definition at line 203 of file eos_tabul.h.

◆ lognb

Tbl* Lorene::Eos_tabul::lognb

protectedinherited

Table of $\log n_b$ .

Definition at line 212 of file eos_tabul.h.

◆ logp

Tbl* Lorene::Eos_tabul::logp

protectedinherited

Table of $\log p$ .

Definition at line 206 of file eos_tabul.h.

◆ name

char Lorene::Eos::name[100]

protectedinherited

EOS name.

Definition at line 215 of file eos.h.

◆ tablename

string Lorene::Eos_tabul::tablename

protectedinherited

Name of the file containing the tabulated data.

Definition at line 192 of file eos_tabul.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Friends

Detailed Description

Constructor & Destructor Documentation

◆ Eos_SLy4() [1/4]

◆ Eos_SLy4() [2/4]

◆ Eos_SLy4() [3/4]

◆ Eos_SLy4() [4/4]

◆ ~Eos_SLy4()

Member Function Documentation

◆ calcule() [1/2]

◆ calcule() [2/2]

◆ csound_square_ent()

◆ csound_square_ent_p()

◆ der_ener_ent() [1/2]

◆ der_ener_ent() [2/2]

◆ der_ener_ent_p()

◆ der_nbar_ent() [1/2]

◆ der_nbar_ent() [2/2]

◆ der_nbar_ent_p()

◆ der_press_ent() [1/2]

◆ der_press_ent() [2/2]

◆ der_press_ent_p()

◆ der_press_nbar_p()

◆ ener_ent() [1/2]

◆ ener_ent() [2/2]

◆ ener_ent_p()

◆ eos_from_file() [1/2]

◆ eos_from_file() [2/2]

◆ get_hmax()

◆ get_hmin()

◆ get_name()

◆ identify()

◆ nbar_ent() [1/2]

◆ nbar_ent() [2/2]

◆ nbar_ent_p()

◆ operator!=()

◆ operator==()

◆ operator>>()

◆ press_ent() [1/2]

◆ press_ent() [2/2]

◆ press_ent_p()

◆ read_table()

◆ sauve()

◆ set_name()

Friends And Related Function Documentation

◆ Eos::eos_from_file

Member Data Documentation

◆ authors

◆ dlpsdlh

◆ dlpsdlnb

◆ hmax

◆ hmin

◆ log_cs2

◆ logh

◆ lognb

◆ logp

◆ name

◆ tablename