29 char et_bin_bhns_extr_um_der_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_bhns_extr_um_der.C,v 1.4 2014/10/13 08:52:55 j_novak Exp $" ;
55 #include "et_bin_bhns_extr.h"
77 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
82 xx_cov.
set(0) = xx + sepa ;
88 xsr_cov = xx_cov / r_bh ;
93 xx_con.
set(0) = xx + sepa ;
99 xsr_con = xx_con / r_bh ;
103 msr = ggrav * mass / r_bh ;
110 lapse_bh2 = 1. / (1.+2.*msr) ;
141 for (
int i=0; i<3; i++) {
142 for (
int j=i; j<3; j++) {
144 - 2*msr() % xsr_con(i) % xsr_con(j) ;
149 tkij_comp = (double(2)/double(3)) *
pow(lapse_bh2, 3.) % msr
160 lapse_bh8 = 1. /
pow(1.+2.*msr, 4.) ;
167 *
pow(2.*msr()+3.*msr()*msr(), 2.) / 3.
168 /
nnn() /
nnn() / r_bh() / r_bh() ;
192 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
197 xx_cov.
set(0) = xx + sepa ;
203 msr = ggrav * mass / r_bh ;
210 tmp = 1. / (1. - 0.25*msr*msr) ;
227 d_beta_comp.
set(0) = 0.5*tmp()%msr()%msr()%xx_cov(0)/r_bh()/r_bh() ;
228 d_beta_comp.
set(1) = 0.5*tmp()%msr()%msr()%xx_cov(1)/r_bh()/r_bh() ;
229 d_beta_comp.
set(2) = 0.5*tmp()%msr()%msr()%xx_cov(2)/r_bh()/r_bh() ;
241 for (
int i=0; i<3; i++) {
242 for (
int j=i; j<3; j++) {
Coord z
z coordinate centered on the grid
Tenseur akcar_comp
Part of the scalar generated by shift_auto and shift_comp , i.e.
Tenseur a_car
Total conformal factor .
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one.
Tensor handling *** DEPRECATED : use class Tensor instead ***.
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 ...
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Tenseur_sym tkij_comp
Part of the extrinsic curvature tensor generated by shift_comp .
Map & mp
Mapping associated with the star.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Cmp pow(const Cmp &, int)
Power .
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
Tenseur nnn
Total lapse function.
Standard units of space, time and mass.
Active physical coordinates and mapping derivatives.
bool kerrschild
Indicator of the background metric: true for the Kerr-Shild metric, false for the conformally flat on...
void change_triad(const Base_vect &new_triad)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Tenseur d_logn_comp
Gradient of logn_comp (Cartesian components with respect to ref_triad )
Coord x
x coordinate centered on the grid
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Coord y
y coordinate centered on the grid
Tenseur d_beta_comp
Gradient of beta_comp (Cartesian components with respect to ref_triad )
virtual void del_deriv() const
Deletes all the derived quantities.
void set_std_base()
Set the standard spectal basis of decomposition for each component.