28char strot_dirac_global_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/strot_dirac_global.C,v 1.13 2014/10/13 08:53:40 j_novak Exp $" ;
81#include "star_rot_dirac.h"
83#include "utilitaires.h"
100 dens.std_spectral_base() ;
127 dens.std_spectral_base() ;
154 phi_kill.set(3).std_spectral_base() ;
162 dens.std_spectral_base() ;
185 dens.std_spectral_base() ;
226 sou_m.std_spectral_base() ;
257 sou_q.std_spectral_base() ;
300 sou_q.std_spectral_base() ;
312 sou_m.std_spectral_base() ;
365 cout <<
"The polar circumferential radius is only well defined\n"
366 <<
"with np = 1!" <<
endl ;
376 for (
int j=0;
j<nt;
j++) {
377 double theta = (+tet)(0, 0,
j, 0) ;
382 for (
int i=0;
i<
nrz;
i++) {
387 rrr.std_spectral_base() ;
393 for (
int j=0;
j<nt;
j++) {
394 double theta = (+tet)(0, 0,
j, 0) ;
397 double grr = gam(1,1).get_spectral_va().val_point_jk(
ls,
xs,
j, 0) ;
398 double grt = gam(1,2).get_spectral_va().val_point_jk(
ls,
xs,
j, 0) ;
399 double gtt = gam(2,2).get_spectral_va().val_point_jk(
ls,
xs,
j, 0) ;
401 double dr =
drrr.get_spectral_va().val_point_jk(
ls,
xs,
j, 0) ;
Active physical coordinates and mapping derivatives.
Time evolution with partial storage (*** under development ***).
Coord tet
coordinate centered on the grid
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates 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.
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
virtual const Sym_tensor & cov() const
Read-only access to the covariant representation.
virtual const Connection & connect() const
Returns the connection.
virtual const Scalar & determinant() const
Returns the determinant.
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_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).
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
const Scalar & dsdt() const
Returns of *this .
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
const Scalar & dsdr() const
Returns of *this .
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
Valeur & set_spectral_va()
Returns va (read/write version)
const Valeur & get_spectral_va() const
Returns va (read only version)
void annule_hard()
Sets the Scalar to zero in a hard way.
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
virtual double mass_g() const
Gravitational mass.
virtual double ellipt() const
Ellipticity e.
virtual double grv3() const
Error on the virial identity GRV3.
virtual double angu_mom() const
Angular momentum.
double * p_grv3
Error on the virial identity GRV3.
virtual double r_circ() const
Circumferential equatorial radius.
double omega
Rotation angular velocity ([f_unit] )
double * p_r_circ
Circumferential equatorial radius.
virtual double tsw() const
Ratio T/W.
Vector j_euler
Momentum density 3-vector with respect to the Eulerian observer.
virtual double rp_circ() const
Circumferential polar radius.
double * p_grv2
Error on the virial identity GRV2.
virtual double mass_b() const
Baryonic mass.
double * p_angu_mom
Angular momentum.
double * p_rp_circ
Circumferential polar radius.
virtual double aplat() const
Flattening r_pole/r_eq.
virtual double grv2() const
Error on the virial identity GRV2.
Scalar ener
Total energy density in the fluid frame.
Scalar logn
Logarithm of the lapse N .
Scalar nn
Lapse function N .
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...
double * p_mass_b
Baryon mass.
Scalar nbar
Baryon density in the fluid frame.
Scalar ener_euler
Total energy density in the Eulerian frame.
const Tbl & xi_surf() const
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate ...
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
double ray_eq() const
Coordinate radius at , [r_unit].
Scalar s_euler
Trace of the stress scalar in the Eulerian frame.
double * p_mass_g
Gravitational mass.
Scalar press
Fluid pressure.
Map & mp
Mapping associated with the star.
int nzet
Number of domains of *mp occupied by the star.
double ray_pole() const
Coordinate radius at [r_unit].
Class intended to describe valence-2 symmetric tensors.
void coef() const
Computes the coeffcients of *this.
Tensor field of valence 1.
Cmp sqrt(const Cmp &)
Square root.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Standard units of space, time and mass.