Class for hot differentially rotating stars in Conformal Flatness Condition
and maximal slicing. More...

#include <hot_star_rot_diff_cfc.h>

Inheritance diagram for Lorene::Hot_star_rot_diff_CFC:

Public Member Functions
	Hot_star_rot_diff_CFC (Map &mp_i, int nzet_i, const Hot_eos &heos_i, int relat_i, int e_type, double const_value_i, int filter, double(omega_frot_i)(double, const Tbl &), double(primfrot_i)(double, const Tbl &), const Tbl &par_frot_i)
	Standard constructor. More...

	Hot_star_rot_diff_CFC (const Hot_star_rot_diff_CFC &)
	Copy constructor. More...

	Hot_star_rot_diff_CFC (Map &mp_i, const Hot_eos &heos_i, FILE fich, double(omega_frot_i)(double, const Tbl &), double(*primfrot_i)(double, const Tbl &))
	Constructor from a file (see `sauve(FILE*)` ). More...

virtual	~Hot_star_rot_diff_CFC ()
	Destructor. More...

void	operator= (const Hot_star_rot_diff_CFC &)
	Assignment to another `Hot_star_rot_diff_CFC`. More...

const Scalar &	get_omega_field () const
	Returns the angular velocity field $\Omega$ . More...

virtual double	get_omega_c () const
	Returns the central value of the rotation angular velocity ([f_unit] ) More...

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

virtual double	tsw () const
	Ratio T/W. More...

virtual void	hydro_euler ()
	Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid frame. More...

void	fait_omega_field (double omeg_min, double omeg_max, double precis, int nitermax)
	Computes $\Omega(r,\theta)$ (member `omega_field` ). More...

void	fait_prim_field ()
	Computes the member `prim_field` from `omega_field` . More...

double	omega_funct (double F) const
	Evaluates $\Omega(F)$ , where F is linked to the components of the fluid 4-velocity by $F = u^t u_\varphi \ .$ . More...

double	prim_funct_omega (double F) const
	Evaluates the primitive of $F\frac{d\Omega(F)}{dF}$ . More...

virtual void	equilibrium (double ent_c, double omega0, double fact_omega, int nzadapt, const Tbl &ent_limit, const Itbl &icontrol, const Tbl &control, double mbar_wanted, double aexp_mass, Tbl &diff, Param *=0x0)
	Computes an equilibrium configuration. More...

int	get_relat () const
	Returns the relativity parameter. More...

bool	is_relativistic () const
	Checks whether the star is computed using a relativistic theory. More...

int	spectral_filter_order () const
	Returns the filtering order. More...

double	get_omega () const
	Returns the rotation angular velocity $\Omega$ . More...

const Scalar &	get_psi4 () const
	Returns the conformal factor $\Psi^4$ . More...

const Scalar &	get_psi2 () const
	Returns $\Psi^2$ . More...

const Scalar &	get_psi () const
	Returns $\Psi$ . More...

const Vector &	get_j_euler () const
	Returns the momentum density 3-vector with respect to the Eulerian observer. More...

const Scalar &	get_v2 () const
	Returns $\gamma_{ij}v^i v^j$ . More...

const Sym_tensor	get_hatA () const
	Returns $\hat{A}^{ij}$ . More...

const Scalar	get_hatA_quad () const
	Returns $\tilde{A}_{ij} A^{ij}$ . More...

virtual double	mass_b () const
	Baryonic mass. More...

virtual double	mass_g () const
	Gravitational mass. More...

virtual double	angu_mom () const
	Angular momentum. More...

virtual double	grv2 () const
	Error on the virial identity GRV2. More...

virtual double	grv3 () const
	Error on the virial identity GRV3. More...

virtual double	aplat () const
	Flattening r_pole/r_eq. More...

virtual double	r_circ () const
	Circumferential equatorial radius. More...

virtual double	rp_circ () const
	Circumferential polar radius. More...

virtual double	ellipt () const
	Ellipticity e. More...

void	update_metric ()
	Computes metric quantities from known potentials. More...

void	equilibrium (double ent_c, double omega0, double fact_omega, int nzadapt, const Tbl &ent_limit, const Itbl &icontrol, const Tbl &control, double mbar_wanted, double aexp_mass, Tbl &diff)
	Computes an equilibrium configuration. More...

void	solve_logn_f (Scalar &ln_f_new) const
	Solution of the `‘matter’' part of the Poisson equation for the lapse for rotating stars in CFC. More...

void	solve_logn_q (Scalar &ln_q_new) const
	Solution of the quadratic part of the Poisson equation for the lapse for rotating stars in CFC. More...

void	solve_psi (Scalar &psi_new)
	Solution of the equations for the conformal factor for rotating stars in CFC. More...

void	solve_shift (Vector &shift_new) const
	Solution of the shift equation for rotating stars in CFC. More...

Map &	set_mp ()
	Read/write of the mapping. More...

void	set_enthalpy (const Scalar &)
	Assignment of the enthalpy field. More...

void	equation_of_state ()
	Computes the proper baryon and energy density, as well as pressure from the enthalpy. More...

virtual void	equilibrium_spher (double ent_c, double precis=1.e-14, const Tbl *pent_limit=0x0)
	Computes a spherical static configuration. More...

const Map &	get_mp () const
	Returns the mapping. More...

int	get_nzet () const
	Returns the number of domains occupied by the star. More...

const Hot_eos &	get_eos () const
	Returns the equation of state. More...

const Scalar &	get_ent () const
	Returns the enthalpy field. More...

const Scalar &	get_entropy () const
	Returns the entropy field. More...

const Scalar &	get_nbar () const
	Returns the proper baryon density. More...

const Scalar &	get_ener () const
	Returns the proper total energy density. More...

const Scalar &	get_press () const
	Returns the fluid pressure. More...

const Scalar &	get_temp () const
	Returns the fluid temperature. More...

const Scalar &	get_ener_euler () const
	Returns the total energy density with respect to the Eulerian observer. More...

const Scalar &	get_s_euler () const
	Returns the trace of the stress tensor in the Eulerian frame. More...

const Scalar &	get_gam_euler () const
	Returns the Lorentz factor between the fluid and Eulerian observers. More...

const Vector &	get_u_euler () const
	Returns the fluid 3-velocity with respect to the Eulerian observer. More...

const Tensor &	get_stress_euler () const
	Returns the spatial part of the stress-energy tensor with respect to the Eulerian observer. More...

const Scalar &	get_logn () const
	Returns the logarithm of the lapse N. More...

const Scalar &	get_nn () const
	Returns the lapse function N. More...

const Vector &	get_beta () const
	Returns the shift vector $\beta^i$ . More...

const Scalar &	get_lnq () const

const Metric &	get_gamma () const
	Returns the 3-metric $\gamma$ . More...

double	ray_eq () const
	Coordinate radius at $\phi=0$ , $\theta=\pi/2$ [r_unit]. More...

double	ray_eq_pis2 () const
	Coordinate radius at $\phi=\pi/2$ , $\theta=\pi/2$ [r_unit]. More...

double	ray_eq_pi () const
	Coordinate radius at $\phi=\pi$ , $\theta=\pi/2$ [r_unit]. More...

double	ray_eq_3pis2 () const
	Coordinate radius at $\phi=3\pi/2$ , $\theta=\pi/2$ [r_unit]. More...

double	ray_pole () const
	Coordinate radius at $\theta=0$ [r_unit]. More...

virtual const Itbl &	l_surf () const
	Description of the stellar surface: returns a 2-D `Itbl` containing the values of the domain index l on the surface at the collocation points in $(\theta', \phi')$ . More...

const Tbl &	xi_surf () const
	Description of the stellar surface: returns a 2-D `Tbl` containing the values of the radial coordinate $\xi$ on the surface at the collocation points in $(\theta', \phi')$ . More...

Protected Member Functions
virtual ostream &	operator>> (ostream &) const
	Operator >> (virtual function called by the operator <<). More...

virtual void	partial_display (ostream &) const
	Printing of some informations, excluding all global quantities. More...

virtual void	del_deriv () const
	Deletes all the derived quantities. More...

void	set_der_0x0 () const
	Sets to `0x0` all the pointers on derived quantities. More...

virtual void	del_hydro_euler ()
	Sets to `ETATNONDEF` (undefined state) the hydrodynamical quantities relative to the Eulerian observer. More...

void	update_entropy ()
	Updates the $\ell=0$ part of the entropy from `entropy_param`. More...

Protected Attributes
double(*	omega_frot )(double, const Tbl &)
	Function $\Omega(F)$ defining the rotation profile. More...

double(*	primfrot )(double, const Tbl &)
	Primitive of the function $F\frac{d\Omega(F)}{dF}$ , which vanishes at the stellar center. More...

Tbl	par_frot
	Parameters of the function $\Omega(F)$ . More...

Scalar	omega_field
	Field $\Omega(r,\theta)$ . More...

double	omega_min
	Minimum value of $\Omega$ . More...

double	omega_max
	Maximum value of $\Omega$ . More...

Scalar	prim_field
	Field $\int_0^F F'\,\frac{d\Omega(F')}{dF'} \, dF'$ . More...

int	relat_type
	Relativistic flag. More...

int	e_type
	Integer defining the type of entropy profile. More...

double	const_value
	Constant value of entropy or temperature. More...

int	spectral_filter
	Spectral exponential filtering order. More...

double	omega
	Rotation angular velocity ([f_unit] ) More...

Scalar	psi
	Conformal factor $\Psi$ . More...

Vector	j_euler
	Momentum density 3-vector with respect to the Eulerian observer. More...

Scalar	v2
	$\gamma_{ij}v^i v^j$ More...

Scalar	psi4
	Conformal factor $\Psi^4$ . More...

Scalar	psi2
	$\Psi^2$ More...

const Metric_flat &	flat
	flat metric $f_{ij}$ (spherical components) More...

Sym_tensor	hatA
	$\hat{A}^{ij}$ More...

Scalar	hatA_quad
	$f_{il}\, f_{jm}\, \hat{A}_{ij} \hat{A}^{lm}$ More...

double *	p_angu_mom
	Angular momentum. More...

double *	p_grv2
	Error on the virial identity GRV2. More...

double *	p_grv3
	Error on the virial identity GRV3. More...

double *	p_tsw
	Ratio T/W. More...

double *	p_r_circ
	Circumferential equatorial radius. More...

double *	p_rp_circ
	Circumferential polar radius. More...

Map &	mp
	Mapping associated with the star. More...

int	nzet
	Number of domains of `*mp` occupied by the star. More...

const Hot_eos &	heos
	Equation of state of the stellar matter. More...

Scalar	ent
	Log-enthalpy. More...

Scalar	nbar
	Baryon density in the fluid frame. More...

Scalar	ener
	Total energy density in the fluid frame. More...

Scalar	press
	Fluid pressure. More...

Scalar	temp
	Temperature field (in MeV) More...

Scalar	entropy
	Entropy per baryon field (in $k_B$) More...

Scalar	ener_euler
	Total energy density in the Eulerian frame. More...

Scalar	s_euler
	Trace of the stress scalar in the Eulerian frame. More...

Scalar	gam_euler
	Lorentz factor between the fluid and Eulerian observers. More...

Vector	u_euler
	Fluid 3-velocity with respect to the Eulerian observer. More...

Sym_tensor	stress_euler
	Spatial part of the stress-energy tensor with respect to the Eulerian observer. More...

Scalar	logn
	Logarithm of the lapse N . More...

Scalar	nn
	Lapse function N . More...

Vector	beta
	Shift vector. More...

Scalar	lnq

Metric	gamma
	3-metric More...

double *	p_ray_eq
	Coordinate radius at $\phi=0$ , $\theta=\pi/2$ . More...

double *	p_ray_eq_pis2
	Coordinate radius at $\phi=\pi/2$ , $\theta=\pi/2$ . More...

double *	p_ray_eq_pi
	Coordinate radius at $\phi=\pi$ , $\theta=\pi/2$ . More...

double *	p_ray_eq_3pis2
	Coordinate radius at $\phi=3\pi/2$ , $\theta=\pi/2$ . More...

double *	p_ray_pole
	Coordinate radius at $\theta=0$ . More...

Itbl *	p_l_surf
	Description of the stellar surface: 2-D `Itbl` containing the values of the domain index l on the surface at the collocation points in $(\theta', \phi')$ . More...

Tbl *	p_xi_surf
	Description of the stellar surface: 2-D `Tbl` containing the values of the radial coordinate $\xi$ on the surface at the collocation points in $(\theta', \phi')$ . More...

double *	p_mass_b
	Baryon mass. More...

double *	p_mass_g
	Gravitational mass. More...

Friends
ostream &	operator<< (ostream &, const Hot_star &)
	Display. More...

Detailed Description

Class for hot differentially rotating stars in Conformal Flatness Condition
and maximal slicing.

()

Definition at line 39 of file hot_star_rot_diff_cfc.h.

Constructor & Destructor Documentation

◆ Hot_star_rot_diff_CFC() [1/3]

Lorene::Hot_star_rot_diff_CFC::Hot_star_rot_diff_CFC	(	Map &	mp_i,
		int	nzet_i,
		const Hot_eos &	heos_i,
		int	relat_i,
		int	e_type,
		double	const_value_i,
		int	filter,
		double(*)(double, const Tbl &)	omega_frot_i,
		double(*)(double, const Tbl &)	primfrot_i,
		const Tbl &	par_frot_i
	)

Standard constructor.

Parameters

mp_i	Mapping on which the star will be defined
nzet_i	Number of domains occupied by the star
heos_i	Hot equation of state of the stellar matter
relat	should be `true` for a relativistic star, `false` for a Newtonian one
e_type_i	Type of entropy profile
const_value_i	Value of constant entropy/temperature
omega_frot_i	Function $\Omega(F)$ defining the rotation profile.
primfrot_i	Primitive of $F\frac{d\Omega(F)}{dF}$ which vanishes at the stellar center
par_frot_i	Parameters of functions `omega_frot_i` and `primfrot_i` , `par_frot_i(0)` being the central value of $\Omega$ .

Definition at line 43 of file hot_star_rot_diff_cfc.C.

References Lorene::Hot_star_rot_CFC::omega, omega_field, omega_max, omega_min, and prim_field.

◆ Hot_star_rot_diff_CFC() [2/3]

Lorene::Hot_star_rot_diff_CFC::Hot_star_rot_diff_CFC ( const Hot_star_rot_diff_CFC & star )

Copy constructor.

Definition at line 69 of file hot_star_rot_diff_cfc.C.

◆ Hot_star_rot_diff_CFC() [3/3]

Lorene::Hot_star_rot_diff_CFC::Hot_star_rot_diff_CFC	(	Map &	mp_i,
		const Hot_eos &	heos_i,
		FILE *	fich,
		double(*)(double, const Tbl &)	omega_frot_i,
		double(*)(double, const Tbl &)	primfrot_i
	)

Constructor from a file (see sauve(FILE*) ).

Parameters

mp_i	Mapping on which the star will be defined
heos_i	Hot equation of state of the stellar matter
fich	input file (must have been created by the function `sauve` )
omega_frot_i	Function $\Omega(F)$ defining the rotation profile.
primfrot_i	Primitive of $F\frac{d\Omega(F)}{dF}$ which vanishes at the stellar center

Definition at line 82 of file hot_star_rot_diff_cfc.C.

References Lorene::contract(), fait_prim_field(), Lorene::Hot_star_rot_CFC::flat, Lorene::fread_be(), Lorene::Map::get_mg(), Lorene::Hot_star_rot_CFC::hatA, Lorene::Hot_star_rot_CFC::hatA_quad, hydro_euler(), Lorene::Hot_star_rot_CFC::j_euler, Lorene::Hot_star::mp, omega_field, omega_max, omega_min, Lorene::Vector::ope_killing_conf(), Lorene::Vector::poisson(), Lorene::Hot_star_rot_CFC::psi2, Lorene::Hot_star_rot_CFC::psi4, and Lorene::Tensor::up_down().

◆ ~Hot_star_rot_diff_CFC()

Lorene::Hot_star_rot_diff_CFC::~Hot_star_rot_diff_CFC ( )

virtual

Destructor.

Definition at line 119 of file hot_star_rot_diff_cfc.C.

Member Function Documentation

◆ angu_mom()

double Lorene::Hot_star_rot_CFC::angu_mom ( ) const

virtualinherited

Angular momentum.

Definition at line 92 of file hot_strot_cfc_global.C.

References Lorene::contract(), Lorene::Metric::cov(), Lorene::Hot_star::gamma, Lorene::Map::get_bvect_spher(), Lorene::Scalar::integrale(), Lorene::Hot_star_rot_CFC::j_euler, Lorene::Hot_star::mp, Lorene::Hot_star_rot_CFC::p_angu_mom, Lorene::Hot_star_rot_CFC::psi2, Lorene::Hot_star_rot_CFC::psi4, Lorene::Vector::set(), Lorene::Scalar::set_etat_zero(), and Lorene::Scalar::std_spectral_base().

◆ aplat()

double Lorene::Hot_star_rot_CFC::aplat ( ) const

virtualinherited

Flattening r_pole/r_eq.

Definition at line 369 of file hot_strot_cfc_global.C.

References Lorene::Hot_star::ray_eq(), and Lorene::Hot_star::ray_pole().

◆ del_deriv()

void Lorene::Hot_star_rot_CFC::del_deriv ( ) const

protectedvirtualinherited

Deletes all the derived quantities.

Reimplemented from Lorene::Hot_star.

Definition at line 201 of file hot_star_rot_cfc.C.

References Lorene::Hot_star::del_deriv(), Lorene::Hot_star_rot_CFC::p_angu_mom, Lorene::Hot_star_rot_CFC::p_grv2, Lorene::Hot_star_rot_CFC::p_grv3, Lorene::Hot_star_rot_CFC::p_r_circ, Lorene::Hot_star_rot_CFC::p_rp_circ, Lorene::Hot_star_rot_CFC::p_tsw, and Lorene::Hot_star_rot_CFC::set_der_0x0().

◆ del_hydro_euler()

void Lorene::Hot_star_rot_CFC::del_hydro_euler ( )

protectedvirtualinherited

Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer.

Reimplemented from Lorene::Hot_star.

Definition at line 229 of file hot_star_rot_cfc.C.

References Lorene::Hot_star_rot_CFC::del_deriv(), Lorene::Hot_star::del_hydro_euler(), Lorene::Hot_star_rot_CFC::j_euler, Lorene::Scalar::set_etat_nondef(), Lorene::Tensor::set_etat_nondef(), and Lorene::Hot_star_rot_CFC::v2.

◆ ellipt()

double Lorene::Hot_star_rot_CFC::ellipt ( ) const

virtualinherited

Ellipticity e.

Defined as $e = \sqrt{1 - \left( \frac{R^c_e}{R^c_p} \right)^2}$ , where $R^c_e$ and $R^c_p$ are, respectively, the equatorial and polar circumferential radius, given by r_circ() and rp_circ().

Definition at line 379 of file hot_strot_cfc_global.C.

References Lorene::Hot_star_rot_CFC::r_circ(), Lorene::Hot_star_rot_CFC::rp_circ(), and Lorene::sqrt().

◆ equation_of_state()

void Lorene::Hot_star::equation_of_state ( )

inherited

Computes the proper baryon and energy density, as well as pressure from the enthalpy.

Definition at line 489 of file hot_star.C.

References Lorene::Scalar::allocate_all(), Lorene::Hot_star::del_deriv(), Lorene::Hot_star::ener, Lorene::Hot_eos::ener_Hs(), Lorene::Hot_star::ent, Lorene::Hot_star::entropy, Lorene::Mg3d::get_grille3d(), Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Hot_star::heos, Lorene::Hot_star::mp, Lorene::Hot_star::nbar, Lorene::Hot_eos::nbar_Hs(), Lorene::Hot_star::nzet, Lorene::Hot_star::press, Lorene::Hot_eos::press_Hs(), Lorene::Mtbl::set(), Lorene::Scalar::set_domain(), Lorene::Mtbl::set_etat_qcq(), Lorene::Tbl::set_etat_qcq(), Lorene::Scalar::set_etat_qcq(), Lorene::Scalar::std_spectral_base(), Lorene::Mtbl::t, Lorene::Hot_star::temp, Lorene::Hot_eos::temp_Hs(), and Lorene::Grille3d::x.

◆ equilibrium() [1/2]

void Lorene::Hot_star_rot_CFC::equilibrium	(	double	ent_c,
		double	omega0,
		double	fact_omega,
		int	nzadapt,
		const Tbl &	ent_limit,
		const Itbl &	icontrol,
		const Tbl &	control,
		double	mbar_wanted,
		double	aexp_mass,
		Tbl &	diff
	)

inherited

Computes an equilibrium configuration.

Definition at line 36 of file hot_strot_cfc_equil.C.

References Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), Lorene::Mg3d::get_type_t(), and Lorene::Hot_star::mp.

◆ equilibrium() [2/2]

void Lorene::Hot_star_rot_diff_CFC::equilibrium	(	double	ent_c,
		double	omega0,
		double	fact_omega,
		int	nzadapt,
		const Tbl &	ent_limit,
		const Itbl &	icontrol,
		const Tbl &	control,
		double	mbar_wanted,
		double	aexp_mass,
		Tbl &	diff,
		Param *	= `0x0`
	)

virtual

Computes an equilibrium configuration.

Parameters

ent_c	[input] Central enthalpy
omega0	[input] Requested central angular velocity (if `fact_omega=1`. )
fact_omega	[input] 1.01 = search for the Keplerian frequency, = otherwise.
nzadapt	[input] Number of (inner) domains where the mapping adaptation to an iso-enthalpy surface should be performed
ent_limit	[input] 1-D `Tbl` of dimension `nzet` which defines the enthalpy at the outer boundary of each domain
icontrol	[input] Set of integer parameters (stored as a 1-D `Itbl` of size 9) to control the iteration: `icontrol(0)` = mer_max : maximum number of steps `icontrol(1)` = mer_rot : step at which the rotation is switched on `icontrol(2)` = mer_change_omega : step at which the rotation velocity is changed to reach the final one `icontrol(3)` = mer_fix_omega : step at which the final rotation velocity must have been reached `icontrol(4)` = mer_mass : the absolute value of `mer_mass` is the step from which the baryon mass is forced to converge, by varying the central enthalpy (`mer_mass` > 0 ) or the angular velocity (`mer_mass` < 0 ) `icontrol(5)` = mermax_poisson : maximum number of steps in `Map_et::poisson` `icontrol(6)` = mer_triax : step at which the 3-D perturbation is switched on `icontrol(7)` = delta_mer_kep : number of steps after `mer_fix_omega` when `omega` starts to be increased by `fact_omega` to search for the Keplerian velocity `icontrol(8)` = mer_diff : step at which differential rotation is switched on
control	[input] Set of parameters (stored as a 1-D `Tbl` of size 7) to control the iteration: `control(0)` = precis : threshold on the enthalpy relative change for ending the computation `control(1)` = omega_ini : initial central angular velocity, switched on only if `mer_rot` < 0 , otherwise 0 is used `control(2)` = relax : relaxation factor in the main iteration `control(3)` = relax_poisson : relaxation factor in `Map_et::poisson` `control(4)` = thres_adapt : threshold on dH/dr for freezing the adaptation of the mapping `control(5)` = ampli_triax : relative amplitude of the 3-D perturbation `control(6)` = precis_adapt : precision for `Map_et::adapt`
mbar_wanted	[input] Requested baryon mass (effective only if `mer_mass` > mer_max )
aexp_mass	[input] Exponent for the increase factor of the central enthalpy to converge to the requested baryon mass
diff	[output] 1-D `Tbl` of size 7 for the storage of some error indicators : `diff(0)` : Relative change in the enthalpy field between two successive steps `diff(1)` : Relative error in the resolution of the Poisson equation for `nuf` `diff(2)` : Relative error in the resolution of the Poisson equation for `nuq` `diff(3)` : Relative error in the resolution of the Poisson equation for `dzeta` `diff(4)` : Relative error in the resolution of the Poisson equation for `tggg` `diff(5)` : Relative error in the resolution of the equation for `shift` (x comp.) `diff(6)` : Relative error in the resolution of the equation for `shift` (y comp.)

Definition at line 40 of file hot_strot_diff_cfc_equil.C.

References Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), Lorene::Mg3d::get_type_t(), and Lorene::Hot_star::mp.

◆ equilibrium_spher()

void Lorene::Hot_star::equilibrium_spher	(	double	ent_c,
		double	precis = `1.e-14`,
		const Tbl *	pent_limit = `0x0`
	)

virtualinherited

Computes a spherical static configuration.

Parameters

ent_c	[input] central value of the enthalpy
precis	[input] threshold in the relative difference between the enthalpy fields of two consecutive steps to stop the iterative procedure (default value: 1.e-14)
ent_limit	[input] : array of enthalpy values to be set at the boundaries between the domains; if set to 0x0 (default), the initial values will be kept.

Definition at line 104 of file hot_star_equil_spher.C.

References Lorene::Map_et::adapt(), Lorene::Param::add_double(), Lorene::Param::add_int(), Lorene::Param::add_int_mod(), Lorene::Param::add_tbl(), Lorene::Scalar::annule(), Lorene::diffrel(), Lorene::Scalar::dsdr(), Lorene::Map_af::dsdr(), Lorene::Hot_star::ener, Lorene::Hot_star::ener_euler, Lorene::Hot_star::ent, Lorene::Hot_star::equation_of_state(), Lorene::exp(), Lorene::Hot_star::gam_euler, Lorene::Hot_star::gamma, Lorene::Map_af::get_alpha(), Lorene::Map_et::get_alpha(), Lorene::Map_af::get_beta(), Lorene::Map_et::get_beta(), Lorene::Map::get_bvect_cart(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Map_af::homothetie(), Lorene::Scalar::integrale(), Lorene::Hot_star::logn, Lorene::Hot_star::mass_b(), Lorene::Hot_star::mass_g(), Lorene::Hot_star::mp, Lorene::Hot_star::nn, Lorene::norme(), Lorene::Hot_star::nzet, Lorene::Map_af::poisson(), Lorene::Hot_star::press, Lorene::Hot_star::s_euler, Lorene::Vector::set(), Lorene::Map_af::set_alpha(), Lorene::Map_af::set_beta(), Lorene::Scalar::set_dzpuis(), Lorene::Cmp::set_etat_qcq(), Lorene::Tensor::set_etat_zero(), Lorene::sqrt(), Lorene::Scalar::std_spectral_base(), Lorene::Hot_star::u_euler, Lorene::Scalar::val_grid_point(), and Lorene::Map::val_r().

◆ fait_omega_field()

void Lorene::Hot_star_rot_diff_CFC::fait_omega_field	(	double	omeg_min,
		double	omeg_max,
		double	precis,
		int	nitermax
	)

Computes $\Omega(r,\theta)$ (member omega_field ).

The computation amounts to solving the equation

$\Omega - \Omega\left({B^2 r^2\sin^2\theta (\Omega - N^\varphi) \over N^2 - B^2 r^2 \sin^2(\Omega-N^\varphi)^2}\right) = 0$

for $\Omega$ .

Parameters

omeg_min	[input] Lower bound of the interval for searching omega
omeg_max	[input] Higher bound of the interval for searching omega
precis	[input] Required precision in the determination of the zero by the secant method
nitermax	[input] Maximum number of iterations in the secant method to compute the zero.

Definition at line 46 of file hot_strot_diff_cfc_faitomeg.C.

References Lorene::Param::add_int(), Lorene::Param::add_scalar(), Lorene::Tensor::annule_domain(), Lorene::Hot_star::beta, Lorene::Scalar::get_etat(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), Lorene::Hot_star::mp, Lorene::Scalar::mult_rsint(), Lorene::Hot_star::nn, omega_field, Lorene::Hot_star_rot_CFC::psi4, and Lorene::Scalar::std_spectral_base().

◆ fait_prim_field()

void Lorene::Hot_star_rot_diff_CFC::fait_prim_field ( )

Computes the member prim_field from omega_field .

Definition at line 176 of file hot_strot_diff_cfc_faitomeg.C.

References Lorene::Scalar::allocate_all(), Lorene::Tensor::annule_domain(), Lorene::Hot_star::beta, Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Hot_star::mp, Lorene::Scalar::mult_rsint(), Lorene::Hot_star::nn, Lorene::Hot_star::nzet, omega_field, par_frot, prim_field, primfrot, Lorene::Hot_star_rot_CFC::psi4, Lorene::Tbl::set(), Lorene::Scalar::set_domain(), Lorene::Scalar::std_spectral_base(), and Lorene::Scalar::val_grid_point().

◆ get_beta()

const Vector& Lorene::Hot_star::get_beta ( ) const

inlineinherited

Returns the shift vector $\beta^i$ .

Definition at line 313 of file hot_star.h.

References Lorene::Hot_star::beta.

◆ get_ener()

const Scalar& Lorene::Hot_star::get_ener ( ) const

inlineinherited

Returns the proper total energy density.

Definition at line 278 of file hot_star.h.

References Lorene::Hot_star::ener.

◆ get_ener_euler()

const Scalar& Lorene::Hot_star::get_ener_euler ( ) const

inlineinherited

Returns the total energy density with respect to the Eulerian observer.

Definition at line 287 of file hot_star.h.

References Lorene::Hot_star::ener_euler.

◆ get_ent()

const Scalar& Lorene::Hot_star::get_ent ( ) const

inlineinherited

Returns the enthalpy field.

Definition at line 269 of file hot_star.h.

References Lorene::Hot_star::ent.

◆ get_entropy()

const Scalar& Lorene::Hot_star::get_entropy ( ) const

inlineinherited

Returns the entropy field.

Definition at line 272 of file hot_star.h.

References Lorene::Hot_star::ent.

◆ get_eos()

const Hot_eos& Lorene::Hot_star::get_eos ( ) const

inlineinherited

Returns the equation of state.

Definition at line 266 of file hot_star.h.

References Lorene::Hot_star::heos.

◆ get_gam_euler()

const Scalar& Lorene::Hot_star::get_gam_euler ( ) const

inlineinherited

Returns the Lorentz factor between the fluid and Eulerian observers.

Definition at line 293 of file hot_star.h.

References Lorene::Hot_star::gam_euler.

◆ get_gamma()

const Metric& Lorene::Hot_star::get_gamma ( ) const

inlineinherited

Returns the 3-metric $\gamma$ .

Definition at line 320 of file hot_star.h.

References Lorene::Hot_star::gamma.

◆ get_hatA()

const Sym_tensor Lorene::Hot_star_rot_CFC::get_hatA ( ) const

inlineinherited

Returns $\hat{A}^{ij}$ .

Definition at line 226 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::hatA.

◆ get_hatA_quad()

const Scalar Lorene::Hot_star_rot_CFC::get_hatA_quad ( ) const

inlineinherited

Returns $\tilde{A}_{ij} A^{ij}$ .

Definition at line 231 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::hatA_quad.

◆ get_j_euler()

const Vector& Lorene::Hot_star_rot_CFC::get_j_euler ( ) const

inlineinherited

Returns the momentum density 3-vector with respect to the Eulerian observer.

Definition at line 213 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::j_euler.

◆ get_logn()

const Scalar& Lorene::Hot_star::get_logn ( ) const

inlineinherited

Returns the logarithm of the lapse N.

In the Newtonian case, this is the Newtonian gravitational potential (in units of $c^2$ ).

Definition at line 307 of file hot_star.h.

References Lorene::Hot_star::logn.

◆ get_mp()

const Map& Lorene::Hot_star::get_mp ( ) const

inlineinherited

Returns the mapping.

Definition at line 260 of file hot_star.h.

References Lorene::Hot_star::mp.

◆ get_nbar()

const Scalar& Lorene::Hot_star::get_nbar ( ) const

inlineinherited

Returns the proper baryon density.

Definition at line 275 of file hot_star.h.

References Lorene::Hot_star::nbar.

◆ get_nn()

const Scalar& Lorene::Hot_star::get_nn ( ) const

inlineinherited

Returns the lapse function N.

Definition at line 310 of file hot_star.h.

References Lorene::Hot_star::nn.

◆ get_nzet()

int Lorene::Hot_star::get_nzet ( ) const

inlineinherited

Returns the number of domains occupied by the star.

Definition at line 263 of file hot_star.h.

References Lorene::Hot_star::nzet.

◆ get_omega()

double Lorene::Hot_star_rot_CFC::get_omega ( ) const

inlineinherited

Returns the rotation angular velocity $\Omega$ .

Definition at line 188 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::omega.

◆ get_omega_c()

double Lorene::Hot_star_rot_diff_CFC::get_omega_c ( ) const

virtual

Returns the central value of the rotation angular velocity ([f_unit] )

Definition at line 255 of file hot_star_rot_diff_cfc.C.

References omega_field, and Lorene::Scalar::val_grid_point().

◆ get_omega_field()

const Scalar& Lorene::Hot_star_rot_diff_CFC::get_omega_field ( ) const

inline

Returns the angular velocity field $\Omega$ .

Definition at line 141 of file hot_star_rot_diff_cfc.h.

References omega_field.

◆ get_press()

const Scalar& Lorene::Hot_star::get_press ( ) const

inlineinherited

Returns the fluid pressure.

Definition at line 281 of file hot_star.h.

References Lorene::Hot_star::press.

◆ get_psi()

const Scalar& Lorene::Hot_star_rot_CFC::get_psi ( ) const

inlineinherited

Returns $\Psi$ .

Definition at line 204 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::psi.

◆ get_psi2()

const Scalar& Lorene::Hot_star_rot_CFC::get_psi2 ( ) const

inlineinherited

Returns $\Psi^2$ .

Definition at line 199 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::psi2.

◆ get_psi4()

const Scalar& Lorene::Hot_star_rot_CFC::get_psi4 ( ) const

inlineinherited

Returns the conformal factor $\Psi^4$ .

Definition at line 194 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::psi4.

◆ get_relat()

int Lorene::Hot_star_rot_CFC::get_relat ( ) const

inlineinherited

Returns the relativity parameter.

Definition at line 177 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::relat_type.

◆ get_s_euler()

const Scalar& Lorene::Hot_star::get_s_euler ( ) const

inlineinherited

Returns the trace of the stress tensor in the Eulerian frame.

Definition at line 290 of file hot_star.h.

References Lorene::Hot_star::s_euler.

◆ get_stress_euler()

const Tensor& Lorene::Hot_star::get_stress_euler ( ) const

inlineinherited

Returns the spatial part of the stress-energy tensor with respect to the Eulerian observer.

Definition at line 301 of file hot_star.h.

References Lorene::Hot_star::stress_euler.

◆ get_temp()

const Scalar& Lorene::Hot_star::get_temp ( ) const

inlineinherited

Returns the fluid temperature.

Definition at line 284 of file hot_star.h.

References Lorene::Hot_star::temp.

◆ get_u_euler()

const Vector& Lorene::Hot_star::get_u_euler ( ) const

inlineinherited

Returns the fluid 3-velocity with respect to the Eulerian observer.

Definition at line 296 of file hot_star.h.

References Lorene::Hot_star::u_euler.

◆ get_v2()

const Scalar& Lorene::Hot_star_rot_CFC::get_v2 ( ) const

inlineinherited

Returns $\gamma_{ij}v^i v^j$ .

Definition at line 218 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::v2.

◆ grv2()

double Lorene::Hot_star_rot_CFC::grv2 ( ) const

virtualinherited

Error on the virial identity GRV2.

Definition at line 153 of file hot_strot_cfc_global.C.

References Lorene::Metric::con(), Lorene::Metric::cov(), Lorene::Scalar::div_r_dzpuis(), Lorene::Scalar::dsdr(), Lorene::Scalar::dsdt(), Lorene::Hot_star::ener_euler, Lorene::Hot_star::gamma, Lorene::Hot_star_rot_CFC::hatA, Lorene::Hot_star::logn, Lorene::Hot_star_rot_CFC::p_grv2, Lorene::Hot_star::press, Lorene::Hot_star_rot_CFC::psi4, Lorene::sqrt(), Lorene::Scalar::std_spectral_base(), and Lorene::Hot_star_rot_CFC::v2.

◆ grv3()

double Lorene::Hot_star_rot_CFC::grv3 ( ) const

virtualinherited

Error on the virial identity GRV3.

Definition at line 223 of file hot_strot_cfc_global.C.

References Lorene::Metric::con(), Lorene::Metric::connect(), Lorene::contract(), Lorene::Scalar::derive_con(), Lorene::Scalar::derive_cov(), Lorene::Hot_star::gamma, Lorene::Connection::get_delta(), Lorene::Hot_star_rot_CFC::hatA_quad, Lorene::Scalar::integrale(), Lorene::Hot_star::logn, Lorene::Hot_star_rot_CFC::p_grv3, Lorene::Hot_star_rot_CFC::psi2, Lorene::Hot_star_rot_CFC::psi4, Lorene::Hot_star::s_euler, and Lorene::Scalar::std_spectral_base().

◆ hydro_euler()

void Lorene::Hot_star_rot_diff_CFC::hydro_euler ( )

virtual

Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid frame.

The calculation is performed starting from the quantities omega_field , ent , ener , press , and a_car ,
which are supposed to be up to date.
From these, the following fields are updated: gam_euler , u_euler , ener_euler , s_euler .

Reimplemented from Lorene::Hot_star_rot_CFC.

Definition at line 38 of file hot_strot_diff_cfc_hydro.C.

References Lorene::Hot_star::beta, Lorene::Metric::con(), Lorene::contract(), Lorene::Metric::cov(), Lorene::Hot_star_rot_CFC::del_deriv(), Lorene::Hot_star::ener, Lorene::Hot_star::ener_euler, Lorene::Vector::exponential_filter_r(), Lorene::Sym_tensor::exponential_filter_r(), Lorene::Scalar::exponential_filter_r(), Lorene::Hot_star::gam_euler, Lorene::Hot_star::gamma, Lorene::Hot_star_rot_CFC::j_euler, Lorene::Scalar::mult_rsint(), Lorene::Hot_star::nn, Lorene::Hot_star::nzet, omega_field, Lorene::Hot_star::press, Lorene::Hot_star::s_euler, Lorene::Vector::set(), Lorene::Scalar::set_etat_zero(), Lorene::Hot_star_rot_CFC::spectral_filter, Lorene::sqrt(), Lorene::Vector::std_spectral_base(), Lorene::Tensor::std_spectral_base(), Lorene::Scalar::std_spectral_base(), Lorene::Hot_star::stress_euler, Lorene::Hot_star::u_euler, and Lorene::Hot_star_rot_CFC::v2.

◆ is_relativistic()

bool Lorene::Hot_star_rot_CFC::is_relativistic ( ) const

inlineinherited

Checks whether the star is computed using a relativistic theory.

Definition at line 180 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::relat_type.

◆ l_surf()

const Itbl & Lorene::Hot_star::l_surf ( ) const

virtualinherited

Description of the stellar surface: returns a 2-D Itbl containing the values of the domain index l on the surface at the collocation points in $(\theta', \phi')$ .

The stellar surface is defined as the location where the enthalpy (member ent) vanishes.

Definition at line 69 of file hot_star_global.C.

References Lorene::Hot_star::ent, Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Scalar::get_spectral_va(), Lorene::Hot_star::mp, Lorene::Hot_star::nzet, Lorene::Hot_star::p_l_surf, and Lorene::Hot_star::p_xi_surf.

◆ mass_b()

double Lorene::Hot_star_rot_CFC::mass_b ( ) const

virtualinherited

Baryonic mass.

Implements Lorene::Hot_star.

Definition at line 44 of file hot_strot_cfc_global.C.

References Lorene::Hot_star::gam_euler, Lorene::Scalar::integrale(), Lorene::Hot_star::nbar, Lorene::Hot_star::p_mass_b, Lorene::Hot_star_rot_CFC::psi2, Lorene::Hot_star_rot_CFC::psi4, and Lorene::Scalar::std_spectral_base().

◆ mass_g()

double Lorene::Hot_star_rot_CFC::mass_g ( ) const

virtualinherited

Gravitational mass.

Implements Lorene::Hot_star.

Definition at line 67 of file hot_strot_cfc_global.C.

References Lorene::Hot_star::beta, Lorene::contract(), Lorene::Metric::cov(), Lorene::Hot_star::ener_euler, Lorene::Hot_star::gamma, Lorene::Scalar::integrale(), Lorene::Hot_star_rot_CFC::j_euler, Lorene::Hot_star::nn, Lorene::Hot_star::p_mass_g, Lorene::Hot_star_rot_CFC::psi2, Lorene::Hot_star_rot_CFC::psi4, Lorene::Hot_star::s_euler, and Lorene::Scalar::std_spectral_base().

◆ omega_funct()

double Lorene::Hot_star_rot_diff_CFC::omega_funct ( double F ) const

Evaluates $\Omega(F)$ , where F is linked to the components of the fluid 4-velocity by

$F = u^t u_\varphi \ .$

.

omega_funct calls frot with the parameters par_frot .

Parameters

F	[input] value of

Returns: value of $\Omega(F))$

Definition at line 243 of file hot_star_rot_diff_cfc.C.

References omega_frot, and par_frot.

◆ operator=()

void Lorene::Hot_star_rot_diff_CFC::operator= ( const Hot_star_rot_diff_CFC & star )

Assignment to another Hot_star_rot_diff_CFC.

Definition at line 128 of file hot_star_rot_diff_cfc.C.

References omega_field, omega_frot, omega_max, omega_min, Lorene::Hot_star_rot_CFC::operator=(), par_frot, prim_field, and primfrot.

◆ operator>>()

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

protectedvirtual

Operator >> (virtual function called by the operator <<).

Reimplemented from Lorene::Hot_star_rot_CFC.

Definition at line 170 of file hot_star_rot_diff_cfc.C.

References Lorene::Hot_star_rot_CFC::angu_mom(), Lorene::Hot_star_rot_CFC::aplat(), Lorene::Hot_star::ener, Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), get_omega_c(), Lorene::Hot_star_rot_CFC::mass_b(), Lorene::Hot_star_rot_CFC::mass_g(), Lorene::max(), Lorene::Hot_star::mp, Lorene::Hot_star::nbar, Lorene::Hot_star::nzet, omega_field, omega_max, omega_min, Lorene::Hot_star_rot_CFC::operator>>(), par_frot, Lorene::pow(), Lorene::Hot_star::press, Lorene::Hot_star::ray_eq(), Lorene::sqrt(), tsw(), and Lorene::Scalar::val_grid_point().

◆ partial_display()

void Lorene::Hot_star_rot_diff_CFC::partial_display ( ostream & ost ) const

protectedvirtual

Printing of some informations, excluding all global quantities.

Definition at line 261 of file hot_star_rot_diff_cfc.C.

References Lorene::Hot_star_rot_CFC::aplat(), Lorene::Hot_star::beta, Lorene::Hot_star::ener, Lorene::Hot_star::ent, Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), get_omega_c(), Lorene::Hot_star::logn, Lorene::Hot_star::mp, Lorene::Hot_star::nbar, Lorene::Hot_star::nn, Lorene::Hot_star::nzet, Lorene::Hot_star::press, Lorene::Hot_star_rot_CFC::psi, Lorene::Hot_star::ray_eq(), Lorene::Hot_star::u_euler, and Lorene::Scalar::val_grid_point().

◆ prim_funct_omega()

double Lorene::Hot_star_rot_diff_CFC::prim_funct_omega ( double F ) const

Evaluates the primitive of $F\frac{d\Omega(F)}{dF}$ .

Parameters

F	[input] value of

Returns: value of $\int_0^F F'\,\frac{d\Omega(F')}{dF'} \, dF'$

Definition at line 249 of file hot_star_rot_diff_cfc.C.

References par_frot, and primfrot.

◆ r_circ()

double Lorene::Hot_star_rot_CFC::r_circ ( ) const

virtualinherited

Circumferential equatorial radius.

Definition at line 282 of file hot_strot_cfc_global.C.

References Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_type_t(), Lorene::Hot_star::mp, Lorene::Hot_star::nzet, and Lorene::Hot_star_rot_CFC::p_r_circ.

◆ ray_eq()

double Lorene::Hot_star::ray_eq ( ) const

inherited

Coordinate radius at $\phi=0$ , $\theta=\pi/2$ [r_unit].

Definition at line 114 of file hot_star_global.C.

References Lorene::Map::get_mg(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), Lorene::Hot_star::mp, and Lorene::Hot_star::p_ray_eq.

◆ ray_eq_3pis2()

double Lorene::Hot_star::ray_eq_3pis2 ( ) const

inherited

Coordinate radius at $\phi=3\pi/2$ , $\theta=\pi/2$ [r_unit].

Definition at line 239 of file hot_star_global.C.

References Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), Lorene::Hot_star::mp, and Lorene::Hot_star::p_ray_eq_3pis2.

◆ ray_eq_pi()

double Lorene::Hot_star::ray_eq_pi ( ) const

inherited

Coordinate radius at $\phi=\pi$ , $\theta=\pi/2$ [r_unit].

Definition at line 192 of file hot_star_global.C.

References Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), Lorene::Hot_star::mp, and Lorene::Hot_star::p_ray_eq_pi.

◆ ray_eq_pis2()

double Lorene::Hot_star::ray_eq_pis2 ( ) const

inherited

Coordinate radius at $\phi=\pi/2$ , $\theta=\pi/2$ [r_unit].

Definition at line 144 of file hot_star_global.C.

References Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), Lorene::Hot_star::mp, and Lorene::Hot_star::p_ray_eq_pis2.

◆ ray_pole()

double Lorene::Hot_star::ray_pole ( ) const

inherited

Coordinate radius at $\theta=0$ [r_unit].

Definition at line 284 of file hot_star_global.C.

References Lorene::Map::get_mg(), Lorene::Mg3d::get_type_t(), Lorene::Hot_star::mp, and Lorene::Hot_star::p_ray_pole.

◆ rp_circ()

double Lorene::Hot_star_rot_CFC::rp_circ ( ) const

virtualinherited

◆ sauve()

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

virtual

Save in a file.

Reimplemented from Lorene::Hot_star_rot_CFC.

Definition at line 153 of file hot_star_rot_diff_cfc.C.

References Lorene::fwrite_be(), omega_field, omega_max, omega_min, par_frot, Lorene::Hot_star_rot_CFC::sauve(), Lorene::Tbl::sauve(), and Lorene::Scalar::sauve().

◆ set_der_0x0()

void Lorene::Hot_star_rot_CFC::set_der_0x0 ( ) const

protectedvirtualinherited

Sets to 0x0 all the pointers on derived quantities.

Reimplemented from Lorene::Hot_star.

Definition at line 217 of file hot_star_rot_cfc.C.

References Lorene::Hot_star_rot_CFC::p_angu_mom, Lorene::Hot_star_rot_CFC::p_grv2, Lorene::Hot_star_rot_CFC::p_grv3, Lorene::Hot_star_rot_CFC::p_r_circ, Lorene::Hot_star_rot_CFC::p_rp_circ, and Lorene::Hot_star_rot_CFC::p_tsw.

◆ set_enthalpy()

void Lorene::Hot_star::set_enthalpy ( const Scalar & ent_i )

inherited

Assignment of the enthalpy field.

Definition at line 397 of file hot_star.C.

References Lorene::Hot_star::del_deriv(), Lorene::Hot_star::ent, and Lorene::Hot_star::equation_of_state().

◆ set_mp()

Map& Lorene::Hot_star::set_mp ( )

inlineinherited

Read/write of the mapping.

Definition at line 227 of file hot_star.h.

References Lorene::Hot_star::mp.

◆ solve_logn_f()

void Lorene::Hot_star_rot_CFC::solve_logn_f ( Scalar & ln_f_new ) const

inherited

Solution of the `‘matter’' part of the Poisson equation for the lapse for rotating stars in CFC.

Definition at line 39 of file hot_strot_cfc_metric.C.

References Lorene::Hot_star::ener_euler, Lorene::Scalar::poisson(), Lorene::Hot_star_rot_CFC::psi4, and Lorene::Hot_star::s_euler.

◆ solve_logn_q()

void Lorene::Hot_star_rot_CFC::solve_logn_q ( Scalar & ln_q_new ) const

inherited

Solution of the quadratic part of the Poisson equation for the lapse for rotating stars in CFC.

Definition at line 53 of file hot_strot_cfc_metric.C.

References Lorene::contract(), Lorene::Scalar::derive_con(), Lorene::Scalar::derive_cov(), Lorene::Hot_star_rot_CFC::flat, Lorene::Hot_star_rot_CFC::hatA_quad, Lorene::log(), Lorene::Hot_star::logn, Lorene::Scalar::poisson(), Lorene::Hot_star_rot_CFC::psi, Lorene::Hot_star_rot_CFC::psi4, Lorene::Scalar::std_spectral_base(), and Lorene::Tensor::up_down().

◆ solve_psi()

void Lorene::Hot_star_rot_CFC::solve_psi ( Scalar & psi_new )

inherited

Solution of the equations for the conformal factor for rotating stars in CFC.

Definition at line 72 of file hot_strot_cfc_metric.C.

References Lorene::contract(), Lorene::Hot_star::ener_euler, Lorene::Hot_star_rot_CFC::flat, Lorene::Scalar::get_etat(), Lorene::Hot_star_rot_CFC::hatA, Lorene::Hot_star_rot_CFC::hatA_quad, Lorene::Hot_star_rot_CFC::j_euler, Lorene::Vector::ope_killing_conf(), Lorene::Vector::poisson(), Lorene::Scalar::poisson(), Lorene::Hot_star_rot_CFC::psi, Lorene::Hot_star_rot_CFC::psi2, Lorene::Hot_star_rot_CFC::psi4, Lorene::Scalar::std_spectral_base(), and Lorene::Tensor::up_down().

◆ solve_shift()

void Lorene::Hot_star_rot_CFC::solve_shift ( Vector & shift_new ) const

inherited

Solution of the shift equation for rotating stars in CFC.

Definition at line 106 of file hot_strot_cfc_metric.C.

References Lorene::Sym_tensor::divergence(), Lorene::Hot_star_rot_CFC::flat, Lorene::Hot_star_rot_CFC::hatA, Lorene::Tensor::inc_dzpuis(), Lorene::Hot_star::nn, Lorene::Vector::poisson(), Lorene::Hot_star_rot_CFC::psi2, Lorene::Hot_star_rot_CFC::psi4, and Lorene::Vector::set().

◆ spectral_filter_order()

int Lorene::Hot_star_rot_CFC::spectral_filter_order ( ) const

inlineinherited

Returns the filtering order.

Definition at line 183 of file hot_star_rot_cfc.h.

References Lorene::Hot_star_rot_CFC::spectral_filter.

◆ tsw()

double Lorene::Hot_star_rot_diff_CFC::tsw ( ) const

virtual

◆ update_entropy()

void Lorene::Hot_star_rot_CFC::update_entropy ( )

protectedinherited

Updates the $\ell=0$ part of the entropy from entropy_param.

Definition at line 291 of file hot_star_rot_cfc.C.

References Lorene::Scalar::annule(), Lorene::Hot_star_rot_CFC::const_value, Lorene::Hot_star_rot_CFC::e_type, Lorene::Hot_star::entropy, Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), Lorene::Hot_star::mp, Lorene::Hot_star::nzet, and Lorene::Scalar::std_spectral_base().

◆ update_metric()

void Lorene::Hot_star_rot_CFC::update_metric ( )

inherited

Computes metric quantities from known potentials.

The calculation is performed starting from psi, logn, hatA, which are supposed to be up to date. From these, the following fields are updated: nnn, psi4, psi2,shift, and hatA_quad.

Definition at line 125 of file hot_strot_cfc_metric.C.

References Lorene::Metric_flat::cov(), Lorene::Hot_star_rot_CFC::del_deriv(), Lorene::exp(), Lorene::Hot_star_rot_CFC::flat, Lorene::Hot_star::gamma, Lorene::Hot_star::logn, Lorene::Hot_star::nn, Lorene::Hot_star_rot_CFC::psi, Lorene::Hot_star_rot_CFC::psi2, Lorene::Hot_star_rot_CFC::psi4, and Lorene::Scalar::std_spectral_base().

◆ xi_surf()

const Tbl & Lorene::Hot_star::xi_surf ( ) const

inherited

Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate $\xi$ on the surface at the collocation points in $(\theta', \phi')$ .

The stellar surface is defined as the location where the enthalpy (member ent) vanishes.

Definition at line 95 of file hot_star_global.C.

References Lorene::Hot_star::l_surf(), Lorene::Hot_star::p_l_surf, and Lorene::Hot_star::p_xi_surf.

Friends And Related Function Documentation

◆ operator<<

ostream& operator<<	(	ostream &	ost,
		const Hot_star &	et
	)

friend

Display.

Definition at line 431 of file hot_star.C.

Member Data Documentation

◆ beta

Vector Lorene::Hot_star::beta

protectedinherited

Shift vector.

Definition at line 133 of file hot_star.h.

◆ const_value

double Lorene::Hot_star_rot_CFC::const_value

protectedinherited

Constant value of entropy or temperature.

Definition at line 59 of file hot_star_rot_cfc.h.

◆ e_type

int Lorene::Hot_star_rot_CFC::e_type

protectedinherited

Integer defining the type of entropy profile.

0: constant entropy. 1: constant temperature

Definition at line 56 of file hot_star_rot_cfc.h.

◆ ener

Scalar Lorene::Hot_star::ener

protectedinherited

Total energy density in the fluid frame.

Definition at line 94 of file hot_star.h.

◆ ener_euler

Scalar Lorene::Hot_star::ener_euler

protectedinherited

Total energy density in the Eulerian frame.

Definition at line 103 of file hot_star.h.

◆ ent

Scalar Lorene::Hot_star::ent

protectedinherited

Log-enthalpy.

Definition at line 91 of file hot_star.h.

◆ entropy

Scalar Lorene::Hot_star::entropy

protectedinherited

Entropy per baryon field (in $k_B$)

Definition at line 99 of file hot_star.h.

◆ flat

const Metric_flat& Lorene::Hot_star_rot_CFC::flat

protectedinherited

flat metric $f_{ij}$ (spherical components)

Definition at line 94 of file hot_star_rot_cfc.h.

◆ gam_euler

Scalar Lorene::Hot_star::gam_euler

protectedinherited

Lorentz factor between the fluid and Eulerian observers.

Definition at line 109 of file hot_star.h.

◆ gamma

Metric Lorene::Hot_star::gamma

protectedinherited

3-metric

Definition at line 140 of file hot_star.h.

◆ hatA

Sym_tensor Lorene::Hot_star_rot_CFC::hatA

protectedinherited

$\hat{A}^{ij}$

Definition at line 96 of file hot_star_rot_cfc.h.

◆ hatA_quad

Scalar Lorene::Hot_star_rot_CFC::hatA_quad

protectedinherited

$f_{il}\, f_{jm}\, \hat{A}_{ij} \hat{A}^{lm}$

Definition at line 97 of file hot_star_rot_cfc.h.

◆ heos

const Hot_eos& Lorene::Hot_star::heos

protectedinherited

Equation of state of the stellar matter.

Definition at line 86 of file hot_star.h.

◆ j_euler

Vector Lorene::Hot_star_rot_CFC::j_euler

protectedinherited

Momentum density 3-vector with respect to the Eulerian observer.

Definition at line 84 of file hot_star_rot_cfc.h.

◆ logn

Scalar Lorene::Hot_star::logn

protectedinherited

Logarithm of the lapse N .

In the Newtonian case, this is the Newtonian gravitational potential (in units of $c^2$ ).

Definition at line 127 of file hot_star.h.

◆ mp

Map& Lorene::Hot_star::mp

protectedinherited

Mapping associated with the star.

Definition at line 81 of file hot_star.h.

◆ nbar

Scalar Lorene::Hot_star::nbar

protectedinherited

Baryon density in the fluid frame.

Definition at line 93 of file hot_star.h.

◆ nn

Scalar Lorene::Hot_star::nn

protectedinherited

Lapse function N .

Definition at line 130 of file hot_star.h.

◆ nzet

int Lorene::Hot_star::nzet

protectedinherited

Number of domains of *mp occupied by the star.

Definition at line 84 of file hot_star.h.

◆ omega

double Lorene::Hot_star_rot_CFC::omega

protectedinherited

Rotation angular velocity ([f_unit] )

Definition at line 69 of file hot_star_rot_cfc.h.

◆ omega_field

Scalar Lorene::Hot_star_rot_diff_CFC::omega_field

protected

Field $\Omega(r,\theta)$ .

Definition at line 73 of file hot_star_rot_diff_cfc.h.

◆ omega_frot

double(* Lorene::Hot_star_rot_diff_CFC::omega_frot) (double, const Tbl &)

protected

Function $\Omega(F)$ defining the rotation profile.

The variable F is linked to the components of the fluid 4-velocity by

$F = u^t u_\varphi \ .$

The first argument of omega_frot must be $F$ ; the second argument contains the parameters, the first of which being necessarily the central value of $\Omega$ .

Definition at line 53 of file hot_star_rot_diff_cfc.h.

◆ omega_max

double Lorene::Hot_star_rot_diff_CFC::omega_max

protected

Maximum value of $\Omega$ .

Definition at line 76 of file hot_star_rot_diff_cfc.h.

◆ omega_min

double Lorene::Hot_star_rot_diff_CFC::omega_min

protected

Minimum value of $\Omega$ .

Definition at line 75 of file hot_star_rot_diff_cfc.h.

◆ p_angu_mom

double* Lorene::Hot_star_rot_CFC::p_angu_mom

mutableprotectedinherited

Angular momentum.

Definition at line 106 of file hot_star_rot_cfc.h.

◆ p_grv2

double* Lorene::Hot_star_rot_CFC::p_grv2

mutableprotectedinherited

Error on the virial identity GRV2.

Definition at line 107 of file hot_star_rot_cfc.h.

◆ p_grv3

double* Lorene::Hot_star_rot_CFC::p_grv3

mutableprotectedinherited

Error on the virial identity GRV3.

Definition at line 108 of file hot_star_rot_cfc.h.

◆ p_l_surf

Itbl* Lorene::Hot_star::p_l_surf

mutableprotectedinherited

Description of the stellar surface: 2-D Itbl containing the values of the domain index l on the surface at the collocation points in $(\theta', \phi')$ .

Definition at line 165 of file hot_star.h.

◆ p_mass_b

double* Lorene::Hot_star::p_mass_b

mutableprotectedinherited

Baryon mass.

Definition at line 173 of file hot_star.h.

◆ p_mass_g

double* Lorene::Hot_star::p_mass_g

mutableprotectedinherited

Gravitational mass.

Definition at line 174 of file hot_star.h.

◆ p_r_circ

double* Lorene::Hot_star_rot_CFC::p_r_circ

mutableprotectedinherited

Circumferential equatorial radius.

Definition at line 110 of file hot_star_rot_cfc.h.

◆ p_ray_eq

double* Lorene::Hot_star::p_ray_eq

mutableprotectedinherited

Coordinate radius at $\phi=0$ , $\theta=\pi/2$ .

Definition at line 147 of file hot_star.h.

◆ p_ray_eq_3pis2

double* Lorene::Hot_star::p_ray_eq_3pis2

mutableprotectedinherited

Coordinate radius at $\phi=3\pi/2$ , $\theta=\pi/2$ .

Definition at line 156 of file hot_star.h.

◆ p_ray_eq_pi

double* Lorene::Hot_star::p_ray_eq_pi

mutableprotectedinherited

Coordinate radius at $\phi=\pi$ , $\theta=\pi/2$ .

Definition at line 153 of file hot_star.h.

◆ p_ray_eq_pis2

double* Lorene::Hot_star::p_ray_eq_pis2

mutableprotectedinherited

Coordinate radius at $\phi=\pi/2$ , $\theta=\pi/2$ .

Definition at line 150 of file hot_star.h.

◆ p_ray_pole

double* Lorene::Hot_star::p_ray_pole

mutableprotectedinherited

Coordinate radius at $\theta=0$ .

Definition at line 159 of file hot_star.h.

◆ p_rp_circ

double* Lorene::Hot_star_rot_CFC::p_rp_circ

mutableprotectedinherited

Circumferential polar radius.

Definition at line 111 of file hot_star_rot_cfc.h.

◆ p_tsw

double* Lorene::Hot_star_rot_CFC::p_tsw

mutableprotectedinherited

Ratio T/W.

Definition at line 109 of file hot_star_rot_cfc.h.

◆ p_xi_surf

Tbl* Lorene::Hot_star::p_xi_surf

mutableprotectedinherited

Description of the stellar surface: 2-D Tbl containing the values of the radial coordinate $\xi$ on the surface at the collocation points in $(\theta', \phi')$ .

Definition at line 171 of file hot_star.h.

◆ par_frot

Tbl Lorene::Hot_star_rot_diff_CFC::par_frot

protected

Parameters of the function $\Omega(F)$ .

To be used as the second argument of functions frot and primfrot . The parameter par_frot(0) must always be the central angular velocity.

Definition at line 70 of file hot_star_rot_diff_cfc.h.

◆ press

Scalar Lorene::Hot_star::press

protectedinherited

Fluid pressure.

Definition at line 95 of file hot_star.h.

◆ prim_field

Scalar Lorene::Hot_star_rot_diff_CFC::prim_field

protected

Field $\int_0^F F'\,\frac{d\Omega(F')}{dF'} \, dF'$ .

Definition at line 79 of file hot_star_rot_diff_cfc.h.

◆ primfrot

double(* Lorene::Hot_star_rot_diff_CFC::primfrot) (double, const Tbl &)

protected

Primitive of the function $F\frac{d\Omega(F)}{dF}$ , which vanishes at the stellar center.

The first argument of primfrot must be $F$ ; the second argument contains the parameters, the first of which being necessarily the central value of $\Omega$ .

Definition at line 61 of file hot_star_rot_diff_cfc.h.

◆ psi

Scalar Lorene::Hot_star_rot_CFC::psi

protectedinherited

Conformal factor $\Psi$ .

Definition at line 74 of file hot_star_rot_cfc.h.

◆ psi2

Scalar Lorene::Hot_star_rot_CFC::psi2

protectedinherited

$\Psi^2$

Definition at line 91 of file hot_star_rot_cfc.h.

◆ psi4

Scalar Lorene::Hot_star_rot_CFC::psi4

protectedinherited

Conformal factor $\Psi^4$ .

Definition at line 90 of file hot_star_rot_cfc.h.

◆ relat_type

int Lorene::Hot_star_rot_CFC::relat_type

protectedinherited

Relativistic flag.

0 : Newtonian theory (not implemented) 3 : CFC with $\hat{A}^{ij}_{TT} = 0$ (see Cordero-Carrion et al. 2009).

Definition at line 50 of file hot_star_rot_cfc.h.

◆ s_euler

Scalar Lorene::Hot_star::s_euler

protectedinherited

Trace of the stress scalar in the Eulerian frame.

Definition at line 106 of file hot_star.h.

◆ spectral_filter

int Lorene::Hot_star_rot_CFC::spectral_filter

protectedinherited

Spectral exponential filtering order.

If 0, no filtering is done (see also Scalar::exponential_filter_r). Filtering is performed only in shells containing matter, i.e. for domain numbers l such that $0 < l <$ nzet .

Definition at line 67 of file hot_star_rot_cfc.h.

◆ stress_euler

Sym_tensor Lorene::Hot_star::stress_euler

protectedinherited

Spatial part of the stress-energy tensor with respect to the Eulerian observer.

Definition at line 117 of file hot_star.h.

◆ temp

Scalar Lorene::Hot_star::temp

protectedinherited

Temperature field (in MeV)

Definition at line 98 of file hot_star.h.

◆ u_euler

Vector Lorene::Hot_star::u_euler

protectedinherited

Fluid 3-velocity with respect to the Eulerian observer.

Definition at line 112 of file hot_star.h.

◆ v2

Scalar Lorene::Hot_star_rot_CFC::v2

protectedinherited

$\gamma_{ij}v^i v^j$

Definition at line 85 of file hot_star_rot_cfc.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ Hot_star_rot_diff_CFC() [1/3]

◆ Hot_star_rot_diff_CFC() [2/3]

◆ Hot_star_rot_diff_CFC() [3/3]

◆ ~Hot_star_rot_diff_CFC()

Member Function Documentation

◆ angu_mom()

◆ aplat()

◆ del_deriv()

◆ del_hydro_euler()

◆ ellipt()

◆ equation_of_state()

◆ equilibrium() [1/2]

◆ equilibrium() [2/2]

◆ equilibrium_spher()

◆ fait_omega_field()

◆ fait_prim_field()

◆ get_beta()

◆ get_ener()

◆ get_ener_euler()

◆ get_ent()

◆ get_entropy()

◆ get_eos()

◆ get_gam_euler()

◆ get_gamma()

◆ get_hatA()

◆ get_hatA_quad()

◆ get_j_euler()

◆ get_logn()

◆ get_mp()

◆ get_nbar()

◆ get_nn()

◆ get_nzet()

◆ get_omega()

◆ get_omega_c()

◆ get_omega_field()

◆ get_press()

◆ get_psi()

◆ get_psi2()

◆ get_psi4()

◆ get_relat()

◆ get_s_euler()

◆ get_stress_euler()

◆ get_temp()

◆ get_u_euler()

◆ get_v2()

◆ grv2()

◆ grv3()

◆ hydro_euler()

◆ is_relativistic()

◆ l_surf()

◆ mass_b()

◆ mass_g()

◆ omega_funct()

◆ operator=()

◆ operator>>()

◆ partial_display()

◆ prim_funct_omega()

◆ r_circ()

◆ ray_eq()

◆ ray_eq_3pis2()

◆ ray_eq_pi()

◆ ray_eq_pis2()

◆ ray_pole()

◆ rp_circ()

◆ sauve()

◆ set_der_0x0()

◆ set_enthalpy()

◆ set_mp()

◆ solve_logn_f()

◆ solve_logn_q()

◆ solve_psi()

◆ solve_shift()

◆ spectral_filter_order()

◆ tsw()