91 double precis = 1.e-15 ;
135 if ( type_t == SYM ) {
136 assert( (type_p == SYM) || (type_p == NONSYM) ) ;
141 double theta = M_PI / 2 ;
149 assert( (type_p == SYM) || (type_p == NONSYM) ) ;
154 double theta = M_PI / 2 ;
183 if ( type_t == SYM ) {
186 double theta = M_PI / 2 ;
187 double phi = M_PI / 2 ;
200 assert( np % 4 == 0 ) ;
209 cout <<
"Etoile::ray_eq_pis2 : the case type_p = " 210 << type_p <<
" is not contemplated yet !" << endl ;
219 double theta = M_PI / 2 ;
220 double phi = M_PI / 2 ;
233 assert( np % 4 == 0 ) ;
242 cout <<
"Etoile::ray_eq_pis2 : the case type_p = " 243 << type_p <<
" is not contemplated yet !" << endl ;
270 if ( type_t == SYM ) {
284 double theta = M_PI / 2 ;
293 cout <<
"Etoile::ray_eq_pi : the case type_t = " << type_t
294 <<
" and type_p = " << type_p << endl ;
295 cout <<
" is not contemplated yet !" << endl ;
313 double theta = M_PI / 2 ;
322 cout <<
"Etoile::ray_eq_pi : the case type_t = " << type_t
323 <<
" and type_p = " << type_p << endl ;
324 cout <<
" is not contemplated yet !" << endl ;
349 if ( type_t == SYM ) {
352 double theta = M_PI / 2 ;
353 double phi = 3. * M_PI / 2 ;
363 assert( np % 4 == 0 ) ;
372 cout <<
"Etoile::ray_eq_3pis2 : the case type_p = " 373 << type_p <<
" is not contemplated yet !" << endl ;
382 double theta = M_PI / 2 ;
383 double phi = 3. * M_PI / 2 ;
393 assert( np % 4 == 0 ) ;
402 cout <<
"Etoile::ray_eq_3pis2 : the case type_p = " 403 << type_p <<
" is not contemplated yet !" << endl ;
426 assert( (type_t == SYM) || (type_t == NONSYM) ) ;
456 if ( type_t == SYM ) {
459 double theta = M_PI / 2 ;
464 cout <<
"Etoile::ray_eq(kk) : the case type_p = " 465 << type_p <<
" is not contemplated yet !" << endl ;
470 double phi = 2. * kk * M_PI / np ;
473 ray_eq_kk =
mp.
val_r(l,xi,theta,phi) ;
478 cout <<
"Etoile::ray_eq(kk) : the case type_p = " 479 << type_p <<
" is not contemplated yet !" << endl ;
488 double theta = M_PI / 2 ;
493 cout <<
"Etoile::ray_eq(kk) : the case type_p = " 494 << type_p <<
" is not contemplated yet !" << endl ;
499 double phi = 2. * kk * M_PI / np ;
502 ray_eq_kk =
mp.
val_r(l,xi,theta,phi) ;
507 cout <<
"Etoile::ray_eq(kk) : the case type_p = " 508 << type_p <<
" is not contemplated yet !" << endl ;
534 "Etoile::mass_b : in the relativistic case, the baryon mass" 536 "computation cannot be performed by the base class Etoile !" 560 "Etoile::mass_g : in the relativistic case, the gravitational mass" 562 "computation cannot be performed by the base class Etoile !" int get_type_p() const
Returns the type of sampling in the direction: SYM : : symmetry with respect to the transformatio...
const Tbl & xi_surf() const
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate...
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
double ray_eq_pis2() const
Coordinate radius at , [r_unit].
virtual double mass_g() const
Gravitational mass.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
double ray_eq() const
Coordinate radius at , [r_unit].
int get_type_t() const
Returns the type of sampling in the direction: SYM : : symmetry with respect to the equatorial pl...
Basic integer array class.
double * p_ray_pole
Coordinate radius at .
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.
Tenseur nbar
Baryon density in the fluid frame.
virtual double mass_b() const
Baryon mass.
Map & mp
Mapping associated with the star.
Tbl * p_xi_surf
Description of the stellar surface: 2-D Tbl containing the values of the radial coordinate on the su...
int get_etat() const
Returns the logical state.
double * p_ray_eq
Coordinate radius at , .
double * p_mass_g
Gravitational mass.
int nzet
Number of domains of *mp occupied by the star.
double * p_ray_eq_pis2
Coordinate radius at , .
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one. ...
double * p_ray_eq_3pis2
Coordinate radius at , .
double ray_pole() const
Coordinate radius at [r_unit].
double ray_eq_pi() const
Coordinate radius at , [r_unit].
Tenseur ent
Log-enthalpy (relativistic case) or specific enthalpy (Newtonian case)
virtual const Itbl & l_surf() const
Description of the stellar surface: returns a 2-D Itbl containing the values of the domain index l o...
double * p_ray_eq_pi
Coordinate radius at , .
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
double ray_eq_3pis2() const
Coordinate radius at , [r_unit].
double * p_mass_b
Baryon mass.
Itbl * p_l_surf
Description of the stellar surface: 2-D Itbl containing the values of the domain index l on the surfa...