28 char 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"
153 phi_kill.
set(3) = 1. ;
189 p_tsw =
new double( tcin / ( mass_p + tcin -
mass_g() ) ) ;
224 sou_m =
sqrt( k_det )*sou_m ;
248 sou_tmp += term_2 + term_3 ;
255 sou_q =
sqrt( k_det )*sou_q ;
259 p_grv2 =
new double(
double(1) + integrale2d(sou_m)/integrale2d(sou_q) ) ;
296 sou_q = sou_q + tmp_1 - tmp_2 ;
316 p_grv3 =
new double( (int_grav + int_mat) / int_mat ) ;
338 int i_b = mg->
get_nr(l_b) - 1 ;
342 double gamma_phi =
gamma.
cov()(3,3).val_grid_point(l_b, k_b, j_b, i_b) ;
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) ;
380 for (
int lz=0; lz<
nzet; lz++) {
382 for (
int i=0; i<nrz; i++) {
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) ;
400 double rr =
mp.
val_r(ls, xs, theta, phi) ;
404 =
sqrt(grr*dr*dr + 2*grt*rr*dr + gtt*rr*rr) ;
const Tensor_sym & get_delta() const
Returns the tensor which defines the connection with respect to the flat one: is the difference bet...
Active physical coordinates and mapping derivatives.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
Coord tet
coordinate centered on the grid
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.
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
void div_r_dzpuis(int ced_mult_r)
Division by r everywhere but with the output flag dzpuis set to ced_mult_r .
const Scalar & dsdt() const
Returns of *this .
double integrale() const
Computes the integral over all space of *this .
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
const Scalar & dsdr() const
Returns of *this .
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
void annule_hard()
Sets the Scalar to zero in a hard way.
void mult_rsint()
Multiplication by everywhere; dzpuis is not changed.
const Valeur & get_spectral_va() const
Returns va (read only version)
Valeur & set_spectral_va()
Returns va (read/write version)
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.
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...
Mtbl_cf * c_cf
Coefficients of the spectral expansion of the function.
void coef() const
Computes the coeffcients of *this.
Tensor field of valence 1.
Scalar & set(int)
Read/write access to a component.
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.