29char blackhole_static_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Black_hole/blackhole_static.C,v 1.3 2014/10/13 08:52:46 j_novak Exp $" ;
57#include "utilitaires.h"
112 for (
int i=1; i<=3; i++)
129 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
141 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
147 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
153 r_are =
r_coord(neumann, first) ;
159 + cc*cc*
pow(mass/r_are/rr,4.)) *
sqrt(r_are) ;
181 for (
int i=1; i<=3; i++) {
182 shift.
set(i) = mass * mass * cc * ll(i) / rr / rr
188 for (
int i=1; i<=3; i++) {
const Scalar r_coord(bool neumann, bool first) const
Expresses the areal radial coordinate by that in spatially isotropic coordinates.
Map & mp
Mapping associated with the black hole.
Scalar lapconf
A function (lapse function * conformal factor) lapconf generated by the black hole.
Vector shift
Shift vector generated by the black hole.
void static_bh(bool neumann, bool first)
Sets the metric quantities to a spherical, static blach-hole analytic solution.
Scalar lapse
Lapse function generated by the black hole.
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
Scalar confo
Conformal factor generated by the black hole.
double mass_bh
Gravitational mass of BH.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Coord r
r coordinate centered on the grid
Tensor field of valence 0 (or component of a tensorial field).
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
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.
Scalar & set(int)
Read/write access to a component.
Cmp sqrt(const Cmp &)
Square root.
Cmp pow(const Cmp &, int)
Power .
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).
Standard units of space, time and mass.