59 #include "et_bin_bhns_extr.h" 107 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
112 xx_con.
set(0) = xx + sepa ;
118 xsr_con = xx_con / r_bh ;
122 msr = ggrav * mass / r_bh ;
126 lapse_bh2 = 1. / (1.+2.*msr) ;
141 for (
int i=0; i<3; i++) {
142 for (
int j=0; j<3; j++) {
143 tmp1.
set(i, j) = -2.*lapse_bh2()%msr()%xsr_con(i)%xsr_con(j) ;
145 tmp2.
set(i, j) = -3.*lapse_bh2()%xsr_con(i)%xsr_con(j)
146 -4.*lapse_bh2()*msr()%xsr_con(i)%xsr_con(j) ;
148 tmp3.
set(i, j) = xsr_con(i)%shift_auto_local(j) ;
161 for (
int i=0; i<3; i++)
162 tmp4.
set() += xsr_con(i) % shift_auto_local(i) ;
175 for (
int i=0; i<3; i++) {
176 for (
int j=i; j<3; j++) {
178 + tmp1dn(i, j) + tmp1dn(j, i)
179 + 2.*lapse_bh2()%msr()/r_bh()%( tmp3(i, j) + tmp3(j, i)
180 + tmp4() % tmp2(i, j) )
181 -
double(2)/double(3) * tmp1(i, j)
182 * (divn() - lapse_bh2() % msr() / r_bh() % tmp4()) ;
185 * (divn() - lapse_bh2() % msr() / r_bh() % tmp4()) ;
197 xx_cov.
set(0) = xx + sepa ;
203 xsr_cov = xx_cov / r_bh ;
211 for (
int i=0; i<3; i++) {
212 for (
int j=0; j<3; j++) {
218 for (
int i=0; i<3; i++) {
219 for (
int j=0; j<3; j++) {
220 tmp5.
set(i, j) = 2.*msr()%xsr_cov(i)%xsr_cov(j) ;
224 for (
int i=0; i<3; i++) {
225 for (
int j=0; j<3; j++) {
226 for (
int k=0; k<3; k++) {
240 for (
int i=0; i<3; i++) {
241 for (
int j=0; j<3; j++) {
242 tmp7.
set() += tmp5(i,j) % tmp6(i,j) ;
252 for (
int i=0; i<3; i++) {
253 for (
int j=0; j<3; j++) {
263 for (
int i=0; i<3; i++) {
264 for (
int j=0; j<3; j++) {
307 for (
int i=0; i<3; i++) {
308 for (
int j=i; j<3; j++) {
325 for (
int i=0; i<3; i++) {
326 for (
int j=0; 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_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...
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.
Map & mp
Mapping associated with the star.
Cmp pow(const Cmp &, int)
Power .
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Tenseur_sym tkij_auto
Part of the extrinsic curvature tensor generated by shift_auto .
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 .
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
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 akcar_auto
Part of the scalar generated by shift_auto , i.e.
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)