30char et_bin_bhns_extr_excurv_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_bhns_extr_excurv.C,v 1.4 2014/10/13 08:52:55 j_novak Exp $" ;
56#include "et_bin_bhns_extr.h"
104 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
109 xx_con.
set(0) = xx + sepa ;
115 xsr_con = xx_con / r_bh ;
119 msr = ggrav * mass / r_bh ;
123 lapse_bh2 = 1. / (1.+2.*msr) ;
138 for (
int i=0; i<3; i++) {
139 for (
int j=0; j<3; j++) {
140 tmp1.
set(i, j) = -2.*lapse_bh2()%msr()%xsr_con(i)%xsr_con(j) ;
142 tmp2.
set(i, j) = -3.*lapse_bh2()%xsr_con(i)%xsr_con(j)
143 -4.*lapse_bh2()*msr()%xsr_con(i)%xsr_con(j) ;
145 tmp3.
set(i, j) = xsr_con(i)%shift_auto_local(j) ;
158 for (
int i=0; i<3; i++)
159 tmp4.
set() += xsr_con(i) % shift_auto_local(i) ;
172 for (
int i=0; i<3; i++) {
173 for (
int j=i; j<3; j++) {
175 + tmp1dn(i, j) + tmp1dn(j, i)
176 + 2.*lapse_bh2()%msr()/r_bh()%( tmp3(i, j) + tmp3(j, i)
177 + tmp4() % tmp2(i, j) )
178 -
double(2)/double(3) * tmp1(i, j)
179 * (divn() - lapse_bh2() % msr() / r_bh() % tmp4()) ;
182 * (divn() - lapse_bh2() % msr() / r_bh() % tmp4()) ;
194 xx_cov.
set(0) = xx + sepa ;
200 xsr_cov = xx_cov / r_bh ;
208 for (
int i=0; i<3; i++) {
209 for (
int j=0; j<3; j++) {
215 for (
int i=0; i<3; i++) {
216 for (
int j=0; j<3; j++) {
217 tmp5.
set(i, j) = 2.*msr()%xsr_cov(i)%xsr_cov(j) ;
221 for (
int i=0; i<3; i++) {
222 for (
int j=0; j<3; j++) {
223 for (
int k=0; k<3; k++) {
237 for (
int i=0; i<3; i++) {
238 for (
int j=0; j<3; j++) {
239 tmp7.
set() += tmp5(i,j) % tmp6(i,j) ;
249 for (
int i=0; i<3; i++) {
250 for (
int j=0; j<3; j++) {
260 for (
int i=0; i<3; i++) {
261 for (
int j=0; j<3; j++) {
304 for (
int i=0; i<3; i++) {
305 for (
int j=i; j<3; j++) {
322 for (
int i=0; i<3; i++) {
323 for (
int j=0; 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 extrinsic_curv_extr(const double &mass, const double &sepa)
Computes tkij_auto and akcar_auto from shift_auto , nnn and a_car .
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
Tenseur shift_auto
Part of the shift vector generated principaly by the star.
Tenseur_sym tkij_auto
Part of the extrinsic curvature tensor generated by shift_auto .
Tenseur akcar_auto
Part of the scalar generated by shift_auto , i.e.
Tenseur nnn
Total lapse function.
Map & mp
Mapping associated with the star.
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,...
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).
const Tenseur & gradient() const
Returns the gradient of *this (Cartesian coordinates)
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.
Cmp pow(const Cmp &, int)
Power .
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Standard units of space, time and mass.