LORENE
|
Public Member Functions | |
Et_magnetisation (Map &mp_i, int nzet_i, bool relat, const Eos &eos_i, bool include_mag=true, bool use_B=true) | |
Standard constructor. | |
Et_magnetisation (const Et_magnetisation &) | |
Copy constructor. | |
Et_magnetisation (Map &mp_i, const Eos &eos_i, FILE *fich) | |
Constructor from a file (see sauve(FILE*) ). | |
virtual | ~Et_magnetisation () |
Destructor. | |
void | operator= (const Et_magnetisation &) |
Assignment to another Et_rot_mag. | |
virtual void | equation_of_state () |
Computes the proper baryon and energy density, as well as pressure from the enthalpy. | |
bool | B_in_eos () const |
Public accessor to the use_B_in_eos flag. | |
bool | use_magnetisation () const |
Public accessor to the include_magnetisation flag. | |
const Scalar & | get_magnetisation () const |
Accessor to the magnetisation scalar field. | |
const Scalar & | get_E_I () const |
Accessor to the interaction energy density. | |
const Vector & | get_J_I () const |
Accessor to the interaction momentum vector. | |
const Sym_tensor & | get_Sij_I () const |
Accessor to the interaction stress tensor. | |
virtual void | sauve (FILE *) const |
Save in a file. | |
virtual ostream & | operator>> (ostream &) const |
Operator >> (virtual function called by the operator <<). | |
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 (ostream *ost=0x0) const |
Error on the virial identity GRV3. | |
virtual double | mom_quad_old () const |
Part of the quadrupole moment. | |
virtual double | mom_quad_Bo () const |
Part of the quadrupole moment. | |
virtual void | magnet_comput (const int adapt_flag, Cmp(*f_j)(const Cmp &x, const double), Param &par_poisson_At, Param &par_poisson_Avect) |
Computes the electromagnetic quantities solving the Maxwell equations (6) and (7) of [Bocquet, Bonazzola, Gourgoulhon and Novak, Astron. | |
virtual void | MHD_comput () |
Computes the electromagnetic part of the stress-energy tensor. | |
void | equilibrium_mag (double ent_c, double omega0, double fact_omega, int nzadapt, const Tbl &ent_limit, const Itbl &icontrol, const Tbl &control, double mbar_wanted, double magmom_wanted, double aexp_mass, Tbl &diff, double Q0, double a_j0, Cmp(*f_j)(const Cmp &x, const double), Cmp(*M_j)(const Cmp &x, const double)) |
Computes an equilibrium configuration. | |
bool | is_conduct () const |
Tells if the star is made of conducting or isolating material. | |
const Cmp & | get_At () const |
Returns the t component of the electromagnetic potential, divided by ![]() | |
const Cmp & | get_Aphi () const |
Returns the ![]() ![]() | |
const Cmp & | get_Bphi () const |
Returns the ![]() | |
const Cmp & | get_jt () const |
Returns the t component of the current 4-vector. | |
const Cmp & | get_jphi () const |
Returns the ![]() | |
const Tenseur & | get_Eem () const |
Returns the electromagnetic energy density in the Eulerian frame. | |
const Tenseur & | get_Jpem () const |
Returns the ![]() | |
const Tenseur & | get_Srrem () const |
Returns the rr-component of the electromagnetic stress 3-tensor, as measured in the Eulerian frame. | |
const Tenseur & | get_Sppem () const |
Returns the ![]() | |
double | get_Q () const |
Returns the requested electric charge in the case of a perfect conductor and the charge/baryon for an isolator. | |
double | get_a_j () const |
Returns the amplitude of the current/charge function. | |
Tenseur | Elec () const |
Computes the electric field spherical components in Lorene's units. | |
Tenseur | Magn () const |
Computes the magnetic field spherical components in Lorene's units. | |
virtual double | tsw () const |
Ratio T/W. | |
double | MagMom () const |
Magnetic Momentum ![]() | |
double | Q_comput () const |
Computed charge deduced from the asymptotic behaviour of At [SI units]. | |
double | Q_int () const |
Computed charge from the integration of charge density over the star (i.e. | |
double | GyroMag () const |
Gyromagnetic ratio ![]() | |
virtual void | magnet_comput_plus (const int adapt_flag, const int initial_j, const Tbl an_j, Cmp(*f_j)(const Cmp &x, const Tbl), const Tbl bn_j, Cmp(*g_j)(const Cmp &x, const Tbl), Cmp(*N_j)(const Cmp &x, const Tbl), Param &par_poisson_At, Param &par_poisson_Avect) |
Computes the electromagnetic quantities solving the Maxwell equations (6) and (7) of [Bocquet, Bonazzola, Gourgoulhon and Novak, Astron. | |
void | equilibrium_mag (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, const double Q0, const double a_j0, Cmp(*f_j)(const Cmp &x, const double), Cmp(*M_j)(const Cmp &x, const double)) |
Computes an equilibrium configuration. | |
void | equilibrium_mag_plus (const Itbl &icontrol, const Tbl &control, Tbl &diff, const int initial_j, const Tbl an_j, Cmp(*f_j)(const Cmp &x, const Tbl), Cmp(*M_j)(const Cmp &x, const Tbl), const Tbl bn_j, Cmp(*g_j)(const Cmp &x, const Tbl), Cmp(*N_j)(const Cmp &x, const Tbl), const double relax_mag) |
Computes an equilibrium configuration. | |
virtual double | get_omega_c () const |
Returns the central value of the rotation angular velocity ( [f_unit] ) | |
const Tenseur & | get_bbb () const |
Returns the metric factor B. | |
const Tenseur & | get_b_car () const |
Returns the square of the metric factor B. | |
const Tenseur & | get_nphi () const |
Returns the metric coefficient ![]() | |
const Tenseur & | get_tnphi () const |
Returns the component ![]() | |
const Tenseur & | get_uuu () const |
Returns the norm of u_euler . | |
const Tenseur & | get_logn () const |
Returns the metric potential ![]() logn_auto . | |
const Tenseur & | get_nuf () const |
Returns the part of the Metric potential ![]() logn generated by the matter terms. | |
const Tenseur & | get_nuq () const |
Returns the Part of the Metric potential ![]() logn generated by the quadratic terms. | |
const Tenseur & | get_dzeta () const |
Returns the Metric potential ![]() beta_auto . | |
const Tenseur & | get_tggg () const |
Returns the Metric potential ![]() | |
const Tenseur & | get_w_shift () const |
Returns the vector ![]() shift , following Shibata's prescription [Prog. | |
const Tenseur & | get_khi_shift () const |
Returns the scalar ![]() shift following Shibata's prescription [Prog. | |
const Tenseur_sym & | get_tkij () const |
Returns the tensor ![]() ![]() | |
const Tenseur & | get_ak_car () const |
Returns the scalar ![]() | |
virtual void | display_poly (ostream &) const |
Display in polytropic units. | |
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 ![]() | |
virtual double | mass_b () const |
Baryon mass. | |
virtual double | r_circ () const |
Circumferential radius. | |
virtual double | area () const |
Surface area. | |
virtual double | mean_radius () const |
Mean radius. | |
virtual double | aplat () const |
Flatening r_pole/r_eq. | |
virtual double | z_eqf () const |
Forward redshift factor at equator. | |
virtual double | z_eqb () const |
Backward redshift factor at equator. | |
virtual double | z_pole () const |
Redshift factor at North pole. | |
virtual double | mom_quad () const |
Quadrupole moment. | |
virtual double | r_isco (ostream *ost=0x0) const |
Circumferential radius of the innermost stable circular orbit (ISCO). | |
virtual double | f_isco () const |
Orbital frequency at the innermost stable circular orbit (ISCO). | |
virtual double | espec_isco () const |
Energy of a particle on the ISCO. | |
virtual double | lspec_isco () const |
Angular momentum of a particle on the ISCO. | |
virtual double | f_eccentric (double ecc, double periast, ostream *ost=0x0) const |
Computation of frequency of eccentric orbits. | |
virtual double | f_eq () const |
Orbital frequency at the equator. | |
virtual void | hydro_euler () |
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid frame. | |
void | update_metric () |
Computes metric coefficients from known potentials. | |
void | fait_shift () |
Computes shift from w_shift and khi_shift according to Shibata's prescription [Prog. | |
void | fait_nphi () |
Computes tnphi and nphi from the Cartesian components of the shift, stored in shift . | |
void | extrinsic_curvature () |
Computes tkij and ak_car from shift , nnn and b_car . | |
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. | |
Map & | set_mp () |
Read/write of the mapping. | |
void | set_enthalpy (const Cmp &) |
Assignment of the enthalpy field. | |
virtual void | equilibrium_spher (double ent_c, double precis=1.e-14, const Tbl *ent_limit=0x0) |
Computes a spherical static configuration. | |
void | equil_spher_regular (double ent_c, double precis=1.e-14) |
Computes a spherical static configuration. | |
virtual void | equil_spher_falloff (double ent_c, double precis=1.e-14) |
Computes a spherical static configuration with the outer boundary condition at a finite radius. | |
const Map & | get_mp () const |
Returns the mapping. | |
int | get_nzet () const |
Returns the number of domains occupied by the star. | |
bool | is_relativistic () const |
Returns true for a relativistic star, false for a Newtonian one. | |
const Eos & | get_eos () const |
Returns the equation of state. | |
const Tenseur & | get_ent () const |
Returns the enthalpy field. | |
const Tenseur & | get_nbar () const |
Returns the proper baryon density. | |
const Tenseur & | get_ener () const |
Returns the proper total energy density. | |
const Tenseur & | get_press () const |
Returns the fluid pressure. | |
const Tenseur & | get_ener_euler () const |
Returns the total energy density with respect to the Eulerian observer. | |
const Tenseur & | get_s_euler () const |
Returns the trace of the stress tensor in the Eulerian frame. | |
const Tenseur & | get_gam_euler () const |
Returns the Lorentz factor between the fluid and Eulerian observers. | |
const Tenseur & | get_u_euler () const |
Returns the fluid 3-velocity with respect to the Eulerian observer. | |
const Tenseur & | get_logn_auto () const |
Returns the logarithm of the part of the lapse N generated principaly by the star. | |
const Tenseur & | get_logn_auto_regu () const |
Returns the regular part of the logarithm of the part of the lapse N generated principaly by the star. | |
const Tenseur & | get_logn_auto_div () const |
Returns the divergent part of the logarithm of the part of the lapse N generated principaly by the star. | |
const Tenseur & | get_d_logn_auto_div () const |
Returns the gradient of logn_auto_div . | |
const Tenseur & | get_beta_auto () const |
Returns the logarithm of the part of the product AN generated principaly by the star. | |
const Tenseur & | get_nnn () const |
Returns the total lapse function N. | |
const Tenseur & | get_shift () const |
Returns the total shift vector ![]() | |
const Tenseur & | get_a_car () const |
Returns the total conformal factor ![]() | |
double | ray_eq () const |
Coordinate radius at ![]() ![]() | |
double | ray_eq (int kk) const |
Coordinate radius at ![]() ![]() | |
double | ray_eq_pis2 () const |
Coordinate radius at ![]() ![]() | |
double | ray_eq_pi () const |
Coordinate radius at ![]() ![]() | |
double | ray_eq_3pis2 () const |
Coordinate radius at ![]() ![]() | |
double | ray_pole () const |
Coordinate radius at ![]() | |
const Tbl & | xi_surf () const |
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate ![]() ![]() | |
Static Public Member Functions | |
static double | lambda_grv2 (const Cmp &sou_m, const Cmp &sou_q) |
Computes the coefficient ![]() | |
Protected Member Functions | |
virtual void | del_deriv () const |
Deletes all the derived quantities. | |
virtual 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 void | partial_display (ostream &) const |
Printing of some informations, excluding all global quantities. | |
Protected Attributes | |
bool | use_B_in_eos |
Flag : true if the value of the magnetic field is used in the Eos. | |
bool | include_magnetisation |
Flag : true if magnetisation terms are included in the equations. | |
Scalar | xmag |
The magnetisation scalar. | |
Scalar | E_I |
Interaction (magnetisation) energy density. | |
Vector | J_I |
Interaction momentum density 3-vector. | |
Sym_tensor | Sij_I |
Interaction stress 3-tensor. | |
Cmp | A_t |
t-component of the elecctromagnetic potential 1-form, divided by ![]() | |
Cmp | A_phi |
![]() ![]() | |
Cmp | B_phi |
![]() | |
Cmp | j_t |
t-component of the current 4-vector | |
Cmp | j_phi |
![]() | |
Tenseur | E_em |
electromagnetic energy density in the Eulerian frame | |
Tenseur | Jp_em |
![]() | |
Tenseur | Srr_em |
rr component of the electromagnetic stress 3-tensor, as measured in the Eulerian frame. (not used and set to 0, should be supressed) | |
Tenseur | Spp_em |
![]() | |
double | Q |
In the case of a perfect conductor, the requated baryonic charge. | |
double | a_j |
Amplitude of the curent/charge function. | |
int | conduc |
Flag: conduc=0->isolator, 1->perfect conductor. | |
double | omega |
Rotation angular velocity ( [f_unit] ) | |
Tenseur | bbb |
Metric factor B. | |
Tenseur | b_car |
Square of the metric factor B. | |
Tenseur | nphi |
Metric coefficient ![]() | |
Tenseur | tnphi |
Component ![]() | |
Tenseur | uuu |
Norm of u_euler . | |
Tenseur & | logn |
Metric potential ![]() logn_auto . | |
Tenseur | nuf |
Part of the Metric potential ![]() logn generated by the matter terms. | |
Tenseur | nuq |
Part of the Metric potential ![]() logn generated by the quadratic terms. | |
Tenseur & | dzeta |
Metric potential ![]() beta_auto . | |
Tenseur | tggg |
Metric potential ![]() | |
Tenseur | w_shift |
Vector ![]() shift , following Shibata's prescription [Prog. | |
Tenseur | khi_shift |
Scalar ![]() shift , following Shibata's prescription [Prog. | |
Tenseur_sym | tkij |
Tensor ![]() ![]() | |
Tenseur | ak_car |
Scalar ![]() | |
Cmp | ssjm1_nuf |
Effective source at the previous step for the resolution of the Poisson equation for nuf by means of Map_et::poisson . | |
Cmp | ssjm1_nuq |
Effective source at the previous step for the resolution of the Poisson equation for nuq by means of Map_et::poisson . | |
Cmp | ssjm1_dzeta |
Effective source at the previous step for the resolution of the Poisson equation for dzeta . | |
Cmp | ssjm1_tggg |
Effective source at the previous step for the resolution of the Poisson equation for tggg . | |
Cmp | ssjm1_khi |
Effective source at the previous step for the resolution of the Poisson equation for the scalar ![]() Map_et::poisson . | |
Tenseur | ssjm1_wshift |
Effective source at the previous step for the resolution of the vector Poisson equation for ![]() | |
double * | p_angu_mom |
Angular momentum. | |
double * | p_tsw |
Ratio T/W. | |
double * | p_grv2 |
Error on the virial identity GRV2. | |
double * | p_grv3 |
Error on the virial identity GRV3. | |
double * | p_r_circ |
Circumferential radius. | |
double * | p_area |
Surface area. | |
double * | p_aplat |
Flatening r_pole/r_eq. | |
double * | p_z_eqf |
Forward redshift factor at equator. | |
double * | p_z_eqb |
Backward redshift factor at equator. | |
double * | p_z_pole |
Redshift factor at North pole. | |
double * | p_mom_quad |
Quadrupole moment. | |
double * | p_mom_quad_old |
Part of the quadrupole moment. | |
double * | p_mom_quad_Bo |
Part of the quadrupole moment. | |
double * | p_r_isco |
Circumferential radius of the ISCO. | |
double * | p_f_isco |
Orbital frequency of the ISCO. | |
double * | p_espec_isco |
Specific energy of a particle on the ISCO. | |
double * | p_lspec_isco |
Specific angular momentum of a particle on the ISCO. | |
double * | p_f_eq |
Orbital frequency at the equator. | |
Map & | mp |
Mapping associated with the star. | |
int | nzet |
Number of domains of *mp occupied by the star. | |
bool | relativistic |
Indicator of relativity: true for a relativistic star, false for a Newtonian one. | |
double | unsurc2 |
![]() unsurc2=1 for a relativistic star, 0 for a Newtonian one. | |
int | k_div |
Index of regularity of the gravitational potential logn_auto . | |
const Eos & | eos |
Equation of state of the stellar matter. | |
Tenseur | ent |
Log-enthalpy (relativistic case) or specific enthalpy (Newtonian case) | |
Tenseur | nbar |
Baryon density in the fluid frame. | |
Tenseur | ener |
Total energy density in the fluid frame. | |
Tenseur | press |
Fluid pressure. | |
Tenseur | ener_euler |
Total energy density in the Eulerian frame. | |
Tenseur | s_euler |
Trace of the stress tensor in the Eulerian frame. | |
Tenseur | gam_euler |
Lorentz factor between the fluid and Eulerian observers. | |
Tenseur | u_euler |
Fluid 3-velocity with respect to the Eulerian observer. | |
Tenseur | logn_auto |
Total of the logarithm of the part of the lapse N generated principaly by the star. | |
Tenseur | logn_auto_regu |
Regular part of the logarithm of the part of the lapse N generated principaly by the star. | |
Tenseur | logn_auto_div |
Divergent part (if k_div!=0 ) of the logarithm of the part of the lapse N generated principaly by the star. | |
Tenseur | d_logn_auto_div |
Gradient of logn_auto_div (if k_div!=0 ) | |
Tenseur | beta_auto |
Logarithm of the part of the product AN generated principaly by by the star. | |
Tenseur | nnn |
Total lapse function. | |
Tenseur | shift |
Total shift vector. | |
Tenseur | a_car |
Total conformal factor ![]() | |
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 ![]() ![]() | |
double * | p_mass_b |
Baryon mass. | |
double * | p_mass_g |
Gravitational mass. | |
Definition at line 610 of file et_rot_mag.h.
Lorene::Et_magnetisation::Et_magnetisation | ( | Map & | mp_i, |
int | nzet_i, | ||
bool | relat, | ||
const Eos & | eos_i, | ||
bool | include_mag = true , |
||
bool | use_B = true |
||
) |
Standard constructor.
Definition at line 87 of file et_magnetisation.C.
References E_I, include_magnetisation, J_I, Lorene::Tensor::set_etat_zero(), Sij_I, use_B_in_eos, and xmag.
Lorene::Et_magnetisation::Et_magnetisation | ( | const Et_magnetisation & | et | ) |
Copy constructor.
Definition at line 154 of file et_magnetisation.C.
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 ) |
Definition at line 120 of file et_magnetisation.C.
References E_I, Lorene::Map::get_bvect_spher(), Lorene::Map::get_mg(), include_magnetisation, J_I, Lorene::Etoile::mp, Sij_I, use_B_in_eos, and xmag.
|
virtual |
Destructor.
Definition at line 168 of file et_magnetisation.C.
|
virtual |
Angular momentum.
Reimplemented from Lorene::Et_rot_mag.
Definition at line 582 of file et_magnetisation_comp.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::bbb, Lorene::Etoile::ener_euler, Lorene::Cmp::integrale(), J_I, Lorene::Et_rot_mag::Jp_em, Lorene::Cmp::mult_r(), Lorene::Etoile::nbar, Lorene::Etoile_rot::p_angu_mom, Lorene::Etoile::press, Lorene::Etoile::relativistic, Lorene::Cmp::std_base_scal(), Lorene::Etoile_rot::uuu, and Lorene::Cmp::va.
|
virtualinherited |
Flatening r_pole/r_eq.
Definition at line 496 of file et_rot_global.C.
References Lorene::Etoile_rot::p_aplat, Lorene::Etoile::ray_eq(), and Lorene::Etoile::ray_pole().
|
virtualinherited |
Surface area.
Definition at line 420 of file et_rot_global.C.
References Lorene::Valeur::annule_hard(), Lorene::Valeur::c_cf, Lorene::Valeur::dsdt(), Lorene::Etoile::get_a_car(), Lorene::Mg3d::get_angu(), Lorene::Etoile_rot::get_bbb(), Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Etoile_rot::l_surf(), Lorene::Etoile::mp, Lorene::Valeur::mult_st(), Lorene::Etoile_rot::p_area, Lorene::Valeur::set(), Lorene::sqrt(), Lorene::Valeur::std_base_scal(), Lorene::Cmp::va, Lorene::Valeur::val_point_jk(), Lorene::Map_radial::val_r_jk(), and Lorene::Etoile::xi_surf().
|
inline |
Public accessor to the use_B_in_eos
flag.
Definition at line 671 of file et_rot_mag.h.
References use_B_in_eos.
|
protectedvirtualinherited |
Deletes all the derived quantities.
Reimplemented from Lorene::Etoile_rot.
Definition at line 258 of file et_rot_mag.C.
References Lorene::Etoile_rot::del_deriv(), and Lorene::Et_rot_mag::set_der_0x0().
|
protectedvirtualinherited |
Sets to ETATNONDEF
(undefined state) the hydrodynamical quantities relative to the Eulerian observer.
Reimplemented from Lorene::Etoile_rot.
Definition at line 273 of file et_rot_mag.C.
References Lorene::Et_rot_mag::del_deriv(), and Lorene::Etoile_rot::del_hydro_euler().
Display in polytropic units.
Reimplemented in Lorene::Et_rot_diff.
Definition at line 690 of file etoile_rot.C.
References Lorene::Etoile_rot::angu_mom(), Lorene::Etoile::ener, Lorene::Etoile::eos, Lorene::Eos_poly::get_gam(), Lorene::Eos_poly::get_kap(), Lorene::Etoile_rot::get_omega_c(), Lorene::Etoile_rot::mass_b(), Lorene::Etoile_rot::mass_g(), Lorene::Etoile::nbar, Lorene::pow(), Lorene::Etoile_rot::r_circ(), Lorene::Etoile::ray_eq(), and Lorene::sqrt().
|
inherited |
Computes the electric field spherical components in Lorene's units.
Definition at line 151 of file et_rot_mag_global.C.
References Lorene::Etoile::a_car, Lorene::Et_rot_mag::A_phi, Lorene::Et_rot_mag::A_t, Lorene::Cmp::dsdr(), Lorene::Map::get_bvect_spher(), Lorene::Etoile::mp, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Tenseur::set(), Lorene::Valeur::set_base(), Lorene::Tenseur::set_etat_qcq(), Lorene::sqrt(), Lorene::Cmp::srdsdt(), and Lorene::Cmp::va.
|
virtual |
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
Reimplemented from Lorene::Etoile.
Definition at line 188 of file et_magnetisation.C.
References Lorene::Etoile::a_car, Lorene::Param::add_double_mod(), Lorene::Cmp::allocate_all(), Lorene::Scalar::allocate_all(), Lorene::Tensor::annule_domain(), Lorene::Tbl::annule_hard(), Lorene::Valeur::c, Lorene::Et_rot_mag::del_deriv(), Lorene::Etoile::ener, Lorene::Eos_mag::ener_ent_p(), Lorene::Etoile::ent, Lorene::Etoile::eos, Lorene::Scalar::exponential_filter_r(), Lorene::Etoile::gam_euler, Lorene::Tbl::get_etat(), Lorene::Mg3d::get_grille3d(), Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), include_magnetisation, Lorene::Eos_mag::mag_ent_p(), Lorene::Et_rot_mag::Magn(), Lorene::Etoile::mp, Lorene::Etoile::nbar, Lorene::Eos_mag::nbar_ent_p(), Lorene::Etoile::nzet, Lorene::Etoile::press, Lorene::Eos_mag::press_ent_p(), Lorene::Tenseur::set(), Lorene::Cmp::set(), Lorene::Mtbl::set(), Lorene::Scalar::set_domain(), Lorene::Valeur::set_etat_c_qcq(), Lorene::Cmp::set_etat_qcq(), Lorene::Mtbl::set_etat_qcq(), Lorene::Tbl::set_etat_qcq(), Lorene::Tenseur::set_etat_qcq(), Lorene::Cmp::set_etat_zero(), Lorene::Scalar::set_etat_zero(), Lorene::Tbl::set_etat_zero(), Lorene::Scalar::set_grid_point(), Lorene::Tenseur::set_std_base(), Lorene::sqrt(), Lorene::Cmp::std_base_scal(), Lorene::Scalar::std_spectral_base(), Lorene::Mtbl::t, use_B_in_eos, Lorene::Cmp::va, Lorene::Grille3d::x, and xmag.
Computes a spherical static configuration with the outer boundary condition at a finite radius.
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) |
Definition at line 57 of file etoile_eqsph_falloff.C.
References Lorene::Etoile::a_car, Lorene::Tenseur::annule(), Lorene::Etoile::beta_auto, Lorene::diffrel(), Lorene::Cmp::dsdr(), Lorene::Map_af::dsdr(), Lorene::Etoile::ener, Lorene::Etoile::ener_euler, Lorene::Etoile::ent, Lorene::Etoile::equation_of_state(), Lorene::exp(), Lorene::Etoile::gam_euler, Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Map_af::homothetie(), Lorene::Etoile::logn_auto, Lorene::Etoile::mass_b(), Lorene::Etoile::mass_g(), Lorene::Etoile::mp, Lorene::Etoile::nbar, Lorene::Etoile::nnn, Lorene::norme(), Lorene::Etoile::nzet, Lorene::Etoile::press, Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Tenseur::set(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_std_base(), Lorene::Etoile::shift, Lorene::sqrt(), Lorene::Etoile::u_euler, Lorene::Etoile::unsurc2, and Lorene::Map::val_r().
Computes a spherical static configuration.
The sources for Poisson equations are regularized by extracting analytical diverging parts.
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) |
Definition at line 115 of file et_equil_spher_regu.C.
References Lorene::Etoile::a_car, Lorene::Tenseur::annule(), Lorene::Etoile::beta_auto, Lorene::Etoile::d_logn_auto_div, Lorene::Eos::der_ener_ent_p(), Lorene::Eos::der_nbar_ent_p(), Lorene::diffrel(), Lorene::Map_af::dsdr(), Lorene::Etoile::ener, Lorene::Etoile::ener_euler, Lorene::Etoile::ent, Lorene::Etoile::eos, Lorene::Etoile::equation_of_state(), Lorene::exp(), Lorene::Etoile::gam_euler, Lorene::Map::get_bvect_spher(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Tenseur::gradient_spher(), Lorene::Map_af::homothetie(), Lorene::Etoile::k_div, Lorene::Etoile::logn_auto, Lorene::Etoile::logn_auto_div, Lorene::Etoile::logn_auto_regu, Lorene::Etoile::mass_b(), Lorene::Etoile::mass_g(), Lorene::Etoile::mp, Lorene::Etoile::nbar, Lorene::Etoile::nnn, Lorene::norme(), Lorene::Etoile::nzet, Lorene::Map_af::poisson(), Lorene::Map_af::poisson_regular(), Lorene::Etoile::press, Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Tenseur::set(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_std_base(), Lorene::Tenseur::set_triad(), Lorene::Etoile::shift, Lorene::sqrt(), Lorene::Cmp::std_base_scal(), Lorene::Etoile::u_euler, Lorene::Etoile::unsurc2, and Lorene::Map::val_r().
|
virtualinherited |
Computes an equilibrium configuration.
ent_c | [input] Central enthalpy |
omega0 | [input] Requested angular velocity (if fact_omega=1 . ) |
fact_omega | [input] 1.01 = search for the Keplerian frequency,
|
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 in Lorene::Et_rot_diff.
Definition at line 147 of file et_rot_equilibrium.C.
References Lorene::Etoile::a_car, Lorene::abs(), Lorene::Map::adapt(), Lorene::Param::add_cmp_mod(), Lorene::Param::add_double(), Lorene::Param::add_double_mod(), Lorene::Param::add_int(), Lorene::Param::add_int_mod(), Lorene::Param::add_tbl(), Lorene::Param::add_tenseur_mod(), Lorene::Etoile_rot::ak_car, Lorene::Cmp::annule(), Lorene::Etoile_rot::bbb, Lorene::Valeur::c_cf, Lorene::Tenseur::change_triad(), Lorene::Map::cmp_zero(), Lorene::Valeur::coef(), Lorene::cos(), Lorene::diffrel(), Lorene::Etoile_rot::dzeta, Lorene::Etoile::ener_euler, Lorene::Etoile::ent, Lorene::Etoile::equation_of_state(), Lorene::Etoile_rot::fait_nphi(), Lorene::flat_scalar_prod(), Lorene::Etoile::gam_euler, Lorene::Map::get_bvect_cart(), Lorene::Tenseur::get_etat(), Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Mg3d::get_type_t(), Lorene::Tenseur::gradient_spher(), Lorene::Etoile_rot::grv2(), Lorene::Map_et::homothetie(), Lorene::Etoile_rot::hydro_euler(), Lorene::Etoile_rot::khi_shift, Lorene::log(), Lorene::log10(), Lorene::Etoile_rot::logn, Lorene::Etoile_rot::mass_b(), Lorene::Etoile_rot::mass_g(), Lorene::Etoile::mp, Lorene::Cmp::mult_rsint(), Lorene::Etoile::nbar, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile_rot::nuf, Lorene::Etoile_rot::nuq, Lorene::Etoile::nzet, Lorene::Etoile_rot::omega, Lorene::Etoile_rot::partial_display(), Lorene::Map::phi, Lorene::Map::poisson2d(), Lorene::Tenseur::poisson_vect(), Lorene::pow(), Lorene::Etoile::press, Lorene::Etoile::ray_eq(), Lorene::Etoile::ray_pole(), Lorene::Map::reevaluate(), Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Tenseur::set(), Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_std_base(), Lorene::Etoile::shift, Lorene::Map::sint, Lorene::sqrt(), Lorene::Etoile_rot::ssjm1_khi, Lorene::Etoile_rot::ssjm1_nuf, Lorene::Etoile_rot::ssjm1_nuq, Lorene::Etoile_rot::ssjm1_tggg, Lorene::Etoile_rot::ssjm1_wshift, Lorene::Cmp::std_base_scal(), Lorene::Etoile_rot::tggg, Lorene::Etoile_rot::tkij, Lorene::Etoile::u_euler, Lorene::Etoile_rot::update_metric(), Lorene::Etoile_rot::uuu, Lorene::Cmp::va, and Lorene::Etoile_rot::w_shift.
|
inherited |
Computes an equilibrium configuration.
ent_c | [input] Central enthalpy |
omega0 | [input] Requested angular velocity (if fact_omega=1 . ) |
fact_omega | [input] 1.01 = search for the Keplerian frequency,
|
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 1 for the storage of some error indicators :
|
Q0 | [input] Requested electric charge for the case of a perfect conductor. Charge per baryon for the case of an isolator. |
a_j0 | [input] Amplitude for the current/charge coupling function |
f_j | [input] current or charge coupling function (see Bocquet et al. 1995). |
M_j | [input] primitive (null for zero) of current/charge coupling function (see Bocquet et al. 1995) used for the first integral of stationary motion. |
Definition at line 114 of file et_rot_mag_equil.C.
References Lorene::Etoile::a_car, Lorene::Et_rot_mag::a_j, Lorene::Et_rot_mag::A_phi, Lorene::Et_rot_mag::A_t, Lorene::abs(), Lorene::Map::adapt(), Lorene::Param::add_cmp_mod(), Lorene::Param::add_double(), Lorene::Param::add_double_mod(), Lorene::Param::add_int(), Lorene::Param::add_int_mod(), Lorene::Param::add_tbl(), Lorene::Param::add_tenseur_mod(), Lorene::Etoile_rot::ak_car, Lorene::Cmp::annule(), Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::bbb, Lorene::Tenseur::change_triad(), Lorene::Map::cmp_zero(), Lorene::Valeur::coef_i(), Lorene::diffrel(), Lorene::Etoile_rot::dzeta, Lorene::Et_rot_mag::E_em, Lorene::Etoile::ener_euler, Lorene::Etoile::ent, Lorene::Etoile::equation_of_state(), Lorene::Etoile_rot::fait_nphi(), Lorene::flat_scalar_prod(), Lorene::Etoile::gam_euler, Lorene::Map::get_bvect_cart(), Lorene::Map::get_bvect_spher(), Lorene::Cmp::get_etat(), Lorene::Tenseur::get_etat(), Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Tenseur::gradient_spher(), Lorene::Et_rot_mag::grv2(), Lorene::Map_et::homothetie(), Lorene::Etoile_rot::hydro_euler(), Lorene::Et_rot_mag::is_conduct(), Lorene::Et_rot_mag::Jp_em, Lorene::Etoile_rot::khi_shift, Lorene::log(), Lorene::log10(), Lorene::Etoile_rot::logn, Lorene::Et_rot_mag::magnet_comput(), Lorene::Etoile_rot::mass_b(), Lorene::Et_rot_mag::mass_g(), Lorene::Et_rot_mag::MHD_comput(), Lorene::Etoile::mp, Lorene::Cmp::mult_rsint(), Lorene::Etoile::nbar, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile_rot::nuf, Lorene::Etoile_rot::nuq, Lorene::Etoile::nzet, Lorene::Etoile_rot::omega, Lorene::Map::poisson2d(), Lorene::Tenseur::poisson_vect(), Lorene::pow(), Lorene::Etoile::press, Lorene::Et_rot_mag::Q, Lorene::Etoile::ray_eq(), Lorene::Etoile::ray_pole(), Lorene::Map::reevaluate(), Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Tenseur::set(), Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), Lorene::Tenseur::set_etat_qcq(), Lorene::Cmp::set_etat_zero(), Lorene::Tenseur::set_etat_zero(), Lorene::Tenseur::set_std_base(), Lorene::Etoile::shift, Lorene::sqrt(), Lorene::Etoile_rot::ssjm1_khi, Lorene::Etoile_rot::ssjm1_nuf, Lorene::Etoile_rot::ssjm1_nuq, Lorene::Etoile_rot::ssjm1_tggg, Lorene::Etoile_rot::ssjm1_wshift, Lorene::Cmp::std_base_scal(), Lorene::Etoile_rot::tggg, Lorene::Etoile_rot::tkij, Lorene::Etoile::u_euler, Lorene::Etoile_rot::update_metric(), Lorene::Etoile_rot::uuu, Lorene::Cmp::va, and Lorene::Etoile_rot::w_shift.
void Lorene::Et_magnetisation::equilibrium_mag | ( | double | ent_c, |
double | omega0, | ||
double | fact_omega, | ||
int | nzadapt, | ||
const Tbl & | ent_limit, | ||
const Itbl & | icontrol, | ||
const Tbl & | control, | ||
double | mbar_wanted, | ||
double | magmom_wanted, | ||
double | aexp_mass, | ||
Tbl & | diff, | ||
double | Q0, | ||
double | a_j0, | ||
Cmp(*)(const Cmp &x, const double) | f_j, | ||
Cmp(*)(const Cmp &x, const double) | M_j | ||
) |
Computes an equilibrium configuration.
ent_c | [input] Central enthalpy |
omega0 | [input] Requested angular velocity (if fact_omega=1 . ) |
fact_omega | [input] 1.01 = search for the Keplerian frequency,
|
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 11) to control the iteration:
|
control | [input] Set of parameters (stored as a 1-D Tbl of size 8) 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 1 for the storage of some error indicators :
|
Q0 | [input] Requested electric charge for the case of a perfect conductor. Charge per baryon for the case of an isolator. |
a_j0 | [input] Amplitude for the current/charge coupling function |
f_j | [input] current or charge coupling function (see Bocquet et al. 1995). |
M_j | [input] primitive (null for zero) of current/charge coupling function (see Bocquet et al. 1995) used for the first integral of stationary motion. |
Definition at line 367 of file et_magnetisation.C.
References Lorene::Etoile::a_car, Lorene::Et_rot_mag::a_j, Lorene::Et_rot_mag::A_phi, Lorene::Et_rot_mag::A_t, Lorene::abs(), Lorene::Map::adapt(), Lorene::Param::add_cmp_mod(), Lorene::Param::add_double(), Lorene::Param::add_double_mod(), Lorene::Param::add_int(), Lorene::Param::add_int_mod(), Lorene::Param::add_tbl(), Lorene::Param::add_tenseur_mod(), Lorene::Etoile_rot::ak_car, Lorene::Cmp::annule(), Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::bbb, Lorene::Tenseur::change_triad(), Lorene::Map::cmp_zero(), Lorene::Valeur::coef_i(), Lorene::diffrel(), Lorene::Etoile_rot::dzeta, Lorene::Et_rot_mag::E_em, E_I, Lorene::Etoile::ener_euler, Lorene::Etoile::ent, equation_of_state(), Lorene::Etoile_rot::fait_nphi(), Lorene::flat_scalar_prod(), Lorene::Etoile::gam_euler, Lorene::Map::get_bvect_cart(), Lorene::Map::get_bvect_spher(), Lorene::Cmp::get_etat(), Lorene::Tenseur::get_etat(), Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Mg3d::get_type_t(), Lorene::Tenseur::gradient_spher(), grv2(), Lorene::Map_et::homothetie(), Lorene::Etoile_rot::hydro_euler(), Lorene::Et_rot_mag::is_conduct(), J_I, Lorene::Et_rot_mag::Jp_em, Lorene::Etoile_rot::khi_shift, Lorene::log(), Lorene::log10(), Lorene::Etoile_rot::logn, Lorene::Et_rot_mag::MagMom(), magnet_comput(), Lorene::Etoile_rot::mass_b(), mass_g(), MHD_comput(), Lorene::Etoile::mp, Lorene::Cmp::mult_rsint(), Lorene::Etoile::nbar, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile_rot::nuf, Lorene::Etoile_rot::nuq, Lorene::Etoile::nzet, Lorene::Etoile_rot::omega, Lorene::Map::poisson2d(), Lorene::Tenseur::poisson_vect(), Lorene::pow(), Lorene::Etoile::press, Lorene::Et_rot_mag::Q, Lorene::Etoile::ray_eq(), Lorene::Etoile::ray_pole(), Lorene::Map::reevaluate(), Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Tenseur::set(), Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), Lorene::Tenseur::set_etat_qcq(), Lorene::Cmp::set_etat_zero(), Lorene::Tenseur::set_etat_zero(), Lorene::Tenseur::set_std_base(), Lorene::Etoile::shift, Sij_I, Lorene::sqrt(), Lorene::Etoile_rot::ssjm1_khi, Lorene::Etoile_rot::ssjm1_nuf, Lorene::Etoile_rot::ssjm1_nuq, Lorene::Etoile_rot::ssjm1_tggg, Lorene::Etoile_rot::ssjm1_wshift, Lorene::Cmp::std_base_scal(), Lorene::Etoile_rot::tggg, Lorene::Etoile_rot::tkij, Lorene::Etoile::u_euler, Lorene::Etoile_rot::update_metric(), Lorene::Etoile_rot::uuu, Lorene::Cmp::va, and Lorene::Etoile_rot::w_shift.
|
inherited |
Computes an equilibrium configuration.
ent_c | [input] Central enthalpy |
omega0 | [input] Requested angular velocity (if fact_omega=1 . ) |
fact_omega | [input] 1.01 = search for the Keplerian frequency,
|
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 1 for the storage of some error indicators :
|
Q0 | [input] Requested electric charge for the case of a perfect conductor. Charge per baryon for the case of an isolator. |
a_j0 | [input] Amplitude for the current/charge coupling function |
f_j | [input] current or charge coupling function (see Bocquet et al. 1995). |
M_j | [input] primitive (null for zero) of current/charge coupling function (see Bocquet et al. 1995) used for the first integral of stationary motion. |
Definition at line 62 of file et_rot_mag_equil_plus.C.
References Lorene::Et_rot_mag::a_j, Lorene::Et_rot_mag::A_phi, Lorene::Et_rot_mag::A_t, Lorene::abs(), Lorene::Param::add_cmp_mod(), Lorene::Param::add_double(), Lorene::Param::add_int(), Lorene::Param::add_int_mod(), Lorene::Param::add_tenseur_mod(), Lorene::Etoile::equation_of_state(), Lorene::Etoile_rot::hydro_euler(), Lorene::Et_rot_mag::j_phi, Lorene::Et_rot_mag::magnet_comput_plus(), Lorene::max(), Lorene::Et_rot_mag::MHD_comput(), Lorene::min(), Lorene::Etoile::mp, Lorene::Cmp::multipole_spectrum(), Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), Lorene::Cmp::set_etat_zero(), Lorene::Etoile_rot::ssjm1_khi, Lorene::Etoile_rot::ssjm1_wshift, Lorene::Cmp::std_base_scal(), and Lorene::Etoile_rot::update_metric().
|
virtualinherited |
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 87 of file etoile_equil_spher.C.
References Lorene::Etoile::a_car, Lorene::Map_et::adapt(), Lorene::Param::add_double(), Lorene::Param::add_int(), Lorene::Param::add_int_mod(), Lorene::Param::add_tbl(), Lorene::Tenseur::annule(), Lorene::Etoile::beta_auto, Lorene::diffrel(), Lorene::Cmp::dsdr(), Lorene::Map_af::dsdr(), Lorene::Etoile::ener, Lorene::Etoile::ener_euler, Lorene::Etoile::ent, Lorene::Etoile::equation_of_state(), Lorene::exp(), Lorene::Etoile::gam_euler, Lorene::Map_af::get_alpha(), Lorene::Map_et::get_alpha(), Lorene::Map_af::get_beta(), Lorene::Map_et::get_beta(), Lorene::Etoile::get_ent(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Etoile::get_press(), Lorene::Map_af::homothetie(), Lorene::Etoile::logn_auto, Lorene::Etoile::mass_b(), Lorene::Etoile::mass_g(), Lorene::Etoile::mp, Lorene::Etoile::nbar, Lorene::Etoile::nnn, Lorene::norme(), Lorene::Etoile::nzet, Lorene::Map_af::poisson(), Lorene::Etoile::press, Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Tenseur::set(), Lorene::Map_af::set_alpha(), Lorene::Map_af::set_beta(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_std_base(), Lorene::Etoile::shift, Lorene::sqrt(), Lorene::Etoile::u_euler, Lorene::Etoile::unsurc2, and Lorene::Map::val_r().
|
virtualinherited |
Energy of a particle on the ISCO.
Definition at line 301 of file et_rot_isco.C.
References Lorene::Etoile_rot::p_espec_isco, and Lorene::Etoile_rot::r_isco().
|
inherited |
Computes tkij
and ak_car
from shift
, nnn
and b_car
.
Definition at line 57 of file et_rot_extr_curv.C.
References Lorene::Etoile_rot::ak_car, Lorene::Etoile_rot::b_car, Lorene::contract(), Lorene::Cmp::get_etat(), Lorene::Tenseur::get_etat(), Lorene::Map::get_mg(), Lorene::Tenseur::gradient(), Lorene::Etoile::mp, Lorene::Cmp::mult_rsint(), Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Tenseur::set(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_etat_zero(), Lorene::Etoile::shift, Lorene::Etoile_rot::tkij, and Lorene::Cmp::va.
|
virtualinherited |
Computation of frequency of eccentric orbits.
ecc | eccentricity of the orbit |
periasrt | periastron of the orbit |
ost | output stream to give details of the computation; if set to 0x0 [default value], no details will be given. |
Definition at line 78 of file et_rot_f_eccentric.C.
References Lorene::Param::add_cmp(), Lorene::Param::add_int(), Lorene::Cmp::annule(), Lorene::Etoile_rot::bbb, Lorene::Cmp::dsdr(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), Lorene::Etoile::mp, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile::nzet, Lorene::Etoile_rot::p_f_isco, Lorene::Etoile_rot::p_r_isco, Lorene::Map::r, Lorene::Etoile::ray_eq(), Lorene::sqrt(), Lorene::Cmp::std_base_scal(), Lorene::Cmp::va, Lorene::Valeur::val_point(), and Lorene::Map::val_r().
|
virtualinherited |
Orbital frequency at the equator.
Definition at line 319 of file et_rot_isco.C.
References Lorene::Etoile_rot::p_f_eq, and Lorene::Etoile_rot::r_isco().
|
virtualinherited |
Orbital frequency at the innermost stable circular orbit (ISCO).
Definition at line 267 of file et_rot_isco.C.
References Lorene::Etoile_rot::p_f_isco, and Lorene::Etoile_rot::r_isco().
|
inherited |
Computes tnphi
and nphi
from the Cartesian components of the shift, stored in shift
.
Definition at line 781 of file etoile_rot.C.
References Lorene::Map::comp_p_from_cartesian(), Lorene::Tenseur::get_etat(), Lorene::Etoile::mp, Lorene::Etoile_rot::nphi, Lorene::Tenseur::set(), Lorene::Tenseur::set_etat_qcq(), Lorene::Etoile::shift, and Lorene::Etoile_rot::tnphi.
|
inherited |
Computes shift
from w_shift
and khi_shift
according to Shibata's prescription [Prog.
Theor. Phys. 101 , 1199 (1999)] :
Definition at line 748 of file etoile_rot.C.
References Lorene::Tenseur::dec2_dzpuis(), Lorene::Tenseur::dec_dzpuis(), Lorene::Tenseur::get_etat(), Lorene::Tenseur::get_triad(), Lorene::Tenseur::gradient(), Lorene::Etoile_rot::khi_shift, Lorene::Tenseur::set(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_triad(), Lorene::Etoile::shift, Lorene::skxk(), and Lorene::Etoile_rot::w_shift.
Returns the total conformal factor
Definition at line 733 of file etoile.h.
References Lorene::Etoile::a_car.
|
inlineinherited |
Returns the amplitude of the current/charge function.
Definition at line 284 of file et_rot_mag.h.
References Lorene::Et_rot_mag::a_j.
Returns the scalar
For axisymmetric stars, this quantity is related to the derivatives of
In particular it is related to the quantities
Definition at line 1799 of file etoile.h.
References Lorene::Etoile_rot::ak_car.
Returns the
Definition at line 251 of file et_rot_mag.h.
References Lorene::Et_rot_mag::A_phi.
Returns the t component of the electromagnetic potential, divided by
Definition at line 246 of file et_rot_mag.h.
References Lorene::Et_rot_mag::A_t.
Returns the square of the metric factor B.
Definition at line 1715 of file etoile.h.
References Lorene::Etoile_rot::b_car.
Returns the metric factor B.
Definition at line 1712 of file etoile.h.
References Lorene::Etoile_rot::bbb.
Returns the logarithm of the part of the product AN generated principaly by the star.
Definition at line 724 of file etoile.h.
References Lorene::Etoile::beta_auto.
Returns the
Definition at line 254 of file et_rot_mag.h.
References Lorene::Et_rot_mag::B_phi.
Returns the gradient of logn_auto_div
.
Definition at line 719 of file etoile.h.
References Lorene::Etoile::d_logn_auto_div.
Returns the Metric potential beta_auto
.
Definition at line 1742 of file etoile.h.
References Lorene::Etoile_rot::dzeta.
Accessor to the interaction energy density.
Definition at line 680 of file et_rot_mag.h.
References E_I.
Returns the electromagnetic energy density in the Eulerian frame.
Definition at line 260 of file et_rot_mag.h.
References Lorene::Et_rot_mag::E_em.
Returns the proper total energy density.
Definition at line 679 of file etoile.h.
References Lorene::Etoile::ener.
Returns the total energy density with respect to the Eulerian observer.
Definition at line 685 of file etoile.h.
References Lorene::Etoile::ener_euler.
Returns the enthalpy field.
Definition at line 673 of file etoile.h.
References Lorene::Etoile::ent.
Returns the equation of state.
Definition at line 670 of file etoile.h.
References Lorene::Etoile::eos.
Returns the Lorentz factor between the fluid and Eulerian observers.
Definition at line 691 of file etoile.h.
References Lorene::Etoile::gam_euler.
Accessor to the interaction momentum vector.
Definition at line 683 of file et_rot_mag.h.
References J_I.
Returns the
Definition at line 265 of file et_rot_mag.h.
References Lorene::Et_rot_mag::Jp_em.
Returns the
Definition at line 258 of file et_rot_mag.h.
References Lorene::Et_rot_mag::j_phi.
Returns the t component of the current 4-vector.
Definition at line 256 of file et_rot_mag.h.
References Lorene::Et_rot_mag::j_t.
Returns the scalar shift
following Shibata's prescription [Prog.
Theor. Phys. 101 , 1199 (1999)] :
NB: w_shift
contains the components of mp
.
Definition at line 1773 of file etoile.h.
References Lorene::Etoile_rot::khi_shift.
Returns the metric potential logn_auto
.
Definition at line 1729 of file etoile.h.
References Lorene::Etoile_rot::logn.
Returns the logarithm of the part of the lapse N generated principaly by the star.
In the Newtonian case, this is the Newtonian gravitational potential (in units of
Definition at line 701 of file etoile.h.
References Lorene::Etoile::logn_auto.
Returns the divergent part of the logarithm of the part of the lapse N generated principaly by the star.
In the Newtonian case, this is the diverging part of the Newtonian gravitational potential (in units of
Definition at line 715 of file etoile.h.
References Lorene::Etoile::logn_auto_div.
Returns the regular part of the logarithm of the part of the lapse N generated principaly by the star.
In the Newtonian case, this is the Newtonian gravitational potential (in units of
Definition at line 708 of file etoile.h.
References Lorene::Etoile::logn_auto_regu.
Accessor to the magnetisation scalar field.
Definition at line 677 of file et_rot_mag.h.
References xmag.
Returns the proper baryon density.
Definition at line 676 of file etoile.h.
References Lorene::Etoile::nbar.
Returns the total lapse function N.
Definition at line 727 of file etoile.h.
References Lorene::Etoile::nnn.
Returns the metric coefficient
Definition at line 1718 of file etoile.h.
References Lorene::Etoile_rot::nphi.
Returns the part of the Metric potential logn
generated by the matter terms.
Definition at line 1734 of file etoile.h.
References Lorene::Etoile_rot::nuf.
Returns the Part of the Metric potential logn
generated by the quadratic terms.
Definition at line 1739 of file etoile.h.
References Lorene::Etoile_rot::nuq.
|
inlineinherited |
Returns the number of domains occupied by the star.
Definition at line 662 of file etoile.h.
References Lorene::Etoile::nzet.
|
virtualinherited |
Returns the central value of the rotation angular velocity ([f_unit] )
Reimplemented in Lorene::Et_rot_diff.
Definition at line 680 of file etoile_rot.C.
References Lorene::Etoile_rot::omega.
Returns the fluid pressure.
Definition at line 682 of file etoile.h.
References Lorene::Etoile::press.
|
inlineinherited |
Returns the requested electric charge in the case of a perfect conductor and the charge/baryon for an isolator.
Definition at line 282 of file et_rot_mag.h.
References Lorene::Et_rot_mag::Q.
Returns the trace of the stress tensor in the Eulerian frame.
Definition at line 688 of file etoile.h.
References Lorene::Etoile::s_euler.
Returns the total shift vector
Definition at line 730 of file etoile.h.
References Lorene::Etoile::shift.
|
inline |
Accessor to the interaction stress tensor.
Definition at line 686 of file et_rot_mag.h.
References Sij_I.
Returns the
Definition at line 276 of file et_rot_mag.h.
References Lorene::Et_rot_mag::Spp_em.
Returns the rr-component of the electromagnetic stress 3-tensor, as measured in the Eulerian frame.
(not used and always equal to 0, should be supressed)
Definition at line 271 of file et_rot_mag.h.
References Lorene::Et_rot_mag::Srr_em.
Returns the Metric potential
Definition at line 1745 of file etoile.h.
References Lorene::Etoile_rot::tggg.
|
inlineinherited |
Returns the tensor
tkij
contains the Cartesian components of
Definition at line 1780 of file etoile.h.
References Lorene::Etoile_rot::tkij.
Returns the component
Definition at line 1723 of file etoile.h.
References Lorene::Etoile_rot::tnphi.
Returns the fluid 3-velocity with respect to the Eulerian observer.
Definition at line 694 of file etoile.h.
References Lorene::Etoile::u_euler.
Returns the norm of u_euler
.
Definition at line 1726 of file etoile.h.
References Lorene::Etoile_rot::uuu.
Returns the vector shift
, following Shibata's prescription [Prog.
Theor. Phys. 101 , 1199 (1999)] :
NB: w_shift
contains the components of mp
.
Definition at line 1759 of file etoile.h.
References Lorene::Etoile_rot::w_shift.
|
virtual |
Error on the virial identity GRV2.
Reimplemented from Lorene::Et_rot_mag.
Definition at line 613 of file et_magnetisation_comp.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::ak_car, Lorene::Etoile_rot::b_car, Lorene::Etoile::ener_euler, Lorene::flat_scalar_prod(), Lorene::Tenseur::gradient_spher(), Lorene::Etoile_rot::lambda_grv2(), Lorene::Etoile_rot::logn, Lorene::Etoile_rot::p_grv2, Lorene::Etoile::press, Sij_I, Lorene::Et_rot_mag::Spp_em, and Lorene::Etoile_rot::uuu.
Error on the virial identity GRV3.
The error is computed as the integral defined by Eq. (43) of [Gourgoulhon and Bonazzola, Class. Quantum Grav. 11 , 443 (1994)] divided by the integral of the matter terms.
ost | output stream to give details of the computation; if set to 0x0 [default value], no details will be given. |
Reimplemented from Lorene::Et_rot_mag.
Definition at line 642 of file et_magnetisation_comp.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::ak_car, Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::bbb, Lorene::Cmp::dsdr(), Lorene::Etoile_rot::dzeta, Lorene::flat_scalar_prod(), Lorene::Cmp::get_dzpuis(), Lorene::Cmp::get_etat(), Lorene::Tenseur::gradient_spher(), Lorene::log(), Lorene::Etoile_rot::logn, Lorene::Etoile::mp, Lorene::Valeur::mult_ct(), Lorene::Etoile::nbar, Lorene::Etoile_rot::p_grv3, Lorene::Etoile::press, Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Cmp::set_dzpuis(), Lorene::Tenseur::set_std_base(), Sij_I, Lorene::Et_rot_mag::Spp_em, Lorene::Cmp::srdsdt(), Lorene::Valeur::ssint(), Lorene::Cmp::std_base_scal(), Lorene::Valeur::sx(), Lorene::Etoile_rot::uuu, Lorene::Cmp::va, and Lorene::Map_radial::xsr.
|
inherited |
Gyromagnetic ratio
Definition at line 268 of file et_rot_mag_global.C.
References Lorene::Et_rot_mag::angu_mom(), Lorene::Et_rot_mag::MagMom(), Lorene::Et_rot_mag::mass_g(), and Lorene::Et_rot_mag::Q_comput().
|
virtualinherited |
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid frame.
The calculation is performed starting from the quantities 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::Etoile.
Reimplemented in Lorene::Et_rot_bifluid, and Lorene::Et_rot_diff.
Definition at line 83 of file et_rot_hydro.C.
References Lorene::Etoile::a_car, Lorene::Tenseur::annule(), Lorene::Etoile_rot::b_car, Lorene::Tenseur::change_triad(), Lorene::Etoile_rot::del_deriv(), Lorene::Etoile::ener, Lorene::Etoile::ener_euler, Lorene::Etoile::gam_euler, Lorene::Map::get_bvect_cart(), Lorene::Map::get_bvect_spher(), Lorene::Tenseur::get_etat(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), Lorene::Etoile::mp, Lorene::Etoile::nnn, Lorene::Etoile_rot::omega, Lorene::Etoile::press, Lorene::Etoile::s_euler, Lorene::Tenseur::set(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_std_base(), Lorene::Tenseur::set_triad(), Lorene::Etoile::shift, Lorene::sqrt(), Lorene::Etoile::u_euler, Lorene::Etoile::unsurc2, Lorene::Etoile_rot::uuu, Lorene::Map::x, and Lorene::Map::y.
|
inlineinherited |
Tells if the star is made of conducting or isolating material.
Definition at line 241 of file et_rot_mag.h.
References Lorene::Et_rot_mag::conduc.
|
inlineinherited |
Returns true
for a relativistic star, false
for a Newtonian one.
Definition at line 667 of file etoile.h.
References Lorene::Etoile::relativistic.
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 from Lorene::Etoile.
Reimplemented in Lorene::Et_rot_bifluid.
Definition at line 95 of file et_rot_global.C.
References Lorene::Etoile::ent, Lorene::Map::get_mg(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Etoile::mp, Lorene::Etoile::nzet, Lorene::Etoile::p_l_surf, and Lorene::Etoile::p_xi_surf.
Computes the coefficient
in order that the total source does not contain any monopolar term, i.e. in order that
where
Then, by construction, the new source
sou_m | [input] matter source term ![]() |
sou_q | [input] quadratic source term ![]() |
Definition at line 79 of file et_rot_lambda_grv2.C.
References Lorene::Valeur::c, Lorene::Cmp::check_dzpuis(), Lorene::Valeur::coef_i(), Lorene::Map_radial::dxdr, Lorene::Map_af::get_alpha(), Lorene::Map_af::get_beta(), Lorene::Cmp::get_etat(), Lorene::Tbl::get_etat(), Lorene::Valeur::get_etat(), Lorene::Mg3d::get_grille3d(), Lorene::Map::get_mg(), Lorene::Cmp::get_mp(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Mg3d::get_type_r(), Lorene::Map_af::set_alpha(), Lorene::Map_af::set_beta(), Lorene::Mtbl::t, Lorene::Tbl::t, Lorene::Cmp::va, Lorene::Map::val_r(), Lorene::Grille3d::x, and Lorene::Map_radial::xsr.
|
virtualinherited |
Angular momentum of a particle on the ISCO.
Definition at line 284 of file et_rot_isco.C.
References Lorene::Etoile_rot::p_lspec_isco, and Lorene::Etoile_rot::r_isco().
|
inherited |
Magnetic Momentum
Definition at line 192 of file et_rot_mag_global.C.
References Lorene::Et_rot_mag::A_phi, Lorene::Cmp::asymptot(), Lorene::Cmp::get_etat(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Etoile::mp, and Lorene::pow().
|
inherited |
Computes the magnetic field spherical components in Lorene's units.
Definition at line 170 of file et_rot_mag_global.C.
References Lorene::Etoile::a_car, Lorene::Et_rot_mag::A_phi, Lorene::Et_rot_mag::B_phi, Lorene::Etoile_rot::bbb, Lorene::Cmp::dsdr(), Lorene::Map::get_bvect_spher(), Lorene::Etoile::mp, Lorene::Tenseur::set(), Lorene::Valeur::set_base(), Lorene::Tenseur::set_etat_qcq(), Lorene::sqrt(), Lorene::Cmp::srdsdt(), and Lorene::Cmp::va.
|
virtual |
Computes the electromagnetic quantities solving the Maxwell equations (6) and (7) of [Bocquet, Bonazzola, Gourgoulhon and Novak, Astron.
Astrophys. 301 , 757 (1995)]. In the case of a perfect conductor, le electromagnetic potential may have a discontinuous derivative across star's surface.
conduc | [input] flag: 0 for an isolator, 1 for a perfect conductor |
adapt_flag | [input] flag: if 0 the mapping is NOT adapted to star's surface |
f_j | [input] current or charge coupling function (see Bocquet et al. 1995). |
par_poisson_At | [input] parameters for controlling the solution of the Poisson equation for At potential (see file et_rot_mag_equil.C) |
par_poisson_Avect | [input] parameters for controlling the solution of vector Poisson equation for magnetic potential (see file et_rot_mag_equil.C) |
Definition at line 97 of file et_magnetisation_comp.C.
References Lorene::Etoile::a_car, Lorene::Et_rot_mag::a_j, Lorene::Et_rot_mag::A_phi, Lorene::Et_rot_mag::A_t, Lorene::Cmp::allocate_all(), Lorene::Cmp::annule(), Lorene::Cmp::asymptot(), Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::bbb, Lorene::Coord::c, Lorene::Tenseur::change_triad(), Lorene::cos(), Lorene::Cmp::dec_dzpuis(), Lorene::Cmp::div_r(), Lorene::Cmp::dsdr(), Lorene::Etoile::ener, Lorene::Scalar::exponential_filter_r(), Lorene::Scalar::exponential_filter_ylm(), Lorene::Coord::fait(), Lorene::flat_scalar_prod_desal(), Lorene::Map::get_bvect_cart(), Lorene::Map::get_bvect_spher(), Lorene::Cmp::get_dzpuis(), Lorene::Tenseur::get_etat(), get_magnetisation(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Tenseur::gradient_spher(), Lorene::Cmp::inc2_dzpuis(), Lorene::Et_rot_mag::j_phi, Lorene::Et_rot_mag::j_t, Lorene::Etoile_rot::l_surf(), Lorene::Cmp::laplacien(), Lorene::log(), Lorene::Etoile_rot::logn, Lorene::Etoile::mp, Lorene::Valeur::mult_ct(), Lorene::Cmp::mult_rsint(), Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile::nzet, Lorene::Etoile_rot::omega, Lorene::Cmp::poisson(), Lorene::Tenseur::poisson_vect(), Lorene::pow(), Lorene::Etoile::press, Lorene::Et_rot_mag::Q, Lorene::Map::r, Lorene::Tenseur::set(), Lorene::Tbl::set(), Lorene::Cmp::set(), Lorene::Tbl::set_etat_qcq(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_std_base(), Lorene::Cmp::srdsdt(), Lorene::Valeur::ssint(), Lorene::Cmp::std_base_scal(), Lorene::Tbl::t, Lorene::Map::tet, Lorene::Etoile_rot::tnphi, Lorene::Cmp::va, Lorene::Valeur::val_point_jk(), Lorene::Map_radial::val_r_jk(), and Lorene::Etoile::xi_surf().
|
virtualinherited |
Computes the electromagnetic quantities solving the Maxwell equations (6) and (7) of [Bocquet, Bonazzola, Gourgoulhon and Novak, Astron.
Astrophys. 301 , 757 (1995)]. In the case of a perfect conductor, le electromagnetic potential may have a discontinuous derivative across star's surface.
adapt_flag | [input] flag: if 0 the mapping is NOT adapted to star's surface |
initial_j | [input] flag: initial current for the iteration: 0= no current, 1=dipolar-like current , 2= quadrupolar-like current |
a_j0 | [input] amplitude of the non-force free current |
f_j | [input] current coupling function (non-FF part) (see Bocquet et al. 1995). |
b_j0 | [input] amplitude of the force free current |
g_j | [input] current coupling function (FF-part) |
N_j | [input] current coupling function (FF-part) |
par_poisson_At | [input] parameters for controlling the solution of the Poisson equation for At potential (see file et_rot_mag_equil.C) |
par_poisson_Avect | [input] parameters for controlling the solution of vector Poisson equation for magnetic potential (see file et_rot_mag_equil.C) |
Definition at line 111 of file et_rot_mag_mag_plus.C.
References Lorene::Etoile::a_car, Lorene::Et_rot_mag::A_phi, Lorene::Et_rot_mag::A_t, Lorene::abs(), Lorene::Cmp::allocate_all(), Lorene::Cmp::annule(), Lorene::Cmp::asymptot(), Lorene::Etoile_rot::b_car, Lorene::Et_rot_mag::B_phi, Lorene::Etoile_rot::bbb, Lorene::Coord::c, Lorene::Tenseur::change_triad(), Lorene::cos(), Lorene::Cmp::dec_dzpuis(), Lorene::Cmp::div_r(), Lorene::Cmp::dsdr(), Lorene::Etoile::ener, Lorene::Coord::fait(), Lorene::flat_scalar_prod_desal(), Lorene::Map::get_bvect_cart(), Lorene::Map::get_bvect_spher(), Lorene::Tbl::get_dim(), Lorene::Cmp::get_dzpuis(), Lorene::Tenseur::get_etat(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Tenseur::gradient_spher(), Lorene::Cmp::inc2_dzpuis(), Lorene::Et_rot_mag::is_conduct(), Lorene::Et_rot_mag::j_phi, Lorene::Et_rot_mag::j_t, Lorene::Etoile_rot::l_surf(), Lorene::Cmp::laplacien(), Lorene::log(), Lorene::Etoile_rot::logn, Lorene::max(), Lorene::Etoile::mp, Lorene::Cmp::mult_cost(), Lorene::Valeur::mult_ct(), Lorene::Cmp::mult_r(), Lorene::Cmp::mult_rsint(), Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile::nzet, Lorene::Etoile_rot::omega, Lorene::Cmp::poisson(), Lorene::Tenseur::poisson_vect(), Lorene::pow(), Lorene::Etoile::press, Lorene::Et_rot_mag::Q, Lorene::Map::r, Lorene::Tenseur::set(), Lorene::Tbl::set(), Lorene::Cmp::set(), Lorene::Cmp::set_dzpuis(), Lorene::Tbl::set_etat_qcq(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_std_base(), Lorene::Cmp::srdsdt(), Lorene::Valeur::ssint(), Lorene::Cmp::std_base_scal(), Lorene::Tbl::t, Lorene::Map::tet, Lorene::Etoile_rot::tnphi, Lorene::Cmp::va, Lorene::Valeur::val_point_jk(), Lorene::Map_radial::val_r_jk(), and Lorene::Etoile::xi_surf().
|
virtualinherited |
Baryon mass.
Reimplemented from Lorene::Etoile.
Reimplemented in Lorene::Et_rot_bifluid.
Definition at line 125 of file et_rot_global.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::bbb, Lorene::Etoile::gam_euler, Lorene::Tenseur::get_etat(), Lorene::Cmp::integrale(), Lorene::Etoile::nbar, Lorene::Etoile::p_mass_b, Lorene::Etoile::relativistic, and Lorene::Cmp::std_base_scal().
|
virtual |
Gravitational mass.
Reimplemented from Lorene::Et_rot_mag.
Definition at line 538 of file et_magnetisation_comp.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::bbb, Lorene::Et_rot_mag::E_em, E_I, Lorene::Etoile::ener_euler, J_I, Lorene::Et_rot_mag::Jp_em, Lorene::Etoile_rot::mass_b(), Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile::p_mass_g, Lorene::Etoile::press, Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Tenseur::set_std_base(), Sij_I, Lorene::Et_rot_mag::Spp_em, Lorene::Etoile_rot::tnphi, and Lorene::Etoile_rot::uuu.
|
virtualinherited |
Mean radius.
Definition at line 472 of file et_rot_global.C.
References Lorene::Etoile_rot::area(), and Lorene::sqrt().
|
virtual |
Computes the electromagnetic part of the stress-energy tensor.
Reimplemented from Lorene::Et_rot_mag.
Definition at line 461 of file et_magnetisation_comp.C.
References Lorene::Etoile::a_car, Lorene::Et_rot_mag::A_phi, Lorene::Et_rot_mag::A_t, Lorene::Etoile_rot::b_car, Lorene::Vector::change_triad(), Lorene::Tensor::down(), Lorene::Et_rot_mag::E_em, E_I, Lorene::Et_rot_mag::Elec(), Lorene::flat_scalar_prod(), Lorene::flat_scalar_prod_desal(), Lorene::Etoile::gam_euler, Lorene::Map::get_bvect_cart(), Lorene::Map::get_bvect_spher(), Lorene::Tenseur::get_etat(), get_magnetisation(), Lorene::Tenseur::gradient_spher(), J_I, Lorene::Et_rot_mag::Jp_em, Lorene::Et_rot_mag::Magn(), Lorene::Etoile::mp, Lorene::Cmp::mult_rsint(), Lorene::Etoile::nnn, Lorene::Tenseur::set(), Lorene::Tensor::set(), Lorene::Vector::set(), Lorene::Scalar::set_dzpuis(), Sij_I, Lorene::Et_rot_mag::Spp_em, Lorene::sqrt(), Lorene::Et_rot_mag::Srr_em, Lorene::Etoile_rot::tnphi, and Lorene::Etoile::u_euler.
|
virtualinherited |
Quadrupole moment.
The quadrupole moment Q is defined according to Eq. (11) of [Pappas and Apostolatos, Physical Review Letters 108, 231104 (2012)]. This is a corrected version of the quadrupole moment defined by [Salgado, Bonazzola, Gourgoulhon and Haensel, Astron. Astrophys. 291 , 155 (1994)]. Following this definition,
Reimplemented in Lorene::Et_rot_bifluid.
Definition at line 600 of file et_rot_global.C.
References Lorene::Etoile_rot::mass_g(), Lorene::Etoile_rot::mom_quad_Bo(), Lorene::Etoile_rot::mom_quad_old(), Lorene::Etoile_rot::p_mom_quad, Lorene::pow(), and Lorene::Etoile::relativistic.
|
virtual |
Part of the quadrupole moment.
Reimplemented from Lorene::Etoile_rot.
Definition at line 828 of file et_magnetisation_comp.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::bbb, Lorene::Cmp::integrale(), Lorene::Cmp::mult_rsint(), Lorene::Etoile::nnn, Lorene::Etoile_rot::p_mom_quad_Bo, Lorene::Etoile::press, Sij_I, and Lorene::Cmp::std_base_scal().
|
virtual |
Part of the quadrupole moment.
This term
Reimplemented from Lorene::Et_rot_mag.
Definition at line 755 of file et_magnetisation_comp.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::ak_car, Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::bbb, Lorene::Cmp::check_dzpuis(), Lorene::Et_rot_mag::E_em, E_I, Lorene::Etoile::ener_euler, Lorene::flat_scalar_prod(), Lorene::Cmp::get_etat(), Lorene::Tenseur::gradient_spher(), Lorene::Cmp::inc2_dzpuis(), Lorene::log(), Lorene::Etoile_rot::logn, Lorene::Etoile::mp, Lorene::Valeur::mult_ct(), Lorene::Cmp::mult_r(), Lorene::Etoile::nbar, Lorene::Etoile_rot::p_mom_quad_old, Lorene::Etoile::relativistic, Lorene::Etoile::s_euler, Lorene::Tenseur::set(), Lorene::Tenseur::set_std_base(), Sij_I, Lorene::Et_rot_mag::Spp_em, and Lorene::Cmp::va.
void Lorene::Et_magnetisation::operator= | ( | const Et_magnetisation & | et | ) |
Assignment to another Et_rot_mag.
Definition at line 174 of file et_magnetisation.C.
References E_I, include_magnetisation, J_I, Lorene::Et_rot_mag::operator=(), Sij_I, use_B_in_eos, and xmag.
Operator >> (virtual function called by the operator <<).
Reimplemented from Lorene::Et_rot_mag.
Definition at line 341 of file et_magnetisation.C.
References include_magnetisation, Lorene::max(), Lorene::maxabs(), Lorene::Et_rot_mag::operator>>(), use_B_in_eos, and xmag.
Printing of some informations, excluding all global quantities.
Reimplemented in Lorene::Et_rot_bifluid.
Definition at line 630 of file etoile_rot.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::aplat(), Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::dzeta, Lorene::Etoile::ener, Lorene::Etoile::ent, Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Etoile_rot::get_omega_c(), Lorene::Etoile_rot::logn, Lorene::Etoile::mp, Lorene::Etoile::nbar, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile::nzet, Lorene::Etoile::press, Lorene::Etoile::ray_eq(), and Lorene::Etoile_rot::uuu.
|
inherited |
Computed charge deduced from the asymptotic behaviour of At [SI units].
Definition at line 217 of file et_rot_mag_global.C.
References Lorene::Et_rot_mag::A_t, Lorene::Cmp::asymptot(), Lorene::Cmp::get_etat(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), Lorene::Etoile::mp, and Lorene::pow().
|
inherited |
Computed charge from the integration of charge density over the star (i.e.
without surface charge) [SI units].
Definition at line 238 of file et_rot_mag_global.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::bbb, Lorene::Tenseur::get_etat(), Lorene::Cmp::integrale(), Lorene::Et_rot_mag::j_t, Lorene::Etoile::nbar, Lorene::Etoile::nnn, Lorene::pow(), Lorene::Etoile::relativistic, and Lorene::Cmp::std_base_scal().
|
virtualinherited |
Circumferential radius.
Definition at line 390 of file et_rot_global.C.
References Lorene::Etoile_rot::bbb, Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_t(), Lorene::Etoile::mp, Lorene::Etoile::nzet, Lorene::Etoile_rot::p_r_circ, and Lorene::Etoile::ray_eq().
Circumferential radius of the innermost stable circular orbit (ISCO).
ost | output stream to give details of the computation; if set to 0x0 [default value], no details will be given. |
Definition at line 84 of file et_rot_isco.C.
References Lorene::Param::add_cmp(), Lorene::Param::add_int(), Lorene::Cmp::annule(), Lorene::Etoile_rot::bbb, Lorene::Cmp::dsdr(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), Lorene::Etoile::mp, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile::nzet, Lorene::Etoile_rot::p_espec_isco, Lorene::Etoile_rot::p_f_eq, Lorene::Etoile_rot::p_f_isco, Lorene::Etoile_rot::p_lspec_isco, Lorene::Etoile_rot::p_r_isco, Lorene::Map::r, Lorene::Etoile::ray_eq(), Lorene::sqrt(), Lorene::Cmp::std_base_scal(), Lorene::Cmp::va, Lorene::Valeur::val_point(), Lorene::Map::val_r(), and Lorene::zerosec().
|
inherited |
Coordinate radius at
Definition at line 120 of file etoile_global.C.
References Lorene::Map::get_mg(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_p(), Lorene::Mg3d::get_type_t(), Lorene::Etoile::l_surf(), Lorene::Etoile::mp, Lorene::Etoile::p_ray_eq, Lorene::Map::val_r(), and Lorene::Etoile::xi_surf().
Coordinate radius at
Definition at line 440 of file etoile_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::Etoile::l_surf(), Lorene::Etoile::mp, Lorene::Map::val_r(), and Lorene::Etoile::xi_surf().
|
inherited |
Coordinate radius at
Definition at line 335 of file etoile_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::Etoile::l_surf(), Lorene::Etoile::mp, Lorene::Etoile::p_ray_eq_3pis2, Lorene::Etoile::ray_eq_pis2(), Lorene::Map::val_r(), and Lorene::Etoile::xi_surf().
|
inherited |
Coordinate radius at
Definition at line 256 of file etoile_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::Etoile::l_surf(), Lorene::Etoile::mp, Lorene::Etoile::p_ray_eq_pi, Lorene::Etoile::ray_eq(), Lorene::Map::val_r(), and Lorene::Etoile::xi_surf().
|
inherited |
Coordinate radius at
Definition at line 169 of file etoile_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::Etoile::l_surf(), Lorene::Etoile::mp, Lorene::Etoile::p_ray_eq_pis2, Lorene::Map::val_r(), and Lorene::Etoile::xi_surf().
|
inherited |
Coordinate radius at
Definition at line 415 of file etoile_global.C.
References Lorene::Map::get_mg(), Lorene::Mg3d::get_type_t(), Lorene::Etoile::l_surf(), Lorene::Etoile::mp, Lorene::Etoile::p_ray_pole, Lorene::Map::val_r(), and Lorene::Etoile::xi_surf().
Save in a file.
Reimplemented from Lorene::Et_rot_mag.
Definition at line 321 of file et_magnetisation.C.
References E_I, include_magnetisation, J_I, Lorene::Et_rot_mag::sauve(), Lorene::Scalar::sauve(), Lorene::Tensor::sauve(), Lorene::Tensor_sym::sauve(), Sij_I, use_B_in_eos, and xmag.
|
protectedvirtualinherited |
Sets to 0x0
all the pointers on derived quantities.
Reimplemented from Lorene::Etoile_rot.
Definition at line 267 of file et_rot_mag.C.
References Lorene::Etoile_rot::set_der_0x0().
Assignment of the enthalpy field.
Definition at line 465 of file etoile.C.
References Lorene::Etoile::del_deriv(), Lorene::Etoile::ent, and Lorene::Etoile::equation_of_state().
|
inlineinherited |
|
virtualinherited |
Ratio T/W.
Reimplemented from Lorene::Etoile_rot.
Definition at line 344 of file et_rot_mag_global.C.
References Lorene::Etoile::a_car, Lorene::Et_rot_mag::angu_mom(), Lorene::Etoile_rot::bbb, Lorene::Etoile::ener, Lorene::Etoile::gam_euler, Lorene::Cmp::integrale(), Lorene::Etoile_rot::logn, Lorene::Et_rot_mag::mass_g(), Lorene::Etoile::nbar, Lorene::Etoile_rot::omega, Lorene::Etoile_rot::p_tsw, Lorene::Etoile::relativistic, and Lorene::Cmp::std_base_scal().
|
inherited |
Computes metric coefficients from known potentials.
The calculation is performed starting from the quantities logn
, dzeta
, tggg
and shift
, which are supposed to be up to date.
From these, the following fields are updated: nnn
, a_car
, bbb
and b_car
.
Definition at line 69 of file et_rot_upmetr.C.
References Lorene::Etoile::a_car, Lorene::Etoile_rot::b_car, Lorene::Etoile_rot::bbb, Lorene::Etoile_rot::del_deriv(), Lorene::Etoile_rot::dzeta, Lorene::exp(), Lorene::Etoile_rot::extrinsic_curvature(), Lorene::Etoile_rot::logn, Lorene::Etoile::nnn, Lorene::Tenseur::set(), Lorene::Tenseur::set_etat_qcq(), Lorene::Tenseur::set_std_base(), Lorene::Etoile_rot::tggg, and Lorene::Etoile::unsurc2.
|
inline |
Public accessor to the include_magnetisation
flag.
Definition at line 674 of file et_rot_mag.h.
References include_magnetisation.
Description of the stellar surface: returns a 2-D Tbl
containing the values of the radial coordinate
The stellar surface is defined as the location where the enthalpy (member ent
) vanishes.
Definition at line 101 of file etoile_global.C.
References Lorene::Etoile::l_surf(), Lorene::Etoile::p_l_surf, and Lorene::Etoile::p_xi_surf.
|
virtualinherited |
Backward redshift factor at equator.
Definition at line 545 of file et_rot_global.C.
References Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_t(), Lorene::Etoile::mp, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile::nzet, Lorene::Etoile_rot::p_z_eqb, Lorene::Etoile_rot::r_circ(), Lorene::sqrt(), and Lorene::Etoile_rot::uuu.
|
virtualinherited |
Forward redshift factor at equator.
Definition at line 513 of file et_rot_global.C.
References Lorene::Map::get_mg(), Lorene::Mg3d::get_nr(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_type_t(), Lorene::Etoile::mp, Lorene::Etoile::nnn, Lorene::Etoile_rot::nphi, Lorene::Etoile::nzet, Lorene::Etoile_rot::p_z_eqf, Lorene::Etoile_rot::r_circ(), Lorene::sqrt(), and Lorene::Etoile_rot::uuu.
|
virtualinherited |
Redshift factor at North pole.
Definition at line 580 of file et_rot_global.C.
References Lorene::Etoile::nnn, Lorene::Etoile_rot::p_z_pole, and Lorene::Etoile::ray_pole().
|
protectedinherited |
|
protectedinherited |
Amplitude of the curent/charge function.
Definition at line 180 of file et_rot_mag.h.
|
protectedinherited |
Definition at line 155 of file et_rot_mag.h.
|
protectedinherited |
t-component of the elecctromagnetic potential 1-form, divided by
Definition at line 150 of file et_rot_mag.h.
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
Definition at line 157 of file et_rot_mag.h.
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
Flag: conduc=0->isolator, 1->perfect conductor.
Definition at line 181 of file et_rot_mag.h.
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
electromagnetic energy density in the Eulerian frame
Definition at line 161 of file et_rot_mag.h.
|
protected |
Interaction (magnetisation) energy density.
Definition at line 624 of file et_rot_mag.h.
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protected |
Flag : true if magnetisation terms are included in the equations.
Definition at line 620 of file et_rot_mag.h.
|
protected |
Interaction momentum density 3-vector.
Definition at line 627 of file et_rot_mag.h.
|
protectedinherited |
Definition at line 159 of file et_rot_mag.h.
|
protectedinherited |
t-component of the current 4-vector
Definition at line 158 of file et_rot_mag.h.
|
protectedinherited |
Definition at line 167 of file et_rot_mag.h.
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
mutableprotectedinherited |
|
protectedinherited |
|
protectedinherited |
In the case of a perfect conductor, the requated baryonic charge.
For an isolator, the charge/baryon.
Definition at line 179 of file et_rot_mag.h.
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protected |
Interaction stress 3-tensor.
Definition at line 630 of file et_rot_mag.h.
|
protectedinherited |
Definition at line 173 of file et_rot_mag.h.
|
protectedinherited |
rr component of the electromagnetic stress 3-tensor, as measured in the Eulerian frame. (not used and set to 0, should be supressed)
Definition at line 170 of file et_rot_mag.h.
|
protectedinherited |
|
protectedinherited |
Effective source at the previous step for the resolution of the Poisson equation for the scalar Map_et::poisson
.
|
protectedinherited |
Effective source at the previous step for the resolution of the Poisson equation for nuf
by means of Map_et::poisson
.
|
protectedinherited |
Effective source at the previous step for the resolution of the Poisson equation for nuq
by means of Map_et::poisson
.
|
protectedinherited |
|
protectedinherited |
Effective source at the previous step for the resolution of the vector Poisson equation for
mp
)
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protected |
Flag : true if the value of the magnetic field is used in the Eos.
Definition at line 617 of file et_rot_mag.h.
|
protectedinherited |
|
protectedinherited |
|
protected |
The magnetisation scalar.
Definition at line 622 of file et_rot_mag.h.