24char binhor_glob_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Bin_hor/binhor_global.C,v 1.10 2014/10/13 08:52:42 j_novak Exp $" ;
79#include "utilitaires.h"
125 assert ((orientation_un==0) || (orientation_un==M_PI)) ;
127 assert ((orientation_deux==0) || (orientation_deux==M_PI)) ;
128 int same_orient = (orientation_un == orientation_deux) ? 1 : -1 ;
143 for (
int i=1 ; i<=3 ; i++)
144 vecteur_un.
set(i) = -ya_un*tkij_un(1, i)+
145 xa_un * tkij_un(2, i) ;
167 for (
int i=1 ; i<=3 ; i++)
168 vecteur_deux.
set(i) = -ya_deux*tkij_deux(1, i)+
169 xa_deux * tkij_deux(2, i) ;
178 return moment_un+same_orient*moment_deux ;
200 assert ((orientation_un==0) || (orientation_un==M_PI)) ;
202 assert ((orientation_deux==0) || (orientation_deux==M_PI)) ;
203 int same_orient = (orientation_un == orientation_deux) ? 1 : -1 ;
207 double* bornes =
new double [nzones+1] ;
209 for (
int i=nzones-1 ; i>0 ; i--) {
210 bornes[i] = courant ;
214 bornes[nzones] = __infinity ;
235 shift_deux.
set(1).
import(same_orient*beta_deux(1)) ;
236 shift_deux.
set(2).
import(same_orient*beta_deux(2)) ;
239 Vector shift_tot (shift_un+shift_deux) ;
241 shift_tot.
annule(0, nzones-2) ;
252 for (
int lig=1 ; lig<=3 ; lig++)
253 for (
int col=lig ; col<=3 ; col++)
257 for (
int i=1 ; i<=3 ; i++)
258 vecteur_un.
set(i) = k_total(1, i) ;
260 Scalar integrant_un (vecteur_un(1)) ;
263 for (
int i=1 ; i<=3 ; i++)
264 vecteur_deux.
set(i) = k_total(2, i) ;
266 Scalar integrant_deux (vecteur_deux(1)) ;
double omega
Angular velocity.
Single_hor hole1
Black hole one.
Single_hor hole2
Black hole two.
double ang_mom_hor() const
Calculates the angular momentum of the black hole using the formula at the horizon.
double adm_mass() const
Calculates the ADM mass of the system.
double ang_mom_adm() const
Calculates the angular momentum of the black hole.
double komar_mass() const
Calculates the Komar mass of the system using : .
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 .
virtual double val_r(int l, double xi, double theta, double pphi) const
Returns the value of the radial coordinate r for a given in a given domain.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Coord ya
Absolute y coordinate.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
const Metric_flat & flat_met_cart() const
Returns the flat metric associated with the Cartesian coordinates and with components expressed in th...
double get_ori_x() const
Returns the x coordinate of the origin.
double get_rot_phi() const
Returns the angle between the x –axis and X –axis.
Coord xa
Absolute x coordinate.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Flat metric for tensor calculation.
int get_nzone() const
Returns the number of domains.
Tensor field of valence 0 (or component of a tensorial field).
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
Valeur & set_spectral_va()
Returns va (read/write version)
const Valeur & get_spectral_va() const
Returns va (read only version)
void mult_r_ced()
Multiplication by r in the compactified external domain (CED), the dzpuis flag is not changed.
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
Scalar psi_auto
Conformal factor .
double radius
Radius of the horizon in LORENE's units.
Vector beta_auto
Shift function .
const Scalar & get_psi4() const
Conformal factor .
Metric_flat ff
3 metric flat
Scalar psi
Conformal factor .
Map_af & mp
Affine mapping.
Scalar n_auto
Lapse function .
Scalar nn
Lapse function .
Sym_tensor hh
Deviation metric.
Sym_tensor aa
Components of the conformal representation of the traceless part of the extrinsic curvature:
Class intended to describe valence-2 symmetric tensors.
const Valeur & mult_st() const
Returns applied to *this.
const Valeur & mult_cp() const
Returns applied to *this.
const Valeur & mult_sp() const
Returns applied to *this.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
double flux(double radius, const Metric &met) const
Computes the flux of the vector accross a sphere r = const.
Sym_tensor ope_killing_conf(const Metric &gam) const
Computes the conformal Killing operator associated with a given metric.
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Scalar & set(int)
Read/write access to a component.
virtual void dec_dzpuis(int dec=1)
Decreases by dec units the value of dzpuis and changes accordingly the values in the compactified ext...
virtual void change_triad(const Base_vect &new_triad)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
virtual void annule(int l_min, int l_max)
Sets the Tensor to zero in several domains.
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values in the compactified ext...
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
Tensor trace(int ind1, int ind2) const
Trace on two different type indices.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .