26char binary_anashift_xcts_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Binary_xcts/binary_anashift_xcts.C,v 1.3 2014/10/13 08:52:45 j_novak Exp $" ;
48#include "binary_xcts.h"
57 for (
int i=0; i<2; i++) {
89 tmp = - 6 * www / a0 * ( 1 - (mp.
r)*(mp.
r) / (3*a0*a0) ) ;
92 tmp_ext = - 4 * www / mp.
r ;
93 tmp_ext.
annule(0, nzet-1) ;
95 w_beta.
set(2) = tmp + tmp_ext ;
106 tmp = 2 * www / a0 * (mp.
y) * ( 1 - 3 * (mp.
r)*(mp.
r) / (5*a0*a0) ) ;
108 tmp_ext = 0.8 * www * a0*a0 * (mp.
sint) * (mp.
sinp)
110 tmp_ext.
annule(0, nzet-1) ;
112 khi_beta = tmp + tmp_ext ;
122 Tenseur x_d_w_temp (
et[i]->get_mp(),2,CON,
et[i]->get_mp().get_bvect_cart()) ;
125 for (
int j=0; j<3; j++)
126 for (
int k=0; k<3; k++)
127 x_d_w_temp.
set(j,k) = xdw_temp(k+1, j+1) ;
132 Vector xdw (
et[i]->get_mp(), CON,
et[i]->get_mp().get_bvect_cart()) ;
133 for (
int j=0; j<3; j++)
134 xdw.
set(j+1) = x_d_w(j) ;
void analytical_shift()
Sets some analytical template for the shift vector (via the members w_shift and khi_shift of the two ...
double separation() const
Returns the coordinate separation of the two stellar centers [r_unit].
Star_bin_xcts * et[2]
Array of the two stars (to perform loops on the stars): et[0] contains the address of star1 and et[1]...
double omega
Angular velocity with respect to an asymptotically inertial observer.
Base class for coordinate mappings.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Coord y
y coordinate centered on the grid
Coord r
r coordinate centered on the grid
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_nzone() const
Returns the number of domains.
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.
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
Vector & set_beta_auto()
Read/write of .
const Metric & get_flat() const
Return the flat metric defined on the mapping (Spherical components with respect to the mapping of th...
virtual double mass_g() const
Gravitational mass.
const Map & get_mp() const
Returns the mapping.
double ray_eq() const
Coordinate radius at , [r_unit].
int get_nzet() const
Returns the number of domains occupied by the star.
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
void dec_dzpuis()
dzpuis -= 1 ;
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Scalar & set(int)
Read/write access to a 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...
const Tensor & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
Tenseur skxk(const Tenseur &)
Contraction of the last index of (*this) with or , depending on the type of S .
Standard units of space, time and mass.