Excision_surf Class Reference
[Stars and black holes]

Surface where boundary conditions for quantities in the bulk will be calculated It relies on geometrical properties of the associated Spheroid() (*** WARNING! under development***). More...

#include <excision_surf.h>

List of all members.

Public Member Functions

 Excision_surf (const Scalar &h_in, const Metric &gij, const Sym_tensor &Kij2, const Scalar &ppsi, const Scalar &nn, const Vector &beta, double timestep, int int_nos)
 Constructor of an excision surface embedded in a 3-slice (Time_slice ) of 3+1 formalism.
 Excision_surf (const Excision_surf &)
 Copy constructor.
 Excision_surf (FILE *)
 Constructor from a file (see sauve(FILE*) ).
virtual ~Excision_surf ()
 Destructor.
void operator= (const Excision_surf &)
 Assignment to another Excision_surf.
void get_evol_params_from_ID (double alpha, double beta, double gamma, Scalar &Ee, Vector &Jj, Sym_tensor &Ss)
 Computes the parameters for the hyperbolic evolution in set_expa_hyperb(), so that the expansion has a C1 matching with initial data.
void set_expa_parab (double c_theta_lap, double c_theta_fin, Scalar &expa_fin)
 Sets a new value for expansion rescaled over lapse (and its derivative), obtained by parabolic evolution.
void set_expa_hyperb (double alph0, double beta0, double gamma0)
 Sets a new value for expansion rescaled over lapse (and its derivative), obtained by hyperbolic evolution.
const Spheroid & get_sph () const
 Returns the spheroid.
const Scalar & get_conf_fact () const
 Returns the conformal factor associated with the surface.
const Scalar & get_lapse () const
 Returns the lapse function.
const Vector & get_shift () const
 Returns the shift vector field.
const Metric & get_gamij () const
 Returns the symmetric tensor $ gamij $.
const Sym_tensor & get_Kij () const
 returns the 3-d extrinsic curvature $ K_{ij}$
double get_delta_t () const
 Returns the timestep used for evolution.
double get_no_of_steps () const
 Returns the internal number of timesteps for one iteration.
const Scalar & get_expa () const
 Returns the assumed expansion associated to the excised surface at t.
const Scalar & get_dt_expa () const
 Returns the assumed time derivative of the expansion at t, evolved by functions of this class;.
Spheroid & set_sph ()
 Sets a new spheroid from data.
Scalar & set_conf_fact ()
 Sets the value of the conformal factor.
Scalar & set_lapse ()
 Sets the lapse function.
Vector & set_shift ()
 Sets the shift vector field.
Metric & set_gamij ()
 Sets the 3d metric of the TimeSlice.
Sym_tensor & set_Kij ()
 Sets the extrinsic curvature.
double set_delta_t ()
double set_no_of_steps ()
Scalar & set_expa ()
 Sets the expansion function on the surface at time t (considering to protect this function).
Scalar & set_dt_expa ()
 Sets the time derivative of the expansion function on the surface at time t (considering to protect this function).
const Scalar & get_BC_conf_fact_1 (bool isMOTS=false) const
 Source for a Neumann BC on the conformal factor. If boolean isMOTS is false, it is based on expansion value of the spheroid or the value of exppa; it is based on zero expansion if isMOTS is true.
const Scalar & get_BC_lapse_1 (double value) const
const Vector & get_BC_shift_1 (double Omega) const
const Scalar & get_BC_Npsi_1 (double value) const
const Scalar & get_BC_conf_fact_2 (double c_psi_lap, double c_psi_fin, Scalar &expa_fin) const
 Source for the Dirichlet BC on the conformal factor, based on a parabolic driver for the conformal factor.
const Scalar & get_BC_conf_fact_3 (double c_theta_lap, double c_theta_fin, Scalar &expa_fin) const
 Source for the Neumann BC on the conformal factor, based on a parabolic driver for the expansio.
const Scalar & get_BC_conf_fact_4 () const
 Source for the Dirchlet BC on the conformal factor, based on the consistency condition derived from the trace.
const Scalar & get_BC_lapse_2 (double lapse_fin, double c_lapse_lap, double c_lapse_fi) const
 Source for Dirichlet BC on the lapse, based on a parabolic driver towards arbitrary constant value.
const Scalar & get_BC_lapse_3 (Scalar &dttheta, Scalar &Ee, Vector &Jj, Sym_tensor &Sij, bool sph_sym=true) const
 Source for Dirichlet BC on the lapse, based on einstein equations.
const Scalar & get_BC_lapse_4 (Scalar &old_nn, Vector &beta_point, Sym_tensor &strain_tens) const
 Source for Dirichlet BC on the lapse, based on einstein equations (conservation of isotropic gauge).
const Scalar & derive_t_expa (Scalar &Ee, Vector &Jj, Sym_tensor &Sij) const
 Forms the prospective time derivative for the expansion using projected Einstein equations. Does NOT modify the member dt_expa: do it by hand!
const Vector & get_BC_shift_2 (double c_bb_lap, double c_bb_fin, double c_V_lap, double epsilon) const
 Source for a Dirichlet BC on the shift, based on a Parabolic driver; no assumptions are made except a global conformal Killing symmetry.
const Vector & get_BC_shift_3 (Scalar &dtpsi, double c_V_lap, double epsilon) const
 Source for a Dirichlet BC on the shift, based on a Parabolic driver; Radial part is dealt with using a kinematical relation.
const Vector & get_BC_shift_4 (Scalar &dttheta, Scalar &Ee, Vector &Jj, Sym_tensor &Sij, double c_V_lap, double epsilon, bool sph_sym=true) const
 Source for a Dirichlet BC on the shift, based on a Parabolic driver; Radial part is dealt with using projection of Einstein Equations.
const Scalar & get_BC_Npsi_2 (double value, double c_npsi_lap, double c_npsi_fin) const
 Source for the Dirichlet BC on (N*Psi1), based on a parabolic driver.
const Scalar & get_BC_Npsi_3 (double n_0, double beta) const
 Source for the Dirichlet BC on (N*Psi1), with Kerr_Schild-like form for the lapse boundary.
const Scalar & get_BC_Npsi_4 (double Kappa) const
 Source for a Dirichlet BC on (N*Psi1), fixing a constant surface gravity in space and time.
const Scalar & get_BC_Npsi_5 (double Kappa) const
 Source for a Neumann BC on (N*Psi1), fixing a constant surface gravity in space and time.
virtual void sauve (FILE *) const
 Save in a file.

Protected Member Functions

virtual void del_deriv () const
 Deletes all the derived quantities.
void set_der_0x0 () const
 Sets to 0x0 all the pointers on derived quantities.

Protected Attributes

Spheroid sph
 The associated Spheroid object.
Scalar conf_fact
 The value of the conformal factor on the 3-slice.
Scalar lapse
 The lapse defined on the 3 slice.
Vector shift
 The Shift 3-vector on the slice.
Metric gamij
 The 3-d metric on the slice.
Sym_tensor Kij
 The 3-d extrinsic curvature on the slice.
double delta_t
 The time step for evolution in parabolic drivers.
double no_of_steps
 The internal number of timesteps for one iteration.
Scalar expa
 The 2d expansion, directly evolved from the initial excision with Einstein Equations.
Scalar dt_expa
 The time derivative of the expansion, derived from Einstein equations and arbitrary evolution.
Scalar * p_get_BC_conf_fact_1
 Source of Neumann boundary condition on $ \psi $,.
Scalar * p_get_BC_lapse_1
 Source of Dirichlet boundary condition of $ N $.
Vector * p_get_BC_shift_1
 Source of Dirichlet BC for the shift vector $ \beta^{i} $.
Scalar * p_get_BC_Npsi_1
 Source of Neumann boundary condition on $ \psi $.
Scalar * p_get_BC_conf_fact_2
 Source of Neumann boundary condition on $ \psi $,.
Scalar * p_get_BC_conf_fact_3
 Source of Neumann boundary condition on $ \psi $,.
Scalar * p_get_BC_conf_fact_4
 Source of Birichlet boundary condition on $ \psi $,.
Scalar * p_get_BC_lapse_2
 Source of Dirichlet boundary condition of $ N $.
Scalar * p_get_BC_lapse_3
 Source of Dirichlet condtion on $ N $, based on einstein equations.
Scalar * p_get_BC_lapse_4
 Source of Dirichlet condtion on $ N $, based on einstein equations (conservation of isotropic gauge).
Scalar * p_derive_t_expa
 Computation of an updated expansion scalar.
Vector * p_get_BC_shift_2
 Source of Dirichlet BC for the shift vector $ \beta^{i} $.
Vector * p_get_BC_shift_3
 Source of Dirichlet BC for the shift vector $ \beta^{i} $, partly derived from kinematical relation.
Vector * p_get_BC_shift_4
 Source of Dirichlet BC for the shift vector $ \beta^{i} $, partly from projection of Einstein Equations.
Scalar * p_get_BC_Npsi_2
 Source of Dirichlet boundary condition on $ N \psi $.
Scalar * p_get_BC_Npsi_3
 Source of Dirichlet boundary condition on $ N \psi $.
Scalar * p_get_BC_Npsi_4
 Source of Dirichlet boundary condition on $ N \psi $.
Scalar * p_get_BC_Npsi_5
 Source of Neumann boundary condition on $ N \psi $.

Friends

ostream & operator<< (ostream &, const Spheroid &)
 Display.

Detailed Description

Surface where boundary conditions for quantities in the bulk will be calculated It relies on geometrical properties of the associated Spheroid() (*** WARNING! under development***).

Definition at line 42 of file excision_surf.h.

Constructor & Destructor Documentation

Excision_surf::Excision_surf ( const Scalar &  h_in,
const Metric &  gij,
const Sym_tensor &  Kij2,
const Scalar &  ppsi,
const Scalar &  nn,
const Vector &  beta,
double  timestep,
int  int_nos = 1 
)

Constructor of an excision surface embedded in a 3-slice (Time_slice ) of 3+1 formalism.

This is done from the Time_slice data.

Parameters:
h_in : the location of the surface r = h_in (WARNING:must be defined on a mono-domain angular grid)
gij : the 3-metric on the 3slice
Kij : the extrinsic curvature of the 3-slice (covariant representation)
timestep : time interval associated with the parabolic-driven boundary conditions.
int_nos : Number of iterations to be done during timestep.

Definition at line 63 of file excision_surf.C.

References dt_expa, set_der_0x0(), and Scalar::set_etat_zero().

Excision_surf::Excision_surf ( const Excision_surf &  exc_in  ) 

Copy constructor.

Definition at line 89 of file excision_surf.C.

References set_der_0x0().

Excision_surf::Excision_surf ( FILE *   ) 

Constructor from a file (see sauve(FILE*) ).

Excision_surf::~Excision_surf (  )  [virtual]

Destructor.

Definition at line 130 of file excision_surf.C.

References del_deriv().

Member Function Documentation

void Excision_surf::del_deriv (  )  const [protected, virtual]

Deletes all the derived quantities.

Definition at line 138 of file excision_surf.C.

References p_derive_t_expa, p_get_BC_conf_fact_1, p_get_BC_conf_fact_2, p_get_BC_conf_fact_3, p_get_BC_conf_fact_4, p_get_BC_lapse_1, p_get_BC_lapse_2, p_get_BC_lapse_3, p_get_BC_lapse_4, p_get_BC_Npsi_1, p_get_BC_Npsi_2, p_get_BC_Npsi_3, p_get_BC_Npsi_4, p_get_BC_Npsi_5, p_get_BC_shift_1, p_get_BC_shift_2, p_get_BC_shift_3, p_get_BC_shift_4, and set_der_0x0().

const Scalar & Excision_surf::derive_t_expa ( Scalar &  Ee,
Vector &  Jj,
Sym_tensor &  Sij 
) const

Forms the prospective time derivative for the expansion using projected Einstein equations. Does NOT modify the member dt_expa: do it by hand!

Definition at line 782 of file excision_surf.C.

References Scalar::annule_hard(), contract(), delta_t, Scalar::derive_cov(), Spheroid::derive_cov2d(), gamij, Spheroid::get_ll(), Map::get_mg(), Tensor::get_mp(), Spheroid::get_qab(), Spheroid::get_ricci(), Kij, lapse, p_derive_t_expa, Metric::radial_vect(), Scalar::set_grid_point(), Scalar::set_spectral_va(), Spheroid::shear(), shift, sph, Scalar::std_spectral_base(), Spheroid::theta_minus(), Spheroid::theta_plus(), Tensor::up_down(), Scalar::val_grid_point(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_conf_fact_1 ( bool  isMOTS = false  )  const

Source for a Neumann BC on the conformal factor. If boolean isMOTS is false, it is based on expansion value of the spheroid or the value of exppa; it is based on zero expansion if isMOTS is true.

Definition at line 279 of file excision_surf.C.

References Scalar::allocate_all(), Tensor::annule_domain(), conf_fact, contract(), Metric::cov(), Scalar::derive_cov(), Vector::divergence(), Scalar::dsdr(), expa, gamij, Valeur::get_etat(), Map::get_mg(), Tensor::get_mp(), Scalar::get_spectral_va(), Kij, lapse, p_get_BC_conf_fact_1, pow(), Metric::radial_vect(), Scalar::set_etat_zero(), Scalar::set_spectral_va(), sph, Scalar::std_spectral_base(), Tensor::std_spectral_base(), Spheroid::theta_plus(), Tensor::trace(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_conf_fact_2 ( double  c_psi_lap,
double  c_psi_fin,
Scalar &  expa_fin 
) const

Source for the Dirichlet BC on the conformal factor, based on a parabolic driver for the conformal factor.

Definition at line 363 of file excision_surf.C.

References Scalar::allocate_all(), Tensor::annule_domain(), conf_fact, contract(), Metric::cov(), delta_t, Scalar::derive_cov(), Vector::divergence(), gamij, Map::get_mg(), Tensor::get_mp(), Kij, lapse, p_get_BC_conf_fact_2, pow(), Metric::radial_vect(), Scalar::set_spectral_va(), sph, Tensor::std_spectral_base(), Scalar::std_spectral_base(), Spheroid::theta_plus(), Tensor::trace(), Scalar::val_grid_point(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_conf_fact_3 ( double  c_theta_lap,
double  c_theta_fin,
Scalar &  expa_fin 
) const

Source for the Neumann BC on the conformal factor, based on a parabolic driver for the expansio.

Definition at line 451 of file excision_surf.C.

References Scalar::allocate_all(), conf_fact, delta_t, get_BC_conf_fact_1(), Map::get_mg(), Tensor::get_mp(), Scalar::lapang(), lapse, p_get_BC_conf_fact_3, pow(), Scalar::set_spectral_va(), sph, Scalar::std_spectral_base(), Spheroid::theta_plus(), Scalar::val_grid_point(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_conf_fact_4 (  )  const

Source for the Dirchlet BC on the conformal factor, based on the consistency condition derived from the trace.

Definition at line 532 of file excision_surf.C.

References Tensor::annule_domain(), conf_fact, contract(), delta_t, Scalar::derive_cov(), Vector::divergence(), Map::flat_met_spher(), gamij, Map::get_mg(), Tensor::get_mp(), Kij, lapse, p_get_BC_conf_fact_4, Scalar::set_spectral_va(), shift, Scalar::std_spectral_base(), Tensor::trace(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_lapse_2 ( double  lapse_fin,
double  c_lapse_lap,
double  c_lapse_fi 
) const

Source for Dirichlet BC on the lapse, based on a parabolic driver towards arbitrary constant value.

Definition at line 602 of file excision_surf.C.

References delta_t, Scalar::lapang(), lapse, p_get_BC_lapse_2, Scalar::set_spectral_va(), Scalar::std_spectral_base(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_lapse_3 ( Scalar &  dttheta,
Scalar &  Ee,
Vector &  Jj,
Sym_tensor &  Sij,
bool  sph_sym = true 
) const

Source for Dirichlet BC on the lapse, based on einstein equations.

Definition at line 646 of file excision_surf.C.

References Scalar::annule_hard(), Valeur::coef_i(), Metric_flat::con(), Metric::con(), contract(), Spheroid::delta(), Scalar::derive_cov(), Spheroid::derive_cov2d(), Spheroid::derive_cov2dflat(), Map::flat_met_spher(), gamij, Spheroid::get_hsurf(), Spheroid::get_ll(), Map::get_mg(), Tensor::get_mp(), Spheroid::get_qab(), Spheroid::get_ricci(), Kij, lapse, p_get_BC_lapse_3, Scalar::poisson_angu(), Metric::radial_vect(), Tensor::set(), Scalar::set_grid_point(), Scalar::set_spectral_va(), Spheroid::shear(), shift, sph, Spheroid::sqrt_q(), Tensor::std_spectral_base(), Scalar::std_spectral_base(), Spheroid::theta_minus(), Spheroid::theta_plus(), Tensor::up_down(), Scalar::val_grid_point(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_lapse_4 ( Scalar &  old_nn,
Vector &  beta_point,
Sym_tensor &  strain_tens 
) const

Source for Dirichlet BC on the lapse, based on einstein equations (conservation of isotropic gauge).

Definition at line 14 of file excision_lapse_4.C.

References Tensor::annule_domain(), Scalar::annule_hard(), Metric_flat::con(), Metric::con(), conf_fact, contract(), Metric_flat::cov(), delta_t, Scalar::derive_con(), Tensor::derive_cov(), Scalar::derive_cov(), Sym_tensor::derive_lie(), Scalar::derive_lie(), Vector::divergence(), gamij, Map::get_bvect_spher(), Map::get_mg(), Tensor::get_mp(), Scalar::inc_dzpuis(), Kij, lapse, log(), Vector::ope_killing_conf(), p_get_BC_lapse_4, Tensor::set(), Scalar::set_spectral_va(), shift, Scalar::std_spectral_base(), Tensor::std_spectral_base(), Tensor::trace(), Tensor::up_down(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_Npsi_2 ( double  value,
double  c_npsi_lap,
double  c_npsi_fin 
) const

Source for the Dirichlet BC on (N*Psi1), based on a parabolic driver.

Definition at line 1716 of file excision_surf.C.

References delta_t, Scalar::lapang(), lapse, p_get_BC_Npsi_2, Scalar::set_spectral_va(), Scalar::std_spectral_base(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_Npsi_3 ( double  n_0,
double  beta 
) const

Source for the Dirichlet BC on (N*Psi1), with Kerr_Schild-like form for the lapse boundary.

Definition at line 1756 of file excision_surf.C.

References Map::cost, Tensor::get_mp(), lapse, p_get_BC_Npsi_3, Scalar::set_spectral_va(), sqrt(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_Npsi_4 ( double  Kappa  )  const

Source for a Dirichlet BC on (N*Psi1), fixing a constant surface gravity in space and time.

Definition at line 1803 of file excision_surf.C.

References contract(), Scalar::derive_cov(), gamij, Kij, lapse, p_get_BC_Npsi_4, Metric::radial_vect(), Scalar::set_dzpuis(), Scalar::set_spectral_va(), Scalar::std_spectral_base(), and Valeur::ylm().

const Scalar & Excision_surf::get_BC_Npsi_5 ( double  Kappa  )  const

Source for a Neumann BC on (N*Psi1), fixing a constant surface gravity in space and time.

Definition at line 1824 of file excision_surf.C.

References Tensor::annule_domain(), contract(), Tensor::derive_cov(), gamij, Map::get_mg(), Tensor::get_mp(), Kij, lapse, p_get_BC_Npsi_5, pow(), Metric::radial_vect(), Scalar::set_spectral_va(), Scalar::std_spectral_base(), and Valeur::ylm().

const Vector & Excision_surf::get_BC_shift_2 ( double  c_bb_lap,
double  c_bb_fin,
double  c_V_lap,
double  epsilon 
) const

Source for a Dirichlet BC on the shift, based on a Parabolic driver; no assumptions are made except a global conformal Killing symmetry.

Definition at line 887 of file excision_surf.C.

References Scalar::allocate_all(), Tensor::annule_domain(), Metric::con(), contract(), Metric::cov(), delta_t, Tensor::derive_cov(), Tensor::down(), gamij, Spheroid::get_ricci(), Scalar::lapang(), lapse, p_get_BC_shift_2, Metric::radial_vect(), Scalar::set_spectral_va(), shift, sph, Vector::std_spectral_base(), Tensor::std_spectral_base(), Scalar::std_spectral_base(), Tensor::up_down(), Scalar::val_grid_point(), and Valeur::ylm().

const Vector & Excision_surf::get_BC_shift_3 ( Scalar &  dtpsi,
double  c_V_lap,
double  epsilon 
) const

Source for a Dirichlet BC on the shift, based on a Parabolic driver; Radial part is dealt with using a kinematical relation.

Definition at line 1126 of file excision_surf.C.

References Scalar::allocate_all(), Tensor::annule_domain(), Metric::con(), conf_fact, contract(), Metric::cov(), delta_t, Tensor::derive_cov(), Scalar::derive_cov(), Vector::divergence(), Tensor::down(), Map::flat_met_spher(), gamij, Map::get_mg(), Tensor::get_mp(), Spheroid::get_ricci(), lapse, p_get_BC_shift_3, Metric::radial_vect(), shift, sph, Vector::std_spectral_base(), Tensor::std_spectral_base(), Tensor::up_down(), and Scalar::val_grid_point().

const Vector & Excision_surf::get_BC_shift_4 ( Scalar &  dttheta,
Scalar &  Ee,
Vector &  Jj,
Sym_tensor &  Sij,
double  c_V_lap,
double  epsilon,
bool  sph_sym = true 
) const

Source for a Dirichlet BC on the shift, based on a Parabolic driver; Radial part is dealt with using projection of Einstein Equations.

Definition at line 1360 of file excision_surf.C.

References Scalar::allocate_all(), Tensor::annule_domain(), Scalar::annule_hard(), Valeur::coef_i(), Metric_flat::con(), Metric::con(), contract(), Metric::cov(), Spheroid::delta(), delta_t, Tensor::derive_cov(), Scalar::derive_cov(), Spheroid::derive_cov2d(), Spheroid::derive_cov2dflat(), Tensor::down(), Map::flat_met_spher(), gamij, Spheroid::get_hsurf(), Spheroid::get_ll(), Map::get_mg(), Tensor::get_mp(), Spheroid::get_qab(), Spheroid::get_ricci(), Scalar::get_spectral_base(), Kij, lapse, p_get_BC_shift_4, Metric::radial_vect(), Tensor::set(), Scalar::set_grid_point(), Scalar::set_spectral_va(), Spheroid::shear(), shift, sph, Spheroid::sqrt_q(), Vector::std_spectral_base(), Tensor::std_spectral_base(), Scalar::std_spectral_base(), Spheroid::theta_minus(), Spheroid::theta_plus(), Tensor::up_down(), Scalar::val_grid_point(), and Valeur::ylm().

const Scalar& Excision_surf::get_conf_fact (  )  const [inline]

Returns the conformal factor associated with the surface.

Definition at line 169 of file excision_surf.h.

References conf_fact.

double Excision_surf::get_delta_t (  )  const [inline]

Returns the timestep used for evolution.

Definition at line 184 of file excision_surf.h.

References delta_t.

const Scalar& Excision_surf::get_dt_expa (  )  const [inline]

Returns the assumed time derivative of the expansion at t, evolved by functions of this class;.

Definition at line 193 of file excision_surf.h.

References dt_expa.

void Excision_surf::get_evol_params_from_ID ( double  alpha,
double  beta,
double  gamma,
Scalar &  Ee,
Vector &  Jj,
Sym_tensor &  Ss 
)

Computes the parameters for the hyperbolic evolution in set_expa_hyperb(), so that the expansion has a C1 matching with initial data.

Sets also values for expa() and dt_expa() accordingly with initial conditions

Definition at line 194 of file excision_surf.C.

References Valeur::c_cf, derive_t_expa(), max(), set_dt_expa(), set_expa(), Scalar::set_spectral_va(), sph, Scalar::std_spectral_base(), Spheroid::theta_plus(), Mtbl_cf::val_in_bound_jk(), and Valeur::ylm().

const Scalar& Excision_surf::get_expa (  )  const [inline]

Returns the assumed expansion associated to the excised surface at t.

Definition at line 190 of file excision_surf.h.

References expa.

const Metric& Excision_surf::get_gamij (  )  const [inline]

Returns the symmetric tensor $ gamij $.

Definition at line 178 of file excision_surf.h.

References gamij.

const Sym_tensor& Excision_surf::get_Kij (  )  const [inline]

returns the 3-d extrinsic curvature $ K_{ij}$

Definition at line 181 of file excision_surf.h.

References Kij.

const Scalar& Excision_surf::get_lapse (  )  const [inline]

Returns the lapse function.

Definition at line 172 of file excision_surf.h.

References lapse.

double Excision_surf::get_no_of_steps (  )  const [inline]

Returns the internal number of timesteps for one iteration.

Definition at line 187 of file excision_surf.h.

References no_of_steps.

const Vector& Excision_surf::get_shift (  )  const [inline]

Returns the shift vector field.

Definition at line 175 of file excision_surf.h.

References shift.

const Spheroid& Excision_surf::get_sph (  )  const [inline]

Returns the spheroid.

Definition at line 166 of file excision_surf.h.

References sph.

void Excision_surf::operator= ( const Excision_surf &  surf_in  ) 

Assignment to another Excision_surf.

Definition at line 109 of file excision_surf.C.

References conf_fact, del_deriv(), delta_t, dt_expa, expa, gamij, Kij, lapse, no_of_steps, shift, and sph.

void Excision_surf::sauve ( FILE *   )  const [virtual]

Save in a file.

Definition at line 1850 of file excision_surf.C.

Scalar& Excision_surf::set_conf_fact (  )  [inline]

Sets the value of the conformal factor.

Definition at line 199 of file excision_surf.h.

References conf_fact, and del_deriv().

void Excision_surf::set_der_0x0 (  )  const [protected]

Sets to 0x0 all the pointers on derived quantities.

Definition at line 160 of file excision_surf.C.

References p_derive_t_expa, p_get_BC_conf_fact_1, p_get_BC_conf_fact_2, p_get_BC_conf_fact_3, p_get_BC_conf_fact_4, p_get_BC_lapse_1, p_get_BC_lapse_2, p_get_BC_lapse_3, p_get_BC_lapse_4, p_get_BC_Npsi_1, p_get_BC_Npsi_2, p_get_BC_Npsi_3, p_get_BC_Npsi_4, p_get_BC_Npsi_5, p_get_BC_shift_1, p_get_BC_shift_2, p_get_BC_shift_3, and p_get_BC_shift_4.

Scalar& Excision_surf::set_dt_expa (  )  [inline]

Sets the time derivative of the expansion function on the surface at time t (considering to protect this function).

Definition at line 223 of file excision_surf.h.

References del_deriv(), and dt_expa.

Scalar& Excision_surf::set_expa (  )  [inline]

Sets the expansion function on the surface at time t (considering to protect this function).

Definition at line 219 of file excision_surf.h.

References del_deriv(), and expa.

void Excision_surf::set_expa_hyperb ( double  alph0,
double  beta0,
double  gamma0 
)

Sets a new value for expansion rescaled over lapse (and its derivative), obtained by hyperbolic evolution.

Parameters for the hyperbolic driver are determined by the function Excision_surf::get_evol_params_from_ID() so that the expansion stays of regularity $C^{1}$ throughout. All manipulated quantities are 2-dimensional.

Definition at line 98 of file set_expa_evol.C.

References delta_t, dt_expa, expa, Valeur::get_etat(), Map::get_mg(), Tensor::get_mp(), Scalar::get_spectral_va(), Scalar::lapang(), lapse, set_dt_expa(), set_expa(), Scalar::set_spectral_va(), Scalar::val_grid_point(), and Valeur::ylm().

void Excision_surf::set_expa_parab ( double  c_theta_lap,
double  c_theta_fin,
Scalar &  expa_fin 
)

Sets a new value for expansion rescaled over lapse (and its derivative), obtained by parabolic evolution.

Definition at line 19 of file set_expa_evol.C.

References Scalar::annule_hard(), delta_t, expa, Valeur::get_etat(), Map::get_mg(), Tensor::get_mp(), Scalar::get_spectral_va(), Scalar::lapang(), lapse, set_dt_expa(), set_expa(), Scalar::set_spectral_va(), Scalar::std_spectral_base(), Scalar::val_grid_point(), and Valeur::ylm().

Metric& Excision_surf::set_gamij (  )  [inline]

Sets the 3d metric of the TimeSlice.

Definition at line 208 of file excision_surf.h.

References del_deriv(), and gamij.

Sym_tensor& Excision_surf::set_Kij (  )  [inline]

Sets the extrinsic curvature.

Definition at line 211 of file excision_surf.h.

References del_deriv(), and Kij.

Scalar& Excision_surf::set_lapse (  )  [inline]

Sets the lapse function.

Definition at line 202 of file excision_surf.h.

References del_deriv(), and lapse.

Vector& Excision_surf::set_shift (  )  [inline]

Sets the shift vector field.

Definition at line 205 of file excision_surf.h.

References del_deriv(), and shift.

Spheroid& Excision_surf::set_sph (  )  [inline]

Sets a new spheroid from data.

Definition at line 196 of file excision_surf.h.

References del_deriv(), and sph.

Friends And Related Function Documentation

ostream& operator<< ( ostream &  ,
const Spheroid &   
) [friend]

Display.

Member Data Documentation

Scalar Excision_surf::conf_fact [protected]

The value of the conformal factor on the 3-slice.

Definition at line 52 of file excision_surf.h.

double Excision_surf::delta_t [protected]

The time step for evolution in parabolic drivers.

Definition at line 67 of file excision_surf.h.

Scalar Excision_surf::dt_expa [protected]

The time derivative of the expansion, derived from Einstein equations and arbitrary evolution.

Definition at line 76 of file excision_surf.h.

Scalar Excision_surf::expa [protected]

The 2d expansion, directly evolved from the initial excision with Einstein Equations.

Definition at line 73 of file excision_surf.h.

Metric Excision_surf::gamij [protected]

The 3-d metric on the slice.

Definition at line 61 of file excision_surf.h.

Sym_tensor Excision_surf::Kij [protected]

The 3-d extrinsic curvature on the slice.

Definition at line 64 of file excision_surf.h.

Scalar Excision_surf::lapse [protected]

The lapse defined on the 3 slice.

Definition at line 55 of file excision_surf.h.

double Excision_surf::no_of_steps [protected]

The internal number of timesteps for one iteration.

Definition at line 70 of file excision_surf.h.

Scalar* Excision_surf::p_derive_t_expa [mutable, protected]

Computation of an updated expansion scalar.

Definition at line 93 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_conf_fact_1 [mutable, protected]

Source of Neumann boundary condition on $ \psi $,.

Definition at line 83 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_conf_fact_2 [mutable, protected]

Source of Neumann boundary condition on $ \psi $,.

Definition at line 87 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_conf_fact_3 [mutable, protected]

Source of Neumann boundary condition on $ \psi $,.

Definition at line 88 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_conf_fact_4 [mutable, protected]

Source of Birichlet boundary condition on $ \psi $,.

Definition at line 89 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_lapse_1 [mutable, protected]

Source of Dirichlet boundary condition of $ N $.

Definition at line 84 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_lapse_2 [mutable, protected]

Source of Dirichlet boundary condition of $ N $.

Definition at line 90 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_lapse_3 [mutable, protected]

Source of Dirichlet condtion on $ N $, based on einstein equations.

Definition at line 91 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_lapse_4 [mutable, protected]

Source of Dirichlet condtion on $ N $, based on einstein equations (conservation of isotropic gauge).

Definition at line 92 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_Npsi_1 [mutable, protected]

Source of Neumann boundary condition on $ \psi $.

Definition at line 86 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_Npsi_2 [mutable, protected]

Source of Dirichlet boundary condition on $ N \psi $.

Definition at line 97 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_Npsi_3 [mutable, protected]

Source of Dirichlet boundary condition on $ N \psi $.

Definition at line 98 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_Npsi_4 [mutable, protected]

Source of Dirichlet boundary condition on $ N \psi $.

Definition at line 99 of file excision_surf.h.

Scalar* Excision_surf::p_get_BC_Npsi_5 [mutable, protected]

Source of Neumann boundary condition on $ N \psi $.

Definition at line 100 of file excision_surf.h.

Vector* Excision_surf::p_get_BC_shift_1 [mutable, protected]

Source of Dirichlet BC for the shift vector $ \beta^{i} $.

Definition at line 85 of file excision_surf.h.

Vector* Excision_surf::p_get_BC_shift_2 [mutable, protected]

Source of Dirichlet BC for the shift vector $ \beta^{i} $.

Definition at line 94 of file excision_surf.h.

Vector* Excision_surf::p_get_BC_shift_3 [mutable, protected]

Source of Dirichlet BC for the shift vector $ \beta^{i} $, partly derived from kinematical relation.

Definition at line 95 of file excision_surf.h.

Vector* Excision_surf::p_get_BC_shift_4 [mutable, protected]

Source of Dirichlet BC for the shift vector $ \beta^{i} $, partly from projection of Einstein Equations.

Definition at line 96 of file excision_surf.h.

Vector Excision_surf::shift [protected]

The Shift 3-vector on the slice.

Definition at line 58 of file excision_surf.h.

Spheroid Excision_surf::sph [protected]

The associated Spheroid object.

Definition at line 49 of file excision_surf.h.

The documentation for this class was generated from the following files:
Generated on 7 Oct 2014 for LORENE by  doxygen 1.6.1