Class for relativistic rotating stars in Dirac gauge and maximal slicing. More...

#include <star_rot_dirac.h>

Inheritance diagram for Lorene::Star_rot_Dirac:

Public Member Functions
	Star_rot_Dirac (Map &mp_i, int nzet_i, const Eos &eos_i, int filter=0)
	Standard constructor.
	Star_rot_Dirac (const Star_rot_Dirac &)
	Copy constructor.
	Star_rot_Dirac (Map &mp_i, const Eos &eos_i, FILE *fich)
	Constructor from a file (see `sauve(FILE*)` ).
virtual	~Star_rot_Dirac ()
	Destructor.
void	operator= (const Star_rot_Dirac &)
	Assignment to another `Star_rot_Dirac`.
int	spectral_filter_order () const
	Returns the filtering order.
double	get_omega () const
	Returns the rotation angular velocity $\Omega$ .
const Scalar &	get_psi4 () const
	Returns the conformal factor $\Psi^4$ .
const Scalar &	get_psi2 () const
	Returns $\Psi^2$ .
const Scalar &	get_qqq () const
	Returns $Q=\Psi^2 N$ .
const Scalar &	get_ln_psi () const
	Returns $ln(\Psi)$ .
const Vector &	get_j_euler () const
	Returns the momentum density 3-vector with respect to the Eulerian observer.
const Scalar &	get_v2 () const
	Reutrns $\gamma_{ij}v^i v^j$ .
const Metric	get_tgamma () const
	Returns the conformal metric $\tilde{\gamma}_{ij}$ .
const Sym_tensor	get_aa () const
	Returns $A^{ij}$ .
const Sym_tensor	get_taa () const
	Returns $\tilde{A}_{ij}$ .
const Scalar	get_aa_quad () const
	Returns $\tilde{A}_{ij} A^{ij}$ .
const Sym_tensor_trans	get_hh () const
	Returns $h^{ij}$ .
virtual void	sauve (FILE *) const
	Save in a file.
virtual double	mass_b () const
	Baryonic mass.
virtual double	mass_g () const
	Gravitational mass.
virtual double	angu_mom () const
	Angular momentum.
virtual double	grv2 () const
	Error on the virial identity GRV2.
virtual double	grv3 () const
	Error on the virial identity GRV3.
virtual double	tsw () const
	Ratio T/W.
virtual double	aplat () const
	Flattening r_pole/r_eq.
virtual double	r_circ () const
	Circumferential equatorial radius.
virtual double	rp_circ () const
	Circumferential polar radius.
virtual double	ellipt () const
	Ellipticity e.
virtual void	hydro_euler ()
	Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid frame.
void	update_metric ()
	Computes metric quantities from known potentials.
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.
void	solve_logn_f (Scalar &ln_f_new) const
	Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
void	solve_logn_q (Scalar &ln_q_new) const
	Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
void	solve_qqq (Scalar &q_new) const
	Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
void	solve_shift (Vector &shift_new) const
	Solution of the shift equation for rotating stars in Dirac gauge.
void	solve_hij (Sym_tensor_trans &hij_new) const
	Solution of the tensor Poisson equation for rotating stars in Dirac gauge.
Map &	set_mp ()
	Read/write of the mapping.
void	set_enthalpy (const Scalar &)
	Assignment of the enthalpy field.
void	equation_of_state ()
	Computes the proper baryon and energy density, as well as pressure from the enthalpy.
virtual void	equilibrium_spher (double ent_c, double precis=1.e-14, const Tbl *pent_limit=0x0)
	Computes a spherical static configuration.
const Map &	get_mp () const
	Returns the mapping.
int	get_nzet () const
	Returns the number of domains occupied by the star.
const Eos &	get_eos () const
	Returns the equation of state.
const Scalar &	get_ent () const
	Returns the enthalpy field.
const Scalar &	get_nbar () const
	Returns the proper baryon density.
const Scalar &	get_ener () const
	Returns the proper total energy density.
const Scalar &	get_press () const
	Returns the fluid pressure.
const Scalar &	get_ener_euler () const
	Returns the total energy density with respect to the Eulerian observer.
const Scalar &	get_s_euler () const
	Returns the trace of the stress tensor in the Eulerian frame.
const Scalar &	get_gam_euler () const
	Returns the Lorentz factor between the fluid and Eulerian observers.
const Vector &	get_u_euler () const
	Returns the fluid 3-velocity with respect to the Eulerian observer.
const Tensor &	get_stress_euler () const
	Returns the spatial part of the stress-energy tensor with respect to the Eulerian observer.
const Scalar &	get_logn () const
	Returns the logarithm of the lapse N.
const Scalar &	get_nn () const
	Returns the lapse function N.
const Vector &	get_beta () const
	Returns the shift vector $\beta^i$ .
const Scalar &	get_lnq () const
const Metric &	get_gamma () const
	Returns the 3-metric $\gamma$ .
double	ray_eq () const
	Coordinate radius at $\phi=0$ , $\theta=\pi/2$ [r_unit].
double	ray_eq_pis2 () const
	Coordinate radius at $\phi=\pi/2$ , $\theta=\pi/2$ [r_unit].
double	ray_eq_pi () const
	Coordinate radius at $\phi=\pi$ , $\theta=\pi/2$ [r_unit].
double	ray_eq_3pis2 () const
	Coordinate radius at $\phi=3\pi/2$ , $\theta=\pi/2$ [r_unit].
double	ray_pole () const
	Coordinate radius at $\theta=0$ [r_unit].
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')$ .
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')$ .

Protected Member Functions
virtual void	del_deriv () const
	Deletes all the derived quantities.
void	set_der_0x0 () const
	Sets to `0x0` all the pointers on derived quantities.
virtual void	del_hydro_euler ()
	Sets to `ETATNONDEF` (undefined state) the hydrodynamical quantities relative to the Eulerian observer.
virtual ostream &	operator>> (ostream &) const
	Operator >> (virtual function called by the operator <<).

Protected Attributes
int	spectral_filter
	Spectral exponential filtering order.
double	omega
	Rotation angular velocity ([f_unit] ).
Scalar	psi4
	Conformal factor $\Psi^4$ .
Scalar	psi2
	$\Psi^2$
Scalar	qqq
	$Q = \Psi^2 N$
Scalar	ln_psi
	$ln(\Psi)$
Vector	j_euler
	Momentum density 3-vector with respect to the Eulerian observer.
Scalar	v2
	$\gamma_{ij}v^i v^j$
const Metric_flat &	flat
	flat metric $f_{ij}$ (spherical components)
Metric	tgamma
	$\tilde{\gamma}_{ij}$
Sym_tensor	aa
	$A^{ij}$
Sym_tensor	taa
	$\tilde{A}_{ij}$
Scalar	aa_quad
	$\tilde{A}_{ij} A^{ij}$
Sym_tensor_trans	hh
	$h^{ij}$ is defined by $\tilde{\gamma}^{ij}=f^{ij}+h^{ij}$ .
double *	p_angu_mom
	Angular momentum.
double *	p_grv2
	Error on the virial identity GRV2.
double *	p_grv3
	Error on the virial identity GRV3.
double *	p_tsw
	Ratio T/W.
double *	p_r_circ
	Circumferential equatorial radius.
double *	p_rp_circ
	Circumferential polar radius.
Map &	mp
	Mapping associated with the star.
int	nzet
	Number of domains of `*mp` occupied by the star.
const Eos &	eos
	Equation of state of the stellar matter.
Scalar	ent
	Log-enthalpy.
Scalar	nbar
	Baryon density in the fluid frame.
Scalar	ener
	Total energy density in the fluid frame.
Scalar	press
	Fluid pressure.
Scalar	ener_euler
	Total energy density in the Eulerian frame.
Scalar	s_euler
	Trace of the stress scalar in the Eulerian frame.
Scalar	gam_euler
	Lorentz factor between the fluid and Eulerian observers.
Vector	u_euler
	Fluid 3-velocity with respect to the Eulerian observer.
Sym_tensor	stress_euler
	Spatial part of the stress-energy tensor with respect to the Eulerian observer.
Scalar	logn
	Logarithm of the lapse N .
Scalar	nn
	Lapse function N .
Vector	beta
	Shift vector.
Scalar	lnq
Metric	gamma
	3-metric
double *	p_ray_eq
	Coordinate radius at $\phi=0$ , $\theta=\pi/2$ .
double *	p_ray_eq_pis2
	Coordinate radius at $\phi=\pi/2$ , $\theta=\pi/2$ .
double *	p_ray_eq_pi
	Coordinate radius at $\phi=\pi$ , $\theta=\pi/2$ .
double *	p_ray_eq_3pis2
	Coordinate radius at $\phi=3\pi/2$ , $\theta=\pi/2$ .
double *	p_ray_pole
	Coordinate radius at $\theta=0$ .
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')$ .
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')$ .
double *	p_mass_b
	Baryon mass.
double *	p_mass_g
	Gravitational mass.

Detailed Description

Class for relativistic rotating stars in Dirac gauge and maximal slicing.

(*** Under development ***)()

Definition at line 46 of file star_rot_dirac.h.

Constructor & Destructor Documentation

◆ Star_rot_Dirac() [1/3]

Lorene::Star_rot_Dirac::Star_rot_Dirac	(	Map &	mp_i,
		int	nzet_i,
		const Eos &	eos_i,
		int	filter = 0 )

Standard constructor.

Parameters

mp_i	Mapping on which the star will be defined
nzet_i	Number of domains occupied by the star
eos_i	Equation of state of the stellar matter
filter	order for spectral exponential filtering

Definition at line 93 of file star_rot_dirac.C.

References aa, aa_quad, flat, Lorene::flat_met_spher(), Lorene::get_bvect_spher(), hh, j_euler, ln_psi, Lorene::Map(), omega, psi2, psi4, qqq, set_der_0x0(), spectral_filter, Lorene::Star::Star(), taa, tgamma, and v2.

◆ Star_rot_Dirac() [2/3]

Lorene::Star_rot_Dirac::Star_rot_Dirac ( const Star_rot_Dirac & star )

Copy constructor.

Definition at line 137 of file star_rot_dirac.C.

References aa, aa_quad, flat, hh, j_euler, ln_psi, omega, psi2, psi4, qqq, set_der_0x0(), spectral_filter, Lorene::Star::Star(), Star_rot_Dirac(), taa, tgamma, and v2.

◆ Star_rot_Dirac() [3/3]

Lorene::Star_rot_Dirac::Star_rot_Dirac	(	Map &	mp_i,
		const Eos &	eos_i,
		FILE *	fich )

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

Parameters

mp_i	Mapping on which the star will be defined
eos_i	Equation of state of the stellar matter
fich	input file (must have been created by the function `sauve` )

Definition at line 164 of file star_rot_dirac.C.

References aa, aa_quad, Lorene::Star::beta, Lorene::Star::equation_of_state(), flat, Lorene::flat_met_spher(), Lorene::fread_be(), Lorene::get_bvect_spher(), Lorene::get_mg(), hh, hydro_euler(), j_euler, ln_psi, Lorene::Map(), omega, psi2, psi4, qqq, set_der_0x0(), spectral_filter, Lorene::Star::Star(), taa, tgamma, update_metric(), and v2.

◆ ~Star_rot_Dirac()

Lorene::Star_rot_Dirac::~Star_rot_Dirac ( )

virtual

Destructor.

Definition at line 207 of file star_rot_dirac.C.

References del_deriv().

Member Function Documentation

◆ angu_mom()

double Lorene::Star_rot_Dirac::angu_mom ( ) const

virtual

Angular momentum.

Definition at line 143 of file strot_dirac_global.C.

References Lorene::contract(), Lorene::Star::gamma, Lorene::Scalar::integrale(), j_euler, Lorene::Star::mp, Lorene::Scalar::mult_rsint(), p_angu_mom, Lorene::Vector::set(), Lorene::Scalar::set_etat_zero(), Lorene::sqrt(), and Lorene::Scalar::std_spectral_base().

◆ aplat()

double Lorene::Star_rot_Dirac::aplat ( ) const

virtual

Flattening r_pole/r_eq.

Definition at line 418 of file strot_dirac_global.C.

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

◆ del_deriv()

void Lorene::Star_rot_Dirac::del_deriv ( ) const

protectedvirtual

Deletes all the derived quantities.

Reimplemented from Lorene::Star.

Definition at line 218 of file star_rot_dirac.C.

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

◆ del_hydro_euler()

void Lorene::Star_rot_Dirac::del_hydro_euler ( )

protectedvirtual

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

Reimplemented from Lorene::Star.

Definition at line 246 of file star_rot_dirac.C.

References del_deriv(), Lorene::Star::del_hydro_euler(), j_euler, and v2.

◆ ellipt()

double Lorene::Star_rot_Dirac::ellipt ( ) const

virtual

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 428 of file strot_dirac_global.C.

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

◆ equation_of_state()

void Lorene::Star::equation_of_state ( )

inherited

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

Definition at line 462 of file star.C.

References Lorene::Param::add_int(), Lorene::Scalar::allocate_all(), del_deriv(), ener, ent, eos, Lorene::Mg3d::get_grille3d(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), mp, nbar, nzet, press, Lorene::Mtbl::set(), Lorene::Scalar::set_domain(), Lorene::Mtbl::set_etat_qcq(), Lorene::Tbl::set_etat_qcq(), Lorene::Scalar::std_spectral_base(), Lorene::Mtbl::t, and Lorene::Grille3d::x.

◆ equilibrium()

void Lorene::Star_rot_Dirac::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.

Definition at line 95 of file strot_dirac_equilibrium.C.

References Lorene::abs(), Lorene::Tbl::annule_hard(), Lorene::Star::beta, Lorene::contract(), Lorene::diffrel(), Lorene::Star::ent, Lorene::Star::equation_of_state(), flat, Lorene::Star::gam_euler, Lorene::Star::gamma, Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Mg3d::get_type_t(), grv2(), hh, hydro_euler(), Lorene::log(), Lorene::log10(), Lorene::Star::logn, mass_b(), mass_g(), Lorene::Star::mp, Lorene::Star::nn, Lorene::Star::nzet, omega, Lorene::pow(), qqq, Lorene::Tbl::set(), Lorene::Tensor::set_etat_zero(), solve_hij(), solve_logn_f(), solve_logn_q(), solve_qqq(), solve_shift(), Lorene::sqrt(), Lorene::Scalar::std_spectral_base(), Lorene::Star::u_euler, update_metric(), v2, and Lorene::Scalar::val_grid_point().

◆ equilibrium_spher()

void Lorene::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 95 of file 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::diffrel(), Lorene::Map_af::dsdr(), Lorene::Scalar::dsdr(), ener, ener_euler, ent, equation_of_state(), Lorene::exp(), gam_euler, gamma, Lorene::Map_af::get_alpha(), Lorene::Map_et::get_alpha(), Lorene::Map_af::get_beta(), Lorene::Map_et::get_beta(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Map_af::homothetie(), Lorene::Scalar::integrale(), logn, mass_b(), mass_g(), mp, nn, Lorene::norme(), nzet, Lorene::Map_af::poisson(), press, s_euler, Lorene::Tensor::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(), u_euler, and Lorene::Scalar::val_grid_point().

◆ get_aa()

const Sym_tensor Lorene::Star_rot_Dirac::get_aa ( ) const

inline

Returns $A^{ij}$ .

Definition at line 229 of file star_rot_dirac.h.

References aa.

◆ get_aa_quad()

const Scalar Lorene::Star_rot_Dirac::get_aa_quad ( ) const

inline

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

Definition at line 239 of file star_rot_dirac.h.

References aa_quad.

◆ get_beta()

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

inlineinherited

Returns the shift vector $\beta^i$ .

Definition at line 402 of file star.h.

References beta.

◆ get_ener()

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

inlineinherited

Returns the proper total energy density.

Definition at line 370 of file star.h.

References ener.

◆ get_ener_euler()

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

inlineinherited

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

Definition at line 376 of file star.h.

References ener_euler.

◆ get_ent()

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

inlineinherited

Returns the enthalpy field.

Definition at line 364 of file star.h.

References ent.

◆ get_eos()

const Eos & Lorene::Star::get_eos ( ) const

inlineinherited

Returns the equation of state.

Definition at line 361 of file star.h.

References eos.

◆ get_gam_euler()

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

inlineinherited

Returns the Lorentz factor between the fluid and Eulerian observers.

Definition at line 382 of file star.h.

References gam_euler.

◆ get_gamma()

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

inlineinherited

Returns the 3-metric $\gamma$ .

Definition at line 409 of file star.h.

References gamma.

◆ get_hh()

const Sym_tensor_trans Lorene::Star_rot_Dirac::get_hh ( ) const

inline

Returns $h^{ij}$ .

Definition at line 244 of file star_rot_dirac.h.

References hh.

◆ get_j_euler()

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

inline

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

Definition at line 209 of file star_rot_dirac.h.

References j_euler.

◆ get_ln_psi()

const Scalar & Lorene::Star_rot_Dirac::get_ln_psi ( ) const

inline

Returns $ln(\Psi)$ .

Definition at line 199 of file star_rot_dirac.h.

References ln_psi.

◆ get_lnq()

const Scalar & Lorene::Star::get_lnq ( ) const

inlineinherited

Definition at line 406 of file star.h.

◆ get_logn()

const Scalar & Lorene::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 396 of file star.h.

References logn.

◆ get_mp()

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

inlineinherited

Returns the mapping.

Definition at line 355 of file star.h.

References Lorene::Map(), and mp.

◆ get_nbar()

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

inlineinherited

Returns the proper baryon density.

Definition at line 367 of file star.h.

References nbar.

◆ get_nn()

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

inlineinherited

Returns the lapse function N.

Definition at line 399 of file star.h.

References nn.

◆ get_nzet()

int Lorene::Star::get_nzet ( ) const

inlineinherited

Returns the number of domains occupied by the star.

Definition at line 358 of file star.h.

References nzet.

◆ get_omega()

double Lorene::Star_rot_Dirac::get_omega ( ) const

inline

Returns the rotation angular velocity $\Omega$ .

Definition at line 178 of file star_rot_dirac.h.

References omega.

◆ get_press()

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

inlineinherited

Returns the fluid pressure.

Definition at line 373 of file star.h.

References press.

◆ get_psi2()

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

inline

Returns $\Psi^2$ .

Definition at line 189 of file star_rot_dirac.h.

References psi2.

◆ get_psi4()

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

inline

Returns the conformal factor $\Psi^4$ .

Definition at line 184 of file star_rot_dirac.h.

References psi4.

◆ get_qqq()

const Scalar & Lorene::Star_rot_Dirac::get_qqq ( ) const

inline

Returns $Q=\Psi^2 N$ .

Definition at line 194 of file star_rot_dirac.h.

References qqq.

◆ get_s_euler()

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

inlineinherited

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

Definition at line 379 of file star.h.

References s_euler.

◆ get_stress_euler()

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

inlineinherited

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

Definition at line 390 of file star.h.

References stress_euler.

◆ get_taa()

const Sym_tensor Lorene::Star_rot_Dirac::get_taa ( ) const

inline

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

Definition at line 234 of file star_rot_dirac.h.

References taa.

◆ get_tgamma()

const Metric Lorene::Star_rot_Dirac::get_tgamma ( ) const

inline

Returns the conformal metric $\tilde{\gamma}_{ij}$ .

Definition at line 222 of file star_rot_dirac.h.

References tgamma.

◆ get_u_euler()

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

inlineinherited

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

Definition at line 385 of file star.h.

References u_euler.

◆ get_v2()

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

inline

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

Definition at line 214 of file star_rot_dirac.h.

References v2.

◆ grv2()

double Lorene::Star_rot_Dirac::grv2 ( ) const

virtual

Error on the virial identity GRV2.

Definition at line 204 of file strot_dirac_global.C.

References aa, Lorene::Scalar::div_r_dzpuis(), Lorene::Star::ener_euler, Lorene::Star::gamma, Lorene::Star::logn, p_grv2, Lorene::Star::press, Lorene::sqrt(), Lorene::Scalar::std_spectral_base(), taa, and v2.

◆ grv3()

double Lorene::Star_rot_Dirac::grv3 ( ) const

virtual

Error on the virial identity GRV3.

Definition at line 273 of file strot_dirac_global.C.

References aa_quad, Lorene::contract(), Lorene::Star::gamma, Lorene::Scalar::integrale(), Lorene::Star::logn, p_grv3, Lorene::Star::s_euler, Lorene::sqrt(), and Lorene::Scalar::std_spectral_base().

◆ hydro_euler()

void Lorene::Star_rot_Dirac::hydro_euler ( )

virtual

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

More later......

Reimplemented from Lorene::Star.

Reimplemented in Lorene::Star_rot_Dirac_diff.

Definition at line 71 of file strot_dirac_hydro.C.

References Lorene::Star::beta, Lorene::contract(), del_deriv(), Lorene::Star::ener, Lorene::Star::ener_euler, Lorene::Star::gam_euler, Lorene::Star::gamma, j_euler, Lorene::Star::nn, Lorene::Star::nzet, omega, Lorene::Star::press, Lorene::Star::s_euler, spectral_filter, Lorene::sqrt(), Lorene::Star::stress_euler, Lorene::Star::u_euler, and v2.

◆ l_surf()

const Itbl & Lorene::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.

Reimplemented in Lorene::Star_rot.

Definition at line 63 of file star_global.C.

References ent, mp, nzet, p_l_surf, and p_xi_surf.

◆ mass_b()

double Lorene::Star_rot_Dirac::mass_b ( ) const

virtual

Baryonic mass.

Implements Lorene::Star.

Definition at line 94 of file strot_dirac_global.C.

References Lorene::Star::gam_euler, Lorene::Star::gamma, Lorene::Scalar::integrale(), Lorene::Star::nbar, Lorene::Star::p_mass_b, Lorene::sqrt(), and Lorene::Scalar::std_spectral_base().

◆ mass_g()

double Lorene::Star_rot_Dirac::mass_g ( ) const

virtual

Gravitational mass.

Implements Lorene::Star.

Definition at line 117 of file strot_dirac_global.C.

References Lorene::Star::beta, Lorene::contract(), Lorene::Star::ener_euler, Lorene::Star::gamma, Lorene::Scalar::integrale(), j_euler, Lorene::Star::nn, Lorene::Star::p_mass_g, Lorene::Star::s_euler, Lorene::sqrt(), and Lorene::Scalar::std_spectral_base().

◆ operator=()

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

Assignment to another Star_rot_Dirac.

Definition at line 266 of file star_rot_dirac.C.

References aa, aa_quad, del_deriv(), flat, hh, j_euler, ln_psi, omega, Lorene::Star::operator=(), psi2, psi4, qqq, spectral_filter, Star_rot_Dirac(), tgamma, and v2.

◆ operator>>()

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

protectedvirtual

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

Reimplemented from Lorene::Star.

Reimplemented in Lorene::Star_rot_Dirac_diff.

Definition at line 315 of file star_rot_dirac.C.

References angu_mom(), aplat(), ellipt(), grv2(), grv3(), mass_g(), Lorene::Star::mp, omega, Lorene::Star::operator>>(), Lorene::pow(), r_circ(), Lorene::Star::ray_eq(), rp_circ(), spectral_filter, and tsw().

◆ r_circ()

double Lorene::Star_rot_Dirac::r_circ ( ) const

virtual

Circumferential equatorial radius.

Definition at line 331 of file strot_dirac_global.C.

References Lorene::Star::gamma, Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_t(), Lorene::Star::mp, Lorene::Star::nzet, p_r_circ, Lorene::Star::ray_eq(), and Lorene::sqrt().

◆ ray_eq()

double Lorene::Star::ray_eq ( ) const

inherited

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

Definition at line 108 of file star_global.C.

References Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), l_surf(), mp, p_ray_eq, Lorene::phi, and xi_surf().

◆ ray_eq_3pis2()

double Lorene::Star::ray_eq_3pis2 ( ) const

inherited

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

Definition at line 233 of file star_global.C.

References Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), l_surf(), mp, p_ray_eq_3pis2, Lorene::phi, ray_eq_pis2(), and xi_surf().

◆ ray_eq_pi()

double Lorene::Star::ray_eq_pi ( ) const

inherited

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

Definition at line 186 of file star_global.C.

References Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), l_surf(), mp, p_ray_eq_pi, Lorene::phi, ray_eq(), and xi_surf().

◆ ray_eq_pis2()

double Lorene::Star::ray_eq_pis2 ( ) const

inherited

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

Definition at line 138 of file star_global.C.

References Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), l_surf(), mp, p_ray_eq_pis2, Lorene::phi, and xi_surf().

◆ ray_pole()

double Lorene::Star::ray_pole ( ) const

inherited

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

Definition at line 278 of file star_global.C.

References Lorene::Mg3d::get_type_t(), l_surf(), mp, p_ray_pole, Lorene::phi, and xi_surf().

◆ rp_circ()

double Lorene::Star_rot_Dirac::rp_circ ( ) const

virtual

◆ sauve()

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

virtual

Save in a file.

Reimplemented from Lorene::Star.

Reimplemented in Lorene::Star_rot_Dirac_diff.

Definition at line 299 of file star_rot_dirac.C.

References Lorene::Star::beta, Lorene::fwrite_be(), hh, omega, qqq, Lorene::Star::sauve(), and spectral_filter.

◆ set_der_0x0()

void Lorene::Star_rot_Dirac::set_der_0x0 ( ) const

protectedvirtual

Sets to 0x0 all the pointers on derived quantities.

Reimplemented from Lorene::Star.

Definition at line 234 of file star_rot_dirac.C.

References p_angu_mom, p_grv2, p_grv3, p_r_circ, p_rp_circ, and p_tsw.

◆ set_enthalpy()

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

inherited

Assignment of the enthalpy field.

Definition at line 379 of file star.C.

References del_deriv(), ent, and equation_of_state().

◆ set_mp()

Map & Lorene::Star::set_mp ( )

inlineinherited

Read/write of the mapping.

Definition at line 322 of file star.h.

References Lorene::Map(), and mp.

◆ solve_hij()

void Lorene::Star_rot_Dirac::solve_hij ( Sym_tensor_trans & hij_new ) const

Solution of the tensor Poisson equation for rotating stars in Dirac gauge.

Definition at line 78 of file strot_dirac_solvehij.C.

References aa, Lorene::Star::beta, Lorene::contract(), Lorene::Tensor::derive_cov(), Lorene::Tensor_sym::derive_cov(), Lorene::Scalar::derive_lie(), Lorene::Sym_tensor::derive_lie(), hh, Lorene::Scalar::inc_dzpuis(), Lorene::Tensor::inc_dzpuis(), ln_psi, Lorene::Star::mp, Lorene::Star::nn, Lorene::Sym_tensor_trans::poisson(), psi2, psi4, qqq, Lorene::Star::s_euler, Lorene::Tensor::set(), Lorene::Scalar::set_dzpuis(), Lorene::Scalar::set_etat_zero(), Lorene::Star::stress_euler, and tgamma.

◆ solve_logn_f()

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

Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.

Definition at line 62 of file strot_dirac_solvenq.C.

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

◆ solve_logn_q()

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

Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.

Definition at line 76 of file strot_dirac_solvenq.C.

References aa_quad, Lorene::contract(), Lorene::Tensor::derive_cov(), hh, Lorene::Tensor::inc_dzpuis(), ln_psi, Lorene::Star::logn, Lorene::Star::mp, Lorene::Scalar::poisson(), psi4, tgamma, and Lorene::Tensor::up_down().

◆ solve_qqq()

void Lorene::Star_rot_Dirac::solve_qqq ( Scalar & q_new ) const

Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.

Definition at line 105 of file strot_dirac_solvenq.C.

References aa_quad, Lorene::contract(), Lorene::Tensor::derive_cov(), Lorene::Scalar::get_etat(), hh, ln_psi, Lorene::Star::mp, Lorene::Star::nn, Lorene::Scalar::poisson(), psi2, psi4, qqq, Lorene::Star::s_euler, Lorene::Scalar::std_spectral_base(), and tgamma.

◆ solve_shift()

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

Solution of the shift equation for rotating stars in Dirac gauge.

Definition at line 54 of file strot_dirac_solveshift.C.

References aa, Lorene::Star::beta, Lorene::contract(), Lorene::Vector::div_free(), hh, j_euler, ln_psi, Lorene::Star::mp, Lorene::Star::nn, Lorene::Vector_divfree::poisson(), psi4, Lorene::Vector::set(), Lorene::Scalar::set_dzpuis(), and tgamma.

◆ spectral_filter_order()

int Lorene::Star_rot_Dirac::spectral_filter_order ( ) const

inline

Returns the filtering order.

Definition at line 173 of file star_rot_dirac.h.

References spectral_filter.

◆ tsw()

double Lorene::Star_rot_Dirac::tsw ( ) const

virtual

Ratio T/W.

Reimplemented in Lorene::Star_rot_Dirac_diff.

Definition at line 177 of file strot_dirac_global.C.

References angu_mom(), Lorene::Star::ener, Lorene::Star::gam_euler, Lorene::Star::gamma, Lorene::Scalar::integrale(), mass_g(), omega, p_tsw, Lorene::sqrt(), and Lorene::Scalar::std_spectral_base().

◆ update_metric()

void Lorene::Star_rot_Dirac::update_metric ( )

Computes metric quantities from known potentials.

The calculation is performed starting from qqq, logn, shift, hh, which are supposed to be up to date. From these, the following fields are updated: nnn, psi4, psi2, ln_psi, tgamma, aa, taa, and aa_quad.

Definition at line 59 of file strot_dirac_upmetr.C.

References aa, aa_quad, Lorene::Star::beta, Lorene::contract(), del_deriv(), Lorene::exp(), flat, Lorene::Star::gamma, hh, ln_psi, Lorene::log(), Lorene::Star::logn, Lorene::Star::nn, psi2, psi4, qqq, Lorene::sqrt(), taa, and tgamma.

◆ xi_surf()

const Tbl & Lorene::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 89 of file star_global.C.

References l_surf(), p_l_surf, and p_xi_surf.

Member Data Documentation

◆ aa

Sym_tensor Lorene::Star_rot_Dirac::aa

protected

$A^{ij}$

Definition at line 87 of file star_rot_dirac.h.

◆ aa_quad

Scalar Lorene::Star_rot_Dirac::aa_quad

protected

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

Definition at line 89 of file star_rot_dirac.h.

◆ beta

Vector Lorene::Star::beta

protectedinherited

Shift vector.

Definition at line 228 of file star.h.

◆ ener

Scalar Lorene::Star::ener

protectedinherited

Total energy density in the fluid frame.

Definition at line 193 of file star.h.

◆ ener_euler

Scalar Lorene::Star::ener_euler

protectedinherited

Total energy density in the Eulerian frame.

Definition at line 198 of file star.h.

◆ ent

Scalar Lorene::Star::ent

protectedinherited

Log-enthalpy.

Definition at line 190 of file star.h.

◆ eos

const Eos& Lorene::Star::eos

protectedinherited

Equation of state of the stellar matter.

Definition at line 185 of file star.h.

◆ flat

const Metric_flat& Lorene::Star_rot_Dirac::flat

protected

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

Definition at line 84 of file star_rot_dirac.h.

◆ gam_euler

Scalar Lorene::Star::gam_euler

protectedinherited

Lorentz factor between the fluid and Eulerian observers.

Definition at line 204 of file star.h.

◆ gamma

Metric Lorene::Star::gamma

protectedinherited

3-metric

Definition at line 235 of file star.h.

◆ hh

Sym_tensor_trans Lorene::Star_rot_Dirac::hh

protected

$h^{ij}$ is defined by $\tilde{\gamma}^{ij}=f^{ij}+h^{ij}$ .

We impose the Dirac gauge $D_j h^{ij}$ explicitly by defining $h^{ij}$ to be a symmetric transverse tensor.

Definition at line 96 of file star_rot_dirac.h.

◆ j_euler

Vector Lorene::Star_rot_Dirac::j_euler

protected

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

Definition at line 77 of file star_rot_dirac.h.

◆ ln_psi

Scalar Lorene::Star_rot_Dirac::ln_psi

protected

$ln(\Psi)$

Definition at line 67 of file star_rot_dirac.h.

◆ lnq

Scalar Lorene::Star::lnq

protectedinherited

Definition at line 232 of file star.h.

◆ logn

Scalar Lorene::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 222 of file star.h.

◆ mp

Map& Lorene::Star::mp

protectedinherited

Mapping associated with the star.

Definition at line 180 of file star.h.

◆ nbar

Scalar Lorene::Star::nbar

protectedinherited

Baryon density in the fluid frame.

Definition at line 192 of file star.h.

◆ nn

Scalar Lorene::Star::nn

protectedinherited

Lapse function N .

Definition at line 225 of file star.h.

◆ nzet

int Lorene::Star::nzet

protectedinherited

Number of domains of *mp occupied by the star.

Definition at line 183 of file star.h.

◆ omega

double Lorene::Star_rot_Dirac::omega

protected

Rotation angular velocity ([f_unit] ).

Definition at line 59 of file star_rot_dirac.h.

◆ p_angu_mom

double* Lorene::Star_rot_Dirac::p_angu_mom

mutableprotected

Angular momentum.

Definition at line 106 of file star_rot_dirac.h.

◆ p_grv2

double* Lorene::Star_rot_Dirac::p_grv2

mutableprotected

Error on the virial identity GRV2.

Definition at line 107 of file star_rot_dirac.h.

◆ p_grv3

double* Lorene::Star_rot_Dirac::p_grv3

mutableprotected

Error on the virial identity GRV3.

Definition at line 108 of file star_rot_dirac.h.

◆ p_l_surf

Itbl* Lorene::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 260 of file star.h.

◆ p_mass_b

double* Lorene::Star::p_mass_b

mutableprotectedinherited

Baryon mass.

Definition at line 268 of file star.h.

◆ p_mass_g

double* Lorene::Star::p_mass_g

mutableprotectedinherited

Gravitational mass.

Definition at line 269 of file star.h.

◆ p_r_circ

double* Lorene::Star_rot_Dirac::p_r_circ

mutableprotected

Circumferential equatorial radius.

Definition at line 110 of file star_rot_dirac.h.

◆ p_ray_eq

double* Lorene::Star::p_ray_eq

mutableprotectedinherited

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

Definition at line 242 of file star.h.

◆ p_ray_eq_3pis2

double* Lorene::Star::p_ray_eq_3pis2

mutableprotectedinherited

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

Definition at line 251 of file star.h.

◆ p_ray_eq_pi

double* Lorene::Star::p_ray_eq_pi

mutableprotectedinherited

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

Definition at line 248 of file star.h.

◆ p_ray_eq_pis2

double* Lorene::Star::p_ray_eq_pis2

mutableprotectedinherited

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

Definition at line 245 of file star.h.

◆ p_ray_pole

double* Lorene::Star::p_ray_pole

mutableprotectedinherited

Coordinate radius at $\theta=0$ .

Definition at line 254 of file star.h.

◆ p_rp_circ

double* Lorene::Star_rot_Dirac::p_rp_circ

mutableprotected

Circumferential polar radius.

Definition at line 111 of file star_rot_dirac.h.

◆ p_tsw

double* Lorene::Star_rot_Dirac::p_tsw

mutableprotected

Ratio T/W.

Definition at line 109 of file star_rot_dirac.h.

◆ p_xi_surf

Tbl* Lorene::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 266 of file star.h.

◆ press

Scalar Lorene::Star::press

protectedinherited

Fluid pressure.

Definition at line 194 of file star.h.

◆ psi2

Scalar Lorene::Star_rot_Dirac::psi2

protected

$\Psi^2$

Definition at line 65 of file star_rot_dirac.h.

◆ psi4

Scalar Lorene::Star_rot_Dirac::psi4

protected

Conformal factor $\Psi^4$ .

Definition at line 64 of file star_rot_dirac.h.

◆ qqq

Scalar Lorene::Star_rot_Dirac::qqq

protected

$Q = \Psi^2 N$

Definition at line 66 of file star_rot_dirac.h.

◆ s_euler

Scalar Lorene::Star::s_euler

protectedinherited

Trace of the stress scalar in the Eulerian frame.

Definition at line 201 of file star.h.

◆ spectral_filter

int Lorene::Star_rot_Dirac::spectral_filter

protected

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 57 of file star_rot_dirac.h.

◆ stress_euler

Sym_tensor Lorene::Star::stress_euler

protectedinherited

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

Definition at line 212 of file star.h.

◆ taa

Sym_tensor Lorene::Star_rot_Dirac::taa

protected

$\tilde{A}_{ij}$

Definition at line 88 of file star_rot_dirac.h.

◆ tgamma

Metric Lorene::Star_rot_Dirac::tgamma

protected

$\tilde{\gamma}_{ij}$

Definition at line 86 of file star_rot_dirac.h.

◆ u_euler

Vector Lorene::Star::u_euler

protectedinherited

Fluid 3-velocity with respect to the Eulerian observer.

Definition at line 207 of file star.h.

◆ v2

Scalar Lorene::Star_rot_Dirac::v2

protected

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

Definition at line 78 of file star_rot_dirac.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Star_rot_Dirac() [1/3]

◆ Star_rot_Dirac() [2/3]

◆ Star_rot_Dirac() [3/3]

◆ ~Star_rot_Dirac()

Member Function Documentation

◆ angu_mom()

◆ aplat()

◆ del_deriv()

◆ del_hydro_euler()

◆ ellipt()

◆ equation_of_state()

◆ equilibrium()

◆ equilibrium_spher()

◆ get_aa()

◆ get_aa_quad()

◆ get_beta()

◆ get_ener()

◆ get_ener_euler()

◆ get_ent()

◆ get_eos()

◆ get_gam_euler()

◆ get_gamma()

◆ get_hh()

◆ get_j_euler()

◆ get_ln_psi()

◆ get_lnq()

◆ get_logn()

◆ get_mp()

◆ get_nbar()

◆ get_nn()

◆ get_nzet()

◆ get_omega()

◆ get_press()

◆ get_psi2()

◆ get_psi4()

◆ get_qqq()

◆ get_s_euler()

◆ get_stress_euler()

◆ get_taa()

◆ get_tgamma()

◆ get_u_euler()

◆ get_v2()

◆ grv2()

◆ grv3()

◆ hydro_euler()

◆ l_surf()

◆ mass_b()

◆ mass_g()

◆ operator=()

◆ operator>>()

◆ 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_hij()

◆ solve_logn_f()

◆ solve_logn_q()

◆ solve_qqq()

◆ solve_shift()

◆ spectral_filter_order()

◆ tsw()

◆ update_metric()

◆ xi_surf()

Member Data Documentation

◆ aa

◆ aa_quad

◆ beta

◆ ener