|
LORENE
|
Class describing an apparent horizon. More...
#include <spheroid.h>
Public Member Functions | |
| App_hor (const Mg3d &grid_angu, double radius) | |
| Standard constructor. | |
| App_hor (const Scalar &h_in, const Metric &gamij, const Sym_tensor &Kij) | |
Constructor of an apparent horizon embedded in a 3-slice (Time_slice ) of 3+1 formalism. | |
| App_hor (const App_hor &) | |
| Copy constructor. | |
| App_hor (FILE *) | |
Constructor from a file (see sauve(FILE*) ) | |
| virtual | ~App_hor () |
| Destructor. | |
| void | operator= (const App_hor &) |
| Assignment to another App_hor. | |
| bool | check_expansion (double thres=1.e-7) const |
| const Sym_tensor & | lie_derive_shear (const Scalar &bb, const Scalar &lapse) |
| Lie derivative of shear tensor with respect to the evolution vector field. | |
| const Sym_tensor & | lie_derive_theta_plus (const Scalar &bb, const Scalar &lapse) |
| Lie derivative of the null outgoing expansion rate with respect to the evolution vector field. | |
| const Sym_tensor & | lie_derive_theta_minus (const Scalar &bb, const Scalar &lapse) |
| Lie derivative of the null ingoing expansion rate with respect to the evolution vector field. | |
| const Sym_tensor & | lie_derive_q_ab (const Scalar &bb, const Scalar &lapse) |
| Lie derivative of 2-metric with respect to the evolution vector field. | |
| const Scalar & | l_non_affinity (const Scalar &bb, const Scalar &lapse) |
| non-affinity (or surface gravity) with respect to the outgoing null vector field | |
| void | set_ephi (const Scalar &) |
| Assigns the conformal Killing vector field to phi. | |
| const Scalar & | get_hsurf () const |
Returns the field h_surf. | |
| const Metric & | get_qab () const |
Returns the metric ![]() | |
| const Scalar & | get_ricci () const |
Returns the 2-ricci scalar ![]() | |
| const Sym_tensor & | get_hh () const |
Returns the symmetric tensor ![]() | |
| const Sym_tensor & | get_qq () const |
returns the 3-d degenerate 2-metric ![]() | |
| const Tensor & | get_proj () const |
returns the 3-d projector on 2-surface ![]() | |
| const Tensor & | get_jac2d () const |
returns the 2-d jacobian of coordinate transformation ![]() | |
| const Scalar & | get_trk () const |
Returns the trace K on the 2-surface. | |
| const Vector & | get_ll () const |
Returns the vector ![]() | |
| const Vector & | get_ss () const |
Returns the vector ![]() | |
| const Vector & | get_ephi () const |
| Returns the conformal Killing symmetry vector on the 2-surface. | |
| const Sym_tensor & | get_jj () const |
Returns the symmetric tensor ![]() | |
| const Scalar & | get_fff () const |
Returns the normalization scalar F. | |
| const Scalar & | get_ggg () const |
Returns the normalization scalar G. | |
| bool | get_issphere () const |
| Returns the flag saying whether or not the horizon is geometrically round. | |
| Scalar & | set_hsurf () |
Sets the field h_surf. | |
| Metric & | set_qab () |
Sets the modified metric (non degenerated) ![]() | |
| Scalar & | set_ricci () |
Sets the 2-Ricci scalar ![]() | |
| Sym_tensor & | set_qq () |
Sets the degenerated metric ![]() | |
| Tensor & | set_proj () |
Sets the projector ![]() | |
| Sym_tensor & | set_hh () |
Sets the symmetric tensor ![]() | |
| Scalar & | set_trk () |
Sets the trace K on the 2-surface. | |
| Vector & | set_ll () |
Sets the vector ![]() | |
| Vector & | set_ss () |
Sets the vector ![]() | |
| Sym_tensor & | set_jj () |
Sets the symmetric tensor ![]() | |
| Scalar & | set_fff () |
Sets the normalization factor F. | |
| Scalar & | set_ggg () |
Sets the normalization factor G. | |
| bool | set_issphere () |
| Sets the boolean linked to geometrical shape of the horizon. | |
| void | update_from_tslice (const Metric &gamij, const Sym_tensor &Kij) |
| Updates from the 3-slice data. | |
| const Scalar & | sqrt_q () const |
| Computes the normal vector field to the 2-surface. | |
| double | area () const |
| Computes the area of the 2-surface. | |
| double | angu_mom () const |
| Computes the angular momentum with respect to a divergence-free vector field tangent to the 2-surface. | |
| double | mass () const |
Computes the mass as defined from the calculus of angular momentum, done with respect to a divergence free tangent vector field ![]() | |
| double | multipole_mass (const int order) const |
| Computes the mass multipole of a given order for the spheroid, assumed to be spherical. | |
| double | multipole_angu (const int order) const |
| Computes the angular multipole of a given order for the spheroid, assumed to be spherical. | |
| double | epsilon_A_minus_one () const |
| Computation of the refined Penrose parameter for axisymmetric spacetimes, and its difference wrt one. | |
| double | epsilon_P_minus_one () const |
| Computation of the classical Penrose parameter, and its difference wrt one. | |
| const Scalar & | theta_plus () const |
Computes the outgoing null expansion ![]() | |
| const Scalar & | theta_minus () const |
Computes the ingoing null expansion ![]() | |
| const Sym_tensor & | shear () const |
Computes the shear of the 2-surface ![]() | |
| Tensor | derive_cov2dflat (const Tensor &uu) const |
| Computes the round covariant derivative on the spheroid. | |
| const Tensor & | delta () const |
| Computes the delta coefficients for covariant derivative. | |
| Tensor | derive_cov2d (const Tensor &uu) const |
| Computes the total covariant derivative on the spheroid. | |
| virtual void | sauve (FILE *) const |
| Save in a file. | |
Protected Member Functions | |
| virtual void | del_deriv () const |
| Deletes all the derived quantities. | |
| void | set_der_0x0 () const |
Sets to 0x0 all the pointers on derived quantities. | |
Protected Attributes | |
| Scalar | h_surf |
The location of the 2-surface as r = h_surf ![]() | |
| Tensor | jac2d |
| The jacobian of the adaptation of coordinates (contravariant/covariant representation) | |
| Tensor | proj |
| The 3-d projector on the 2-surface (contravariant-covariant form). | |
| Sym_tensor | |
| The 3-d covariant degenerated 2-metric on the surface. | |
| Vector | ss |
| The adapted normal vector field to spheroid in the 3-slice. | |
| Vector | ephi |
The conformal Killing vector field on the 2-surface (set to by default to the axial vector associated with coordinate ![]() | |
| Metric | qab |
| Scalar | ricci |
Induced metric on the 2-surface ![]() | |
| Sym_tensor | hh |
| The ricci scalar on the 2-surface. | |
| Scalar | trk |
| Trace K of the extrinsic curvature of the 3-slice. | |
| Vector | ll |
| Normal-tangent component of the extrinsic curvature of the 3-slice. | |
| Sym_tensor | jj |
| Tangent components of the extrinsic curvature of the 3-slice. | |
| Scalar | fff |
| Normalization function for the outgoing null vector l. | |
| Scalar | ggg |
| Normalization function for the ingoing null vector k. | |
| Scalar | zeta |
| bool | issphere |
| Flag to know whether the horizon is geometrically round or distorted. | |
| Scalar * | p_sqrt_q |
Surface element ![]() | |
| double * | p_area |
| The area of the 2-surface. | |
| double * | p_angu_mom |
| The angular momentum. | |
| double * | p_mass |
| Mass defined from angular momentum. | |
| double * | p_multipole_mass |
| Mass multipole for the spheroid. | |
| double * | p_multipole_angu |
| Angular momentum multipole for the spheroid. | |
| double * | p_epsilon_A_minus_one |
| double * | p_epsilon_P_minus_one |
| Refined Penrose parameter, difference wrt one. | |
| Scalar * | p_theta_plus |
| Classical Penrose parameter, difference wrt one. | |
| Scalar * | p_theta_minus |
| Null ingoing expansion. | |
| Sym_tensor * | p_shear |
| The shear tensor. | |
| Tensor * | p_delta |
Class describing an apparent horizon.
It is derived from the class Spheroid, with the property that the outgoing null expansion 
Definition at line 408 of file spheroid.h.
Constructor of an apparent horizon embedded in a 3-slice (Time_slice ) of 3+1 formalism.
This is done from the Time_slice data.
| h_in | : the location of the surface r = h_in |
| gamij | : the 3-metric on the 3-slice |
| Kij | : the extrinsic curvature of the 3-slice (covariant representation) |
| Lorene::App_hor::App_hor | ( | FILE * | ) |
Constructor from a file (see sauve(FILE*) )
|
virtual |
Destructor.
|
inherited |
Computes the angular momentum with respect to a divergence-free vector field tangent to the 2-surface.
This is defined as
![\[
{\cal J} = \int h^2 L_i \phi^i \sqrt{\det q_{ab}} \sin \theta {\rm d}\theta
{\rm d}\varphi \]](form_747.png)
| phi | : the divergence-free vector field ![]() |
Definition at line 721 of file spheroid.C.
References Lorene::contract(), Lorene::Map::get_bvect_spher(), Lorene::Spheroid::get_ephi(), Lorene::Tensor::get_mp(), Lorene::Spheroid::h_surf, Lorene::Map_af::integrale_surface(), Lorene::Spheroid::jac2d, Lorene::Spheroid::ll, Lorene::Spheroid::p_angu_mom, and Lorene::Spheroid::sqrt_q().
|
inherited |
Computes the area of the 2-surface.
This is defined as
![\[
{\cal A} = \int h^2 \sqrt{\det q_{ab}} \sin \theta {\rm d}\theta
{\rm d}\varphi \]](form_746.png)
Definition at line 708 of file spheroid.C.
References Lorene::Tensor::get_mp(), Lorene::Spheroid::h_surf, Lorene::Map_af::integrale_surface(), Lorene::Spheroid::p_area, and Lorene::Spheroid::sqrt_q().
|
protectedvirtualinherited |
Deletes all the derived quantities.
Definition at line 650 of file spheroid.C.
References Lorene::Spheroid::p_angu_mom, Lorene::Spheroid::p_area, Lorene::Spheroid::p_epsilon_P_minus_one, Lorene::Spheroid::p_mass, Lorene::Spheroid::p_multipole_angu, Lorene::Spheroid::p_multipole_mass, Lorene::Spheroid::p_shear, Lorene::Spheroid::p_sqrt_q, Lorene::Spheroid::p_theta_minus, Lorene::Spheroid::p_theta_plus, and Lorene::Spheroid::set_der_0x0().
Computes the delta coefficients for covariant derivative.
Definition at line 1203 of file spheroid.C.
References Lorene::Metric::con(), Lorene::Metric::cov(), Lorene::Spheroid::derive_cov2dflat(), Lorene::Map::get_bvect_spher(), Lorene::Metric::get_mp(), Lorene::Tensor::set(), Lorene::Tensor::set_etat_zero(), and Lorene::Tensor::set_index_type().
Computes the total covariant derivative on the spheroid.
Definition at line 1240 of file spheroid.C.
References Lorene::contract(), Lorene::Spheroid::delta(), Lorene::Spheroid::derive_cov2dflat(), Lorene::Tensor::get_index_type(), and Lorene::Tensor::get_valence().
Computes the round covariant derivative on the spheroid.
Definition at line 954 of file spheroid.C.
References Lorene::Scalar::div_r_dzpuis(), Lorene::Scalar::div_tant(), Lorene::Map::get_bvect_spher(), Lorene::Tensor::get_index_type(), Lorene::Tensor::get_mp(), Lorene::Tensor::get_n_comp(), Lorene::Tensor::get_triad(), Lorene::Tensor::get_valence(), Lorene::Tensor::indices(), Lorene::Tensor::set(), Lorene::Itbl::set(), Lorene::Tensor_sym::sym_index1(), and Lorene::Tensor_sym::sym_index2().
|
inherited |
Computation of the refined Penrose parameter for axisymmetric spacetimes, and its difference wrt one.
Definition at line 864 of file spheroid.C.
References Lorene::Spheroid::angu_mom(), Lorene::Spheroid::area(), Lorene::Spheroid::mass(), Lorene::pow(), and Lorene::sqrt().
|
inherited |
Computation of the classical Penrose parameter, and its difference wrt one.
To use in replacement of epsilon_A_minus_one when the computed spacetime is not axisymmetric.
Definition at line 876 of file spheroid.C.
References Lorene::Spheroid::angu_mom(), Lorene::Spheroid::area(), Lorene::Spheroid::mass(), Lorene::Spheroid::p_epsilon_P_minus_one, and Lorene::pow().
Returns the conformal Killing symmetry vector on the 2-surface.
Definition at line 253 of file spheroid.h.
References Lorene::Spheroid::ephi.
Returns the normalization scalar F.
Definition at line 259 of file spheroid.h.
References Lorene::Spheroid::fff.
Returns the normalization scalar G.
Definition at line 262 of file spheroid.h.
References Lorene::Spheroid::ggg.
|
inlineinherited |
Returns the symmetric tensor 
Definition at line 232 of file spheroid.h.
References Lorene::Spheroid::hh.
Returns the field h_surf.
Definition at line 223 of file spheroid.h.
References Lorene::Spheroid::h_surf.
|
inlineinherited |
Returns the flag saying whether or not the horizon is geometrically round.
Definition at line 265 of file spheroid.h.
References Lorene::Spheroid::issphere.
returns the 2-d jacobian of coordinate transformation 
Definition at line 241 of file spheroid.h.
References Lorene::Spheroid::jac2d.
|
inlineinherited |
Returns the symmetric tensor 
Definition at line 256 of file spheroid.h.
References Lorene::Spheroid::jj.
returns the 3-d projector on 2-surface 
Definition at line 238 of file spheroid.h.
References Lorene::Spheroid::proj.
Returns the metric 
Definition at line 226 of file spheroid.h.
|
inlineinherited |
returns the 3-d degenerate 2-metric 
Definition at line 235 of file spheroid.h.
References Lorene::Spheroid::qq.
Returns the 2-ricci scalar 
Definition at line 229 of file spheroid.h.
References Lorene::Spheroid::ricci.
Returns the trace K on the 2-surface.
Definition at line 244 of file spheroid.h.
References Lorene::Spheroid::trk.
non-affinity (or surface gravity) with respect to the outgoing null vector field
Lie derivative of 2-metric with respect to the evolution vector field.
Lie derivative of shear tensor with respect to the evolution vector field.
| const Sym_tensor & Lorene::App_hor::lie_derive_theta_minus | ( | const Scalar & | bb, |
| const Scalar & | lapse | ||
| ) |
Lie derivative of the null ingoing expansion rate with respect to the evolution vector field.
| const Sym_tensor & Lorene::App_hor::lie_derive_theta_plus | ( | const Scalar & | bb, |
| const Scalar & | lapse | ||
| ) |
Lie derivative of the null outgoing expansion rate with respect to the evolution vector field.
|
inherited |
Computes the mass as defined from the calculus of angular momentum, done with respect to a divergence free tangent vector field 
Spheroid has to be a real sphere (flag issphere true), of constant radius
![\[R_{s} \]](form_750.png)
. defined as
![\[ M = \frac{1}{2 R_{s}} \sqrt{R_{s}^{4} + 4{\cal J}^{2}} \]](form_751.png)
Definition at line 736 of file spheroid.C.
References Lorene::Spheroid::angu_mom(), Lorene::Spheroid::area(), Lorene::Spheroid::p_mass, and Lorene::sqrt().
Computes the angular multipole of a given order for the spheroid, assumed to be spherical.

Definition at line 801 of file spheroid.C.
References Lorene::Spheroid::area(), Lorene::contract(), Lorene::Map::get_bvect_spher(), Lorene::Spheroid::get_ephi(), Lorene::Tensor::get_mp(), Lorene::Spheroid::h_surf, Lorene::Map_af::integrale_surface(), Lorene::Spheroid::jac2d, Lorene::Spheroid::ll, Lorene::Spheroid::p_multipole_angu, Lorene::Scalar::set_spectral_va(), Lorene::sqrt(), Lorene::Spheroid::sqrt_q(), Lorene::Scalar::std_spectral_base(), and Lorene::Valeur::ylm().
Computes the mass multipole of a given order for the spheroid, assumed to be spherical.
WARNING: For technical reasons, only even orders are supported by the code.
Definition at line 748 of file spheroid.C.
References Lorene::Spheroid::area(), Lorene::Tensor::get_mp(), Lorene::Spheroid::get_ricci(), Lorene::Spheroid::h_surf, Lorene::Map_af::integrale_surface(), Lorene::Spheroid::mass(), Lorene::Spheroid::p_multipole_mass, Lorene::Scalar::set_spectral_va(), Lorene::sqrt(), Lorene::Spheroid::sqrt_q(), Lorene::Scalar::std_spectral_base(), and Lorene::Valeur::ylm().
Save in a file.
Definition at line 1186 of file spheroid.C.
|
protectedinherited |
Sets to 0x0 all the pointers on derived quantities.
Definition at line 666 of file spheroid.C.
References Lorene::Spheroid::p_angu_mom, Lorene::Spheroid::p_area, Lorene::Spheroid::p_epsilon_P_minus_one, Lorene::Spheroid::p_mass, Lorene::Spheroid::p_multipole_angu, Lorene::Spheroid::p_multipole_mass, Lorene::Spheroid::p_shear, Lorene::Spheroid::p_sqrt_q, Lorene::Spheroid::p_theta_minus, and Lorene::Spheroid::p_theta_plus.
Assigns the conformal Killing vector field to phi.
|
inlineinherited |
Sets the normalization factor F.
Definition at line 298 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::fff.
|
inlineinherited |
Sets the normalization factor G.
Definition at line 301 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::ggg.
|
inlineinherited |
Sets the symmetric tensor 
Definition at line 283 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::hh.
|
inlineinherited |
Sets the field h_surf.
Definition at line 268 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::h_surf.
|
inlineinherited |
Sets the boolean linked to geometrical shape of the horizon.
Definition at line 304 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::issphere.
|
inlineinherited |
Sets the symmetric tensor 
Definition at line 295 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::jj.
|
inlineinherited |
Sets the vector 
Definition at line 289 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::ll.
|
inlineinherited |
Sets the projector 
Definition at line 280 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::proj.
|
inlineinherited |
Sets the modified metric (non degenerated) 
Definition at line 271 of file spheroid.h.
References Lorene::Spheroid::del_deriv().
|
inlineinherited |
Sets the degenerated metric 
Definition at line 277 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::qq.
|
inlineinherited |
Sets the 2-Ricci scalar 
Definition at line 274 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::ricci.
|
inlineinherited |
Sets the vector 
Definition at line 292 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::ss.
|
inlineinherited |
Sets the trace K on the 2-surface.
Definition at line 286 of file spheroid.h.
References Lorene::Spheroid::del_deriv(), and Lorene::Spheroid::trk.
|
inherited |
Computes the shear of the 2-surface 
Definition at line 925 of file spheroid.C.
References Lorene::Metric::cov(), Lorene::Spheroid::fff, Lorene::Spheroid::hh, Lorene::Spheroid::jj, Lorene::Spheroid::p_shear, Lorene::Tensor::std_spectral_base(), and Lorene::Tensor::trace().
Computes the normal vector field to the 2-surface.
Computes the square root of the determinant of 
Definition at line 693 of file spheroid.C.
References Lorene::Spheroid::get_qq(), Lorene::Spheroid::p_sqrt_q, Lorene::sqrt(), and Lorene::Scalar::std_spectral_base().
Computes the ingoing null expansion 
Definition at line 909 of file spheroid.C.
References Lorene::Spheroid::ggg, Lorene::Spheroid::hh, Lorene::Spheroid::jj, Lorene::Spheroid::p_theta_minus, Lorene::Scalar::std_spectral_base(), and Lorene::Tensor::trace().
Computes the outgoing null expansion 
Definition at line 889 of file spheroid.C.
References Lorene::Spheroid::fff, Lorene::Spheroid::hh, Lorene::Spheroid::jj, Lorene::Spheroid::p_theta_plus, Lorene::Scalar::set_spectral_va(), Lorene::Scalar::std_spectral_base(), Lorene::Tensor::trace(), and Lorene::Valeur::ylm().
|
inherited |
Updates from the 3-slice data.
|
protectedinherited |
The conformal Killing vector field on the 2-surface (set to by default to the axial vector associated with coordinate 
Definition at line 113 of file spheroid.h.
|
protectedinherited |
Normalization function for the outgoing null vector l.
Definition at line 138 of file spheroid.h.
|
protectedinherited |
Normalization function for the ingoing null vector k.
Definition at line 143 of file spheroid.h.
|
protectedinherited |
The location of the 2-surface as r = h_surf 
Definition at line 91 of file spheroid.h.
|
protectedinherited |
The ricci scalar on the 2-surface.
Extrinsic curvature of the 2-surface in the 3-slice. 
Definition at line 122 of file spheroid.h.
|
protectedinherited |
Flag to know whether the horizon is geometrically round or distorted.
Definition at line 151 of file spheroid.h.
|
protectedinherited |
The jacobian of the adaptation of coordinates (contravariant/covariant representation)
Definition at line 96 of file spheroid.h.
|
protectedinherited |
Tangent components of the extrinsic curvature of the 3-slice.

Definition at line 134 of file spheroid.h.
|
protectedinherited |
Normal-tangent component of the extrinsic curvature of the 3-slice.

Definition at line 129 of file spheroid.h.
|
mutableprotectedinherited |
The angular momentum.
Definition at line 159 of file spheroid.h.
|
mutableprotectedinherited |
The area of the 2-surface.
Definition at line 158 of file spheroid.h.
|
mutableprotectedinherited |
Definition at line 168 of file spheroid.h.
|
mutableprotectedinherited |
Definition at line 163 of file spheroid.h.
|
mutableprotectedinherited |
Refined Penrose parameter, difference wrt one.
Definition at line 164 of file spheroid.h.
|
mutableprotectedinherited |
Mass defined from angular momentum.
Definition at line 160 of file spheroid.h.
|
mutableprotectedinherited |
Angular momentum multipole for the spheroid.
Definition at line 162 of file spheroid.h.
|
mutableprotectedinherited |
Mass multipole for the spheroid.
Definition at line 161 of file spheroid.h.
|
mutableprotectedinherited |
The shear tensor.
Definition at line 167 of file spheroid.h.
|
mutableprotectedinherited |
Surface element 
Definition at line 157 of file spheroid.h.
|
mutableprotectedinherited |
Null ingoing expansion.
Definition at line 166 of file spheroid.h.
|
mutableprotectedinherited |
Classical Penrose parameter, difference wrt one.
Null outgoing expansion
Definition at line 165 of file spheroid.h.
|
protectedinherited |
The 3-d projector on the 2-surface (contravariant-covariant form).
Definition at line 100 of file spheroid.h.
|
protectedinherited |
Definition at line 115 of file spheroid.h.
|
protectedinherited |
The 3-d covariant degenerated 2-metric on the surface.
Definition at line 104 of file spheroid.h.
|
protectedinherited |
Induced metric on the 2-surface 
Definition at line 117 of file spheroid.h.
|
protectedinherited |
The adapted normal vector field to spheroid in the 3-slice.
Definition at line 108 of file spheroid.h.
|
protectedinherited |
Trace K of the extrinsic curvature of the 3-slice.
Definition at line 124 of file spheroid.h.
|
protectedinherited |
Definition at line 147 of file spheroid.h.