29char et_bin_bhns_extr_kinema_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_bhns_extr_kinema.C,v 1.3 2014/10/13 08:52:55 j_novak Exp $" ;
49#include "et_bin_bhns_extr.h"
77 bsn.
set(1) = - omega * xa ;
99 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
104 xx_cov.
set(0) = xx + sepa ;
110 xsr_cov = xx_cov / r_bh ;
114 msr = ggrav * mass / r_bh ;
127 for (
int i=0; i<3; i++)
128 tmp1.
set() += xsr_cov(i) %
bsn(i) ;
132 Tenseur tmp2 = 2.*msr % tmp1 % tmp1 ;
135 for (
int i=0; i<3; i++)
138 tmp2 =
a_car % tmp2 ;
146 cout <<
"BH-NS binary system should be relativistic !!!" << endl ;
176 bsn.
set(1) = - omega * xa ;
202 for (
int i=0; i<3; i++)
213 cout <<
"BH-NS binary system should be relativistic !!!" << endl ;
Active physical coordinates and mapping derivatives.
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 bsn
3-vector shift, divided by N , of the rotating coordinates, .
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
virtual void del_deriv() const
Deletes all the derived quantities.
Tenseur pot_centri
Centrifugal potential.
Tenseur nnn
Total lapse function.
Map & mp
Mapping associated with the star.
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one.
Tenseur shift
Total shift vector.
Tenseur a_car
Total conformal factor .
Coord y
y coordinate centered on the grid
Coord ya
Absolute y coordinate.
Coord z
z coordinate centered on the grid
Coord x
x coordinate centered on the grid
Coord xa
Absolute x coordinate.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_nzone() const
Returns the number of domains.
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 annule(int l)
Sets the Tenseur to zero in a given domain.
void set_std_base()
Set the standard spectal basis of decomposition for each component.
Cmp sqrt(const Cmp &)
Square root.
Cmp pow(const Cmp &, int)
Power .
Cmp log(const Cmp &)
Neperian logarithm.
Standard units of space, time and mass.