30char star_equil_spher_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/star_equil_spher.C,v 1.15 2014/10/13 08:53:39 j_novak Exp $" ;
165 cout <<
"-----------------------------------------------" <<
endl ;
301 double*
bornes =
new double[nz+1] ;
304 for(
int l=0;
l<nz;
l++) {
320 for(
int l=0;
l<nz;
l++) {
366 source.std_spectral_base() ;
382 a_car = qq * qq / (
nn *
nn ) ;
414 a_car = qq * qq / (
nn *
nn ) ;
419 for (
int i=1;
i<=3;
i++){
436 <<
"Characteristics of the star obtained by Etoile::equilibrium_spher : "
438 <<
"-----------------------------------------------------------------"
442 cout <<
"Coordinate radius : " <<
ray / km <<
" km" <<
endl ;
448 cout <<
"Circumferential radius R : " <<
rcirc/km <<
" km" <<
endl ;
450 cout <<
"Gravitational mass M : " <<
mass_g()/msol <<
" Mo" <<
endl ;
460 source.std_spectral_base() ;
468 - 0.5 *
tmp.dsdr() *
tmp.dsdr() )
471 source.std_spectral_base() ;
478 cout <<
"Virial theorem GRV3 : " <<
endl ;
481 cout <<
" relative error : " << grv3 <<
endl ;
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Time evolution with partial storage (*** under development ***).
Radial mapping of rather general form.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
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_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Tensor field of valence 0 (or component of a tensorial field).
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
const Scalar & dsdr() const
Returns of *this .
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
Scalar ener
Total energy density in the fluid frame.
virtual double mass_g() const =0
Gravitational mass.
Scalar logn
Logarithm of the lapse N .
Scalar nn
Lapse function N .
virtual void equilibrium_spher(double ent_c, double precis=1.e-14, const Tbl *pent_limit=0x0)
Computes a spherical static configuration.
Scalar ener_euler
Total energy density in the Eulerian frame.
void equation_of_state()
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
Scalar s_euler
Trace of the stress scalar in the Eulerian frame.
Scalar press
Fluid pressure.
Vector u_euler
Fluid 3-velocity with respect to the Eulerian observer.
Map & mp
Mapping associated with the star.
int nzet
Number of domains of *mp occupied by the star.
virtual double mass_b() const =0
Baryon mass.
Class intended to describe valence-2 symmetric tensors.
Scalar & set(int)
Read/write access to a component.
Cmp sqrt(const Cmp &)
Square root.
Cmp exp(const Cmp &)
Exponential.
Tbl diffrel(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (norme version).
Tbl norme(const Cmp &)
Sums of the absolute values of all the values of the Cmp in each domain.
Standard units of space, time and mass.