single_hor.C
/*
* Methods of class single_hor
*
* (see file isol_hor.h for documentation).
*
*/
/*
* Copyright (c) 2004 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
*
*/
char single_hor_C[] = "$Header: /cvsroot/Lorene/C++/Source/Isol_hor/single_hor.C,v 1.3 2014/10/13 08:53:01 j_novak Exp $" ;
/*
* $Id: single_hor.C,v 1.3 2014/10/13 08:53:01 j_novak Exp $
* $Log: single_hor.C,v $
* Revision 1.3 2014/10/13 08:53:01 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.2 2014/10/06 15:13:11 j_novak
* Modified #include directives to use c++ syntax.
*
* Revision 1.1 2007/04/13 15:28:35 f_limousin
* Lots of improvements, generalisation to an arbitrary state of
* rotation, implementation of the spatial metric given by Samaya.
*
*
* $Header: /cvsroot/Lorene/C++/Source/Isol_hor/single_hor.C,v 1.3 2014/10/13 08:53:01 j_novak Exp $
*
*/
// C headers
#include <cstdlib>
#include <cassert>
// Lorene headers
#include "param.h"
#include "utilitaires.h"
#include "time_slice.h"
#include "isol_hor.h"
#include "tensor.h"
#include "metric.h"
#include "evolution.h"
//#include "graphique.h"
//--------------//
// Constructors //
//--------------//
// Standard constructor
// --------------------
namespace Lorene {
Single_hor::Single_hor(Map_af& mpi) :
mp(mpi), nz(mpi.get_mg()->get_nzone()), radius ((mpi.get_alpha())[0]),
omega(0),regul(0),
n_auto(mpi), n_comp(mpi), nn(mpi),
psi_auto(mpi), psi_comp(mpi), psi(mpi),
dn(mpi, COV, mpi.get_bvect_spher()),
dpsi(mpi, COV, mpi.get_bvect_spher()),
beta_auto(mpi, CON, mpi.get_bvect_spher()),
beta_comp(mpi, CON, mpi.get_bvect_spher()),
beta(mpi, CON, mpi.get_bvect_spher()),
aa_auto(mpi, CON, mpi.get_bvect_spher()),
aa_comp(mpi, CON, mpi.get_bvect_spher()),
aa(mpi, CON, mpi.get_bvect_spher()),
tgam(mpi.flat_met_spher()),
ff(mpi.flat_met_spher()),
hh(mpi, CON, mpi.get_bvect_spher()),
gamt_point(mpi, CON, mpi.get_bvect_spher()),
trK(mpi), trK_point(mpi), decouple(mpi){
hh.set_etat_zero() ;
set_der_0x0() ;
}
// Copy constructor
// ----------------
Single_hor::Single_hor(const Single_hor& singlehor_in)
: mp(singlehor_in.mp),
nz(singlehor_in.nz),
radius(singlehor_in.radius),
omega(singlehor_in.omega),
regul(singlehor_in.regul),
n_auto(singlehor_in.n_auto),
n_comp(singlehor_in.n_comp),
nn(singlehor_in.nn),
psi_auto(singlehor_in.psi_auto),
psi_comp(singlehor_in.psi_comp),
psi(singlehor_in.psi),
dn(singlehor_in.dn),
dpsi(singlehor_in.dpsi),
beta_auto(singlehor_in.beta_auto),
beta_comp(singlehor_in.beta_comp),
beta(singlehor_in.beta),
aa_auto(singlehor_in.aa_auto),
aa_comp(singlehor_in.aa_comp),
aa(singlehor_in.aa),
tgam(singlehor_in.tgam),
ff(singlehor_in.ff),
hh(singlehor_in.hh),
gamt_point(singlehor_in.gamt_point),
trK(singlehor_in.trK),
trK_point(singlehor_in.trK_point),
decouple(singlehor_in.decouple){
set_der_0x0() ;
}
// Constructor from a file
// -----------------------
Single_hor::Single_hor(Map_af& mpi, FILE* fich)
: mp(mpi), nz(mpi.get_mg()->get_nzone()), radius ((mpi.get_alpha())[0]),
omega(0),regul(0),
n_auto(mpi, *(mpi.get_mg()), fich), n_comp(mpi),
nn(mpi),
psi_auto(mpi, *(mpi.get_mg()), fich), psi_comp(mpi),
psi(mpi),
dn(mpi, COV, mpi.get_bvect_spher()),
dpsi(mpi, COV, mpi.get_bvect_spher()),
beta_auto(mpi, mpi.get_bvect_spher(), fich),
beta_comp(mpi, CON, mpi.get_bvect_spher()),
beta(mpi, CON, mpi.get_bvect_spher()),
aa_auto(mpi, CON, mpi.get_bvect_spher()),
aa_comp(mpi, CON, mpi.get_bvect_spher()),
aa(mpi, CON, mpi.get_bvect_spher()),
tgam(mpi.flat_met_spher()),
ff(mpi.flat_met_spher()),
hh(mpi, CON, mpi.get_bvect_spher()),
gamt_point(mpi, CON, mpi.get_bvect_spher()),
trK(mpi), trK_point(mpi), decouple(mpi){
fread_be(&omega, sizeof(double), 1, fich) ;
// tgam, gamt_point, trK, trK_point
Sym_tensor met_file (mp, mp.get_bvect_spher(), fich) ;
tgam = met_file ;
Sym_tensor gamt_point_file (mp, mp.get_bvect_spher(), fich) ;
gamt_point = gamt_point_file ;
Scalar trK_file (mp, *(mp.get_mg()), fich) ;
trK = trK_file ;
Scalar trK_point_file (mp, *(mp.get_mg()), fich) ;
trK_point = trK_point_file ;
set_der_0x0() ;
hh = tgam.con() - ff.con() ;
}
//--------------//
// Destructor //
//--------------//
Single_hor::~Single_hor(){
Single_hor::del_deriv() ;
}
//-----------------------//
// Mutators / assignment //
//-----------------------//
void Single_hor::operator=(const Single_hor& singlehor_in) {
mp = singlehor_in.mp ;
nz = singlehor_in.nz ;
radius = singlehor_in.radius ;
omega = singlehor_in.omega ;
regul = singlehor_in.regul ;
n_auto = singlehor_in.n_auto ;
n_comp = singlehor_in.n_comp ;
nn = singlehor_in.nn ;
psi_auto = singlehor_in.psi_auto ;
psi_comp = singlehor_in.psi_comp ;
psi = singlehor_in.psi ;
dn = singlehor_in.dn ;
dpsi = singlehor_in.dpsi ;
beta_auto = singlehor_in.beta_auto ;
beta_comp = singlehor_in.beta_comp ;
beta = singlehor_in.beta ;
aa_auto = singlehor_in.aa_auto ;
aa_comp = singlehor_in.aa_comp ;
aa = singlehor_in.aa ;
tgam = singlehor_in.tgam ;
ff = singlehor_in.ff ;
hh = singlehor_in.hh ;
gamt_point = singlehor_in.gamt_point ;
trK = singlehor_in.trK ;
trK_point = singlehor_in.trK_point ;
decouple = singlehor_in.decouple ;
}
//---------------------//
// Memory management //
//---------------------//
void Single_hor::del_deriv() const {
if (p_psi4 != 0x0) delete p_psi4 ;
if (p_gam != 0x0) delete p_gam ;
if (p_k_dd != 0x0) delete p_k_dd ;
set_der_0x0() ;
}
void Single_hor::set_der_0x0() const {
p_psi4 = 0x0 ;
p_gam = 0x0 ;
p_k_dd = 0x0 ;
}
//--------------------------//
// Save in a file //
//--------------------------//
void Single_hor::sauve(FILE* fich) const {
n_auto.sauve(fich) ;
psi_auto.sauve(fich) ;
beta_auto.sauve(fich) ;
fwrite_be (&omega, sizeof(double), 1, fich) ;
tgam.con().sauve(fich) ;
gamt_point.sauve(fich) ;
trK.sauve(fich) ;
trK_point.sauve(fich) ;
}
// Accessors
// ---------
const Scalar& Single_hor::get_n_auto() const {
return n_auto ;
}
const Scalar& Single_hor::get_n_comp() const {
return n_comp ;
}
const Scalar& Single_hor::get_nn() const {
return nn ;
}
const Scalar& Single_hor::get_psi_auto() const {
return psi_auto ;
}
const Scalar& Single_hor::get_psi_comp() const {
return psi_comp ;
}
const Scalar& Single_hor::get_psi() const {
return psi ;
}
const Scalar& Single_hor::get_psi4() const {
if (p_psi4 == 0x0) {
p_psi4 = new Scalar( pow( psi, 4.) ) ;
p_psi4->std_spectral_base() ;
}
return *p_psi4 ;
}
const Vector& Single_hor::get_dn() const {
return dn ;
}
const Vector& Single_hor::get_dpsi() const {
return dpsi ;
}
const Vector& Single_hor::get_beta_auto() const {
return beta_auto ;
}
const Vector& Single_hor::get_beta_comp() const {
return beta_comp ;
}
const Vector& Single_hor::get_beta() const {
return beta ;
}
const Sym_tensor& Single_hor::get_aa_auto() const {
return aa_auto ;
}
const Sym_tensor& Single_hor::get_aa_comp() const {
return aa_comp ;
}
const Sym_tensor& Single_hor::get_aa() const {
return aa ;
}
const Metric& Single_hor::get_gam() const {
if (p_gam == 0x0) {
p_gam = new Metric( get_psi4()*tgam.cov() ) ;
}
return *p_gam ;
}
const Sym_tensor& Single_hor::get_k_dd() const {
if (p_k_dd == 0x0) {
Sym_tensor temp (aa/get_psi4()+1./3.*trK*get_gam().con()) ;
p_k_dd = new Sym_tensor( temp.up_down(get_gam()) ) ;
p_k_dd->std_spectral_base() ;
}
return *p_k_dd ;
}
void Single_hor::set_psi_auto(const Scalar& psi_in) {
psi_auto = psi_in ;
}
void Single_hor::set_n_auto(const Scalar& n_in) {
n_auto = n_in ;
}
void Single_hor::set_beta_auto(const Scalar& beta_in) {
beta_auto = beta_in ;
}
void Single_hor::set_aa_auto(const Scalar& aa_in) {
aa_auto = aa_in ;
}
void Single_hor::set_aa_comp(const Scalar& aa_in) {
aa_comp = aa_in ;
}
void Single_hor::set_aa(const Scalar& aa_in) {
aa = aa_in ;
}
// Import the lapse from the companion (Bhole case)
void Single_hor::n_comp_import(const Single_hor& comp) {
Scalar temp (mp) ;
temp.import(comp.n_auto) ;
temp.std_spectral_base() ;
n_comp = temp ;
nn = temp + n_auto ;
Vector dn_comp (mp, COV, mp.get_bvect_cart()) ;
dn_comp.set_etat_qcq() ;
Vector auxi (comp.n_auto.derive_cov(comp.ff)) ;
auxi.dec_dzpuis(2) ;
auxi.change_triad(auxi.get_mp().get_bvect_cart()) ;
for (int i=1 ; i<=3 ; i++){
if (auxi(i).get_etat() != ETATZERO)
auxi.set(i).raccord(3) ;
}
auxi.change_triad(mp.get_bvect_cart()) ;
assert ( *(auxi.get_triad()) == *(dn_comp.get_triad())) ;
for (int i=1 ; i<=3 ; i++){
dn_comp.set(i).import(auxi(i)) ;
dn_comp.set(i).set_spectral_va().set_base(auxi(i).get_spectral_va().
get_base()) ;
}
dn_comp.inc_dzpuis(2) ;
dn_comp.change_triad(mp.get_bvect_spher()) ;
dn = n_auto.derive_cov(ff) + dn_comp ;
}
// Import the conformal factor from the companion (Bhole case)
void Single_hor::psi_comp_import(const Single_hor& comp) {
Scalar temp (mp) ;
temp.import(comp.psi_auto) ;
temp.std_spectral_base() ;
psi_comp = temp ;
psi = temp + psi_auto ;
Vector dpsi_comp (mp, COV, mp.get_bvect_cart()) ;
dpsi_comp.set_etat_qcq() ;
Vector auxi (comp.psi_auto.derive_cov(comp.ff)) ;
auxi.dec_dzpuis(2) ;
auxi.change_triad(auxi.get_mp().get_bvect_cart()) ;
for (int i=1 ; i<=3 ; i++){
if (auxi(i).get_etat() != ETATZERO)
auxi.set(i).raccord(3) ;
}
auxi.change_triad(mp.get_bvect_cart()) ;
assert ( *(auxi.get_triad()) == *(dpsi_comp.get_triad())) ;
for (int i=1 ; i<=3 ; i++){
dpsi_comp.set(i).import(auxi(i)) ;
dpsi_comp.set(i).set_spectral_va().set_base(auxi(i).get_spectral_va().
get_base()) ;
}
dpsi_comp.inc_dzpuis(2) ;
dpsi_comp.change_triad(mp.get_bvect_spher()) ;
/*
Vector dpsi_comp_zec (psi_comp().derive_cov(ff)) ;
for (int i=1 ; i<=3 ; i++)
for (int l=nz-1 ; l<=nz-1 ; l++) {
if (dpsi_comp.set(i).get_etat() == ETATQCQ)
dpsi_comp.set(i).set_domain(l) = dpsi_comp_zec(i).domain(l) ;
}
*/
dpsi = psi_auto.derive_cov(ff) + dpsi_comp ;
}
void Single_hor::beta_comp_import(const Single_hor& comp) {
Vector tmp_vect (mp, CON, mp.get_bvect_cart()) ;
Vector shift_comp (comp.beta_auto) ;
shift_comp.change_triad(comp.mp.get_bvect_cart()) ;
shift_comp.change_triad(mp.get_bvect_cart()) ;
assert (*(shift_comp.get_triad()) == *(tmp_vect.get_triad())) ;
tmp_vect.set(1).import(shift_comp(1)) ;
tmp_vect.set(2).import(shift_comp(2)) ;
tmp_vect.set(3).import(shift_comp(3)) ;
tmp_vect.std_spectral_base() ;
tmp_vect.change_triad(mp.get_bvect_spher()) ;
beta_comp = tmp_vect ;
beta = beta_auto + beta_comp ;
}
//Initialisation to Schwartzchild
void Single_hor::init_bhole () {
Scalar auxi(mp) ;
// Initialisation of the lapse different of zero on the horizon
// at the first step
auxi = 0.5 - 0.5/mp.r ;
auxi.annule(0, 0);
auxi.set_dzpuis(0) ;
Scalar temp(mp) ;
temp = auxi;
temp.std_spectral_base() ;
temp.raccord(1) ;
n_auto = temp ;
temp = 0.5 ;
temp.std_spectral_base() ;
n_comp = temp ;
nn = n_auto ;
auxi = 0.5 + radius/mp.r ;
auxi.annule(0, 0);
auxi.set_dzpuis(0) ;
temp = auxi;
temp.std_spectral_base() ;
temp.raccord(1) ;
psi_auto = temp ;
temp = 0.5 ;
temp.std_spectral_base() ;
psi_comp = temp ;
psi = psi_auto + psi_comp ;
dn = nn.derive_cov(ff) ;
dpsi = psi.derive_cov(ff) ;
Vector temp_vect1(mp, CON, mp.get_bvect_spher()) ;
temp_vect1.set(1) = 0.0/mp.r/mp.r ;
temp_vect1.set(2) = 0. ;
temp_vect1.set(3) = 0. ;
temp_vect1.std_spectral_base() ;
Vector temp_vect2(mp, CON, mp.get_bvect_spher()) ;
temp_vect2.set_etat_zero() ;
beta_auto = temp_vect1 ;
beta_comp = temp_vect2 ;
beta = temp_vect1 ;
}
void Single_hor::init_met_trK() {
Metric flat (mp.flat_met_spher()) ;
tgam = flat ;
gamt_point.set_etat_zero() ;
trK.set_etat_zero() ;
trK_point.set_etat_zero() ;
}
void Single_hor::init_bhole_seul () {
Scalar auxi(mp) ;
auxi = (1-radius/mp.r)/(1+radius/mp.r) ;
auxi.annule(0, 0);
auxi.set_outer_boundary((*mp.get_mg()).get_nzone(), 1.) ;
auxi.set_dzpuis(0) ;
Scalar temp(mp) ;
temp = auxi;
temp.std_spectral_base() ;
temp.raccord(1) ;
n_auto = temp;
temp.set_etat_zero() ;
n_comp = temp ;
nn = temp ;
auxi = 1 + radius/mp.r ;
auxi.annule(0, 0);
auxi.set_outer_boundary((*mp.get_mg()).get_nzone(), 1.) ;
auxi.set_dzpuis(0) ;
temp = auxi;
temp.std_spectral_base() ;
temp.raccord(1) ;
psi_auto = temp ;
temp.set_etat_zero() ;
psi_comp = temp ;
psi = temp ;
dn = nn.derive_cov(ff) ;
dpsi = psi.derive_cov(ff) ;
Vector temp_vect(mp, CON, mp.get_bvect_spher()) ;
temp_vect.set_etat_zero() ;
beta_auto = temp_vect ;
beta_comp = temp_vect ;
beta = temp_vect ;
}
}
Lorene::Map_af
Affine radial mapping.
Definition map.h:2027
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::Scalar::raccord
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
Definition scalar_raccord.C:62
Lorene::Scalar::derive_cov
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this.
Definition scalar_deriv.C:390
Lorene::Scalar::std_spectral_base
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
Definition scalar.C:784
Lorene::Scalar::set_etat_zero
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Definition scalar.C:324
Lorene::Scalar::annule
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
Definition scalar.C:391
Lorene::Scalar::set_spectral_va
Valeur & set_spectral_va()
Returns va (read/write version).
Definition scalar.h:604
Lorene::Scalar::set_dzpuis
void set_dzpuis(int)
Modifies the dzpuis flag.
Definition scalar.C:808
Lorene::Scalar::set_outer_boundary
void set_outer_boundary(int l, double x)
Sets the value of the Scalar at the outer boundary of a given domain.
Definition scalar_manip.C:315
Lorene::Scalar::import
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
Definition scalar_import.C:68
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::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::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::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::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::init_bhole
void init_bhole()
Sets the values of the fields to :
Definition single_hor.C:484
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::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::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::init_bhole_seul
void init_bhole_seul()
Initiates for a single black hole.
Definition single_hor.C:548
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::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::gamt_point
Sym_tensor gamt_point
Time derivative of the 3-metric tilde.
Definition isol_hor.h:986
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::p_gam
Metric * p_gam
Spatial metric .
Definition isol_hor.h:953
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::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::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::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
Class intended to describe valence-2 symmetric tensors.
Definition sym_tensor.h:223
Lorene::Valeur::set_base
void set_base(const Base_val &)
Sets the bases for spectral expansions (member base ).
Definition valeur.C:810
Lorene::Vector
Tensor field of valence 1.
Definition vector.h:188
Lorene::Vector::std_spectral_base
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Definition vector.C:316
Lorene::Vector::change_triad
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Definition vector_change_triad.C:75
Lorene::Vector::set
Scalar & set(int)
Read/write access to a component.
Definition vector.C:296
Lorene::pow
Cmp pow(const Cmp &, int)
Power .
Definition cmp_math.C:348
Lorene::fread_be
int fread_be(int *aa, int size, int nb, FILE *fich)
Reads integer(s) from a binary file according to the big endian convention.
Definition fread_be.C:69
Lorene::fwrite_be
int fwrite_be(const int *aa, int size, int nb, FILE *fich)
Writes integer(s) into a binary file according to the big endian convention.
Definition fwrite_be.C:70
Lorene::Tensor::get_mp
const Map & get_mp() const
Returns the mapping.
Definition tensor.h:861
Lorene::Tensor::dec_dzpuis
virtual void dec_dzpuis(int dec=1)
Decreases by dec units the value of dzpuis and changes accordingly the values in the compactified ext...
Definition tensor.C:808
Lorene::Tensor::up_down
Tensor up_down(const Metric &gam) const
Computes a new tensor by raising or lowering all the indices of *this .
Definition tensor_calculus.C:305
Lorene::Tensor::inc_dzpuis
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values in the compactified ext...
Definition tensor.C:816
Lorene::Tensor::get_triad
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Definition tensor.h:866
Lorene::Tensor::set_etat_qcq
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Definition tensor.C:481
Lorene::Tensor::set_etat_zero
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Definition tensor.C:497
Lorene
Lorene prototypes.
Definition app_hor.h:64
Lorene::get_bvect_spher
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
Definition map.h:783
Lorene::get_mg
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Definition map.h:765
Lorene::flat_met_spher
const Metric_flat & flat_met_spher() const
Returns the flat metric associated with the spherical coordinates and with components expressed in th...
Definition map.C:321