Class for relativistic differentially rotating stars in Dirac gauge and maximal slicing. More...
#include <star_rot_dirac_diff.h>
Public Member Functions | |
| Star_rot_Dirac_diff (Map &mp_i, int nzet_i, const Eos &eos_i, double(*frot_i)(double, const Tbl &), double(*primfrot_i)(double, const Tbl &), const Tbl &par_frot_i) | |
| Standard constructor. | |
| Star_rot_Dirac_diff (const Star_rot_Dirac_diff &) | |
| Copy constructor. | |
| Star_rot_Dirac_diff (Map &mp_i, const Eos &eos_i, FILE *fich, double(*frot_i)(double, const Tbl &), double(*primfrot_i)(double, const Tbl &)) | |
Constructor from a file (see sauve(FILE*) ). | |
| virtual | ~Star_rot_Dirac_diff () |
| Destructor. | |
| void | operator= (const Star_rot_Dirac_diff &) |
Assignment to another Star_rot_Dirac_diff. | |
| const Scalar & | get_omega_field () const |
Returns the angular velocity field . | |
| virtual double | get_omega_c () const |
Returns the central value of the rotation angular velocity ([f_unit] ). | |
| virtual void | sauve (FILE *) const |
| Save in a file. | |
| virtual double | tsw () const |
| Ratio T/W. | |
| virtual void | hydro_euler () |
| Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid frame. | |
| void | fait_omega_field (double omeg_min, double omeg_max, double precis, int nitermax) |
Computes (member omega_field ). | |
| void | fait_prim_field () |
Computes the member prim_field from omega_field . | |
| double | funct_omega (double omeg) const |
Evaluates , where F is the function defining the rotation profile. | |
| double | prim_funct_omega (double omeg) const |
Evaluates the primitive of , where F is the function defining the rotation profile. | |
| 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) |
| Computes an equilibrium configuration. | |
| int | spectral_filter_order () const |
| Returns the filtering order. | |
| double | get_omega () const |
Returns the rotation angular velocity . | |
| const Scalar & | get_psi4 () const |
Returns the conformal factor . | |
| const Scalar & | get_psi2 () const |
Returns . | |
| const Scalar & | get_qqq () const |
Returns . | |
| const Scalar & | get_ln_psi () const |
Returns . | |
| const Vector & | get_j_euler () const |
| Returns the momentum density 3-vector with respect to the Eulerian observer. | |
| const Scalar & | get_v2 () const |
Reutrns . | |
| const Metric | get_tgamma () const |
Returns the conformal metric . | |
| const Sym_tensor | get_aa () const |
Returns . | |
| const Sym_tensor | get_taa () const |
Returns . | |
| const Scalar | get_aa_quad () const |
Returns . | |
| const Sym_tensor_trans | get_hh () const |
Returns . | |
| 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 | 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. | |
| void | update_metric () |
| Computes metric quantities from known potentials. | |
| 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 . | |
| const Scalar & | get_lnq () const |
| const Metric & | get_gamma () const |
Returns the 3-metric . | |
| double | ray_eq () const |
Coordinate radius at , [r_unit]. | |
| double | ray_eq_pis2 () const |
Coordinate radius at , [r_unit]. | |
| double | ray_eq_pi () const |
Coordinate radius at , [r_unit]. | |
| double | ray_eq_3pis2 () const |
Coordinate radius at , [r_unit]. | |
| double | ray_pole () const |
Coordinate radius at [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 . | |
| const Tbl & | xi_surf () const |
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate on the surface at the collocation points in . | |
Protected Member Functions | |
| virtual ostream & | operator>> (ostream &) const |
| Operator >> (virtual function called by the operator <<). | |
| 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. | |
Protected Attributes | |
| double(* | frot )(double, const Tbl &) |
Function defining the rotation profile. | |
| double(* | primfrot )(double, const Tbl &) |
Primitive of the function , which vanishes at the stellar center. | |
| Tbl | par_frot |
Parameters of the function . | |
| Scalar | omega_field |
Field . | |
| double | omega_min |
Minimum value of . | |
| double | omega_max |
Maximum value of . | |
| Scalar | prim_field |
Field . | |
| int | spectral_filter |
| Spectral exponential filtering order. | |
| double | omega |
Rotation angular velocity ([f_unit] ). | |
| Scalar | psi4 |
Conformal factor . | |
| Scalar | psi2 |
| |
| Scalar | qqq |
| |
| Scalar | ln_psi |
| |
| Vector | j_euler |
| Momentum density 3-vector with respect to the Eulerian observer. | |
| Scalar | v2 |
| |
| const Metric_flat & | flat |
flat metric (spherical components) | |
| Metric | tgamma |
| |
| Sym_tensor | aa |
| |
| Sym_tensor | taa |
| |
| Scalar | aa_quad |
| |
| Sym_tensor_trans | hh |
is defined by . | |
| 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 , . | |
| double * | p_ray_eq_pis2 |
Coordinate radius at , . | |
| double * | p_ray_eq_pi |
Coordinate radius at , . | |
| double * | p_ray_eq_3pis2 |
Coordinate radius at , . | |
| double * | p_ray_pole |
Coordinate radius at . | |
| 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 . | |
| Tbl * | p_xi_surf |
Description of the stellar surface: 2-D Tbl containing the values of the radial coordinate on the surface at the collocation points in . | |
| double * | p_mass_b |
| Baryon mass. | |
| double * | p_mass_g |
| Gravitational mass. | |
Friends | |
| ostream & | operator<< (ostream &, const Star &) |
| Display. | |
| ostream & | operator<< (ostream &, const Star &) |
| Display. | |
Class for relativistic differentially rotating stars in Dirac gauge and maximal slicing.
(*** Under development ***) ()
Definition at line 48 of file star_rot_dirac_diff.h.
| Star_rot_Dirac_diff::Star_rot_Dirac_diff | ( | Map & | mp_i, | |
| int | nzet_i, | |||
| const Eos & | eos_i, | |||
| double(*)(double, const Tbl &) | frot_i, | |||
| double(*)(double, const Tbl &) | primfrot_i, | |||
| const Tbl & | par_frot_i | |||
| ) |
Standard constructor.
| mp_i | Mapping on which the star will be defined | |
| nzet_i | Number of domains occupied by the star | |
| relat | should be true for a relativistic star, false for a Newtonian one | |
| eos_i | Equation of state of the stellar matter | |
| frot_i | Function defining the rotation profile. | |
| primfrot_i | Primitive of which vanishes at the stellar center | |
| par_frot_i | Parameters of functions frot_i and primfrot_i , par_frot_i(0) being the central value of . |
Definition at line 52 of file star_rot_dirac_diff.C.
References omega_field, omega_max, omega_min, and prim_field.
| Star_rot_Dirac_diff::Star_rot_Dirac_diff | ( | const Star_rot_Dirac_diff & | star | ) |
Copy constructor.
Definition at line 75 of file star_rot_dirac_diff.C.
| Star_rot_Dirac_diff::Star_rot_Dirac_diff | ( | Map & | mp_i, | |
| const Eos & | eos_i, | |||
| FILE * | fich, | |||
| double(*)(double, const Tbl &) | frot_i, | |||
| double(*)(double, const Tbl &) | primfrot_i | |||
| ) |
Constructor from a file (see sauve(FILE*) ).
| 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 ) | |
| frot_i | Function defining the rotation profile. | |
| primfrot_i | Primitive of which vanishes at the stellar center |
Definition at line 89 of file star_rot_dirac_diff.C.
References fait_prim_field(), fread_be(), Map::get_mg(), Star::mp, omega_field, omega_max, and omega_min.
| Star_rot_Dirac_diff::~Star_rot_Dirac_diff | ( | ) | [virtual] |
Destructor.
Definition at line 115 of file star_rot_dirac_diff.C.
| double Star_rot_Dirac::angu_mom | ( | ) | const [virtual, inherited] |
Angular momentum.
Definition at line 139 of file strot_dirac_global.C.
References contract(), Metric::cov(), Metric::determinant(), Star::gamma, Map::get_bvect_spher(), Scalar::integrale(), Star_rot_Dirac::j_euler, Star::mp, Star_rot_Dirac::p_angu_mom, sqrt(), and Scalar::std_spectral_base().
| double Star_rot_Dirac::aplat | ( | ) | const [virtual, inherited] |
Flattening r_pole/r_eq.
Definition at line 414 of file strot_dirac_global.C.
References Star::ray_eq(), and Star::ray_pole().
| void Star_rot_Dirac::del_deriv | ( | ) | const [protected, virtual, inherited] |
Deletes all the derived quantities.
Reimplemented from Star.
Definition at line 214 of file star_rot_dirac.C.
References Star_rot_Dirac::p_angu_mom, Star_rot_Dirac::p_grv2, Star_rot_Dirac::p_grv3, Star_rot_Dirac::p_r_circ, Star_rot_Dirac::p_rp_circ, Star_rot_Dirac::p_tsw, and Star_rot_Dirac::set_der_0x0().
| void Star_rot_Dirac::del_hydro_euler | ( | ) | [protected, virtual, inherited] |
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer.
Reimplemented from Star.
Definition at line 242 of file star_rot_dirac.C.
References Star_rot_Dirac::del_deriv(), Star_rot_Dirac::j_euler, Scalar::set_etat_nondef(), Tensor::set_etat_nondef(), and Star_rot_Dirac::v2.
| double Star_rot_Dirac::ellipt | ( | ) | const [virtual, inherited] |
Ellipticity e.
Defined as
, where
and
are, respectively, the equatorial and polar circumferential radius, given by r_circ() and rp_circ().
Definition at line 424 of file strot_dirac_global.C.
References Star_rot_Dirac::r_circ(), Star_rot_Dirac::rp_circ(), and sqrt().
| void Star::equation_of_state | ( | ) | [inherited] |
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
Reimplemented in Gravastar.
Definition at line 458 of file star.C.
References Param::add_int(), Scalar::allocate_all(), Star::del_deriv(), Star::ener, Eos::ener_ent(), Star::ent, Star::eos, Mg3d::get_grille3d(), Map::get_mg(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Star::mp, Star::nbar, Eos::nbar_ent(), Star::nzet, Star::press, Eos::press_ent(), Mtbl::set(), Scalar::set_domain(), Scalar::set_etat_qcq(), Tbl::set_etat_qcq(), Mtbl::set_etat_qcq(), Scalar::std_spectral_base(), Mtbl::t, and Grille3d::x.
| void Star_rot_Dirac_diff::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 | |||
| ) | [virtual] |
Computes an equilibrium configuration.
| 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, 1. = 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 8) to control the iteration:
| |
| control | [input] Set of parameters (stored as a 1-D Tbl of size 7) to control the iteration:
| |
| 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 :
|
Reimplemented from Star_rot_Dirac.
Definition at line 46 of file strot_dirac_diff_equil.C.
References Scalar::annule(), Star::beta, contract(), Metric::cov(), diffrel(), Star::ent, Star::equation_of_state(), fait_omega_field(), Star_rot_Dirac::flat, Star::gam_euler, Star::gamma, Map::get_bvect_spher(), Map::get_mg(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Mg3d::get_type_t(), Star_rot_Dirac::grv2(), Star_rot_Dirac::hh, Map::homothetie(), hydro_euler(), log(), log10(), Star::logn, Star::mp, Scalar::mult_rsint(), Star::nn, Star::nzet, Star_rot_Dirac::omega, omega_field, par_frot, prim_field, Star_rot_Dirac::qqq, Star::ray_eq(), Vector::set(), Tbl::set(), Tbl::set_etat_qcq(), Scalar::set_etat_zero(), Star_rot_Dirac::solve_hij(), Star_rot_Dirac::solve_logn_f(), Star_rot_Dirac::solve_logn_q(), Star_rot_Dirac::solve_qqq(), Star_rot_Dirac::solve_shift(), sqrt(), Scalar::std_spectral_base(), Star::u_euler, Star_rot_Dirac::update_metric(), Star_rot_Dirac::v2, and Scalar::val_grid_point().
| void Star::equilibrium_spher | ( | double | ent_c, | |
| double | precis = 1.e-14, |
|||
| const Tbl * | pent_limit = 0x0 | |||
| ) | [virtual, inherited] |
Computes a spherical static configuration.
| 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 91 of file star_equil_spher.C.
References Map_et::adapt(), Param::add_double(), Param::add_int(), Param::add_int_mod(), Param::add_tbl(), Scalar::annule(), diffrel(), Scalar::dsdr(), Map_af::dsdr(), Star::ener, Star::ener_euler, Star::ent, Star::equation_of_state(), exp(), Star::gam_euler, Star::gamma, Map_et::get_alpha(), Map_af::get_alpha(), Map_et::get_beta(), Map_af::get_beta(), Map::get_bvect_cart(), Map::get_mg(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Map_af::homothetie(), Scalar::integrale(), Star::logn, Star::mass_b(), Star::mass_g(), Star::mp, Star::nn, norme(), Star::nzet, Map_af::poisson(), Star::press, Star::s_euler, Vector::set(), Map_af::set_alpha(), Map_af::set_beta(), Scalar::set_dzpuis(), Cmp::set_etat_qcq(), Scalar::set_etat_zero(), sqrt(), Scalar::std_spectral_base(), Star::u_euler, Scalar::val_grid_point(), and Map::val_r().
| void Star_rot_Dirac_diff::fait_omega_field | ( | double | omeg_min, | |
| double | omeg_max, | |||
| double | precis, | |||
| int | nitermax | |||
| ) |
Computes
(member omega_field ).
The computation amounts to solving the equation
for
.
| 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 51 of file strot_dirac_diff_faitomeg.C.
References Param::add_int(), Param::add_scalar(), Param::add_star(), Scalar::allocate_all(), Scalar::annule(), Star::beta, Metric::cov(), Scalar::div_rsint(), Scalar::domain(), fait_prim_field(), Star::gamma, Scalar::get_etat(), Map::get_mg(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), max(), min(), Star::mp, Scalar::mult_rsint(), Star::nzet, omega_field, omega_max, omega_min, Star_rot_Dirac::psi4, Star_rot_Dirac::qqq, Tbl::set(), Scalar::set_domain(), Scalar::std_spectral_base(), and zerosec().
| void Star_rot_Dirac_diff::fait_prim_field | ( | ) |
Computes the member prim_field from omega_field .
Definition at line 175 of file strot_dirac_diff_faitomeg.C.
References Scalar::allocate_all(), Map::get_mg(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Star::mp, Star::nzet, omega_field, par_frot, prim_field, primfrot, Tbl::set(), Scalar::set_domain(), Scalar::std_spectral_base(), and Scalar::val_grid_point().
| double Star_rot_Dirac_diff::funct_omega | ( | double | omeg | ) | const |
Evaluates
, where F is the function defining the rotation profile.
This function is linked to the components of the fluid 4-velocity by
funct_omega calls frot with the parameters par_frot .
| omeg | [input] value of |
Definition at line 249 of file star_rot_dirac_diff.C.
| const Sym_tensor Star_rot_Dirac::get_aa | ( | ) | const [inline, inherited] |
| const Scalar Star_rot_Dirac::get_aa_quad | ( | ) | const [inline, inherited] |
| const Vector& Star::get_beta | ( | ) | const [inline, inherited] |
| const Scalar& Star::get_ener | ( | ) | const [inline, inherited] |
Returns the proper total energy density.
Definition at line 366 of file star.h.
References Star::ener.
| const Scalar& Star::get_ener_euler | ( | ) | const [inline, inherited] |
Returns the total energy density with respect to the Eulerian observer.
Definition at line 372 of file star.h.
References Star::ener_euler.
| const Scalar& Star::get_ent | ( | ) | const [inline, inherited] |
| const Eos& Star::get_eos | ( | ) | const [inline, inherited] |
| const Scalar& Star::get_gam_euler | ( | ) | const [inline, inherited] |
Returns the Lorentz factor between the fluid and Eulerian observers.
Definition at line 378 of file star.h.
References Star::gam_euler.
| const Metric& Star::get_gamma | ( | ) | const [inline, inherited] |
| const Sym_tensor_trans Star_rot_Dirac::get_hh | ( | ) | const [inline, inherited] |
| const Vector& Star_rot_Dirac::get_j_euler | ( | ) | const [inline, inherited] |
Returns the momentum density 3-vector with respect to the Eulerian observer.
Definition at line 208 of file star_rot_dirac.h.
References Star_rot_Dirac::j_euler.
| const Scalar& Star_rot_Dirac::get_ln_psi | ( | ) | const [inline, inherited] |
| const Scalar& Star::get_logn | ( | ) | const [inline, inherited] |
Returns the logarithm of the lapse N.
In the Newtonian case, this is the Newtonian gravitational potential (in units of
).
Definition at line 392 of file star.h.
References Star::logn.
| const Map& Star::get_mp | ( | ) | const [inline, inherited] |
| const Scalar& Star::get_nbar | ( | ) | const [inline, inherited] |
| const Scalar& Star::get_nn | ( | ) | const [inline, inherited] |
| int Star::get_nzet | ( | ) | const [inline, inherited] |
Returns the number of domains occupied by the star.
Definition at line 354 of file star.h.
References Star::nzet.
| double Star_rot_Dirac::get_omega | ( | ) | const [inline, inherited] |
Returns the rotation angular velocity
.
Definition at line 177 of file star_rot_dirac.h.
References Star_rot_Dirac::omega.
| double Star_rot_Dirac_diff::get_omega_c | ( | ) | const [virtual] |
Returns the central value of the rotation angular velocity ([f_unit] ).
Definition at line 261 of file star_rot_dirac_diff.C.
References omega_field, and Scalar::val_grid_point().
| const Scalar& Star_rot_Dirac_diff::get_omega_field | ( | ) | const [inline] |
Returns the angular velocity field
.
Definition at line 148 of file star_rot_dirac_diff.h.
References omega_field.
| const Scalar& Star::get_press | ( | ) | const [inline, inherited] |
| const Scalar& Star_rot_Dirac::get_psi2 | ( | ) | const [inline, inherited] |
| const Scalar& Star_rot_Dirac::get_psi4 | ( | ) | const [inline, inherited] |
Returns the conformal factor
.
Definition at line 183 of file star_rot_dirac.h.
References Star_rot_Dirac::psi4.
| const Scalar& Star_rot_Dirac::get_qqq | ( | ) | const [inline, inherited] |
| const Scalar& Star::get_s_euler | ( | ) | const [inline, inherited] |
Returns the trace of the stress tensor in the Eulerian frame.
Definition at line 375 of file star.h.
References Star::s_euler.
| const Tensor& Star::get_stress_euler | ( | ) | const [inline, inherited] |
Returns the spatial part of the stress-energy tensor with respect to the Eulerian observer.
Definition at line 386 of file star.h.
References Star::stress_euler.
| const Sym_tensor Star_rot_Dirac::get_taa | ( | ) | const [inline, inherited] |
| const Metric Star_rot_Dirac::get_tgamma | ( | ) | const [inline, inherited] |
Returns the conformal metric
.
Definition at line 221 of file star_rot_dirac.h.
References Star_rot_Dirac::tgamma.
| const Vector& Star::get_u_euler | ( | ) | const [inline, inherited] |
Returns the fluid 3-velocity with respect to the Eulerian observer.
Definition at line 381 of file star.h.
References Star::u_euler.
| const Scalar& Star_rot_Dirac::get_v2 | ( | ) | const [inline, inherited] |
| double Star_rot_Dirac::grv2 | ( | ) | const [virtual, inherited] |
Error on the virial identity GRV2.
Definition at line 200 of file strot_dirac_global.C.
References Star_rot_Dirac::aa, Metric::con(), Metric::cov(), Scalar::div_r_dzpuis(), Scalar::dsdr(), Scalar::dsdt(), Star::ener_euler, Star::gamma, Star::logn, Star_rot_Dirac::p_grv2, Star::press, sqrt(), Scalar::std_spectral_base(), Star_rot_Dirac::taa, and Star_rot_Dirac::v2.
| double Star_rot_Dirac::grv3 | ( | ) | const [virtual, inherited] |
Error on the virial identity GRV3.
Definition at line 269 of file strot_dirac_global.C.
References Star_rot_Dirac::aa_quad, Metric::con(), Metric::connect(), contract(), Scalar::derive_con(), Scalar::derive_cov(), Metric::determinant(), Star::gamma, Connection::get_delta(), Scalar::integrale(), Star::logn, Star_rot_Dirac::p_grv3, Star::s_euler, sqrt(), and Scalar::std_spectral_base().
| void Star_rot_Dirac_diff::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 Star_rot_Dirac.
Definition at line 43 of file strot_dirac_diff_hydro.C.
References Star::beta, Metric::con(), contract(), Metric::cov(), Star_rot_Dirac::del_deriv(), Star::ener, Star::ener_euler, Star::gam_euler, Star::gamma, Star_rot_Dirac::j_euler, Scalar::mult_rsint(), Star::nn, omega_field, Star::press, Star::s_euler, Vector::set(), Scalar::set_etat_zero(), sqrt(), Tensor::std_spectral_base(), Vector::std_spectral_base(), Scalar::std_spectral_base(), Star::stress_euler, Star::u_euler, and Star_rot_Dirac::v2.
| const Itbl & Star::l_surf | ( | ) | const [virtual, inherited] |
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
.
The stellar surface is defined as the location where the enthalpy (member ent) vanishes.
Reimplemented in Star_rot.
Definition at line 59 of file star_global.C.
References Star::ent, Map::get_mg(), Mg3d::get_np(), Mg3d::get_nt(), Scalar::get_spectral_va(), Star::mp, Star::nzet, Star::p_l_surf, and Star::p_xi_surf.
| double Star_rot_Dirac::mass_b | ( | ) | const [virtual, inherited] |
Baryonic mass.
Implements Star.
Definition at line 90 of file strot_dirac_global.C.
References Metric::determinant(), Star::gam_euler, Star::gamma, Scalar::integrale(), Star::nbar, Star::p_mass_b, sqrt(), and Scalar::std_spectral_base().
| double Star_rot_Dirac::mass_g | ( | ) | const [virtual, inherited] |
Gravitational mass.
Implements Star.
Definition at line 113 of file strot_dirac_global.C.
References Star::beta, contract(), Metric::cov(), Metric::determinant(), Star::ener_euler, Star::gamma, Scalar::integrale(), Star_rot_Dirac::j_euler, Star::nn, Star::p_mass_g, Star::s_euler, sqrt(), and Scalar::std_spectral_base().
| void Star_rot_Dirac_diff::operator= | ( | const Star_rot_Dirac_diff & | star | ) |
Assignment to another Star_rot_Dirac_diff.
Reimplemented from Star_rot_Dirac.
Definition at line 125 of file star_rot_dirac_diff.C.
References frot, omega_field, omega_max, omega_min, par_frot, prim_field, and primfrot.
| ostream & Star_rot_Dirac_diff::operator>> | ( | ostream & | ost | ) | const [protected, virtual] |
Operator >> (virtual function called by the operator <<).
Reimplemented from Star_rot_Dirac.
Definition at line 169 of file star_rot_dirac_diff.C.
References Star_rot_Dirac::angu_mom(), Star_rot_Dirac::aplat(), Star_rot_Dirac::ellipt(), Map::get_mg(), Mg3d::get_nr(), Mg3d::get_nt(), Star_rot_Dirac::grv2(), Star_rot_Dirac::grv3(), Star_rot_Dirac::mass_g(), Star::mp, Star::nzet, Star_rot_Dirac::omega, omega_field, omega_max, omega_min, par_frot, pow(), Star_rot_Dirac::r_circ(), Star::ray_eq(), Star_rot_Dirac::rp_circ(), tsw(), and Scalar::val_grid_point().
| double Star_rot_Dirac_diff::prim_funct_omega | ( | double | omeg | ) | const |
Evaluates the primitive of
, where F is the function defining the rotation profile.
| omeg | [input] value of |
Definition at line 255 of file star_rot_dirac_diff.C.
| double Star_rot_Dirac::r_circ | ( | ) | const [virtual, inherited] |
Circumferential equatorial radius.
Definition at line 327 of file strot_dirac_global.C.
References Metric::cov(), Star::gamma, Map::get_mg(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_type_t(), Star::mp, Star::nzet, Star_rot_Dirac::p_r_circ, Star::ray_eq(), and sqrt().
| double Star::ray_eq | ( | ) | const [inherited] |
Coordinate radius at
,
[r_unit].
Definition at line 104 of file star_global.C.
References Map::get_mg(), Mg3d::get_nt(), Mg3d::get_type_p(), Mg3d::get_type_t(), Star::l_surf(), Star::mp, Star::p_ray_eq, Map::val_r(), and Star::xi_surf().
| double Star::ray_eq_3pis2 | ( | ) | const [inherited] |
Coordinate radius at
,
[r_unit].
Definition at line 229 of file star_global.C.
References Map::get_mg(), Mg3d::get_np(), Mg3d::get_nt(), Mg3d::get_type_p(), Mg3d::get_type_t(), Star::l_surf(), Star::mp, Star::p_ray_eq_3pis2, Star::ray_eq_pis2(), Map::val_r(), and Star::xi_surf().
| double Star::ray_eq_pi | ( | ) | const [inherited] |
Coordinate radius at
,
[r_unit].
Definition at line 182 of file star_global.C.
References Map::get_mg(), Mg3d::get_np(), Mg3d::get_nt(), Mg3d::get_type_p(), Mg3d::get_type_t(), Star::l_surf(), Star::mp, Star::p_ray_eq_pi, Star::ray_eq(), Map::val_r(), and Star::xi_surf().
| double Star::ray_eq_pis2 | ( | ) | const [inherited] |
Coordinate radius at
,
[r_unit].
Definition at line 134 of file star_global.C.
References Map::get_mg(), Mg3d::get_np(), Mg3d::get_nt(), Mg3d::get_type_p(), Mg3d::get_type_t(), Star::l_surf(), Star::mp, Star::p_ray_eq_pis2, Map::val_r(), and Star::xi_surf().
| double Star::ray_pole | ( | ) | const [inherited] |
Coordinate radius at
[r_unit].
Definition at line 274 of file star_global.C.
References Map::get_mg(), Mg3d::get_type_t(), Star::l_surf(), Star::mp, Star::p_ray_pole, Map::val_r(), and Star::xi_surf().
| double Star_rot_Dirac::rp_circ | ( | ) | const [virtual, inherited] |
Circumferential polar radius.
Definition at line 353 of file strot_dirac_global.C.
References Scalar::annule(), Scalar::annule_hard(), Valeur::c_cf, Valeur::coef(), Metric::cov(), Scalar::dsdt(), Star::gamma, Map::get_mg(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Star::l_surf(), Star::mp, Star::nzet, Star_rot_Dirac::p_rp_circ, Scalar::set_grid_point(), Scalar::set_spectral_va(), sqrt(), Scalar::std_spectral_base(), Map::tet, Valeur::val_point_jk(), Map::val_r(), and Star::xi_surf().
| void Star_rot_Dirac_diff::sauve | ( | FILE * | fich | ) | const [virtual] |
Save in a file.
Reimplemented from Star_rot_Dirac.
Definition at line 150 of file star_rot_dirac_diff.C.
References fwrite_be(), omega_field, omega_max, omega_min, par_frot, Scalar::sauve(), and Tbl::sauve().
| void Star_rot_Dirac::set_der_0x0 | ( | ) | const [protected, virtual, inherited] |
Sets to 0x0 all the pointers on derived quantities.
Reimplemented from Star.
Definition at line 230 of file star_rot_dirac.C.
References Star_rot_Dirac::p_angu_mom, Star_rot_Dirac::p_grv2, Star_rot_Dirac::p_grv3, Star_rot_Dirac::p_r_circ, Star_rot_Dirac::p_rp_circ, and Star_rot_Dirac::p_tsw.
| void Star::set_enthalpy | ( | const Scalar & | ent_i | ) | [inherited] |
Assignment of the enthalpy field.
Definition at line 375 of file star.C.
References Star::del_deriv(), Star::ent, and Star::equation_of_state().
| Map& Star::set_mp | ( | ) | [inline, inherited] |
| void Star_rot_Dirac::solve_hij | ( | Sym_tensor_trans & | hij_new | ) | const [inherited] |
Solution of the tensor Poisson equation for rotating stars in Dirac gauge.
Definition at line 74 of file strot_dirac_solvehij.C.
References Star_rot_Dirac::aa, Star::beta, Metric::con(), contract(), Metric::cov(), Tensor_sym::derive_con(), Scalar::derive_con(), Tensor::derive_cov(), Tensor_sym::derive_cov(), Scalar::derive_cov(), Scalar::derive_lie(), Sym_tensor::derive_lie(), Vector::divergence(), Map::flat_met_spher(), Map::get_bvect_spher(), Map::get_mg(), Mg3d::get_np(), Star_rot_Dirac::hh, Scalar::inc_dzpuis(), Tensor::inc_dzpuis(), Star_rot_Dirac::ln_psi, Star::mp, Star::nn, Vector::ope_killing_conf(), Sym_tensor_trans::poisson(), Star_rot_Dirac::psi2, Star_rot_Dirac::psi4, Star_rot_Dirac::qqq, Star::s_euler, Tensor::set(), Scalar::set_dzpuis(), Scalar::set_etat_zero(), Star::stress_euler, Star_rot_Dirac::tgamma, and Tensor::trace().
| void Star_rot_Dirac::solve_logn_f | ( | Scalar & | ln_f_new | ) | const [inherited] |
Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
Definition at line 58 of file strot_dirac_solvenq.C.
References Star::ener_euler, Scalar::poisson(), Star_rot_Dirac::psi4, and Star::s_euler.
| void Star_rot_Dirac::solve_logn_q | ( | Scalar & | ln_q_new | ) | const [inherited] |
Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
Definition at line 72 of file strot_dirac_solvenq.C.
References Star_rot_Dirac::aa_quad, contract(), Scalar::derive_con(), Tensor::derive_cov(), Scalar::derive_cov(), Map::flat_met_spher(), Map::get_bvect_spher(), Star_rot_Dirac::hh, Tensor::inc_dzpuis(), Star_rot_Dirac::ln_psi, Star::logn, Star::mp, Scalar::poisson(), Star_rot_Dirac::psi4, Star_rot_Dirac::tgamma, and Tensor::up_down().
| void Star_rot_Dirac::solve_qqq | ( | Scalar & | q_new | ) | const [inherited] |
Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
Definition at line 101 of file strot_dirac_solvenq.C.
References Star_rot_Dirac::aa_quad, contract(), Metric::cov(), Scalar::derive_con(), Tensor::derive_cov(), Tensor_sym::derive_cov(), Scalar::derive_cov(), Map::flat_met_spher(), Scalar::get_etat(), Star_rot_Dirac::hh, Tenseur::inc_dzpuis(), Star_rot_Dirac::ln_psi, Star::mp, Star::nn, Scalar::poisson(), Star_rot_Dirac::psi2, Star_rot_Dirac::psi4, Star_rot_Dirac::qqq, Star::s_euler, Scalar::std_spectral_base(), and Star_rot_Dirac::tgamma.
| void Star_rot_Dirac::solve_shift | ( | Vector & | shift_new | ) | const [inherited] |
Solution of the shift equation for rotating stars in Dirac gauge.
Definition at line 50 of file strot_dirac_solveshift.C.
References Star_rot_Dirac::aa, Star::beta, Metric::connect(), contract(), Tensor::derive_cov(), Scalar::derive_cov(), Vector::div_free(), Vector::divergence(), Map::flat_met_spher(), Connection::get_delta(), Star_rot_Dirac::hh, Tensor::inc_dzpuis(), Star_rot_Dirac::j_euler, Star_rot_Dirac::ln_psi, Star::mp, Star::nn, Vector_divfree::poisson(), Star_rot_Dirac::psi4, Vector::set(), Scalar::set_dzpuis(), and Star_rot_Dirac::tgamma.
| int Star_rot_Dirac::spectral_filter_order | ( | ) | const [inline, inherited] |
Returns the filtering order.
Definition at line 172 of file star_rot_dirac.h.
References Star_rot_Dirac::spectral_filter.
| double Star_rot_Dirac_diff::tsw | ( | ) | const [virtual] |
Ratio T/W.
Reimplemented from Star_rot_Dirac.
Definition at line 47 of file strot_dirac_diff_global.C.
References contract(), Metric::cov(), Metric::determinant(), Star::ener, Star::gam_euler, Star::gamma, Map::get_bvect_spher(), Scalar::integrale(), Star_rot_Dirac::j_euler, Star_rot_Dirac::mass_g(), Star::mp, omega_field, Star_rot_Dirac::p_tsw, sqrt(), and Scalar::std_spectral_base().
| void Star_rot_Dirac::update_metric | ( | ) | [inherited] |
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 55 of file strot_dirac_upmetr.C.
References Star_rot_Dirac::aa, Star_rot_Dirac::aa_quad, Star::beta, Metric::con(), Metric_flat::con(), contract(), Star_rot_Dirac::del_deriv(), Sym_tensor::derive_lie(), exp(), Star_rot_Dirac::flat, Star::gamma, Star_rot_Dirac::hh, Star_rot_Dirac::ln_psi, log(), Star::logn, Star::nn, Vector::ope_killing_conf(), Star_rot_Dirac::psi2, Star_rot_Dirac::psi4, Star_rot_Dirac::qqq, sqrt(), Scalar::std_spectral_base(), Star_rot_Dirac::taa, Star_rot_Dirac::tgamma, and Tensor::up_down().
| const Tbl & Star::xi_surf | ( | ) | const [inherited] |
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate
on the surface at the collocation points in
.
The stellar surface is defined as the location where the enthalpy (member ent) vanishes.
Definition at line 85 of file star_global.C.
References Star::l_surf(), Star::p_l_surf, and Star::p_xi_surf.
| ostream& operator<< | ( | ostream & | , | |
| const Star & | ||||
| ) | [friend, inherited] |
Display.
| ostream& operator<< | ( | ostream & | ost, | |
| const Star & | et | |||
| ) | [friend] |
Sym_tensor Star_rot_Dirac::aa [protected, inherited] |
Definition at line 86 of file star_rot_dirac.h.
Scalar Star_rot_Dirac::aa_quad [protected, inherited] |
Definition at line 88 of file star_rot_dirac.h.
Vector Star::beta [protected, inherited] |
Scalar Star::ener [protected, inherited] |
Scalar Star::ener_euler [protected, inherited] |
const Metric_flat& Star_rot_Dirac::flat [protected, inherited] |
flat metric
(spherical components)
Definition at line 83 of file star_rot_dirac.h.
double(* Star_rot_Dirac_diff::frot)(double, const Tbl &) [protected] |
Function
defining the rotation profile.
This function is linked to the components of the fluid 4-velocity by
The first argument of frot must be
; the second argument contains the parameters, the first of which being necessarily the central value of
.
Scalar Star::gam_euler [protected, inherited] |
Metric Star::gamma [protected, inherited] |
Sym_tensor_trans Star_rot_Dirac::hh [protected, inherited] |
is defined by
.
We impose the Dirac gauge
explicitly by defining
to be a symmetric transverse tensor.
Definition at line 95 of file star_rot_dirac.h.
Vector Star_rot_Dirac::j_euler [protected, inherited] |
Momentum density 3-vector with respect to the Eulerian observer.
Definition at line 76 of file star_rot_dirac.h.
Scalar Star_rot_Dirac::ln_psi [protected, inherited] |
Definition at line 66 of file star_rot_dirac.h.
Scalar Star::logn [protected, inherited] |
Scalar Star::nbar [protected, inherited] |
int Star::nzet [protected, inherited] |
double Star_rot_Dirac::omega [protected, inherited] |
Rotation angular velocity ([f_unit] ).
Definition at line 58 of file star_rot_dirac.h.
Scalar Star_rot_Dirac_diff::omega_field [protected] |
Field
.
Definition at line 83 of file star_rot_dirac_diff.h.
double Star_rot_Dirac_diff::omega_max [protected] |
Maximum value of
.
Definition at line 86 of file star_rot_dirac_diff.h.
double Star_rot_Dirac_diff::omega_min [protected] |
Minimum value of
.
Definition at line 85 of file star_rot_dirac_diff.h.
double* Star_rot_Dirac::p_angu_mom [mutable, protected, inherited] |
Angular momentum.
Definition at line 105 of file star_rot_dirac.h.
double* Star_rot_Dirac::p_grv2 [mutable, protected, inherited] |
Error on the virial identity GRV2.
Definition at line 106 of file star_rot_dirac.h.
double* Star_rot_Dirac::p_grv3 [mutable, protected, inherited] |
Error on the virial identity GRV3.
Definition at line 107 of file star_rot_dirac.h.
Itbl* Star::p_l_surf [mutable, protected, inherited] |
double* Star::p_mass_b [mutable, protected, inherited] |
double* Star::p_mass_g [mutable, protected, inherited] |
double* Star_rot_Dirac::p_r_circ [mutable, protected, inherited] |
Circumferential equatorial radius.
Definition at line 109 of file star_rot_dirac.h.
double* Star::p_ray_eq [mutable, protected, inherited] |
double* Star::p_ray_eq_3pis2 [mutable, protected, inherited] |
double* Star::p_ray_eq_pi [mutable, protected, inherited] |
double* Star::p_ray_eq_pis2 [mutable, protected, inherited] |
double* Star::p_ray_pole [mutable, protected, inherited] |
double* Star_rot_Dirac::p_rp_circ [mutable, protected, inherited] |
Circumferential polar radius.
Definition at line 110 of file star_rot_dirac.h.
double* Star_rot_Dirac::p_tsw [mutable, protected, inherited] |
Ratio T/W.
Definition at line 108 of file star_rot_dirac.h.
Tbl* Star::p_xi_surf [mutable, protected, inherited] |
Tbl Star_rot_Dirac_diff::par_frot [protected] |
Parameters of the function
.
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 80 of file star_rot_dirac_diff.h.
Scalar Star::press [protected, inherited] |
Scalar Star_rot_Dirac_diff::prim_field [protected] |
Field
.
Definition at line 89 of file star_rot_dirac_diff.h.
double(* Star_rot_Dirac_diff::primfrot)(double, const Tbl &) [protected] |
Primitive of the function
, which vanishes at the stellar center.
The first argument of primfrot must be
; the second argument contains the parameters, the first of which being necessarily the central value of
.
Scalar Star_rot_Dirac::psi2 [protected, inherited] |
Definition at line 64 of file star_rot_dirac.h.
Scalar Star_rot_Dirac::psi4 [protected, inherited] |
Conformal factor
.
Definition at line 63 of file star_rot_dirac.h.
Scalar Star_rot_Dirac::qqq [protected, inherited] |
Definition at line 65 of file star_rot_dirac.h.
Scalar Star::s_euler [protected, inherited] |
int Star_rot_Dirac::spectral_filter [protected, inherited] |
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
nzet .
Definition at line 56 of file star_rot_dirac.h.
Sym_tensor Star::stress_euler [protected, inherited] |
Sym_tensor Star_rot_Dirac::taa [protected, inherited] |
Definition at line 87 of file star_rot_dirac.h.
Metric Star_rot_Dirac::tgamma [protected, inherited] |
Definition at line 85 of file star_rot_dirac.h.
Vector Star::u_euler [protected, inherited] |
Scalar Star_rot_Dirac::v2 [protected, inherited] |
Definition at line 77 of file star_rot_dirac.h.
1.6.1