52 #include "et_bin_bhns_extr.h" 80 bsn.
set(1) = - omega * xa ;
102 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
107 xx_cov.
set(0) = xx + sepa ;
113 xsr_cov = xx_cov / r_bh ;
117 msr = ggrav * mass / r_bh ;
130 for (
int i=0; i<3; i++)
131 tmp1.
set() += xsr_cov(i) %
bsn(i) ;
135 Tenseur tmp2 = 2.*msr % tmp1 % tmp1 ;
138 for (
int i=0; i<3; i++)
141 tmp2 =
a_car % tmp2 ;
149 cout <<
"BH-NS binary system should be relativistic !!!" << endl ;
179 bsn.
set(1) = - omega * xa ;
205 for (
int i=0; i<3; i++)
216 cout <<
"BH-NS binary system should be relativistic !!!" << endl ;
Coord xa
Absolute x coordinate.
Cmp log(const Cmp &)
Neperian logarithm.
void annule(int l)
Sets the Tenseur to zero in a given domain.
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
Cmp sqrt(const Cmp &)
Square root.
void set_std_base()
Set the standard spectal basis of decomposition for each component.
Tenseur pot_centri
Centrifugal potential.
Standard units of space, time and mass.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tenseur nnn
Total lapse function.
void kinematics_extr(double omega, const double &mass, const double &sepa)
Computes the quantities bsn and pot_centri in the Kerr-Schild background metric or in the conformally...
bool kerrschild
Indicator of the background metric: true for the Kerr-Shild metric, false for the conformally flat on...
Tenseur shift
Total shift vector.
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
virtual void del_deriv() const
Deletes all the derived quantities.
Tenseur bsn
3-vector shift, divided by N , of the rotating coordinates, .
Map & mp
Mapping associated with the star.
int get_nzone() const
Returns the number of domains.
Cmp pow(const Cmp &, int)
Power .
Active physical coordinates and mapping derivatives.
Tenseur a_car
Total conformal factor .
Coord ya
Absolute y coordinate.
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one. ...
Coord y
y coordinate centered on the grid
Coord x
x coordinate centered on the grid
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Coord z
z coordinate centered on the grid
Tensor handling *** DEPRECATED : use class Tensor instead ***.