Lorene::Isol_hole Class Reference

Class to compute quasistationary single black hole spacetimes in vacuum. More...

#include <isol_hole.h>

List of all members.

Public Member Functions

 Isol_hole (const Map &mp_i, double Omega_i, bool NorKappa_i, Scalar NoK_i, bool isCF_i=false)
 Standard constructor.
 Isol_hole (const Isol_hole &)
 Copy constructor.
 Isol_hole (const Map &mp_i, double Omega_i, bool NorKappa_i, Scalar NoK_i, bool isCF_i, FILE *fich)
 Constructor from a file (see sauve(FILE* )).
virtual ~Isol_hole ()
 Destructor.
void operator= (const Isol_hole &)
 Assignment to another Isol_hole.
void compute_stat_metric (double precis, double relax, int mer_max, int mer_max2, bool isvoid=true)
 Computes a quasi-stationary 3-slice from the chosen parameters.
void secmembre_kerr (Sym_tensor &source_hh)
 Computes the rhs of hyperbolic equation for conformal metric assuming statioarity; WARNING; up to now, we are only able to handle void spacetimes.
const Mapget_mp () const
 Returns the mapping.
double get_Omega () const
 Returns the rotation rate.
const Scalarget_boundN () const
 Returns the boundary value used for the lapse (if it is the one used).
const Scalarget_Kappa () const
 Returns the surface gravity value at the boundary (if it is the one used).
const Scalarget_lapse () const
 Returns the lapse function N.
const Scalarget_conf_fact () const
 Returns the conformal factor.
const Vectorget_shift () const
 Returns the shift vector $\beta^i$.
const Sym_tensorget_hij () const
 Returns the deviation tensor $ h^{ij} $.
const Sym_tensorget_hatA () const
 Returns the rescaled tracefree extrinsic curvature $\hat{A}^{ij}$.
virtual void sauve (FILE *) const
 Save in a file.
void Einstein_errors ()
 Prints out errors in Einstein equations for the data obtained.
Spheroid hor ()
 Spheroid at the excised boundary associated with the black hole MOTS on the slice.
double adm_mass ()
 Computation of the ADM mass of the BH spacetime.
double komar_angmom ()
 Computation of the Komar angular momentum w.r.t.
double virial_residue ()
 Computation of the Virial residual, as difference at infinity between the ADM mass and the Komar integral associated to the mass.

Protected Member Functions

virtual void del_deriv () const
 Deletes all the derived quantities.
virtual void set_der_0x0 () const
 Sets to 0x0 all the pointers on derived quantities.

Protected Attributes

const Mapmp
 Mapping associated with the star.
double Omega
bool NorKappa
 Rotation rate of the horizon in the azimuthal direction.
Scalar boundNoK
 Indicates if the boundary value for the lapse or the surface gravity is used.
bool isCF
 Stores the boundary value of the lapse or surface gravity.
Scalar lapse
 Lapse function.
Scalar conf_fact
Vector shift
 Shift vector.
Sym_tensor hij
 Deviation tensor( non-flat part of the conformal 3-metric on the slice; see Bonazzola et al.
Sym_tensor hatA
 Rescaled tracefree extrinsic curvature tensor: see Cordero et al.
Spheroidp_hor
 Computation of the spheroid associated with the black hole horizon.
double * p_adm_mass
 Computation of the ADM mass of the BH spacetime.
double * p_komar_angmom
 Computation of the Komar angular momentum w.r.t.
double * p_virial_residue
 Computation of the Virial residual, as difference at infinity between the ADM mass and the Komar integral associated to the mass.

Detailed Description

Class to compute quasistationary single black hole spacetimes in vacuum.

It takes as arguments:

The main goal of this class is the computation of metric data from those parameters, as well as some global quantities related to the physical characteristics of the spacetime. Those metric data are given on a spacelike 3-slice, assuming a global timelike Killing field and using adapted coordinates.

Definition at line 50 of file isol_hole.h.


Constructor & Destructor Documentation

Lorene::Isol_hole::Isol_hole ( const Map mp_i,
double  Omega_i,
bool  NorKappa_i,
Scalar  NoK_i,
bool  isCF_i = false 
)

Standard constructor.

Parameters:
mp_i Mapping on which the black hole slice will be defined
Omega_i rotation rate of the horizon
NorKappa_i FALSE: use of boundary value for the lapse TRUE use of boundary value for the surface gravity (Warning! not implemented yet!)
NoK_i value either for the lapse or surface gravity.
isCF_i FALSE: Full GR 3+1 equations to be verified TRUE: IWM approximation used in determination of the spacetime geometry;

Definition at line 51 of file isol_hole.C.

References boundNoK, Lorene::Tensor::get_mp(), hatA, hij, lapse, set_der_0x0(), Lorene::Tensor::set_etat_zero(), and shift.

Lorene::Isol_hole::Isol_hole ( const Isol_hole ih  ) 

Copy constructor.

Definition at line 83 of file isol_hole.C.

References set_der_0x0().

Lorene::Isol_hole::Isol_hole ( const Map mp_i,
double  Omega_i,
bool  NorKappa_i,
Scalar  NoK_i,
bool  isCF_i,
FILE *  fich 
)

Constructor from a file (see sauve(FILE* )).

Parameters:
mp_i Mapping on which the star will be defined
Omega_i rotation rate of the horizon
NorKappa_i FALSE: use of boundary value for the lapse TRUE use of boundary value for the surface gravity (Warning! not implemented yet!)
NoK_i value either for the lapse or surface gravity.
isCF_i FALSE: Full GR 3+1 equations to be verified TRUE: IWM approximation used in determination of the spacetime geometry;
fich input file (must have been created by the function sauve)

Definition at line 101 of file isol_hole.C.

References set_der_0x0().

Lorene::Isol_hole::~Isol_hole (  )  [virtual]

Destructor.

Definition at line 125 of file isol_hole.C.

References del_deriv().


Member Function Documentation

double Lorene::Isol_hole::adm_mass (  ) 
void Lorene::Isol_hole::compute_stat_metric ( double  precis,
double  relax,
int  mer_max,
int  mer_max2,
bool  isvoid = true 
)

Computes a quasi-stationary 3-slice from the chosen parameters.

Parameters:
precis [input] threshold in the relative difference between the radial shift for two consecutive steps to stop the iterative procedure (default value: 1.e-11)

Definition at line 50 of file Isol_hole/isol_hole_compute_metric.C.

References Lorene::Scalar::allocate_all(), Lorene::Tensor::annule(), Lorene::Tensor::annule_domain(), Lorene::Scalar::annule_hard(), boundNoK, Lorene::Valeur::c_cf, Lorene::Valeur::coef(), Lorene::Metric::con(), Lorene::Metric_flat::con(), Lorene::contract(), Lorene::Metric::cov(), Lorene::Scalar::dec_dzpuis(), Lorene::Tensor::dec_dzpuis(), Lorene::Scalar::derive_con(), Lorene::Tensor::derive_con(), Lorene::Tensor_sym::derive_con(), Lorene::Tensor::derive_cov(), Lorene::Scalar::derive_cov(), Lorene::Tensor_sym::derive_cov(), Lorene::Sym_tensor::derive_lie(), Lorene::Metric::determinant(), Lorene::Sym_tensor::divergence(), Lorene::Vector::divergence(), Lorene::Scalar::dsdr(), Lorene::Vector::eta(), Lorene::Scalar::get_etat(), Lorene::Scalar::get_spectral_va(), hatA, hij, Lorene::Tensor::inc_dzpuis(), Lorene::Scalar::inc_dzpuis(), isCF, Lorene::Scalar::laplacian(), lapse, Lorene::log(), Lorene::max(), Lorene::maxabs(), mp, Lorene::Vector::mu(), Lorene::Scalar::mult_rsint(), NorKappa, secmembre_kerr(), Lorene::Vector::set(), Lorene::Tensor::set(), Lorene::Scalar::set_dzpuis(), Lorene::Scalar::set_etat_qcq(), Lorene::Scalar::set_grid_point(), Lorene::Scalar::set_spectral_va(), shift, Lorene::Scalar::sol_elliptic_boundary(), Lorene::sqrt(), Lorene::Vector::std_spectral_base(), Lorene::Scalar::std_spectral_base(), Lorene::Tensor::trace(), Lorene::Tensor::up(), Lorene::Tensor::up_down(), Lorene::Scalar::val_grid_point(), and Lorene::Valeur::ylm().

void Lorene::Isol_hole::del_deriv (  )  const [protected, virtual]

Deletes all the derived quantities.

Definition at line 136 of file isol_hole.C.

References p_adm_mass, p_hor, p_komar_angmom, p_virial_residue, and set_der_0x0().

void Lorene::Isol_hole::Einstein_errors (  ) 
const Scalar& Lorene::Isol_hole::get_boundN (  )  const [inline]

Returns the boundary value used for the lapse (if it is the one used).

Definition at line 202 of file isol_hole.h.

References boundNoK, and NorKappa.

const Scalar& Lorene::Isol_hole::get_conf_fact (  )  const [inline]

Returns the conformal factor.

Definition at line 222 of file isol_hole.h.

const Sym_tensor& Lorene::Isol_hole::get_hatA (  )  const [inline]

Returns the rescaled tracefree extrinsic curvature $\hat{A}^{ij}$.

Definition at line 231 of file isol_hole.h.

References hatA.

const Sym_tensor& Lorene::Isol_hole::get_hij (  )  const [inline]

Returns the deviation tensor $ h^{ij} $.

Definition at line 228 of file isol_hole.h.

References hij.

const Scalar& Lorene::Isol_hole::get_Kappa (  )  const [inline]

Returns the surface gravity value at the boundary (if it is the one used).

Definition at line 210 of file isol_hole.h.

References boundNoK, and NorKappa.

const Scalar& Lorene::Isol_hole::get_lapse (  )  const [inline]

Returns the lapse function N.

Definition at line 219 of file isol_hole.h.

References lapse.

const Map& Lorene::Isol_hole::get_mp (  )  const [inline]

Returns the mapping.

Definition at line 196 of file isol_hole.h.

References mp.

double Lorene::Isol_hole::get_Omega (  )  const [inline]

Returns the rotation rate.

Definition at line 199 of file isol_hole.h.

const Vector& Lorene::Isol_hole::get_shift (  )  const [inline]

Returns the shift vector $\beta^i$.

Definition at line 225 of file isol_hole.h.

References shift.

Spheroid Lorene::Isol_hole::hor (  ) 

Spheroid at the excised boundary associated with the black hole MOTS on the slice.

Set by default at the position $ r=1 $.

Definition at line 247 of file isol_hole.C.

References Lorene::Metric_flat::con(), hatA, hij, mp, Lorene::pow(), Lorene::Tensor::std_spectral_base(), Lorene::Scalar::std_spectral_base(), Lorene::Tensor::up_down(), and Lorene::Scalar::val_grid_point().

double Lorene::Isol_hole::komar_angmom (  ) 
void Lorene::Isol_hole::operator= ( const Isol_hole ih  ) 

Assignment to another Isol_hole.

Definition at line 160 of file isol_hole.C.

References boundNoK, del_deriv(), hatA, hij, isCF, lapse, mp, NorKappa, and shift.

void Lorene::Isol_hole::sauve ( FILE *  fich  )  const [virtual]

Save in a file.

Definition at line 183 of file isol_hole.C.

References hatA, hij, lapse, Lorene::Tensor_sym::sauve(), Lorene::Tensor::sauve(), Lorene::Scalar::sauve(), and shift.

void Lorene::Isol_hole::secmembre_kerr ( Sym_tensor source_hh  ) 
void Lorene::Isol_hole::set_der_0x0 (  )  const [protected, virtual]

Sets to 0x0 all the pointers on derived quantities.

Definition at line 145 of file isol_hole.C.

References p_adm_mass, p_hor, p_komar_angmom, and p_virial_residue.

double Lorene::Isol_hole::virial_residue (  ) 

Computation of the Virial residual, as difference at infinity between the ADM mass and the Komar integral associated to the mass.

Definition at line 334 of file isol_hole.C.

References Lorene::Scalar::asymptot(), Lorene::Map::get_mg(), Lorene::Mg3d::get_nzone(), lapse, and mp.


Member Data Documentation

Indicates if the boundary value for the lapse or the surface gravity is used.

Definition at line 65 of file isol_hole.h.

Rescaled tracefree extrinsic curvature tensor: see Cordero et al.

( * 2009).

Definition at line 92 of file isol_hole.h.

Deviation tensor( non-flat part of the conformal 3-metric on the slice; see Bonazzola et al.

(2004)).

Definition at line 87 of file isol_hole.h.

bool Lorene::Isol_hole::isCF [protected]

Stores the boundary value of the lapse or surface gravity.

Indicates if the CF approximation is used.

Definition at line 69 of file isol_hole.h.

Lapse function.

Definition at line 76 of file isol_hole.h.

const Map& Lorene::Isol_hole::mp [protected]

Mapping associated with the star.

Definition at line 55 of file isol_hole.h.

bool Lorene::Isol_hole::NorKappa [protected]

Rotation rate of the horizon in the azimuthal direction.

Definition at line 61 of file isol_hole.h.

double* Lorene::Isol_hole::p_adm_mass [mutable, protected]

Computation of the ADM mass of the BH spacetime.

Definition at line 103 of file isol_hole.h.

Spheroid* Lorene::Isol_hole::p_hor [mutable, protected]

Computation of the spheroid associated with the black hole horizon.

Definition at line 101 of file isol_hole.h.

double* Lorene::Isol_hole::p_komar_angmom [mutable, protected]

Computation of the Komar angular momentum w.r.t.

assumed rotational symmetry

Definition at line 108 of file isol_hole.h.

double* Lorene::Isol_hole::p_virial_residue [mutable, protected]

Computation of the Virial residual, as difference at infinity between the ADM mass and the Komar integral associated to the mass.

Definition at line 113 of file isol_hole.h.

Shift vector.

Definition at line 82 of file isol_hole.h.


The documentation for this class was generated from the following files:

Generated on 7 Dec 2019 for LORENE by  doxygen 1.6.1