63 #include "hole_bhns.h" 65 #include "utilitaires.h" 78 if (p_mass_irr_bhns == 0x0) {
86 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
91 double mirr =
sqrt(a_ah/16./M_PI) / ggrav ;
93 p_mass_irr_bhns =
new double( mirr ) ;
97 return *p_mass_irr_bhns ;
106 const double& theta_i,
const int& nrk_phi,
107 const int& nrk_theta)
const {
115 double mass = ggrav *
mass_bh ;
136 ll.set(1) = st % cp ;
137 ll.
set(2) = st % sp ;
141 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
145 cout <<
"Not yet prepared!!!" << endl ;
159 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
171 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
177 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
191 nrk_phi, nrk_theta) ;
192 killing_spher.std_spectral_base() ;
199 killing_spher.std_spectral_base() ;
204 killing.set(1) = (killing_spher(2) * ct * cp - killing_spher(3) * sp)
206 killing.set(2) = (killing_spher(2) * ct * sp + killing_spher(3) * cp)
208 killing.set(3) = - killing_spher(2) * st / radius_ah ;
230 *
sqrt(1. - 2.*mass/r_are/rr + cc*cc*
pow(mass/r_are/rr,4.))
231 * (ll(1)*killing(1) + ll(2)*killing(2) + ll(3)*killing(3))
236 source_surf = source_1 + source_2 ;
244 double spin_angmom = 0.5 * spin / qpig ;
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).
double integrale_surface(const Cmp &ci, double rayon) const
Performs the surface integration of ci on the sphere of radius rayon .
Cmp sqrt(const Cmp &)
Square root.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
double mass_bh
Gravitational mass of BH.
Standard units of space, time and mass.
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.
double * p_spin_am_bhns
Irreducible mass of BH.
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
virtual double val_r(int l, double xi, double theta, double pphi) const =0
Returns the value of the radial coordinate r for a given in a given domain.
Vector killing_vect(const Tbl &xi_i, const double &phi_i, const double &theta_i, const int &nrk_phi, const int &nrk_theta) const
Compute the Killing vector of a black hole normalized so that its affine length is 2 M_PI...
Cmp pow(const Cmp &, int)
Power .
Scalar confo_tot
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.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Sym_tensor taij_tot_rs
Part of the extrinsic curvature tensor from the numerical computation.
virtual double mass_irr_bhns() const
Irreducible mass of the black hole.
bool bc_lapconf_fs
true for the first type BC for the lapconf function, false for the second type BH ...
bool bc_lapconf_nd
true for the Neumann type BC for the lapconf function, false for the Dirichlet type BH ...
const Scalar r_coord(bool neumann, bool first) const
Expresses the areal radial coordinate by that in spatially isotropic coordinates. ...
double spin_am_bhns(const Tbl &xi_i, const double &phi_i, const double &theta_i, const int &nrk_phi, const int &nrk_theta) const
Spin angular momentum of the black hole.
virtual void dec_dzpuis(int dec=1)
Decreases by dec units the value of dzpuis and changes accordingly the values of the Scalar in the co...
Scalar lapconf_tot
Total lapconf function.
Coord r
r coordinate centered on the grid