isol_hor.h
/*
* Definition of Lorene class Isol_Hor
*
*/
/*
* Copyright (c) 2004-2005 Jose Luis Jaramillo
* Francois Limousin
*
* This file is part of LORENE.
*
* LORENE is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License version 2
* as published by the Free Software Foundation.
*
* LORENE is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with LORENE; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*/
#ifndef __Isol_hor_H_
#define __Isol_hor_H_
/*
* $Id: isol_hor.h,v 1.54 2014/10/13 08:52:35 j_novak Exp $
* $Log: isol_hor.h,v $
* Revision 1.54 2014/10/13 08:52:35 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.53 2007/08/22 16:10:37 f_limousin
* Correction of many errors in binhor_hh.C
*
* Revision 1.52 2007/04/13 15:29:44 f_limousin
* Lots of improvements, generalisation to an arbitrary state of
* rotation, implementation of the spatial metric given by Samaya.
*
* Revision 1.51 2006/08/01 14:36:25 f_limousin
* New argument for the functions vv_bound_cart_bin( )
*
* Revision 1.50 2006/06/29 08:58:57 f_limousin
* Boundary conditions in argument of write_global()
*
* Revision 1.49 2006/06/28 13:36:26 f_limousin
* Convergence to a given irreductible mass
*
* Revision 1.48 2006/05/24 16:53:51 f_limousin
* Funtion import(Bin_hor&)
*
* Revision 1.47 2006/02/22 16:55:42 f_limousin
* boundary_nn_Neu_kk(int nn = 1)
*
* Revision 1.46 2006/02/22 16:20:54 jl_jaramillo
* 2nc correction Valeur boundary_nn_Neu_kk(int nn)
*
* Revision 1.45 2006/02/20 16:48:14 jl_jaramillo
* corrections numerical viscosity
*
* Revision 1.44 2006/02/20 15:01:37 jl_jaramillo
* function for computing generalised Conformal THin Sandwich eqautions
*
* Revision 1.43 2006/01/18 12:29:18 jl_jaramillo
* new method (chi,theta) for the spherically symmetric case
*
* Revision 1.42 2006/01/16 17:15:34 jl_jaramillo
* function for solving the spherical case
*
* Revision 1.41 2005/11/02 16:10:38 jl_jaramillo
* change in boundary_nn_Dir_lapl
*
* Revision 1.40 2005/10/24 16:45:16 jl_jaramillo
* Cook boundary condition
*
* Revision 1.39 2005/10/23 12:32:00 f_limousin
* Pure Dirichlet boundary condition for Psi.
*
* Revision 1.38 2005/09/13 18:33:17 f_limousin
* New function vv_bound_cart_bin(double) for computing binaries with
* berlin condition for the shift vector.
* Suppress all the symy and asymy in the importations.
*
* Revision 1.37 2005/07/11 08:19:54 f_limousin
* New function axi_break() to compute the departure to axisymmetry.
*
* Revision 1.36 2005/06/09 07:56:24 f_limousin
* Implement a new function sol_elliptic_boundary() and
* Vector::poisson_boundary(...) which solve the vectorial poisson
* equation (method 6) with an inner boundary condition.
*
* Revision 1.35 2005/05/12 14:48:55 f_limousin
* New boundary condition for the lapse : boundary_nn_lapl().
*
* Revision 1.34 2005/04/29 14:05:03 f_limousin
* Important changes : manage the dependances between quantities (for
* instance psi and psi4). New function write_global(ost).
*
* Revision 1.33 2005/04/15 11:20:39 jl_jaramillo
* Function adapt_hor(double c_min, double c_max) for adapting a given surface
* to the excised surface
*
* Revision 1.32 2005/04/08 12:15:38 f_limousin
* Function set_psi(). And dependance in phi.
*
* Revision 1.31 2005/04/03 19:48:38 f_limousin
* Implementation of set_psi(psi_in).
*
* Revision 1.30 2005/04/02 15:50:08 f_limousin
* New data member nz (number of zones). Delete ww.
*
* Revision 1.29 2005/03/31 09:48:04 f_limousin
* New functions compute_ww(..) and aa_kerr_ww() and new data ww.
*
* Revision 1.28 2005/03/28 19:45:41 f_limousin
* Implement Isol_hor::aa_kerr_perturb(...) and new member aa_quad_evol.
*
* Revision 1.27 2005/03/24 16:50:40 f_limousin
* Add parameters solve_shift and solve_psi in par_isol.d and in function
* init_dat(...). Implement Isolhor::kerr_perturb().
*
* Revision 1.26 2005/03/10 16:57:01 f_limousin
* Improve the convergence of the code coal_bh.
*
* Revision 1.25 2005/03/10 10:19:42 f_limousin
* Add the regularisation of the shift in the case of a single black hole
* and lapse zero on the horizon.
*
* Revision 1.24 2005/03/09 10:28:37 f_limousin
* Delete functions init_data_b_neumann(...) and init_data_berlin(...)
* --> New parameter solve_lapse in the function init_data(...).
* New function update_aa().
*
* Revision 1.23 2005/03/06 16:59:33 f_limousin
* New function Isol_hor::aa() (the one belonging to the class
* Time_slice_conf need to compute the time derivative of hh and thus
* cannot work in the class Isol_hor).
*
* Revision 1.22 2005/03/04 09:39:31 f_limousin
* Implement the constructor from a file, operator>>, operator<<
* and function sauve in the class Bin_hor.
*
* Revision 1.21 2005/03/03 10:25:16 f_limousin
* In the class Isol_hor :
* - Add the boost in x and z-direction (members boost_x and boost_z,
* and functions get_boost_x(), set_boost_x(double))
* - Add function area_hor()
* - Put the boundary conditions for the lapse, psi and beta in
* the parameter file.
* In the class bin_hor :
* - Introduce function to compute global quantities as ADM mass,
* Komar mass and angular momentum.
*
* Revision 1.20 2005/02/24 17:22:53 f_limousin
* Suppression of the function beta_bound_cart().
* The boundary conditions for psi, N and beta are now some parameters
* in par_init.D and par_coal.d.
*
* Revision 1.19 2005/02/07 10:30:09 f_limousin
* Add the regularisation in the case N=0 on the horizon.
*
* Revision 1.18 2004/12/31 15:33:37 f_limousin
* Change the constructor from a file and the standard constructor.
*
* Revision 1.17 2004/12/29 16:30:00 f_limousin
* Improve comments for doxygen
*
* Revision 1.16 2004/12/29 16:10:25 f_limousin
* Add the new class Bin_hor.
*
* Revision 1.14 2004/11/24 19:32:05 jl_jaramillo
* Method for initial data with Berlin boundary conditions
*
* Revision 1.13 2004/11/18 10:53:03 jl_jaramillo
* Declarations for Berlin boundary conditions
*
* Revision 1.12 2004/11/05 11:01:13 f_limousin
* Delete arguments ener_dens, mom_dens and trace_stress in all functions
* source_nn, source_psi, source_beta, init_data. Delete also
* argument partial_save in function save.
*
* Revision 1.11 2004/11/05 10:11:23 f_limousin
* The member Metric met_gamt replace Sym_tensor gamt.
*
* Revision 1.10 2004/11/03 17:15:46 f_limousin
* Change the standart constructor. Add 4 memebers : trK, trK_point,
* gamt and gamt_point.
* Add also a constructor from a file.
*
* Revision 1.9 2004/11/02 17:42:33 f_limousin
* New method sauve(...) to save in a binary file.
*
* Revision 1.8 2004/11/02 16:15:12 f_limousin
* Add new argument ang_vel in function init_dat(...).
*
* Revision 1.7 2004/10/29 15:46:14 jl_jaramillo
* Remove 2 members, add ADM angular momentum and change name
* of functions.
*
* Revision 1.6 2004/10/01 16:51:16 f_limousin
* Pure Dirichlet boundary condition added
*
* Revision 1.5 2004/09/28 16:03:58 f_limousin
* Add parameter niter in the parameter file par_hor.d. It appears in
* argument of the function init_data_schwarz(...).
*
* Revision 1.4 2004/09/17 13:35:25 f_limousin
* Introduction of relaxation in init_data_schwarz
*
* Revision 1.3 2004/09/16 08:36:57 f_limousin
* New boundary conditions for lapse and psi.
*
* Revision 1.2 2004/09/09 17:04:27 jl_jaramillo
* Elimination of _ih
*
*
*
* $Header: /cvsroot/Lorene/C++/Include/isol_hor.h,v 1.54 2014/10/13 08:52:35 j_novak Exp $
*
*/
#include "time_slice.h"
#include "proto.h"
#include "headcpp.h"
#include "cmp.h"
#include "evolution.h"
namespace Lorene {
class Sym_tensor_trans ;
class Sym_tensor ;
class Vector ;
class Scalar ;
class Metric ;
class Metric_flat ;
class Base_vect ;
class Map ;
class Tbl ;
class Time_slice ;
class Time_slice_conf ;
//----------------------------//
// class Isol_Hor //
//----------------------------//
class Isol_hor : public Time_slice_conf {
// Data :
// -----
protected:
Map_af& mp ;
int nz ;
double radius ;
double omega ;
double boost_x ;
double boost_z ;
double regul ;
mutable Evolution_std<Scalar> n_auto_evol ;
mutable Evolution_std<Scalar> n_comp_evol ;
mutable Evolution_std<Scalar> psi_auto_evol ;
mutable Evolution_std<Scalar> psi_comp_evol ;
mutable Evolution_std<Vector> dn_evol ;
mutable Evolution_std<Vector> dpsi_evol ;
mutable Evolution_std<Vector> beta_auto_evol ;
mutable Evolution_std<Vector> beta_comp_evol ;
mutable Evolution_std<Sym_tensor> aa_auto_evol ;
mutable Evolution_std<Sym_tensor> aa_comp_evol ;
mutable Evolution_std<Sym_tensor> aa_nn ;
mutable Evolution_std<Scalar> aa_quad_evol ;
Metric met_gamt ;
Sym_tensor gamt_point ;
Scalar trK ;
Scalar trK_point ;
Scalar decouple ;
// Constructors - Destructor
// -------------------------
public:
Isol_hor(Map_af& mpi, int depth_in = 3) ;
Isol_hor(Map_af& mpi, const Scalar& lapse_in, const Scalar& psi_in,
const Vector& shift_in, const Sym_tensor& aa_in,
const Metric& gamt, const Sym_tensor& gamt_point,
const Scalar& trK, const Scalar& trK_point,
const Metric_flat& ff_in, int depth_in = 3) ;
Isol_hor(const Isol_hor& ) ;
Isol_hor (Map_af& mp, FILE* fich,
bool partial_read, int depth_in = 3) ;
virtual ~Isol_hor() ;
// Mutators / assignment
// ---------------------
public:
void operator=(const Isol_hor&) ;
public:
const Map_af& get_mp() const {return mp;} ;
Map_af& set_mp() {return mp; } ;
double get_radius() const {return radius;} ;
void set_radius(double rad) {radius = rad ;} ;
double get_omega() const {return omega ;} ;
void set_omega(double ome) {omega = ome ;} ;
double get_boost_x() const {return boost_x ;} ;
void set_boost_x(double bo) {boost_x = bo ;} ;
double get_boost_z() const {return boost_z ;} ;
void set_boost_z(double bo) {boost_z = bo ;} ;
// Accessors
// ---------
public:
virtual const Scalar& n_auto() const ;
virtual const Scalar& n_comp() const ;
virtual const Scalar& psi_auto() const ;
virtual const Scalar& psi_comp() const ;
virtual const Vector& dnn() const ;
virtual const Vector& dpsi() const ;
virtual const Vector& beta_auto() const ;
virtual const Vector& beta_comp() const ;
virtual const Sym_tensor& aa_auto() const ;
virtual const Sym_tensor& aa_comp() const ;
virtual const Scalar& aa_quad() const ;
virtual const Metric& tgam() const {return met_gamt ;}
const Scalar get_decouple() const {return decouple ;}
public:
void n_comp (const Isol_hor& comp) ;
void psi_comp (const Isol_hor& comp) ;
void beta_comp (const Isol_hor& comp) ;
double viriel_seul () const ;
void init_bhole () ;
void init_met_trK() ;
void init_bhole_seul () ;
void set_psi(const Scalar& psi_in) ;
void set_nn(const Scalar& nn_in) ;
void set_gamt(const Metric& gam_tilde) ;
// Physical parameters
//--------------------
public:
const Vector radial_vect_hor() const ;
const Vector tradial_vect_hor() const ;
const Scalar b_tilde() const ;
const Scalar darea_hor() const ;
double area_hor() const ;
double radius_hor() const ;
double ang_mom_hor() const ;
double mass_hor() const ;
double kappa_hor() const ;
double omega_hor() const ;
double ang_mom_adm() const ;
Scalar expansion() const ;
//Computational methods
//---------------------
public:
/* function to compute initial data for a single black hole
* @param bound_nn boundary condition for the lapse
* @param lim_nn value of the boundary condition for the lapse
* @param bound_psi boundary condition for \f$ \Psi \f$
* @param bound_beta boundary condition for the shift
* @param solve_lapse do we solve the equation for the lapse ?
* @param precis precision for the convergence
* @param relax relaxation
* @param niter number of iterations
*/
void init_data(int bound_nn, double lim_nn, int bound_psi, int bound_beta,
int solve_lapse, int solve_psi, int solve_shift,
double precis = 1.e-12,
double relax_nn = 0.5, double relax_psi = 0.5,
double relax_beta = 0.5, int niter = 100) ;
void init_data_loop(int bound_nn, double lim_nn, int bound_psi,
int bound_beta, int solve_lapse, int solve_psi,
int solve_shift, double precis= 1.e-12,
double precis_loop= 1.e-12,
double relax_nn = 1., double relax_psi= 1.,
double relax_beta = 1., double relax_loop = 1.,
int niter = 100) ;
void init_data_spher(int bound_nn, double lim_nn, int bound_psi,
int bound_beta, int solve_lapse, int solve_psi,
int solve_shift, double precis = 1.e-12,
double relax = 1., int niter = 100) ;
void init_data_alt(int bound_nn, double lim_nn, int bound_psi,
int bound_beta, int solve_lapse, int solve_psi,
int solve_shift, double precis = 1.e-12,
double relax = 1., int niter = 100) ;
void init_data_CTS_gen(int bound_nn, double lim_nn, int bound_psi, int bound_beta,
int solve_lapse, int solve_psi, int solve_shift,
double precis = 1.e-12, double relax_nn = 1.,
double relax_psi = 1., double relax_beta = 1.,
int niter = 100, double a = 1., double zeta = 4.) ;
//Sources
//-------
const Scalar source_psi() const ;
const Scalar source_nn() const ;
const Vector source_beta() const ;
const Scalar source_b_tilde() const ;
const Vector source_vector_b() const ;
// BOUNDARY CONDITIONS
//--------------------
const Valeur boundary_psi_Dir_evol() const ;
const Valeur boundary_psi_Neu_evol() const ;
const Valeur boundary_psi_Dir_spat() const ;
const Valeur boundary_psi_Neu_spat() const ;
const Valeur boundary_psi_app_hor() const ;
const Valeur boundary_psi_Dir() const ;
const Valeur boundary_nn_Dir_kk() const ;
const Valeur boundary_nn_Neu_kk(int nn = 1) const ;
const Valeur boundary_nn_Neu_Cook() const ;
const Valeur boundary_nn_Dir_eff(double aa) const ;
const Valeur boundary_nn_Dir_lapl(int mer = 1) const ;
const Valeur boundary_nn_Neu_eff(double aa) const ;
const Valeur boundary_nn_Dir(double aa) const ;
const Valeur boundary_beta_r() const ;
const Valeur boundary_beta_theta() const ;
const Valeur boundary_beta_phi(double om) const ;
const Valeur boundary_beta_x(double om) const ;
const Valeur boundary_beta_y(double om) const ;
const Valeur boundary_beta_z() const ;
const Valeur beta_boost_x() const ;
const Valeur beta_boost_z() const ;
const Vector vv_bound_cart(double om) const ;
const Vector vv_bound_cart_bin(double om, int hole = 0) const ;
const Valeur boundary_vv_x(double om) const ;
const Valeur boundary_vv_y(double om) const ;
const Valeur boundary_vv_z(double om) const ;
const Valeur boundary_vv_x_bin(double om, int hole = 0) const ;
const Valeur boundary_vv_y_bin(double om, int hole = 0) const ;
const Valeur boundary_vv_z_bin(double om, int hole = 0) const ;
const Valeur boundary_b_tilde_Neu() const ;
const Valeur boundary_b_tilde_Dir() const ;
void update_aa() ;
double regularisation (const Vector& shift_auto, const Vector& shift_comp,
double ang_vel) ;
double regularise_one() ;
void met_kerr_perturb() ;
/* Perturbation of Kerr using \f$ A^{ij}A_{ij} \f$ from
* equation (14) of Dain (2002).
* @param mm mass of the Kerr black hole metric.
* @param aa rotation parameter of the Kerr black hole metric.
*/
void aa_kerr_ww(double mm, double aa) ;
double axi_break() const ;
/* Calculation of the outermost trapped surface and adaptation
* of all necessary quantities
*/
void adapt_hor(double c_min, double c_max) ;
// Outputs
// -------
protected:
virtual ostream& operator>>(ostream& ) const ;
public :
virtual void sauve(FILE* fich, bool partial_save) const ;
friend class Bin_hor ;
};
class Single_hor {
// Data :
// -----
protected:
Map_af& mp ;
int nz ;
double radius ;
double omega ;
double regul ;
Scalar n_auto ;
Scalar n_comp ;
Scalar nn ;
Scalar psi_auto ;
Scalar psi_comp ;
Scalar psi ;
mutable Scalar* p_psi4 ;
Vector dn ;
Vector dpsi ;
Vector beta_auto ;
Vector beta_comp ;
Vector beta ;
mutable Metric* p_gam ;
Sym_tensor aa_auto ;
Sym_tensor aa_comp ;
Sym_tensor aa ;
mutable Sym_tensor* p_k_dd ;
Metric tgam ;
Metric_flat ff ;
Sym_tensor hh ;
Sym_tensor gamt_point ;
Scalar trK ;
Scalar trK_point ;
Scalar decouple ;
// Constructors - Destructor
// -------------------------
public:
Single_hor(Map_af& mpi) ;
Single_hor(const Single_hor& ) ;
Single_hor (Map_af& mp, FILE* fich) ;
virtual ~Single_hor() ;
// Mutators / assignment
// ---------------------
public:
void operator=(const Single_hor&) ;
public:
const Map_af& get_mp() const {return mp;} ;
Map_af& set_mp() {return mp; } ;
double get_radius() const {return radius;} ;
void set_radius(double rad) {radius = rad ;} ;
double get_omega() const {return omega ;} ;
void set_omega(double ome) {omega = ome ;} ;
// Memory management
// -----------------
protected:
void del_deriv() const ;
void set_der_0x0() const ;
// Accessors
// ---------
public:
const Scalar& get_n_auto() const ;
const Scalar& get_n_comp() const ;
const Scalar& get_nn() const ;
const Scalar& get_psi_auto() const ;
const Scalar& get_psi_comp() const ;
const Scalar& get_psi() const ;
const Scalar& get_psi4() const ;
const Vector& get_dn() const ;
const Vector& get_dpsi() const ;
const Vector& get_beta_auto() const ;
const Vector& get_beta_comp() const ;
const Vector& get_beta() const ;
const Sym_tensor& get_aa_auto() const ;
const Sym_tensor& get_aa_comp() const ;
const Sym_tensor& get_aa() const ;
const Metric& get_tgam() const {return tgam ;}
const Metric& get_gam() const ;
const Sym_tensor& get_k_dd() const ;
const Scalar get_decouple() const {return decouple ;}
public:
void n_comp_import (const Single_hor& comp) ;
void psi_comp_import (const Single_hor& comp) ;
void beta_comp_import (const Single_hor& comp) ;
double viriel_seul () const ;
void init_bhole () ;
void init_met_trK() ;
void init_bhole_seul () ;
void set_psi_auto(const Scalar& psi_in) ;
void set_n_auto(const Scalar& nn_in) ;
void set_beta_auto(const Scalar& shift_in) ;
void set_aa_auto(const Scalar& aa_auto_in) ;
void set_aa_comp(const Scalar& aa_comp_in) ;
void set_aa(const Scalar& aa_in) ;
// Physical parameters
//--------------------
public:
const Scalar b_tilde() const ;
const Scalar darea_hor() const ;
double area_hor() const ;
double radius_hor() const ;
double ang_mom_hor() const ;
double mass_hor() const ;
double kappa_hor() const ;
double omega_hor() const ;
double ang_mom_adm() const ;
Scalar expansion() const ;
// BOUNDARY CONDITIONS
//--------------------
const Valeur boundary_psi_app_hor() const ;
const Valeur boundary_nn_Dir(double aa) const ;
const Valeur boundary_nn_Neu(double aa) const ;
const Valeur boundary_beta_x(double om_orb, double om_loc) const ;
const Valeur boundary_beta_y(double om_orb, double om_loc) const ;
const Valeur boundary_beta_z() const ;
double regularisation (const Vector& shift_auto, const Vector& shift_comp,
double ang_vel) ;
double regularise_one() ;
public :
virtual void sauve(FILE* fich) const ;
friend class Bin_hor ;
};
class Bin_hor {
// data :
private:
// lThe two black holes
Single_hor hole1 ;
Single_hor hole2 ;
Single_hor* holes[2] ;
double omega ;
public:
Bin_hor(Map_af& mp1, Map_af& mp2) ;
Bin_hor(const Bin_hor& ) ;
Bin_hor (Map_af& mp1, Map_af& mp2, FILE* fich) ;
virtual ~Bin_hor() ;
public :
void sauve(FILE* fich) const ;
void write_global(ostream&, double lim_nn, int bound_nn,
int bound_psi, int bound_beta, double alpha) const ;
public:
void operator=(const Bin_hor&) ;
Single_hor& set(int i)
{ assert( (i==1) || (i==2) );
return *holes[i-1] ;} ;
void set_omega(double ome) {omega = ome ;
hole1.set_omega (ome) ;
hole2.set_omega (ome) ;} ;
public:
*/ const Single_hor& operator()(int i) const
{ assert( (i==1) || (i==2) );
return *holes[i-1] ;} ;
double get_omega() const {return omega; } ;
void init_bin_hor() ;
double viriel() const ;
void extrinsic_curvature () ;
void decouple () ;
public:
void set_statiques (double precis, double relax, int bound_nn,
double lim_nn, int bound_psi) ;
double coal (double ang_vel, double relax, int nb_om,
int nb_it, int bound_nn, double lim_nn,
int bound_psi, int bound_beta, double omega_eff,
double alpha,
ostream& fich_iteration, ostream& fich_correction,
ostream& fich_viriel, ostream& fich_kss,
int step, int search_mass, double mass_irr,
const int sortie = 0) ;
void solve_lapse (double precis, double relax, int bound_nn,
double lim_nn) ;
void solve_psi (double precis, double relax, int bound_psi) ;
void solve_shift (double precis, double relax, int bound_beta,
double omega_eff) ;
void import_bh (const Bin_hor& bin) ;
double adm_mass() const ;
double komar_mass() const ;
double ang_mom_hor() const ;
double ang_mom_adm() const ;
double proper_distance(const int nr = 65) const ;
Sym_tensor hh_Samaya_hole1() ;
Sym_tensor hh_Samaya_hole2() ;
void set_hh_Samaya() ;
} ;
}
#endif
Lorene::Base_vect
Vectorial bases (triads) with respect to which the tensorial components are defined.
Definition base_vect.h:105
Lorene::Bin_hor
Definition isol_hor.h:1337
Lorene::Bin_hor::viriel
double viriel() const
Computes the viriel error, that is the difference between the ADM and the Komar masses,...
Definition binhor_viriel.C:71
Lorene::Bin_hor::solve_psi
void solve_psi(double precis, double relax, int bound_psi)
Solves the equation for the conformal factor : The fields are the total values excpet those with subs...
Definition binhor_equations.C:301
Lorene::Bin_hor::omega
double omega
Angular velocity.
Definition isol_hor.h:1348
Lorene::Bin_hor::coal
double coal(double ang_vel, double relax, int nb_om, int nb_it, int bound_nn, double lim_nn, int bound_psi, int bound_beta, double omega_eff, double alpha, ostream &fich_iteration, ostream &fich_correction, ostream &fich_viriel, ostream &fich_kss, int step, int search_mass, double mass_irr, const int sortie=0)
Solves the equation for a particular angular velocity, the systeme being initialized to Misner-Lindqu...
Definition binhor_coal.C:133
Lorene::Bin_hor::hole1
Single_hor hole1
Black hole one.
Definition isol_hor.h:1342
Lorene::Bin_hor::holes
Single_hor * holes[2]
Array on the black holes.
Definition isol_hor.h:1346
Lorene::Bin_hor::hole2
Single_hor hole2
Black hole two.
Definition isol_hor.h:1343
Lorene::Bin_hor::ang_mom_hor
double ang_mom_hor() const
Calculates the angular momentum of the black hole using the formula at the horizon.
Definition binhor_global.C:118
Lorene::Bin_hor::init_bin_hor
void init_bin_hor()
Initialisation of the system.
Definition bin_hor.C:161
Lorene::Bin_hor::adm_mass
double adm_mass() const
Calculates the ADM mass of the system.
Definition binhor_global.C:83
Lorene::Bin_hor::solve_shift
void solve_shift(double precis, double relax, int bound_beta, double omega_eff)
Solves the equation for the shift, using the Oohara-Nakarmure scheme : The fields are the total value...
Definition binhor_equations.C:453
Lorene::Bin_hor::import_bh
void import_bh(const Bin_hor &bin)
Function to initialize a Bin_hor from a solution computed with a smaller number of colocation points.
Lorene::Bin_hor::ang_mom_adm
double ang_mom_adm() const
Calculates the angular momentum of the black hole.
Definition binhor_global.C:184
Lorene::Bin_hor::sauve
void sauve(FILE *fich) const
Total or partial saves in a binary file.
Definition bin_hor.C:151
Lorene::Bin_hor::komar_mass
double komar_mass() const
Calculates the Komar mass of the system using : .
Definition binhor_global.C:101
Lorene::Bin_hor::decouple
void decouple()
Calculates decouple which is used to obtain tkij_auto and tkij_comp.
Definition binhor_kij.C:473
Lorene::Bin_hor::operator=
void operator=(const Bin_hor &)
Affectation operator.
Definition bin_hor.C:139
Lorene::Bin_hor::Bin_hor
Bin_hor(Map_af &mp1, Map_af &mp2)
Standard constructor.
Definition bin_hor.C:97
Lorene::Bin_hor::set_statiques
void set_statiques(double precis, double relax, int bound_nn, double lim_nn, int bound_psi)
Initialize the systeme to Misner Lindquist solution, that is solving for N and in the case .
Definition binhor_coal.C:96
Lorene::Bin_hor::set_hh_Samaya
void set_hh_Samaya()
Calculation of the Post-Newtonian correction to .
Definition binhor_hh.C:786
Lorene::Bin_hor::hh_Samaya_hole2
Sym_tensor hh_Samaya_hole2()
Calculation of the hole2 part of the Post-Newtonian correction to .
Definition binhor_hh.C:453
Lorene::Bin_hor::solve_lapse
void solve_lapse(double precis, double relax, int bound_nn, double lim_nn)
Solves the equation for the lapse : The fields are the total values except those with subscript ,...
Definition binhor_equations.C:124
Lorene::Bin_hor::~Bin_hor
virtual ~Bin_hor()
Destructor.
Definition bin_hor.C:132
Lorene::Bin_hor::set
Single_hor & set(int i)
Read/write of a component of the system.
Definition isol_hor.h:1403
Lorene::Bin_hor::set_omega
void set_omega(double ome)
Sets the orbital velocity to ome.
Definition isol_hor.h:1409
Lorene::Bin_hor::proper_distance
double proper_distance(const int nr=65) const
Calculation of the proper distance between the two spheres of inversion, along the x axis.
Lorene::Bin_hor::operator()
const Single_hor & operator()(int i) const
Read only of a component of the system.
Definition isol_hor.h:1417
Lorene::Bin_hor::write_global
void write_global(ostream &, double lim_nn, int bound_nn, int bound_psi, int bound_beta, double alpha) const
Write global quantities in a formatted file.
Definition bin_hor.C:178
Lorene::Bin_hor::extrinsic_curvature
void extrinsic_curvature()
Calculation of the extrinsic curvature tensor.
Definition binhor_kij.C:88
Lorene::Bin_hor::hh_Samaya_hole1
Sym_tensor hh_Samaya_hole1()
Calculation of the hole1 part of the Post-Newtonian correction to .
Definition binhor_hh.C:63
Lorene::Bin_hor::get_omega
double get_omega() const
Returns the angular velocity.
Definition isol_hor.h:1422
Lorene::Evolution_std
Time evolution with partial storage (*** under development ***).
Definition evolution.h:371
Lorene::Isol_hor
Spacelike time-slice of an Isolated Horizon in a 3+1 spacetime with conformal decomposition.
Definition isol_hor.h:254
Lorene::Isol_hor::dnn
virtual const Vector & dnn() const
Covariant derivative of the lapse function at the current time step jtime.
Definition isol_hor.C:479
Lorene::Isol_hor::expansion
Scalar expansion() const
Expansion of the outgoing null normal ( ).
Definition phys_param.C:263
Lorene::Isol_hor::met_kerr_perturb
void met_kerr_perturb()
Initialisation of the metric tilde from equation (15) of Dain (2002).
Definition isol_hor.C:895
Lorene::Isol_hor::get_mp
const Map_af & get_mp() const
Returns the mapping (readonly).
Definition isol_hor.h:418
Lorene::Isol_hor::set_omega
void set_omega(double ome)
Sets the angular velocity to ome .
Definition isol_hor.h:440
Lorene::Isol_hor::aa_quad
virtual const Scalar & aa_quad() const
Conformal representation .
Definition isol_hor.C:884
Lorene::Isol_hor::aa_nn
Evolution_std< Sym_tensor > aa_nn
Values at successive time steps of the components .
Definition isol_hor.h:320
Lorene::Isol_hor::tgam
virtual const Metric & tgam() const
Conformal metric Returns the value at the current time step (jtime ).
Definition isol_hor.h:514
Lorene::Isol_hor::init_bhole
void init_bhole()
Sets the values of the fields to :
Definition isol_hor.C:733
Lorene::Isol_hor::boundary_vv_z_bin
const Valeur boundary_vv_z_bin(double om, int hole=0) const
Component z of boundary value of .
Definition bound_hor.C:1625
Lorene::Isol_hor::boundary_vv_y
const Valeur boundary_vv_y(double om) const
Component y of boundary value of .
Definition bound_hor.C:1521
Lorene::Isol_hor::gamt_point
Sym_tensor gamt_point
Time derivative of the 3-metric tilde.
Definition isol_hor.h:329
Lorene::Isol_hor::boundary_psi_Dir_evol
const Valeur boundary_psi_Dir_evol() const
Dirichlet boundary condition for (evolution).
Definition bound_hor.C:174
Lorene::Isol_hor::boundary_psi_Dir_spat
const Valeur boundary_psi_Dir_spat() const
Dirichlet boundary condition for (spatial).
Definition bound_hor.C:231
Lorene::Isol_hor::source_beta
const Vector source_beta() const
Source for .
Definition sources_hor.C:190
Lorene::Isol_hor::get_omega
double get_omega() const
Returns the angular velocity.
Definition isol_hor.h:436
Lorene::Isol_hor::regularise_one
double regularise_one()
Corrects the shift in the innermost shell, so that it remains and that equals zero on the horizon.
Definition regularisation.C:183
Lorene::Isol_hor::aa_quad_evol
Evolution_std< Scalar > aa_quad_evol
Values at successive time steps of the components .
Definition isol_hor.h:323
Lorene::Isol_hor::omega
double omega
Angular velocity in LORENE's units.
Definition isol_hor.h:269
Lorene::Isol_hor::boundary_nn_Neu_eff
const Valeur boundary_nn_Neu_eff(double aa) const
Neumann boundary condition on nn (effectif) .
Definition bound_hor.C:622
Lorene::Isol_hor::regul
double regul
Intensity of the correction on the shift vector.
Definition isol_hor.h:278
Lorene::Isol_hor::boundary_vv_x
const Valeur boundary_vv_x(double om) const
Component x of boundary value of .
Definition bound_hor.C:1493
Lorene::Isol_hor::boundary_beta_y
const Valeur boundary_beta_y(double om) const
Component y of boundary value of .
Definition bound_hor.C:1071
Lorene::Isol_hor::boundary_psi_app_hor
const Valeur boundary_psi_app_hor() const
Neumann boundary condition for (spatial).
Definition bound_hor.C:297
Lorene::Isol_hor::n_comp
virtual const Scalar & n_comp() const
Lapse function at the current time step jtime.
Definition isol_hor.C:461
Lorene::Isol_hor::source_vector_b
const Vector source_vector_b() const
Source for vector_b.
Lorene::Isol_hor::trK
Scalar trK
Trace of the extrinsic curvature.
Definition isol_hor.h:332
Lorene::Isol_hor::operator>>
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator<<).
Definition isol_hor.C:378
Lorene::Isol_hor::boundary_beta_x
const Valeur boundary_beta_x(double om) const
Component x of boundary value of .
Definition bound_hor.C:1021
Lorene::Isol_hor::darea_hor
const Scalar darea_hor() const
Element of area of the horizon.
Definition phys_param.C:146
Lorene::Isol_hor::met_gamt
Metric met_gamt
3 metric tilde
Definition isol_hor.h:326
Lorene::Isol_hor::boundary_nn_Neu_Cook
const Valeur boundary_nn_Neu_Cook() const
Neumann boundary condition for N using Cook's boundary condition.
Definition bound_hor.C:489
Lorene::Isol_hor::boundary_b_tilde_Neu
const Valeur boundary_b_tilde_Neu() const
Neumann boundary condition for b_tilde.
Definition bound_hor.C:1169
Lorene::Isol_hor::b_tilde
const Scalar b_tilde() const
Radial component of the shift with respect to the conformal metric.
Definition phys_param.C:136
Lorene::Isol_hor::source_psi
const Scalar source_psi() const
Source for .
Definition sources_hor.C:108
Lorene::Isol_hor::dpsi_evol
Evolution_std< Vector > dpsi_evol
Values at successive time steps of the covariant derivative of the conformal factor .
Definition isol_hor.h:298
Lorene::Isol_hor::beta_boost_x
const Valeur beta_boost_x() const
Boundary value for a boost in x-direction.
Definition bound_hor.C:1144
Lorene::Isol_hor::psi_auto_evol
Evolution_std< Scalar > psi_auto_evol
Values at successive time steps of the conformal factor .
Definition isol_hor.h:287
Lorene::Isol_hor::boundary_beta_r
const Valeur boundary_beta_r() const
Component r of boundary value of .
Definition bound_hor.C:897
Lorene::Isol_hor::nz
int nz
Number of zones.
Definition isol_hor.h:263
Lorene::Isol_hor::boundary_vv_z
const Valeur boundary_vv_z(double om) const
Component z of boundary value of .
Definition bound_hor.C:1547
Lorene::Isol_hor::boundary_vv_y_bin
const Valeur boundary_vv_y_bin(double om, int hole=0) const
Component y of boundary value of .
Definition bound_hor.C:1599
Lorene::Isol_hor::ang_mom_hor
double ang_mom_hor() const
Angular momentum (modulo).
Definition phys_param.C:178
Lorene::Isol_hor::set_boost_z
void set_boost_z(double bo)
Sets the boost velocity in z-direction to bo .
Definition isol_hor.h:458
Lorene::Isol_hor::aa_comp_evol
Evolution_std< Sym_tensor > aa_comp_evol
Values at successive time steps of the components of the conformal representation of the traceless p...
Definition isol_hor.h:316
Lorene::Isol_hor::set_psi
void set_psi(const Scalar &psi_in)
Sets the conformal factor relating the physical metric to the conformal one: .
Definition isol_hor.C:515
Lorene::Isol_hor::viriel_seul
double viriel_seul() const
Computes the viriel error, that is the difference between the ADM and the Komar masses,...
Lorene::Isol_hor::beta_comp
virtual const Vector & beta_comp() const
Shift function at the current time step jtime.
Definition isol_hor.C:497
Lorene::Isol_hor::init_met_trK
void init_met_trK()
Sets the 3-metric tilde to the flat metric and gamt_point, trK and trK_point to zero.
Definition isol_hor.C:785
Lorene::Isol_hor::trK_point
Scalar trK_point
Time derivative of the trace of the extrinsic curvature.
Definition isol_hor.h:335
Lorene::Isol_hor::beta_auto_evol
Evolution_std< Vector > beta_auto_evol
Values at successive time steps of the shift function .
Definition isol_hor.h:301
Lorene::Isol_hor::dpsi
virtual const Vector & dpsi() const
Covariant derivative with respect to the flat metric of the conformal factor at the current time ste...
Definition isol_hor.C:485
Lorene::Isol_hor::set_nn
void set_nn(const Scalar &nn_in)
Sets the lapse.
Definition isol_hor.C:540
Lorene::Isol_hor::boundary_nn_Dir_lapl
const Valeur boundary_nn_Dir_lapl(int mer=1) const
Dirichlet boundary condition for N fixing the divergence of the connection form .
Definition bound_hor.C:693
Lorene::Isol_hor::boost_z
double boost_z
Boost velocity in z-direction.
Definition isol_hor.h:275
Lorene::Isol_hor::boundary_beta_theta
const Valeur boundary_beta_theta() const
Component theta of boundary value of .
Definition bound_hor.C:938
Lorene::Isol_hor::psi_auto
virtual const Scalar & psi_auto() const
Conformal factor at the current time step jtime.
Definition isol_hor.C:467
Lorene::Isol_hor::beta_boost_z
const Valeur beta_boost_z() const
Boundary value for a boost in z-direction.
Definition bound_hor.C:1155
Lorene::Isol_hor::boundary_vv_x_bin
const Valeur boundary_vv_x_bin(double om, int hole=0) const
Component x of boundary value of .
Definition bound_hor.C:1572
Lorene::Isol_hor::ang_mom_adm
double ang_mom_adm() const
ADM angular Momentum.
Definition phys_param.C:248
Lorene::Isol_hor::aa_comp
virtual const Sym_tensor & aa_comp() const
Conformal representation of the traceless part of the extrinsic curvature: Returns the value at the ...
Definition isol_hor.C:509
Lorene::Isol_hor::set_radius
void set_radius(double rad)
Sets the radius of the horizon to rad .
Definition isol_hor.h:431
Lorene::Isol_hor::vv_bound_cart
const Vector vv_bound_cart(double om) const
Vector for boundary conditions in cartesian.
Definition bound_hor.C:1251
Lorene::Isol_hor::decouple
Scalar decouple
Function used to construct from the total .
Definition isol_hor.h:345
Lorene::Isol_hor::kappa_hor
double kappa_hor() const
Surface gravity.
Definition phys_param.C:218
Lorene::Isol_hor::update_aa
void update_aa()
Conformal representation of the traceless part of the extrinsic curvature: .
Definition isol_hor.C:842
Lorene::Isol_hor::boundary_beta_z
const Valeur boundary_beta_z() const
Component z of boundary value of .
Definition bound_hor.C:1121
Lorene::Isol_hor::n_auto
virtual const Scalar & n_auto() const
Lapse function at the current time step jtime.
Definition isol_hor.C:455
Lorene::Isol_hor::~Isol_hor
virtual ~Isol_hor()
Destructor.
Definition isol_hor.C:336
Lorene::Isol_hor::radius
double radius
Radius of the horizon in LORENE's units.
Definition isol_hor.h:266
Lorene::Isol_hor::aa_auto
virtual const Sym_tensor & aa_auto() const
Conformal representation of the traceless part of the extrinsic curvature: Returns the value at the ...
Definition isol_hor.C:503
Lorene::Isol_hor::get_decouple
const Scalar get_decouple() const
Returns the function used to construct tkij_auto from tkij_tot .
Definition isol_hor.h:519
Lorene::Isol_hor::boundary_nn_Dir_kk
const Valeur boundary_nn_Dir_kk() const
Dirichlet boundary condition for N using the extrinsic curvature.
Definition bound_hor.C:371
Lorene::Isol_hor::mass_hor
double mass_hor() const
Mass computed at the horizon.
Definition phys_param.C:206
Lorene::Isol_hor::init_bhole_seul
void init_bhole_seul()
Initiates for a single black hole.
Definition isol_hor.C:797
Lorene::Isol_hor::aa_auto_evol
Evolution_std< Sym_tensor > aa_auto_evol
Values at successive time steps of the components of the conformal representation of the traceless p...
Definition isol_hor.h:310
Lorene::Isol_hor::set_gamt
void set_gamt(const Metric &gam_tilde)
Sets the conformal metric to gam_tilde.
Definition isol_hor.C:550
Lorene::Isol_hor::radius_hor
double radius_hor() const
Radius of the horizon.
Definition phys_param.C:167
Lorene::Isol_hor::get_boost_x
double get_boost_x() const
Returns the boost velocity in x-direction.
Definition isol_hor.h:445
Lorene::Isol_hor::tradial_vect_hor
const Vector tradial_vect_hor() const
Vector radial normal tilde.
Definition phys_param.C:116
Lorene::Isol_hor::mp
Map_af & mp
Affine mapping.
Definition isol_hor.h:260
Lorene::Isol_hor::boundary_psi_Neu_spat
const Valeur boundary_psi_Neu_spat() const
Neumann boundary condition for (spatial).
Definition bound_hor.C:267
Lorene::Isol_hor::boundary_beta_phi
const Valeur boundary_beta_phi(double om) const
Component phi of boundary value of .
Definition bound_hor.C:978
Lorene::Isol_hor::psi_comp
virtual const Scalar & psi_comp() const
Conformal factor at the current time step jtime.
Definition isol_hor.C:473
Lorene::Isol_hor::radial_vect_hor
const Vector radial_vect_hor() const
Vector radial normal.
Definition phys_param.C:95
Lorene::Isol_hor::axi_break
double axi_break() const
Breaking of the axial symmetry on the horizon.
Definition isol_hor.C:1087
Lorene::Isol_hor::get_boost_z
double get_boost_z() const
Returns the boost velocity in z-direction.
Definition isol_hor.h:454
Lorene::Isol_hor::n_auto_evol
Evolution_std< Scalar > n_auto_evol
Values at successive time steps of the lapse function .
Definition isol_hor.h:281
Lorene::Isol_hor::n_comp_evol
Evolution_std< Scalar > n_comp_evol
Values at successive time steps of the lapse function .
Definition isol_hor.h:284
Lorene::Isol_hor::set_boost_x
void set_boost_x(double bo)
Sets the boost velocity in x-direction to bo .
Definition isol_hor.h:449
Lorene::Isol_hor::boundary_nn_Dir_eff
const Valeur boundary_nn_Dir_eff(double aa) const
Dirichlet boundary condition for N (effectif) .
Definition bound_hor.C:586
Lorene::Isol_hor::dn_evol
Evolution_std< Vector > dn_evol
Values at successive time steps of the covariant derivative of the lapse with respect to the flat met...
Definition isol_hor.h:294
Lorene::Isol_hor::Isol_hor
Isol_hor(Map_af &mpi, int depth_in=3)
Standard constructor.
Definition isol_hor.C:176
Lorene::Isol_hor::regularisation
double regularisation(const Vector &shift_auto, const Vector &shift_comp, double ang_vel)
Corrects shift_auto in such a way that the total is equal to zero in the horizon,...
Definition regularisation.C:88
Lorene::Isol_hor::boundary_nn_Dir
const Valeur boundary_nn_Dir(double aa) const
Dirichlet boundary condition .
Definition bound_hor.C:647
Lorene::Isol_hor::boundary_nn_Neu_kk
const Valeur boundary_nn_Neu_kk(int nn=1) const
Neumann boundary condition for N using the extrinsic curvature.
Definition bound_hor.C:420
Lorene::Isol_hor::boost_x
double boost_x
Boost velocity in x-direction.
Definition isol_hor.h:272
Lorene::Isol_hor::boundary_psi_Dir
const Valeur boundary_psi_Dir() const
Dirichlet boundary condition for (spatial).
Definition bound_hor.C:324
Lorene::Isol_hor::operator=
void operator=(const Isol_hor &)
Assignment to another Isol_hor.
Definition isol_hor.C:343
Lorene::Isol_hor::beta_comp_evol
Evolution_std< Vector > beta_comp_evol
Values at successive time steps of the shift function .
Definition isol_hor.h:304
Lorene::Isol_hor::area_hor
double area_hor() const
Area of the horizon.
Definition phys_param.C:157
Lorene::Isol_hor::boundary_b_tilde_Dir
const Valeur boundary_b_tilde_Dir() const
Dirichlet boundary condition for b_tilde.
Definition bound_hor.C:1212
Lorene::Isol_hor::boundary_psi_Neu_evol
const Valeur boundary_psi_Neu_evol() const
Neumann boundary condition for (evolution).
Definition bound_hor.C:203
Lorene::Isol_hor::vv_bound_cart_bin
const Vector vv_bound_cart_bin(double om, int hole=0) const
Vector for boundary conditions in cartesian for binary systems.
Definition bound_hor.C:1379
Lorene::Isol_hor::source_nn
const Scalar source_nn() const
Source for N.
Definition sources_hor.C:145
Lorene::Isol_hor::beta_auto
virtual const Vector & beta_auto() const
Shift function at the current time step jtime.
Definition isol_hor.C:491
Lorene::Isol_hor::sauve
virtual void sauve(FILE *fich, bool partial_save) const
Total or partial saves in a binary file.
Definition isol_hor.C:401
Lorene::Isol_hor::set_mp
Map_af & set_mp()
Read/write of the mapping.
Definition isol_hor.h:421
Lorene::Isol_hor::omega_hor
double omega_hor() const
Orbital velocity.
Definition phys_param.C:233
Lorene::Isol_hor::get_radius
double get_radius() const
Returns the radius of the horizon.
Definition isol_hor.h:426
Lorene::Isol_hor::psi_comp_evol
Evolution_std< Scalar > psi_comp_evol
Values at successive time steps of the lapse function .
Definition isol_hor.h:290
Lorene::Isol_hor::source_b_tilde
const Scalar source_b_tilde() const
Source for b_tilde.
Lorene::Map_af
Affine radial mapping.
Definition map.h:2027
Lorene::Metric_flat
Flat metric for tensor calculation.
Definition metric.h:261
Lorene::Metric
Metric for tensor calculation.
Definition metric.h:90
Lorene::Scalar
Tensor field of valence 0 (or component of a tensorial field).
Definition scalar.h:387
Lorene::Single_hor
Binary black holes system.
Definition isol_hor.h:894
Lorene::Single_hor::radius_hor
double radius_hor() const
Radius of the horizon.
Definition single_param.C:97
Lorene::Single_hor::dpsi
Vector dpsi
Covariant derivative of the conformal factor .
Definition isol_hor.h:941
Lorene::Single_hor::psi_auto
Scalar psi_auto
Conformal factor .
Definition isol_hor.h:924
Lorene::Single_hor::radius
double radius
Radius of the horizon in LORENE's units.
Definition isol_hor.h:906
Lorene::Single_hor::get_gam
const Metric & get_gam() const
metric
Definition single_hor.C:339
Lorene::Single_hor::trK_point
Scalar trK_point
Time derivative of the trace of the extrinsic curvature.
Definition isol_hor.h:992
Lorene::Single_hor::darea_hor
const Scalar darea_hor() const
Element of area of the horizon.
Definition single_param.C:77
Lorene::Single_hor::boundary_psi_app_hor
const Valeur boundary_psi_app_hor() const
Neumann boundary condition for .
Definition single_bound.C:69
Lorene::Single_hor::boundary_beta_z
const Valeur boundary_beta_z() const
Component z of boundary value of .
Definition single_bound.C:260
Lorene::Single_hor::boundary_beta_y
const Valeur boundary_beta_y(double om_orb, double om_loc) const
Component y of boundary value of .
Definition single_bound.C:217
Lorene::Single_hor::decouple
Scalar decouple
Function used to construct from the total .
Definition isol_hor.h:1002
Lorene::Single_hor::beta_auto
Vector beta_auto
Shift function .
Definition isol_hor.h:944
Lorene::Single_hor::beta_comp_import
void beta_comp_import(const Single_hor &comp)
Imports the part of due to the companion hole comp.
Definition single_hor.C:465
Lorene::Single_hor::set_psi_auto
void set_psi_auto(const Scalar &psi_in)
Sets the conformal factor relating the physical metric to the conformal one: .
Definition single_hor.C:362
Lorene::Single_hor::mass_hor
double mass_hor() const
Mass computed at the horizon.
Definition single_param.C:135
Lorene::Single_hor::n_comp
Scalar n_comp
Lapse function .
Definition isol_hor.h:918
Lorene::Single_hor::get_k_dd
const Sym_tensor & get_k_dd() const
k_dd
Definition single_hor.C:348
Lorene::Single_hor::operator=
void operator=(const Single_hor &)
Assignment to another Single_hor.
Definition single_hor.C:186
Lorene::Single_hor::set_omega
void set_omega(double ome)
Sets the angular velocity to ome .
Definition isol_hor.h:1060
Lorene::Single_hor::aa_auto
Sym_tensor aa_auto
Components of the conformal representation of the traceless part of the extrinsic curvature:
Definition isol_hor.h:959
Lorene::Single_hor::get_psi_comp
const Scalar & get_psi_comp() const
Conformal factor .
Definition single_hor.C:280
Lorene::Single_hor::set_beta_auto
void set_beta_auto(const Scalar &shift_in)
Sets the shift.
Definition single_hor.C:370
Lorene::Single_hor::get_psi4
const Scalar & get_psi4() const
Conformal factor .
Definition single_hor.C:288
Lorene::Single_hor::set_n_auto
void set_n_auto(const Scalar &nn_in)
Sets the lapse.
Definition single_hor.C:366
Lorene::Single_hor::get_n_auto
const Scalar & get_n_auto() const
Lapse function .
Definition single_hor.C:261
Lorene::Single_hor::get_n_comp
const Scalar & get_n_comp() const
Lapse function .
Definition single_hor.C:266
Lorene::Single_hor::ff
Metric_flat ff
3 metric flat
Definition isol_hor.h:980
Lorene::Single_hor::psi
Scalar psi
Conformal factor .
Definition isol_hor.h:930
Lorene::Single_hor::mp
Map_af & mp
Affine mapping.
Definition isol_hor.h:900
Lorene::Single_hor::get_dpsi
const Vector & get_dpsi() const
Covariant derivative with respect to the flat metric of the conformal factor .
Definition single_hor.C:304
Lorene::Single_hor::set_radius
void set_radius(double rad)
Sets the radius of the horizon to rad .
Definition isol_hor.h:1051
Lorene::Single_hor::get_beta_comp
const Vector & get_beta_comp() const
Shift function .
Definition single_hor.C:314
Lorene::Single_hor::init_met_trK
void init_met_trK()
Sets the 3-metric tilde to the flat metric and gamt_point, trK and trK_point to zero.
Definition single_hor.C:536
Lorene::Single_hor::area_hor
double area_hor() const
Area of the horizon.
Definition single_param.C:88
Lorene::Single_hor::init_bhole
void init_bhole()
Sets the values of the fields to :
Definition single_hor.C:484
Lorene::Single_hor::expansion
Scalar expansion() const
Expansion of the outgoing null normal ( ).
Definition single_param.C:190
Lorene::Single_hor::ang_mom_adm
double ang_mom_adm() const
ADM angular Momentum.
Definition single_param.C:174
Lorene::Single_hor::get_nn
const Scalar & get_nn() const
Lapse function .
Definition single_hor.C:270
Lorene::Single_hor::p_psi4
Scalar * p_psi4
Conformal factor .
Definition isol_hor.h:933
Lorene::Single_hor::dn
Vector dn
Covariant derivative of the lapse with respect to the flat metric .
Definition isol_hor.h:937
Lorene::Single_hor::viriel_seul
double viriel_seul() const
Computes the viriel error, that is the difference between the ADM and the Komar masses,...
Definition binhor_viriel.C:57
Lorene::Single_hor::Single_hor
Single_hor(Map_af &mpi)
Standard constructor.
Definition single_hor.C:70
Lorene::Single_hor::set_der_0x0
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
Definition single_hor.C:231
Lorene::Single_hor::get_beta_auto
const Vector & get_beta_auto() const
Shift function .
Definition single_hor.C:309
Lorene::Single_hor::boundary_nn_Neu
const Valeur boundary_nn_Neu(double aa) const
Neumann boundary condition on nn .
Definition single_bound.C:106
Lorene::Single_hor::get_aa
const Sym_tensor & get_aa() const
Conformal representation of the traceless part of the extrinsic curvature:
Definition single_hor.C:334
Lorene::Single_hor::get_beta
const Vector & get_beta() const
Shift function .
Definition single_hor.C:319
Lorene::Single_hor::set_aa_comp
void set_aa_comp(const Scalar &aa_comp_in)
Sets aa_comp.
Definition single_hor.C:378
Lorene::Single_hor::nz
int nz
Number of zones.
Definition isol_hor.h:903
Lorene::Single_hor::regularisation
double regularisation(const Vector &shift_auto, const Vector &shift_comp, double ang_vel)
Corrects shift_auto in such a way that the total is equal to zero in the horizon,...
Definition single_regul.C:59
Lorene::Single_hor::init_bhole_seul
void init_bhole_seul()
Initiates for a single black hole.
Definition single_hor.C:548
Lorene::Single_hor::kappa_hor
double kappa_hor() const
Surface gravity.
Definition single_param.C:146
Lorene::Single_hor::aa_comp
Sym_tensor aa_comp
Components of the conformal representation of the traceless part of the extrinsic curvature:
Definition isol_hor.h:965
Lorene::Single_hor::get_mp
const Map_af & get_mp() const
Returns the mapping (readonly).
Definition isol_hor.h:1038
Lorene::Single_hor::del_deriv
void del_deriv() const
Deletes all the derived quantities.
Definition single_hor.C:221
Lorene::Single_hor::set_aa
void set_aa(const Scalar &aa_in)
Sets aa.
Definition single_hor.C:382
Lorene::Single_hor::regul
double regul
Intensity of the correction on the shift vector.
Definition isol_hor.h:912
Lorene::Single_hor::beta_comp
Vector beta_comp
Shift function .
Definition isol_hor.h:947
Lorene::Single_hor::ang_mom_hor
double ang_mom_hor() const
Angular momentum (modulo).
Definition single_param.C:107
Lorene::Single_hor::boundary_beta_x
const Valeur boundary_beta_x(double om_orb, double om_loc) const
Component x of boundary value of .
Definition single_bound.C:174
Lorene::Single_hor::gamt_point
Sym_tensor gamt_point
Time derivative of the 3-metric tilde.
Definition isol_hor.h:986
Lorene::Single_hor::set_mp
Map_af & set_mp()
Read/write of the mapping.
Definition isol_hor.h:1041
Lorene::Single_hor::omega_hor
double omega_hor() const
Orbital velocity.
Definition single_param.C:160
Lorene::Single_hor::get_aa_comp
const Sym_tensor & get_aa_comp() const
Conformal representation of the traceless part of the extrinsic curvature:
Definition single_hor.C:329
Lorene::Single_hor::get_dn
const Vector & get_dn() const
Covariant derivative of the lapse function .
Definition single_hor.C:299
Lorene::Single_hor::~Single_hor
virtual ~Single_hor()
Destructor.
Definition single_hor.C:176
Lorene::Single_hor::omega
double omega
Angular velocity in LORENE's units.
Definition isol_hor.h:909
Lorene::Single_hor::sauve
virtual void sauve(FILE *fich) const
Total or partial saves in a binary file.
Definition single_hor.C:244
Lorene::Single_hor::set_aa_auto
void set_aa_auto(const Scalar &aa_auto_in)
Sets aa_auto.
Definition single_hor.C:374
Lorene::Single_hor::get_radius
double get_radius() const
Returns the radius of the horizon.
Definition isol_hor.h:1046
Lorene::Single_hor::p_gam
Metric * p_gam
Spatial metric .
Definition isol_hor.h:953
Lorene::Single_hor::get_decouple
const Scalar get_decouple() const
Returns the function used to construct tkij_auto from tkij_tot .
Definition isol_hor.h:1146
Lorene::Single_hor::n_comp_import
void n_comp_import(const Single_hor &comp)
Imports the part of N due to the companion hole comp .
Definition single_hor.C:390
Lorene::Single_hor::n_auto
Scalar n_auto
Lapse function .
Definition isol_hor.h:915
Lorene::Single_hor::regularise_one
double regularise_one()
Corrects the shift in the innermost shell, so that it remains and that equals zero on the horizon.
Definition single_regul.C:153
Lorene::Single_hor::boundary_nn_Dir
const Valeur boundary_nn_Dir(double aa) const
Dirichlet boundary condition for N .
Definition single_bound.C:135
Lorene::Single_hor::nn
Scalar nn
Lapse function .
Definition isol_hor.h:921
Lorene::Single_hor::get_aa_auto
const Sym_tensor & get_aa_auto() const
Conformal representation of the traceless part of the extrinsic curvature:
Definition single_hor.C:324
Lorene::Single_hor::get_psi
const Scalar & get_psi() const
Conformal factor .
Definition single_hor.C:284
Lorene::Single_hor::p_k_dd
Sym_tensor * p_k_dd
Components of the extrinsic curvature:
Definition isol_hor.h:974
Lorene::Single_hor::get_psi_auto
const Scalar & get_psi_auto() const
Conformal factor .
Definition single_hor.C:275
Lorene::Single_hor::hh
Sym_tensor hh
Deviation metric.
Definition isol_hor.h:983
Lorene::Single_hor::get_omega
double get_omega() const
Returns the angular velocity.
Definition isol_hor.h:1056
Lorene::Single_hor::aa
Sym_tensor aa
Components of the conformal representation of the traceless part of the extrinsic curvature:
Definition isol_hor.h:971
Lorene::Single_hor::beta
Vector beta
Shift function .
Definition isol_hor.h:950
Lorene::Single_hor::get_tgam
const Metric & get_tgam() const
Conformal metric .
Definition isol_hor.h:1133
Lorene::Single_hor::b_tilde
const Scalar b_tilde() const
Radial component of the shift with respect to the conformal metric.
Definition single_param.C:68
Lorene::Single_hor::trK
Scalar trK
Trace of the extrinsic curvature.
Definition isol_hor.h:989
Lorene::Single_hor::psi_comp
Scalar psi_comp
Conformal factor .
Definition isol_hor.h:927
Lorene::Single_hor::psi_comp_import
void psi_comp_import(const Single_hor &comp)
Imports the part of due to the companion hole comp .
Definition single_hor.C:424
Lorene::Single_hor::tgam
Metric tgam
3 metric tilde
Definition isol_hor.h:977
Lorene::Sym_tensor_trans
Transverse symmetric tensors of rank 2.
Definition sym_tensor.h:608
Lorene::Sym_tensor
Class intended to describe valence-2 symmetric tensors.
Definition sym_tensor.h:223
Lorene::Tbl
Basic array class.
Definition tbl.h:161
Lorene::Time_slice_conf
Spacelike time slice of a 3+1 spacetime with conformal decomposition.
Definition time_slice.h:498
Lorene::Time_slice_conf::nn
virtual const Scalar & nn() const
Lapse function N at the current time step (jtime ).
Definition time_slice_conf.C:591
Lorene::Time_slice_conf::aa
virtual Sym_tensor aa() const
Conformal representation of the traceless part of the extrinsic curvature: .
Definition time_slice_conf.C:765
Lorene::Time_slice_conf::Time_slice_conf
Time_slice_conf(const Scalar &lapse_in, const Vector &shift_in, const Metric_flat &ff_in, const Scalar &psi_in, const Sym_tensor &hh_in, const Sym_tensor &hata_in, const Scalar &trk_in, int depth_in=3)
Constructor from conformal decomposition.
Definition time_slice_conf.C:118
Lorene::Time_slice
Spacelike time slice of a 3+1 spacetime.
Definition time_slice.h:173
Lorene::Valeur
Values and coefficients of a (real-value) function.
Definition valeur.h:287
Lorene::Vector
Tensor field of valence 1.
Definition vector.h:188
Lorene
Lorene prototypes.
Definition app_hor.h:64
Lorene::Map
Map(const Mg3d &)
Constructor from a multi-domain 3D grid.
Definition map.C:139