58 #include "blackhole.h" 60 #include "utilitaires.h" 115 for (
int i=1; i<=3; i++)
132 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
144 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
150 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
156 r_are =
r_coord(neumann, first) ;
162 + cc*cc*
pow(mass/r_are/rr,4.)) *
sqrt(r_are) ;
184 for (
int i=1; i<=3; i++) {
185 shift.
set(i) = mass * mass * cc * ll(i) / rr / rr
191 for (
int i=1; i<=3; i++) {
Map & mp
Mapping associated with the black hole.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Cmp sqrt(const Cmp &)
Square root.
double mass_bh
Gravitational mass of BH.
Standard units of space, time and mass.
void static_bh(bool neumann, bool first)
Sets the metric quantities to a spherical, static blach-hole analytic solution.
Tensor field of valence 0 (or component of a tensorial field).
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
Scalar confo
Conformal factor generated by the black hole.
Cmp pow(const Cmp &, int)
Power .
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Scalar lapse
Lapse function generated by the black hole.
Vector shift
Shift vector generated by the black hole.
Scalar & set(int)
Read/write access to a component.
const Scalar r_coord(bool neumann, bool first) const
Expresses the areal radial coordinate by that in spatially isotropic coordinates. ...
Scalar lapconf
A function (lapse function * conformal factor) lapconf generated by the black hole.
Coord r
r coordinate centered on the grid