58 #include "et_bin_bhns_extr.h" 85 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
90 xx_con.
set(0) = xx + sepa ;
96 xsr_con = xx_con / r_bh ;
100 msr = ggrav * mass / r_bh ;
104 lapse_bh = 1. /
sqrt( 1.+2.*msr ) ;
118 shift_comp.
set(0) = -2.*lapse_bh()*lapse_bh()*msr()*xsr_con(0) ;
119 shift_comp.
set(1) = -2.*lapse_bh()*lapse_bh()*msr()*xsr_con(1) ;
120 shift_comp.
set(2) = -2.*lapse_bh()*lapse_bh()*msr()*xsr_con(2) ;
189 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
193 msr = ggrav * mass / r_bh ;
197 lapse_bh = (1.-0.5*msr) / (1.+0.5*msr) ;
205 lappsi = 1. - 0.25*msr*msr ;
Cmp log(const Cmp &)
Neperian logarithm.
Tenseur shift_comp
Part of the shift vector generated principaly by the companion star.
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
Cmp exp(const Cmp &)
Exponential.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Cmp sqrt(const Cmp &)
Square root.
void set_std_base()
Set the standard spectal basis of decomposition for each component.
Tenseur logn_auto_regu
Regular part of the logarithm of the part of the lapse N generated principaly by the star...
Standard units of space, time and mass.
Tenseur nnn
Total lapse function.
void update_metric_extr(const double &mass, const double &sepa)
Computes metric coefficients from known potentials, when the companion is a black hole with the Kerr-...
double unsurc2
: unsurc2=1 for a relativistic star, 0 for a Newtonian one.
bool kerrschild
Indicator of the background metric: true for the Kerr-Shild metric, false for the conformally flat on...
Tenseur d_beta_auto
Gradient of beta_auto (Cartesian components with respect to ref_triad )
Tenseur shift
Total shift vector.
Tenseur shift_auto
Part of the shift vector generated principaly by the star.
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
void change_triad(const Base_vect &new_triad)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
virtual void del_deriv() const
Deletes all the derived quantities.
Tenseur logn_comp
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principaly by ...
Tenseur d_logn_auto
Gradient of logn_auto (Cartesian components with respect to ref_triad )
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Tenseur beta_comp
Part of the logarithm of AN generated principaly by the companion star.
Map & mp
Mapping associated with the star.
Cmp pow(const Cmp &, int)
Power .
Active physical coordinates and mapping derivatives.
void extrinsic_curv_extr(const double &mass, const double &sepa)
Computes tkij_auto and akcar_auto from shift_auto , nnn and a_car .
Tenseur a_car
Total conformal factor .
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one. ...
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
Tenseur logn_auto
Total of the logarithm of the part of the lapse N generated principaly by the star.
Tenseur d_logn_auto_regu
Gradient of logn_auto_regu (Cartesian components with respect to ref_triad )
Coord y
y coordinate centered on the grid
Coord x
x coordinate centered on the grid
Tenseur beta_auto
Logarithm of the part of the product AN generated principaly by by the star.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tenseur d_logn_auto_div
Gradient of logn_auto_div (if k_div!=0 )
Coord z
z coordinate centered on the grid
Tensor handling *** DEPRECATED : use class Tensor instead ***.
const Tenseur & gradient() const
Returns the gradient of *this (Cartesian coordinates)