84 #include "star_rot_dirac.h" 86 #include "utilitaires.h" 155 phi_kill.set(2).set_etat_zero() ;
156 phi_kill.set(3) = 1. ;
157 phi_kill.set(3).std_spectral_base() ;
158 phi_kill.set(3).mult_rsint() ;
192 p_tsw =
new double( tcin / ( mass_p + tcin -
mass_g() ) ) ;
227 sou_m =
sqrt( k_det )*sou_m ;
251 sou_tmp += term_2 + term_3 ;
258 sou_q =
sqrt( k_det )*sou_q ;
262 p_grv2 =
new double(
double(1) + integrale2d(sou_m)/integrale2d(sou_q) ) ;
299 sou_q = sou_q + tmp_1 - tmp_2 ;
319 p_grv3 =
new double( (int_grav + int_mat) / int_mat ) ;
341 int i_b = mg->
get_nr(l_b) - 1 ;
345 double gamma_phi =
gamma.
cov()(3,3).val_grid_point(l_b, k_b, j_b, i_b) ;
368 cout <<
"The polar circumferential radius is only well defined\n" 369 <<
"with np = 1!" << endl ;
379 for (
int j=0; j<nt; j++) {
380 double theta = (+tet)(0, 0, j, 0) ;
383 for (
int lz=0; lz<
nzet; lz++) {
385 for (
int i=0; i<nrz; i++) {
396 for (
int j=0; j<nt; j++) {
397 double theta = (+tet)(0, 0, j, 0) ;
400 double grr = gam(1,1).get_spectral_va().val_point_jk(ls, xs, j, 0) ;
401 double grt = gam(1,2).get_spectral_va().val_point_jk(ls, xs, j, 0) ;
402 double gtt = gam(2,2).get_spectral_va().val_point_jk(ls, xs, j, 0) ;
403 double rr =
mp.
val_r(ls, xs, theta, phi) ;
407 =
sqrt(grr*dr*dr + 2*grt*rr*dr + gtt*rr*rr) ;
virtual double grv3() const
Error on the virial identity GRV3.
double * p_mass_b
Baryon mass.
double * p_angu_mom
Angular momentum.
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
double * p_rp_circ
Circumferential polar radius.
Mtbl_cf * c_cf
Coefficients of the spectral expansion of the function.
double omega
Rotation angular velocity ([f_unit] )
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.
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).
double integrale() const
Computes the integral over all space of *this .
int get_type_t() const
Returns the type of sampling in the direction: SYM : : symmetry with respect to the equatorial pl...
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
virtual double r_circ() const
Circumferential equatorial radius.
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.
virtual double angu_mom() const
Angular momentum.
int nzet
Number of domains of *mp occupied by the star.
double * p_r_circ
Circumferential equatorial radius.
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 tsw() const
Ratio T/W.
virtual double mass_g() const
Gravitational mass.
Scalar press
Fluid pressure.
virtual double ellipt() const
Ellipticity e.
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 grv2() const
Error on the virial identity GRV2.
int get_nzone() const
Returns the number of domains.
double * p_grv3
Error on the virial identity GRV3.
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...
double * p_grv2
Error on the virial identity GRV2.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Scalar logn
Logarithm of the lapse N .
virtual double mass_b() const
Baryonic mass.
Active physical coordinates and mapping derivatives.
virtual const Scalar & determinant() const
Returns the determinant.
virtual double aplat() const
Flattening r_pole/r_eq.
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 .
virtual double rp_circ() const
Circumferential polar radius.
double ray_eq() const
Coordinate radius at , [r_unit].
const Scalar & dsdr() const
Returns of *this .
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.
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)