28char blackhole_global_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Black_hole/blackhole_global.C,v 1.5 2014/10/13 08:52:46 j_novak Exp $" ;
62#include "utilitaires.h"
83 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
88 double mirr =
sqrt(a_ah/16./M_PI) / ggrav ;
115 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
141 ll.
set(1) = st * cp ;
142 ll.
set(2) = st * sp ;
151 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
232 divshf =
shift(1).deriv(1) +
shift(2).deriv(2)
233 +
shift(3).deriv(3) ;
247 source_surf =
confo/rr
267 madm = integ_s / qpig + integ_v / qpig ;
273 cout <<
"Another ADM mass : " << mmm / msol << endl ;
301 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
327 ll.
set(1) = st * cp ;
328 ll.
set(2) = st * sp ;
334 for (
int i=1; i<=3; i++)
341 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
453 source_surf = lldlap ;
467 for (
int i=1; i<=3; i++)
474 - 2.*(dlap(1)%dlcnf(1) + dlap(2)%dlcnf(2) + dlap(3)%dlcnf(3)) ;
483 mkom = integ_s / qpig + integ_v / qpig ;
488 cout <<
"Another Komar mass : " << mmm / msol << endl ;
529 const Tbl& xi_i,
const double& phi_i,
530 const double& theta_i,
const int& nrk_phi,
531 const int& nrk_theta)
const {
554 ll.
set(1) = st % cp ;
555 ll.
set(2) = st % sp ;
559 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
563 cout <<
"Not yet prepared!!!" << endl ;
573 nrk_phi, nrk_theta) ;
576 killing_spher.
set(2) =
confo *
confo * radius_ah * killing_spher(2) ;
577 killing_spher.
set(3) =
confo *
confo * radius_ah * killing_spher(3) ;
585 killing.
set(1) = (killing_spher(2) * ct * cp - killing_spher(3) * sp)
587 killing.
set(2) = (killing_spher(2) * ct * sp + killing_spher(3) * cp)
589 killing.
set(3) = - killing_spher(2) * st / radius_ah ;
596 source_1 = (ll(1) * (
taij(1,1) * killing(1)
597 +
taij(1,2) * killing(2)
598 +
taij(1,3) * killing(3))
599 + ll(2) * (
taij(2,1) * killing(1)
600 +
taij(2,2) * killing(2)
601 +
taij(2,3) * killing(3))
602 + ll(3) * (
taij(3,1) * killing(1)
603 +
taij(3,2) * killing(2)
604 +
taij(3,3) * killing(3)))
610 source_surf = source_1 ;
618 double spin_angmom = 0.5 * spin / qpig ;
645 double integ_bh_s_x ;
646 double integ_bh_s_y ;
647 double integ_bh_s_z ;
649 double radius_ah =
mp.
val_r(1,-1.,M_PI/2.,0.) ;
677 ll.
set(1) = st % cp ;
678 ll.
set(2) = st % sp ;
684 for (
int i=1; i<=3; i++)
691 for (
int i=1; i<=3; i++)
698 for (
int i=1; i<=3; i++)
761 sou_bh_sx = jbh_x(1)%ll(1) + jbh_x(2)%ll(2) + jbh_x(3)%ll(3) ;
782 sou_bh_sy = jbh_y(1)%ll(1) + jbh_y(2)%ll(2) + jbh_y(3)%ll(3) ;
803 sou_bh_sz = jbh_z(1)%ll(1) + jbh_z(2)%ll(2) + jbh_z(3)%ll(3) ;
virtual double mass_kom() const
Komar mass.
double spin_am_bh(bool bclapconf_nd, bool bclapconf_fs, const Tbl &xi_i, const double &phi_i, const double &theta_i, const int &nrk_phi, const int &nrk_theta) const
Spin angular momentum.
Scalar taij_quad
Part of the scalar generated by .
Sym_tensor taij
Trace of the extrinsic curvature.
virtual double mass_irr() const
Irreducible mass of the black hole.
double * p_spin_am_bh
Radius of the apparent horizon.
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.
virtual double rad_ah() const
Radius of the apparent horizon.
Tbl * p_angu_mom_bh
Spin angular momentum.
Vector killing_vect_bh(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.
virtual double mass_adm() const
ADM mass.
double * p_rad_ah
Komar mass.
double * p_mass_kom
ADM mass.
const Tbl & angu_mom_bh() const
Total angular momentum.
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
Scalar confo
Conformal factor generated by the black hole.
double * p_mass_adm
Irreducible mass of the black hole.
double * p_mass_irr
Conformal metric .
double integrale_surface_infini(const Cmp &ci) const
Performs the surface integration of ci at infinity.
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 y
y coordinate centered on the grid
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.
Coord z
z coordinate centered on the grid
Coord x
x coordinate centered on the grid
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.
const Scalar & deriv(int i) const
Returns of *this , where .
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
double integrale() const
Computes the integral over all space of *this .
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
const Scalar & dsdr() const
Returns of *this .
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...
void annule_hard()
Sets the Tbl to zero in a hard way.
double & set(int i)
Read/write of a particular element (index i) (1D case)
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.