27char etoile_global_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/etoile_global.C,v 1.7 2014/10/13 08:52:59 j_novak Exp $" ;
88 double precis = 1.e-15 ;
132 if ( type_t == SYM ) {
133 assert( (type_p == SYM) || (type_p == NONSYM) ) ;
138 double theta = M_PI / 2 ;
146 assert( (type_p == SYM) || (type_p == NONSYM) ) ;
151 double theta = M_PI / 2 ;
180 if ( type_t == SYM ) {
183 double theta = M_PI / 2 ;
184 double phi = M_PI / 2 ;
197 assert( np % 4 == 0 ) ;
206 cout <<
"Etoile::ray_eq_pis2 : the case type_p = "
207 << type_p <<
" is not contemplated yet !" << endl ;
216 double theta = M_PI / 2 ;
217 double phi = M_PI / 2 ;
230 assert( np % 4 == 0 ) ;
239 cout <<
"Etoile::ray_eq_pis2 : the case type_p = "
240 << type_p <<
" is not contemplated yet !" << endl ;
267 if ( type_t == SYM ) {
281 double theta = M_PI / 2 ;
290 cout <<
"Etoile::ray_eq_pi : the case type_t = " << type_t
291 <<
" and type_p = " << type_p << endl ;
292 cout <<
" is not contemplated yet !" << endl ;
310 double theta = M_PI / 2 ;
319 cout <<
"Etoile::ray_eq_pi : the case type_t = " << type_t
320 <<
" and type_p = " << type_p << endl ;
321 cout <<
" is not contemplated yet !" << endl ;
346 if ( type_t == SYM ) {
349 double theta = M_PI / 2 ;
350 double phi = 3. * M_PI / 2 ;
360 assert( np % 4 == 0 ) ;
369 cout <<
"Etoile::ray_eq_3pis2 : the case type_p = "
370 << type_p <<
" is not contemplated yet !" << endl ;
379 double theta = M_PI / 2 ;
380 double phi = 3. * M_PI / 2 ;
390 assert( np % 4 == 0 ) ;
399 cout <<
"Etoile::ray_eq_3pis2 : the case type_p = "
400 << type_p <<
" is not contemplated yet !" << endl ;
423 assert( (type_t == SYM) || (type_t == NONSYM) ) ;
453 if ( type_t == SYM ) {
456 double theta = M_PI / 2 ;
461 cout <<
"Etoile::ray_eq(kk) : the case type_p = "
462 << type_p <<
" is not contemplated yet !" << endl ;
467 double phi = 2. * kk * M_PI / np ;
470 ray_eq_kk =
mp.
val_r(l,xi,theta,phi) ;
475 cout <<
"Etoile::ray_eq(kk) : the case type_p = "
476 << type_p <<
" is not contemplated yet !" << endl ;
485 double theta = M_PI / 2 ;
490 cout <<
"Etoile::ray_eq(kk) : the case type_p = "
491 << type_p <<
" is not contemplated yet !" << endl ;
496 double phi = 2. * kk * M_PI / np ;
499 ray_eq_kk =
mp.
val_r(l,xi,theta,phi) ;
504 cout <<
"Etoile::ray_eq(kk) : the case type_p = "
505 << type_p <<
" is not contemplated yet !" << endl ;
531 "Etoile::mass_b : in the relativistic case, the baryon mass"
533 "computation cannot be performed by the base class Etoile !"
557 "Etoile::mass_g : in the relativistic case, the gravitational mass"
559 "computation cannot be performed by the base class Etoile !"
int nzet
Number of domains of *mp occupied by the star.
double * p_mass_b
Baryon mass.
double * p_mass_g
Gravitational mass.
Tbl * p_xi_surf
Description of the stellar surface: 2-D Tbl containing the values of the radial coordinate on the su...
double ray_eq_pi() const
Coordinate radius at , [r_unit].
double ray_eq() const
Coordinate radius at , [r_unit].
Tenseur nbar
Baryon density in the fluid frame.
Itbl * p_l_surf
Description of the stellar surface: 2-D Itbl containing the values of the domain index l on the surfa...
Map & mp
Mapping associated with the star.
double * p_ray_eq_3pis2
Coordinate radius at , .
double * p_ray_eq_pi
Coordinate radius at , .
virtual double mass_b() const
Baryon mass.
double ray_eq_pis2() const
Coordinate radius at , [r_unit].
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 ...
double ray_eq_3pis2() const
Coordinate radius at , [r_unit].
virtual double mass_g() const
Gravitational mass.
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one.
double * p_ray_eq_pis2
Coordinate radius at , .
Tenseur ent
Log-enthalpy (relativistic case) or specific enthalpy (Newtonian case)
const Tbl & xi_surf() const
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinat...
double * p_ray_eq
Coordinate radius at , .
double ray_pole() const
Coordinate radius at [r_unit].
Basic integer array class.
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.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_type_t() const
Returns the type of sampling in the direction: SYM : : symmetry with respect to the equatorial pl...
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
int get_type_p() const
Returns the type of sampling in the direction: SYM : : symmetry with respect to the transformatio...
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_etat() const
Returns the logical state.