29char et_bin_bhns_extr_upmetr_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_bhns_extr_upmetr.C,v 1.4 2014/10/13 08:52:55 j_novak Exp $" ;
55#include "et_bin_bhns_extr.h"
82 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
87 xx_con.
set(0) = xx + sepa ;
93 xsr_con = xx_con / r_bh ;
97 msr = ggrav * mass / r_bh ;
101 lapse_bh = 1. /
sqrt( 1.+2.*msr ) ;
115 shift_comp.
set(0) = -2.*lapse_bh()*lapse_bh()*msr()*xsr_con(0) ;
116 shift_comp.
set(1) = -2.*lapse_bh()*lapse_bh()*msr()*xsr_con(1) ;
117 shift_comp.
set(2) = -2.*lapse_bh()*lapse_bh()*msr()*xsr_con(2) ;
186 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
190 msr = ggrav * mass / r_bh ;
194 lapse_bh = (1.-0.5*msr) / (1.+0.5*msr) ;
202 lappsi = 1. - 0.25*msr*msr ;
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_extr(const double &mass, const double &sepa)
Computes metric coefficients from known potentials, when the companion is a black hole with the Kerr-...
void extrinsic_curv_extr(const double &mass, const double &sepa)
Computes tkij_auto and akcar_auto from shift_auto , nnn and a_car .
Tenseur shift_comp
Part of the shift vector generated principaly by the companion star.
Tenseur d_beta_auto
Gradient of beta_auto (Cartesian components with respect to ref_triad )
Tenseur d_logn_auto
Gradient of logn_auto (Cartesian components with respect to ref_triad )
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
Tenseur logn_comp
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principaly by ...
Tenseur beta_comp
Part of the logarithm of AN generated principaly by the companion star.
virtual void del_deriv() const
Deletes all the derived quantities.
Tenseur shift_auto
Part of the shift vector generated principaly by the star.
Tenseur d_logn_auto_regu
Gradient of logn_auto_regu (Cartesian components with respect to ref_triad )
Tenseur logn_auto_regu
Regular part of the logarithm of the part of the lapse N generated principaly by the star.
Tenseur nnn
Total lapse function.
Tenseur logn_auto
Total of the logarithm of the part of the lapse N generated principaly by the star.
Map & mp
Mapping associated with the star.
Tenseur d_logn_auto_div
Gradient of logn_auto_div (if k_div!=0 )
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one.
Tenseur shift
Total shift vector.
Tenseur beta_auto
Logarithm of the part of the product AN generated principaly by by the star.
Tenseur a_car
Total conformal factor .
double unsurc2
: unsurc2=1 for a relativistic star, 0 for a Newtonian one.
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 ***.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
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.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Cmp sqrt(const Cmp &)
Square root.
Cmp exp(const Cmp &)
Exponential.
Cmp pow(const Cmp &, int)
Power .
Cmp log(const Cmp &)
Neperian logarithm.
Standard units of space, time and mass.