283 Star(
Map& mp_i,
int nzet_i,
const Eos& eos_i) ;
295 Star(
Map& mp_i,
const Eos& eos_i, FILE* fich) ;
349 const Tbl* pent_limit = 0x0 ) ;
406 const Scalar& get_lnq()
const {
return lnq;} ;
414 virtual void sauve(FILE* )
const ;
421 virtual ostream&
operator>>(ostream& )
const ;
460 virtual double mass_b()
const = 0 ;
463 virtual double mass_g()
const = 0 ;
704 bool irrot,
bool conf_flat) ;
743 Scalar& set_pot_centri() ;
910 bool is_conf_flat()
const {
return conf_flat; } ;
915 virtual void sauve(FILE* )
const ;
919 virtual ostream&
operator>>(ostream& )
const ;
925 virtual double mass_b()
const ;
928 virtual double mass_g()
const ;
931 virtual double xa_barycenter()
const ;
969 void update_metric(
const Star_bin& comp,
double omega) ;
981 double relax,
double omega) ;
988 void update_metric_der_comp(
const Star_bin& comp,
double omega) ;
1000 void kinematics(
double omega,
double x_axe) ;
1012 void extrinsic_curvature(
double omega) ;
1031 void equilibrium(
double ent_c,
int mermax,
int mermax_potvit,
1032 int mermax_poisson,
double relax_poisson,
1033 double relax_potvit,
double thres_adapt,
Tbl& diff,
1050 double velocity_potential(
int mermax,
double precis,
double relax) ;
1063 void relaxation(
const Star_bin& star_prev,
double relax_ent,
1064 double relax_met,
int mer,
int fmer_met) ;
1068 void test_K_Hi()
const ;
1071 void helical(
double omega)
const ;
1291 Scalar& set_pot_centri() ;
1306 Vector& set_beta_auto() ;
1427 virtual void sauve(FILE* )
const ;
1431 virtual ostream&
operator>>(ostream& )
const ;
1437 virtual double mass_b()
const ;
1440 virtual double mass_g()
const ;
1443 virtual double xa_barycenter()
const ;
1511 void kinematics(
double omega,
double x_axe) ;
1521 void extrinsic_curvature() ;
1544 void equilibrium(
double ent_c,
int mermax,
int mermax_potvit,
1545 int mermax_poisson,
double relax_poisson,
1546 double relax_potvit,
double thres_adapt,
1547 const Tbl& fact,
const Tbl* pent_limit,
1563 double velocity_potential(
int mermax,
double precis,
double relax) ;
1576 void relaxation(
const Star_bin_xcts& star_prev,
double relax_ent,
1577 double relax_met,
int mer,
int fmer_met) ;
Sym_tensor hij_comp
Deviation of the inverse conformal metric from the inverse flat metric generated principally by the ...
Vector dcov_phi
Covariant derivative of the logarithm of the conformal factor.
friend ostream & operator<<(ostream &, const Star &)
Display.
Scalar chi
Total function .
double * p_mass_b
Baryon mass.
bool irrotational
true for an irrotational star, false for a corotating one
Metric for tensor calculation.
Vector bsn
3-vector shift, divided by N, of the rotating coordinates, .
const Vector & get_wit_w() const
Returns the spatial projection of the fluid 3-velocity with respect to the co-orbiting observer...
void set_logn_auto(const Scalar &logn_auto_new)
Assignment of a new logn_auto.
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
Vector dcov_chi
Covariant derivative of the function .
const Scalar & get_press() const
Returns the fluid pressure.
Scalar psi0
Scalar potential of the non-translational part of fluid 4-velocity (in the irrotational case) ...
const Scalar & get_ener_euler() const
Returns the total energy density with respect to the Eulerian observer.
Vector wit_w
Spatial projection of the fluid 3-velocity with respect to the co-orbiting observer.
const Vector & get_wit_w() const
Returns the spatial projection of the fluid 3-velocity with respect to the co-orbiting observer...
const Vector & get_beta() const
Returns the shift vector .
const Sym_tensor & get_haij_auto() const
Returns the part of the extrinsic curvature tensor generated by beta_auto.
Vector wit_w
Spatial projection of the fluid 3-velocity with respect to the co-orbiting observer.
Map & mp
Mapping associated with the star.
Scalar ssjm1_lnq
Effective source at the previous step for the resolution of the Poisson equation for lnq_auto...
Sym_tensor hij_auto
Deviation of the inverse conformal metric from the inverse flat metric generated principally by the ...
const Vector & get_u_euler() const
Returns the fluid 3-velocity with respect to the Eulerian observer.
double * p_ray_eq_3pis2
Coordinate radius at , .
const Vector & get_bsn() const
Returns the shift vector, divided by N, of the rotating coordinates, .
double * p_ray_eq_pi
Coordinate radius at , .
const Scalar & get_loggam() const
Returns the logarithm of the Lorentz factor between the fluid and the co-orbiting observer...
const Scalar & get_lnq_comp() const
Returns the part of the vector field generated principally by the companion star.
const Vector & get_beta_comp() const
Returns the part of the shift vector generated principally by the companion ( ). ...
Vector d_psi
Gradient of (in the irrotational case) (Spherical components with respect to the mapping of the star...
Scalar ssjm1_h31
Effective source at the previous step for the resolution of the Poisson equation for h20_auto...
Binary systems in eXtended Conformal Thin Sandwich formulation.
const Eos & eos
Equation of state of the stellar matter.
Scalar ssjm1_h11
Effective source at the previous step for the resolution of the Poisson equation for h00_auto...
Equation of state base class.
Sym_tensor stress_euler
Spatial part of the stress-energy tensor with respect to the Eulerian observer.
Flat metric for tensor calculation.
Tensor field of valence 0 (or component of a tensorial field).
const Scalar & get_ener() const
Returns the proper total energy density.
Base class for coordinate mappings.
const Vector & get_bsn() const
Returns the shift vector, divided by N, of the rotating coordinates, .
bool irrotational
true for an irrotational star, false for a corotating one
const Vector & get_beta_auto() const
Returns the part of the shift vector generated principally by the star ( ).
Scalar ssjm1_h33
Effective source at the previous step for the resolution of the Poisson equation for h22_auto...
Scalar logn_comp
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principally by...
Vector dcov_Psi
Covariant derivative of the conformal factor .
const Scalar & get_logn_comp() const
Returns the part of the lapse logarithm (gravitational potential at the Newtonian limit) generated pr...
virtual void equilibrium_spher(double ent_c, double precis=1.e-14, const Tbl *pent_limit=0x0)
Computes a spherical static configuration.
Basic integer array class.
Scalar loggam
Logarithm of the Lorentz factor between the fluid and the co-orbiting observer.
const Sym_tensor & get_tkij_comp() const
Returns the part of the extrinsic curvature tensor generated by beta_comp.
Scalar kcar_auto
Part of the scalar generated by beta_auto, i.e.
const Metric & get_gamma() const
Returns the 3-metric .
const Map & get_mp() const
Returns the mapping.
Class for stars in binary system in eXtended Conformal Thin Sandwich formulation. ...
Scalar nbar
Baryon density in the fluid frame.
double * p_xa_barycenter
Absolute coordinate X of the barycenter of the baryon density.
const Scalar get_decouple() const
Returns the function used to construct beta_auto from beta.
Scalar s_euler
Trace of the stress scalar in the Eulerian frame.
Tensor field of valence 1.
const Scalar & get_pot_centri() const
Returns the centrifugal potential.
Scalar logn_auto
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principally by...
double * p_ray_eq
Coordinate radius at , .
Sym_tensor haij_auto
Part of the extrinsic curvature tensor generated by beta_auto.
Scalar ssjm1_chi
Effective source at the previous step for the resolution of the Poisson equation for ...
const Scalar & get_Psi() const
Return the conformal factor .
const Scalar & get_psi4() const
Return the conformal factor .
Scalar ssjm1_h21
Effective source at the previous step for the resolution of the Poisson equation for h10_auto...
const Tbl & xi_surf() const
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate ...
const Scalar & get_loggam() const
Returns the logarithm of the Lorentz factor between the fluid and the co-orbiting observer...
Tbl * p_xi_surf
Description of the stellar surface: 2-D Tbl containing the values of the radial coordinate on the su...
double * p_ray_eq_pis2
Coordinate radius at , .
void operator=(const Star &)
Assignment to another Star.
const Tensor & get_stress_euler() const
Returns the spatial part of the stress-energy tensor with respect to the Eulerian observer...
int nzet
Number of domains of *mp occupied by the star.
const Scalar & get_logn() const
Returns the logarithm of the lapse N.
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
void set_conf_flat(bool confflat)
Write if conformally flat.
Scalar loggam
Logarithm of the Lorentz factor between the fluid and the co-orbiting observer.
Scalar decouple
Function used to construct the part generated by the star from the total .
double ray_eq_3pis2() const
Coordinate radius at , [r_unit].
Scalar ener
Total energy density in the fluid frame.
Scalar pot_centri
Centrifugal potential.
Sym_tensor haij_comp
Part of the extrinsic curvature tensor generated by beta_comp.
Sym_tensor hij
Total deviation of the inverse conformal metric from the inverse flat metric.
const Scalar & get_kcar_auto() const
Returns the part of generated by beta_auto.
Vector ssjm1_wbeta
Effective source at the previous step for wbeta of the vector Poisson equation for wbeta (needed for ...
Scalar Psi_auto
Scalar field generated principally by the star.
const Scalar & get_psi0() const
Returns the non-translational part of the velocity potential.
Vector d_psi
Gradient of (in the irrotational case) (Spherical components with respect to the mapping of the star...
virtual void hydro_euler()
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
Scalar pot_centri
Centrifugal potential.
Scalar press
Fluid pressure.
const Metric & get_gtilde() const
Return the conformal 3-metric .
Vector w_beta
Solution for the vector part of the vector Poisson equation for .
const Vector & get_beta_auto() const
Returns the part of the shift vector generated principally by the star.
Star(Map &mp_i, int nzet_i, const Eos &eos_i)
Standard constructor.
const Scalar & get_Psi_auto() const
Returns the scalar field generated principally by the star.
const Vector & get_dcov_Psi() const
Returns the covariant derivative of the conformal factor .
Scalar chi_comp
Scalar field generated principally by the companion star.
const Metric & get_flat() const
Return the flat metric defined on the mapping (Spherical components with respect to the mapping of th...
bool is_irrotational() const
Returns true for an irrotational motion, false for a corotating one.
double ray_eq_pis2() const
Coordinate radius at , [r_unit].
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the star) ...
Scalar lnq_auto
Scalar field generated principally by the star.
const Vector & get_beta_comp() const
Returns the part of the shift vector generated principally by the star.
const Vector & get_dcon_phi() const
Returns the contravariant derivative of (logarithm of the conformal factor).
double * p_ray_pole
Coordinate radius at .
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...
Scalar lnq_comp
Scalar field generated principally by the companion star.
const Vector & get_dcov_phi() const
Returns the covariant derivative of (logarithm of the conformal factor).
const Metric & get_flat() const
Return the flat metric defined on the mapping (Spherical components with respect to the mapping of th...
Vector beta_comp
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
const Scalar & get_s_euler() const
Returns the trace of the stress tensor in the Eulerian frame.
const Scalar & get_lnq_auto() const
Returns the part of the vector field generated principally by the star.
virtual ~Star()
Destructor.
const Scalar & get_chi_comp() const
Returns the scalar field generated principally by the companion star.
Vector u_euler
Fluid 3-velocity with respect to the Eulerian observer.
Sym_tensor tkij_comp
Part of the extrinsic curvature tensor generated by beta_comp.
double ray_eq_pi() const
Coordinate radius at , [r_unit].
const Scalar & get_chi() const
Return the function .
Scalar chi_auto
Scalar field generated principally by the star.
Vector dcov_logn
Covariant derivative of the total logarithm of the lapse.
Scalar logn
Logarithm of the lapse N .
Class for stars in binary system.
const Sym_tensor & get_tkij_auto() const
Returns the part of the extrinsic curvature tensor generated by beta_auto.
Scalar ssjm1_khi
Effective source at the previous step for the resolution of the Poisson equation for khi...
Scalar psi4
Conformal factor .
Metric gtilde
Conformal metric .
const Scalar & get_Psi_comp() const
Returns the scalar field generated principally by the companion star.
Vector dcon_logn
Contravariant derivative of the total logarithm of the lapse.
Scalar psi0
Scalar potential of the non-translational part of fluid 4-velocity (in the irrotational case) ...
const Vector & get_dcov_logn() const
Returns the covariant derivative of .
double ray_pole() const
Coordinate radius at [r_unit].
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the star) ...
const Scalar & get_hacar_auto() const
Returns the part of generated by beta_auto.
Scalar Psi
Total conformal factor .
const Vector & get_dcov_chi() const
Returns the covariant derivative of .
const Scalar & get_psi4() const
Return the conformal factor .
Scalar hacar_auto
Part of the scalar generated by beta_auto, i.e.
const Scalar & get_psi0() const
Returns the non-translational part of the velocity potential.
const Scalar & get_kcar_comp() const
Returns the part of generated by beta_comp.
virtual void del_deriv() const
Deletes all the derived quantities.
const Scalar & get_nbar() const
Returns the proper baryon density.
virtual void del_hydro_euler()
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer...
Scalar nn
Lapse function N .
int get_nzet() const
Returns the number of domains occupied by the star.
Scalar ssjm1_h22
Effective source at the previous step for the resolution of the Poisson equation for h11_auto...
Vector beta_comp
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
const Sym_tensor & get_haij_comp() const
Returns the part of the extrinsic curvature tensor generated by beta_comp.
const Scalar & get_pot_centri() const
Returns the centrifugal potential.
Vector beta_auto
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
Sym_tensor tkij_auto
Part of the extrinsic curvature tensor generated by beta_auto.
double * p_xa_barycenter
Absolute coordinate X of the barycenter of the baryon density.
const Scalar & get_ent() const
Returns the enthalpy field.
const Scalar & get_hacar_comp() const
Returns the part of generated by beta_comp.
bool is_irrotational() const
Returns true for an irrotational motion, false for a corotating one.
Map & set_mp()
Read/write of the mapping.
bool conf_flat
true if the 3-metric is conformally flat, false for a more general metric.
double ray_eq() const
Coordinate radius at , [r_unit].
Scalar ssjm1_h32
Effective source at the previous step for the resolution of the Poisson equation for h21_auto...
Vector beta_auto
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
Scalar Psi_comp
Scalar field generated principally by the companion star.
Scalar kcar_comp
Part of the scalar generated by beta_auto and beta_comp, i.e.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
Scalar hacar_comp
Part of the scalar generated by beta_auto and beta_comp, i.e.
virtual double mass_b() const =0
Baryon mass.
virtual void sauve(FILE *) const
Save in a file.
const Eos & get_eos() const
Returns the equation of state.
const Vector & get_d_psi() const
Returns the covariant derivative of the velocity potential (Spherical components with respect to the ...
Scalar ssjm1_khi
Effective source at the previous step for the resolution of the Poisson equation for khi...
const Sym_tensor & get_hij() const
Return the total deviation of the inverse conformal metric from the inverse flat metric...
const Scalar & get_nn() const
Returns the lapse function N.
const Scalar & get_gam_euler() const
Returns the Lorentz factor between the fluid and Eulerian observers.
void equation_of_state()
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
const Scalar & get_chi_auto() const
Returns the scalar field generated principally by the star.
Scalar khi
Solution for the scalar part of the vector Poisson equation for .
const Scalar & get_logn_auto() const
Returns the part of the lapse logarithm (gravitational potential at the Newtonian limit) generated pr...
Vector bsn
3-vector shift, divided by N, of the rotating coordinates, .
Scalar ssjm1_psi
Effective source at the previous step for the resolution of the Poisson equation for ...
const Vector & get_d_psi() const
Returns the covariant derivative of the velocity potential (Spherical components with respect to the ...
Scalar ener_euler
Total energy density in the Eulerian frame.
const Sym_tensor & get_hij_comp() const
Return the deviation of the inverse conformal metric from the inverse flat metric generated principa...
Itbl * p_l_surf
Description of the stellar surface: 2-D Itbl containing the values of the domain index l on the surfa...
double * p_mass_g
Gravitational mass.
Scalar psi4
Conformal factor .
void set_enthalpy(const Scalar &)
Assignment of the enthalpy field.
virtual double mass_g() const =0
Gravitational mass.
Class intended to describe valence-2 symmetric tensors.
void set_lnq_auto(const Scalar &lnq_auto_new)
Assignment of a new lnq_auto.
Vector dcon_phi
Contravariant derivative of the logarithm of the conformal factor.
const Sym_tensor & get_hij_auto() const
Return the deviation of the inverse conformal metric from the inverse flat metric principally genera...
const Vector & get_dcon_logn() const
Returns the contravariant derivative of .
Scalar ssjm1_logn
Effective source at the previous step for the resolution of the Poisson equation for logn_auto...