78 for (
int i=0; i<2; i++) {
110 tmp = - 6 * www / a0 * ( 1 - (mp.
r)*(mp.
r) / (3*a0*a0) ) ;
113 tmp_ext = - 4 * www / mp.
r ;
114 tmp_ext.
annule(0, nzet-1) ;
116 w_beta.set(2) = tmp + tmp_ext ;
127 tmp = 2 * www / a0 * (mp.
y) * ( 1 - 3 * (mp.
r)*(mp.
r) / (5*a0*a0) ) ;
129 tmp_ext = 0.8 * www * a0*a0 * (mp.
sint) * (mp.
sinp)
131 tmp_ext.annule(0, nzet-1) ;
133 khi_beta = tmp + tmp_ext ;
144 Tenseur x_d_w_temp (
et[i]->get_mp(),2,CON,
et[i]->get_mp().get_bvect_cart()) ;
146 for (
int j=0; j<3; j++)
147 for (
int k=0; k<3; k++)
148 x_d_w_temp.
set(j,k) = xdw_temp(k+1, j+1) ;
153 Vector xdw (
et[i]->get_mp(), CON,
et[i]->get_mp().get_bvect_cart()) ;
154 for (
int j=0; j<3; j++)
155 xdw.
set(j+1) = x_d_w(j) ;
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
void analytical_shift()
Sets some analytical template for the shift vector (via the members w_shift and khi_shift of the two ...
Standard units of space, time and mass.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
double omega
Angular velocity with respect to an asymptotically inertial observer.
Tensor field of valence 0 (or component of a tensorial field).
Base class for coordinate mappings.
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
const Map & get_mp() const
Returns the mapping.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
virtual void dec_dzpuis(int dec=1)
Decreases by dec units the value of dzpuis and changes accordingly the values in the compactified ext...
double separation() const
Returns the coordinate separation of the two stellar centers [r_unit].
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Tenseur skxk(const Tenseur &)
Contraction of the last index of (*this) with or , depending on the type of S .
const Metric & get_flat() const
Return the flat metric defined on the mapping (Spherical components with respect to the mapping of th...
int get_nzone() const
Returns the number of domains.
Star_bin * et[2]
Array of the two stars (to perform loops on the stars): et[0] contains the address of star1 and et[1]...
const Tensor & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
void dec_dzpuis()
dzpuis -= 1 ;
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
int get_nzet() const
Returns the number of domains occupied by the star.
Coord y
y coordinate centered on the grid
double ray_eq() const
Coordinate radius at , [r_unit].
virtual double mass_g() const
Gravitational mass.
Scalar & set(int)
Read/write access to a component.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Coord r
r coordinate centered on the grid
Vector & set_beta_auto()
Read/write of .