28char tslice_check_einstein_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Time_slice/tslice_check_einstein.C,v 1.10 2014/10/13 08:53:47 j_novak Exp $" ;
75#include "time_slice.h"
99 if (
verb)
comment =
"Check of the Hamiltonian constraint" ;
128 if (
verb)
comment =
"Check of the momentum constraint" ;
148 int nz =
dyn_lhs.get_mp().get_mg()->get_nzone() ;
157 const Scalar* ener = 0x0 ;
176 if (
new_e)
delete ener ;
181 + qpig * (*matter))) ;
190 if (
verb)
comment =
"Check of the dynamical equations" ;
201 for (
int i=1;
i<=3;
i++) {
202 for (
int j=1;
j<=
i;
j++) {
Time evolution with partial storage (*** under development ***).
TyT time_derive(int j, int n=2) const
Computes the time derivative at time step j by means of a n-th order scheme, from the values at steps...
int position(int j) const
Gives the position in the arrays step, the_time and val corresponding to the time step j.
Basic integer array class.
Tensor field of valence 0 (or component of a tensorial field).
Class intended to describe valence-2 symmetric tensors.
virtual const Sym_tensor & k_dd() const
Extrinsic curvature tensor (covariant components ) at the current time step (jtime )
int jtime
Time step index of the latest slice.
Evolution_std< Sym_tensor > k_dd_evol
Values at successive time steps of the covariant components of the extrinsic curvature tensor .
virtual const Scalar & trk() const
Trace K of the extrinsic curvature at the current time step (jtime )
virtual const Scalar & nn() const
Lapse function N at the current time step (jtime )
Tbl check_hamiltonian_constraint(const Scalar *energy_density=0x0, ostream &ost=cout, bool verb=true) const
Checks the level at which the hamiltonian constraint is verified.
int scheme_order
Order of the finite-differences scheme for the computation of time derivatives.
virtual const Vector & beta() const
shift vector at the current time step (jtime )
const Metric & gam() const
Induced metric at the current time step (jtime )
virtual const Sym_tensor & k_uu() const
Extrinsic curvature tensor (contravariant components ) at the current time step (jtime )
Tbl check_dynamical_equations(const Sym_tensor *strain_tensor=0x0, const Scalar *energy_density=0x0, ostream &ost=cout, bool verb=true) const
Checks the level at which the dynamical equations are verified.
Tbl check_momentum_constraint(const Vector *momentum_density=0x0, ostream &ost=cout, bool verb=true) const
Checks the level at which the momentum constraints are verified.
Tensor field of valence 1.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Tbl maxabs(const Tensor &aa, const char *comment=0x0, ostream &ost=cout, bool verb=true)
Maxima in each domain of the absolute values of the tensor components.
Standard units of space, time and mass.