31 #include "star_rot_cfc.h" 33 #include "utilitaires.h" 101 phi_kill.set(2).set_etat_zero() ;
102 phi_kill.set(3) = 1. ;
103 phi_kill.set(3).std_spectral_base() ;
104 phi_kill.set(3).mult_rsint() ;
138 p_tsw =
new double( tcin / ( mass_p + tcin -
mass_g() ) ) ;
173 sou_m =
sqrt( k_det )*sou_m ;
198 sou_tmp += term_2 + term_3 ;
205 sou_q =
sqrt( k_det )*sou_q ;
209 p_grv2 =
new double(
double(1) + integrale2d(sou_m)/integrale2d(sou_q) ) ;
247 sou_q = sou_q + psi6*(tmp_1 - tmp_2) ;
261 sou_m = psi6 * sou_m ;
265 double int_mat = sou_m.integrale() ;
267 p_grv3 =
new double( (int_grav + int_mat) / int_mat ) ;
289 int i_b = mg->
get_nr(l_b) - 1 ;
293 double gamma_phi =
gamma.
cov()(3,3).val_grid_point(l_b, k_b, j_b, i_b) ;
316 cout <<
"The polar circumferential radius is only well defined\n" 317 <<
"with np = 1!" << endl ;
327 for (
int j=0; j<nt; j++) {
328 double theta = (+tet)(0, 0, j, 0) ;
331 for (
int lz=0; lz<
nzet; lz++) {
333 for (
int i=0; i<nrz; i++) {
344 for (
int j=0; j<nt; j++) {
345 double theta = (+tet)(0, 0, j, 0) ;
348 double grr = gam(1,1).get_spectral_va().val_point_jk(ls, xs, j, 0) ;
349 double grt = gam(1,2).get_spectral_va().val_point_jk(ls, xs, j, 0) ;
350 double gtt = gam(2,2).get_spectral_va().val_point_jk(ls, xs, j, 0) ;
351 double rr =
mp.
val_r(ls, xs, theta, phi) ;
355 =
sqrt(grr*dr*dr + 2*grt*rr*dr + gtt*rr*rr) ;
double * p_grv2
Error on the virial identity GRV2.
double * p_mass_b
Baryon mass.
virtual double tsw() const
Ratio T/W.
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
Mtbl_cf * c_cf
Coefficients of the spectral expansion of the function.
Map & mp
Mapping associated with the star.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
void coef() const
Computes the coeffcients of *this.
Cmp sqrt(const Cmp &)
Square root.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
virtual double grv3() const
Error on the virial identity GRV3.
Standard units of space, time and mass.
const Scalar & dsdt() const
Returns of *this .
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
virtual double angu_mom() const
Angular momentum.
double integrale() const
Computes the integral over all space of *this .
virtual double r_circ() const
Circumferential equatorial radius.
double * p_rp_circ
Circumferential polar radius.
int get_type_t() const
Returns the type of sampling in the direction: SYM : : symmetry with respect to the equatorial pl...
Scalar psi4
Conformal factor .
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
Scalar nbar
Baryon density in the fluid frame.
Scalar s_euler
Trace of the stress scalar in the Eulerian frame.
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
Tensor field of valence 1.
const Tensor_sym & get_delta() const
Returns the tensor which defines the connection with respect to the flat one: is the difference bet...
Coord tet
coordinate centered on the grid
const Tbl & xi_surf() const
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate ...
void annule_hard()
Sets the Scalar to zero in a hard way.
int nzet
Number of domains of *mp occupied by the star.
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
virtual double val_r(int l, double xi, double theta, double pphi) const =0
Returns the value of the radial coordinate r for a given in a given domain.
Scalar ener
Total energy density in the fluid frame.
virtual double aplat() const
Flattening r_pole/r_eq.
double omega
Rotation angular velocity ([f_unit] )
Scalar press
Fluid pressure.
virtual const Connection & connect() const
Returns the connection.
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...
virtual double rp_circ() const
Circumferential polar radius.
double * p_angu_mom
Angular momentum.
double * p_r_circ
Circumferential equatorial radius.
int get_nzone() const
Returns the number of domains.
virtual double mass_b() const
Baryonic mass.
double val_point_jk(int l, double x, int j, int k) const
Computes the value of the field represented by *this at an arbitrary point in , but collocation point...
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Scalar logn
Logarithm of the lapse N .
Active physical coordinates and mapping derivatives.
double ray_pole() const
Coordinate radius at [r_unit].
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
virtual const Sym_tensor & cov() const
Read-only access to the covariant representation.
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
Scalar nn
Lapse function N .
double ray_eq() const
Coordinate radius at , [r_unit].
const Scalar & dsdr() const
Returns of *this .
double * p_grv3
Error on the virial identity GRV3.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Valeur & set_spectral_va()
Returns va (read/write version)
Scalar & set(int)
Read/write access to a component.
virtual double grv2() const
Error on the virial identity GRV2.
void div_r_dzpuis(int ced_mult_r)
Division by r everywhere but with the output flag dzpuis set to ced_mult_r .
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
Vector j_euler
Momentum density 3-vector with respect to the Eulerian observer.
Scalar ener_euler
Total energy density in the Eulerian frame.
double * p_mass_g
Gravitational mass.
Class intended to describe valence-2 symmetric tensors.
const Valeur & get_spectral_va() const
Returns va (read only version)
virtual double mass_g() const
Gravitational mass.
virtual double ellipt() const
Ellipticity e.