Classes
class	Lorene::Bhole
	Black hole. More...
class	Lorene::Bhole_binaire
	Binary black holes system. More...
class	Lorene::Bin_bhns
	Class for computing a black hole - neutron star binary system with comparable mass (). More...
class	Lorene::Bin_bhns_extr
	Class for computing a Black hole - Neutron star binary system with an extreme mass ratio. More...
class	Lorene::Binary
	Binary systems. More...
class	Lorene::Binary_xcts
	Binary systems in eXtended Conformal Thin Sandwich formulation. More...
class	Lorene::Black_hole
	Base class for black holes. More...
class	Lorene::Et_bin_bhns_extr
	Class for a neutron star in black hole - neutron star binary systems. More...
class	Lorene::Et_rot_bifluid
	Class for two-fluid rotating relativistic stars. More...
class	Lorene::Et_rot_diff
	Class for differentially rotating stars. More...
class	Lorene::Et_rot_mag
	Class for magnetized (isolator or perfect conductor), rigidly rotating stars. More...
class	Lorene::Etoile
	Base class for stars * DEPRECATED : use class `Star` instead *. More...
class	Lorene::Etoile_bin
	Class for stars in binary system. More...
class	Lorene::Etoile_rot
	Class for isolated rotating stars * DEPRECATED : use class `Star_rot` instead *. More...
class	Lorene::Excision_hor
	Surface where boundary conditions for quantities in the bulk will be calculated It relies on geometrical properties of the associated Spheroid() (*** WARNING! More...
class	Lorene::Excision_surf
	Surface where boundary conditions for quantities in the bulk will be calculated It relies on geometrical properties of the associated Spheroid() (*** WARNING! More...
class	Lorene::Gravastar
	Class for perfect fluid rotating gravastar. More...
class	Lorene::Hole_bhns
	Class for black holes in black hole-neutron star binaries. More...
class	Lorene::Single_hor
	Binary black holes system. More...
class	Lorene::Spheroid
	Spheroidal 2-surfaces embedded in a time-slice of the 3+1 formalism. More...
class	Lorene::Star
	Base class for stars. More...
class	Lorene::Star_bin
	Class for stars in binary system. More...
class	Lorene::Star_bin_xcts
	Class for stars in binary system in eXtended Conformal Thin Sandwich formulation. More...
class	Lorene::Star_bhns
	Class for stars in black hole-neutron star binaries. More...
class	Lorene::Star_rot
	Class for isolated rotating stars. More...
class	Lorene::Star_rot_Dirac
	Class for relativistic rotating stars in Dirac gauge and maximal slicing. More...
class	Lorene::Star_rot_Dirac_diff
	Class for relativistic differentially rotating stars in Dirac gauge and maximal slicing. More...

Functions
bool	Lorene::ah_finder (const Metric &gamma, const Sym_tensor &k_dd_in, Valeur &h, Scalar &ex_fcn, double a_axis, double b_axis, double c_axis, bool verbose=true, bool print=false, double precis=1.e-8, double precis_exp=1.e-6, int it_max=200, int it_relax=200, double relax_fac=1.)
	Class for apparent horizon (under heavy development).

Detailed Description

Function Documentation

◆ ah_finder()

bool Lorene::ah_finder	(	const Metric &	gamma,
		const Sym_tensor &	k_dd_in,
		Valeur &	h,
		Scalar &	ex_fcn,
		double	a_axis,
		double	b_axis,
		double	c_axis,
		bool	verbose = true,
		bool	print = false,
		double	precis = 1.e-8,
		double	precis_exp = 1.e-6,
		int	it_max = 200,
		int	it_relax = 200,
		double	relax_fac = 1. )

Class for apparent horizon (under heavy development).

Apparent horizon finder.()

Find the apparent horizon (AH) on a spatial slice for given 3-metric $\gamma_{ij}$ and extrinsic curvature $K_{ij}$ . Method: We solve the apparent-horizon equation $\Theta=0$ (where $\Theta$ is the expansion function of the outgoing null geodesic) as a generalized angular Poisson equation ( $\Delta_{\theta\phi}u+\lambda u=\sigma$ ) in the form:

$\‍[ \Delta_{\theta\phi}h - 2h = \Psi^4 |DF| h^2\Theta + \Delta_{\theta\phi}h - 2h \‍]$

where $h=h(\theta,\phi)$ (h) is the 2-surface of the AH, $\Psi$ is the conformal factor, and $|DF|:= (\gamma^{ij}D_iFD_jF)^{1/2}$ (where $F$ is the level set function $F:=r-h(\theta,\phi)$ and $D_i$ is the covariant derivative w.r.t. $\gamma_{ij}$ ).
We solve the above equation iteratively with the right hand side of the equation as a source term.

Parameters

gamma	: [input] the 3-metric $\gamma_{ij}$ w.r.t. which the AH is to be found.
k_dd_in	: [input] the extrinsic curvature $K_{ij}$ .
Valeur	h : [output] the 2-surface of the apparent horizon
Scalar	ex_fcn : [output] the expansion function defined from the level set function $F:= r - h(\theta,\phi)$
a_axis	: [input] the initial guess for is a triaxial ellipsoidal surface with a_axis the axis along the x-axis of the Cartesian grid
b_axis	: [input] axis along the y-axis (cf a_axis)
c_axis	: [input] axis along the z-axis (cf a_axis)
bool	print : [input] screen printout during iterations? (default: false)
precis	: [input] threshold in the relative difference between the 2-surface of two consecutive steps to stop the iterative procedure (default value: 1.e-8)
precis_exp	: [input] maximum error of the expansion function evaluated on the 2-surface (should be zero by definition) after the iteration is stopped (default value: 1.e-6)
it_max	: [input] maximum number of steps (default value: 200)
it_relax	: [input] step at which relaxation is switched on (default value: 200)
relax_fac	: [input] relaxation factor (default value: 1 for no relaxation)

Returns: bool ah_flag : a flag to indicate whether an apparent horizon is found

Definition at line 91 of file app_hor_finder.C.

Classes

Functions

Detailed Description

Function Documentation

◆ ah_finder()