58 #include "et_bin_bhns_extr.h" 80 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
85 xx_cov.
set(0) = xx + sepa ;
91 xsr_cov = xx_cov / r_bh ;
96 xx_con.
set(0) = xx + sepa ;
102 xsr_con = xx_con / r_bh ;
106 msr = ggrav * mass / r_bh ;
113 lapse_bh2 = 1. / (1.+2.*msr) ;
144 for (
int i=0; i<3; i++) {
145 for (
int j=i; j<3; j++) {
147 - 2*msr() % xsr_con(i) % xsr_con(j) ;
152 tkij_comp = (double(2)/double(3)) *
pow(lapse_bh2, 3.) % msr
163 lapse_bh8 = 1. /
pow(1.+2.*msr, 4.) ;
170 *
pow(2.*msr()+3.*msr()*msr(), 2.) / 3.
171 /
nnn() /
nnn() / r_bh() / r_bh() ;
195 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
200 xx_cov.
set(0) = xx + sepa ;
206 msr = ggrav * mass / r_bh ;
213 tmp = 1. / (1. - 0.25*msr*msr) ;
230 d_beta_comp.
set(0) = 0.5*tmp()%msr()%msr()%xx_cov(0)/r_bh()/r_bh() ;
231 d_beta_comp.
set(1) = 0.5*tmp()%msr()%msr()%xx_cov(1)/r_bh()/r_bh() ;
232 d_beta_comp.
set(2) = 0.5*tmp()%msr()%msr()%xx_cov(2)/r_bh()/r_bh() ;
244 for (
int i=0; i<3; i++) {
245 for (
int j=i; j<3; j++) {
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
void set_std_base()
Set the standard spectal basis of decomposition for each component.
Standard units of space, time and mass.
Tenseur nnn
Total lapse function.
bool kerrschild
Indicator of the background metric: true for the Kerr-Shild metric, false for the conformally flat on...
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Tenseur d_beta_comp
Gradient of beta_comp (Cartesian components with respect to ref_triad )
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.
void update_metric_der_comp_extr(const double &mass, const double &sepa)
Computes the derivative of metric functions related to the companion black hole with the Kerr-Schild ...
Tenseur akcar_comp
Part of the scalar generated by shift_auto and shift_comp , i.e.
Map & mp
Mapping associated with the star.
Cmp pow(const Cmp &, int)
Power .
Active physical coordinates and mapping derivatives.
Tenseur a_car
Total conformal factor .
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).
Tenseur d_logn_comp
Gradient of logn_comp (Cartesian components with respect to ref_triad )
Coord z
z coordinate centered on the grid
Tenseur_sym tkij_comp
Part of the extrinsic curvature tensor generated by shift_comp .
Tensor handling *** DEPRECATED : use class Tensor instead ***.