29char 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++) {
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 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.
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
Tenseur d_logn_comp
Gradient of logn_comp (Cartesian components with respect to ref_triad )
virtual void del_deriv() const
Deletes all the derived quantities.
Tenseur_sym tkij_comp
Part of the extrinsic curvature tensor generated by shift_comp .
Tenseur d_beta_comp
Gradient of beta_comp (Cartesian components with respect to ref_triad )
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 a_car
Total conformal factor .
Coord y
y coordinate centered on the grid
Coord z
z coordinate centered on the grid
Coord x
x coordinate centered on the grid
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 set_std_base()
Set the standard spectal basis of decomposition for each component.
void change_triad(const Base_vect &new_triad)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Cmp pow(const Cmp &, int)
Power .
Standard units of space, time and mass.