28char hole_bhns_global_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Hole_bhns/hole_bhns_global.C,v 1.5 2014/10/13 08:53:00 j_novak Exp $" ;
62#include "utilitaires.h"
75 if (p_mass_irr_bhns == 0x0) {
83 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
88 double mirr =
sqrt(a_ah/16./M_PI) / ggrav ;
90 p_mass_irr_bhns =
new double( mirr ) ;
94 return *p_mass_irr_bhns ;
103 const double& theta_i,
const int& nrk_phi,
104 const int& nrk_theta)
const {
112 double mass = ggrav *
mass_bh ;
133 ll.
set(1) = st % cp ;
134 ll.
set(2) = st % sp ;
138 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
142 cout <<
"Not yet prepared!!!" << endl ;
156 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
168 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
174 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
188 nrk_phi, nrk_theta) ;
201 killing.
set(1) = (killing_spher(2) * ct * cp - killing_spher(3) * sp)
203 killing.
set(2) = (killing_spher(2) * ct * sp + killing_spher(3) * cp)
205 killing.
set(3) = - killing_spher(2) * st / radius_ah ;
227 *
sqrt(1. - 2.*mass/r_are/rr + cc*cc*
pow(mass/r_are/rr,4.))
228 * (ll(1)*killing(1) + ll(2)*killing(2) + ll(3)*killing(3))
233 source_surf = source_1 + source_2 ;
241 double spin_angmom = 0.5 * spin / qpig ;
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.
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
double mass_bh
Gravitational mass of BH.
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.
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.
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.
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
Scalar lapconf_tot
Total lapconf function.
Scalar confo_tot
Total conformal factor.
double * p_spin_am_bhns
Irreducible mass of BH.
double integrale_surface(const Cmp &ci, double rayon) const
Performs the surface integration of ci on the sphere of radius rayon .
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
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
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.
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.
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...
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.