LORENE
|
Class describing an apparent horizon. More...
#include <spheroid.h>
Public Member Functions | |
App_hor (const Mg3d &grid_angu, double radius) | |
Standard constructor. More... | |
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. More... | |
App_hor (const App_hor &) | |
Copy constructor. More... | |
App_hor (FILE *) | |
Constructor from a file (see sauve(FILE*) ) More... | |
virtual | ~App_hor () |
Destructor. More... | |
void | operator= (const App_hor &) |
Assignment to another App_hor. More... | |
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. More... | |
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. More... | |
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. More... | |
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. More... | |
const Scalar & | l_non_affinity (const Scalar &bb, const Scalar &lapse) |
non-affinity (or surface gravity) with respect to the outgoing null vector field More... | |
void | set_ephi (const Scalar &) |
Assigns the conformal Killing vector field to phi. More... | |
const Scalar & | get_hsurf () const |
Returns the field h_surf . More... | |
const Metric & | get_qab () const |
Returns the metric . More... | |
const Scalar & | get_ricci () const |
Returns the 2-ricci scalar . More... | |
const Sym_tensor & | get_hh () const |
Returns the symmetric tensor . More... | |
const Sym_tensor & | get_qq () const |
returns the 3-d degenerate 2-metric More... | |
const Tensor & | get_proj () const |
returns the 3-d projector on 2-surface More... | |
const Tensor & | get_jac2d () const |
returns the 2-d jacobian of coordinate transformation More... | |
const Scalar & | get_trk () const |
Returns the trace K on the 2-surface. More... | |
const Vector & | get_ll () const |
Returns the vector . More... | |
const Vector & | get_ss () const |
Returns the vector . More... | |
const Vector & | get_ephi () const |
Returns the conformal Killing symmetry vector on the 2-surface. More... | |
const Sym_tensor & | get_jj () const |
Returns the symmetric tensor . More... | |
const Scalar & | get_fff () const |
Returns the normalization scalar F . More... | |
const Scalar & | get_ggg () const |
Returns the normalization scalar G . More... | |
bool | get_issphere () const |
Returns the flag saying whether or not the horizon is geometrically round. More... | |
Scalar & | set_hsurf () |
Sets the field h_surf . More... | |
Metric & | set_qab () |
Sets the modified metric (non degenerated) . More... | |
Scalar & | set_ricci () |
Sets the 2-Ricci scalar . More... | |
Sym_tensor & | set_qq () |
Sets the degenerated metric . More... | |
Tensor & | set_proj () |
Sets the projector . More... | |
Sym_tensor & | set_hh () |
Sets the symmetric tensor . More... | |
Scalar & | set_trk () |
Sets the trace K on the 2-surface. More... | |
Vector & | set_ll () |
Sets the vector . More... | |
Vector & | set_ss () |
Sets the vector . More... | |
Sym_tensor & | set_jj () |
Sets the symmetric tensor . More... | |
Scalar & | set_fff () |
Sets the normalization factor F . More... | |
Scalar & | set_ggg () |
Sets the normalization factor G . More... | |
bool | set_issphere () |
Sets the boolean linked to geometrical shape of the horizon. More... | |
void | update_from_tslice (const Metric &gamij, const Sym_tensor &Kij) |
Updates from the 3-slice data. More... | |
const Scalar & | sqrt_q () const |
Computes the normal vector field to the 2-surface. More... | |
double | area () const |
Computes the area of the 2-surface. More... | |
double | angu_mom () const |
Computes the angular momentum with respect to a divergence-free vector field tangent to the 2-surface. More... | |
double | mass () const |
Computes the mass as defined from the calculus of angular momentum, done with respect to a divergence free tangent vector field . More... | |
double | multipole_mass (const int order) const |
Computes the mass multipole of a given order for the spheroid, assumed to be spherical. More... | |
double | multipole_angu (const int order) const |
Computes the angular multipole of a given order for the spheroid, assumed to be spherical. More... | |
double | epsilon_A_minus_one () const |
Computation of the refined Penrose parameter for axisymmetric spacetimes, and its difference wrt one. More... | |
double | epsilon_P_minus_one () const |
Computation of the classical Penrose parameter, and its difference wrt one. More... | |
const Scalar & | theta_plus () const |
Computes the outgoing null expansion . More... | |
const Scalar & | theta_minus () const |
Computes the ingoing null expansion . More... | |
const Sym_tensor & | shear () const |
Computes the shear of the 2-surface . More... | |
Tensor | derive_cov2dflat (const Tensor &uu) const |
Computes the round covariant derivative on the spheroid. More... | |
const Tensor & | delta () const |
Computes the delta coefficients for covariant derivative. More... | |
Tensor | derive_cov2d (const Tensor &uu) const |
Computes the total covariant derivative on the spheroid. More... | |
virtual void | sauve (FILE *) const |
Save in a file. More... | |
Protected Member Functions | |
virtual void | del_deriv () const |
Deletes all the derived quantities. More... | |
void | set_der_0x0 () const |
Sets to 0x0 all the pointers on derived quantities. More... | |
Protected Attributes | |
Scalar | h_surf |
The location of the 2-surface as r = h_surf . More... | |
Tensor | jac2d |
The jacobian of the adaptation of coordinates (contravariant/covariant representation) More... | |
Tensor | proj |
The 3-d projector on the 2-surface (contravariant-covariant form). More... | |
Sym_tensor | |
The 3-d covariant degenerated 2-metric on the surface. More... | |
Vector | ss |
The adapted normal vector field to spheroid in the 3-slice. More... | |
Vector | ephi |
The conformal Killing vector field on the 2-surface (set to by default to the axial vector associated with coordinate ) More... | |
Metric | qab |
Scalar | ricci |
Induced metric on the 2-surface . More... | |
Sym_tensor | hh |
The ricci scalar on the 2-surface. More... | |
Scalar | trk |
Trace K of the extrinsic curvature of the 3-slice. More... | |
Vector | ll |
Normal-tangent component of the extrinsic curvature of the 3-slice. More... | |
Sym_tensor | jj |
Tangent components of the extrinsic curvature of the 3-slice. More... | |
Scalar | fff |
Normalization function for the outgoing null vector l. More... | |
Scalar | ggg |
Normalization function for the ingoing null vector k. More... | |
Scalar | zeta |
bool | issphere |
Flag to know whether the horizon is geometrically round or distorted. More... | |
Scalar * | p_sqrt_q |
Surface element . More... | |
double * | p_area |
The area of the 2-surface. More... | |
double * | p_angu_mom |
The angular momentum. More... | |
double * | p_mass |
Mass defined from angular momentum. More... | |
double * | p_multipole_mass |
Mass multipole for the spheroid. More... | |
double * | p_multipole_angu |
Angular momentum multipole for the spheroid. More... | |
double * | p_epsilon_A_minus_one |
double * | p_epsilon_P_minus_one |
Refined Penrose parameter, difference wrt one. More... | |
Scalar * | p_theta_plus |
Classical Penrose parameter, difference wrt one. More... | |
Scalar * | p_theta_minus |
Null ingoing expansion. More... | |
Sym_tensor * | p_shear |
The shear tensor. More... | |
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.
Lorene::App_hor::App_hor | ( | const Mg3d & | grid_angu, |
double | radius | ||
) |
Standard constructor.
Lorene::App_hor::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.
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 | ( | const App_hor & | ) |
Copy constructor.
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
phi | : the divergence-free vector field |
Definition at line 724 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
Definition at line 711 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 653 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().
|
inherited |
Computes the delta coefficients for covariant derivative.
Definition at line 1206 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(), and Lorene::Tensor::set_index_type().
Computes the total covariant derivative on the spheroid.
Definition at line 1243 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 957 of file spheroid.C.
References 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::Itbl::set(), Lorene::Tensor::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 867 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 879 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().
|
inlineinherited |
Returns the conformal Killing symmetry vector on the 2-surface.
Definition at line 253 of file spheroid.h.
References Lorene::Spheroid::ephi.
|
inlineinherited |
Returns the normalization scalar F
.
Definition at line 259 of file spheroid.h.
References Lorene::Spheroid::fff.
|
inlineinherited |
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.
|
inlineinherited |
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.
|
inlineinherited |
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.
|
inlineinherited |
|
inlineinherited |
returns the 3-d projector on 2-surface
Definition at line 238 of file spheroid.h.
References Lorene::Spheroid::proj.
|
inlineinherited |
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.
|
inlineinherited |
Returns the 2-ricci scalar .
Definition at line 229 of file spheroid.h.
References Lorene::Spheroid::ricci.
|
inlineinherited |
|
inlineinherited |
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
const Sym_tensor& Lorene::App_hor::lie_derive_q_ab | ( | const Scalar & | bb, |
const Scalar & | lapse | ||
) |
Lie derivative of 2-metric with respect to the evolution vector field.
const Sym_tensor& Lorene::App_hor::lie_derive_shear | ( | const Scalar & | bb, |
const Scalar & | lapse | ||
) |
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
. defined as
Definition at line 739 of file spheroid.C.
References Lorene::Spheroid::angu_mom(), Lorene::Spheroid::area(), Lorene::Spheroid::p_mass, and Lorene::sqrt().
|
inherited |
Computes the angular multipole of a given order for the spheroid, assumed to be spherical.
is a divergence free tangent vector field. WARNING:order has to be strictly higher than zero (no topological defects here...), and an odd number for technical reasons.
Definition at line 804 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().
|
inherited |
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 751 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().
|
virtualinherited |
Save in a file.
Definition at line 1189 of file spheroid.C.
|
protectedinherited |
Sets to 0x0
all the pointers on derived quantities.
Definition at line 669 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.
|
inherited |
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 928 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().
|
inherited |
Computes the normal vector field to the 2-surface.
Computes the square root of the determinant of .
Definition at line 696 of file spheroid.C.
References Lorene::Spheroid::get_qq(), Lorene::Spheroid::p_sqrt_q, Lorene::sqrt(), and Lorene::Scalar::std_spectral_base().
|
inherited |
Computes the ingoing null expansion .
Definition at line 912 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().
|
inherited |
Computes the outgoing null expansion .
Definition at line 892 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. (covariant representation)
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.
(covariant representation)
Definition at line 134 of file spheroid.h.
|
protectedinherited |
Normal-tangent component of the extrinsic curvature of the 3-slice.
(covariant representation)
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 |
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 |
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.