28char tslice_dirac_max_solve_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Time_slice/tslice_dirac_max_solve.C,v 1.18 2014/10/13 08:53:48 j_novak Exp $" ;
108#include "time_slice.h"
110#include "graphique.h"
123 const Map& map =
npsi().get_mp() ;
160 "Absolute error in the resolution of the equation for N") ;
176 const Map& map =
psi().get_mp() ;
184 const Vector& dpsi =
psi().derive_cov(ff) ;
211 "Absolute error in the resolution of the equation for Psi") ;
230 -
contract(
hh(), 0, 1,
beta().derive_cov(ff).derive_cov(ff), 1, 2)
231 - 0.3333333333333333*
contract(
hh(), 1,
beta().divergence(ff).derive_cov(ff), 0) ;
234 source_beta +=
sym_tmp.divergence(ff) ;
245 + 0.3333333333333333 * (
beta_new.divergence(ff)).derive_con(ff) ;
248 "Absolute error in the resolution of the equation for beta") ;
Time evolution with partial storage (*** under development ***).
Base class for coordinate mappings.
Tensor field of valence 0 (or component of a tensorial field).
Scalar poisson() const
Solves the scalar Poisson equation with *this as a source.
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values of the Scalar in the co...
Class intended to describe valence-2 symmetric tensors.
Symmetric tensors (with respect to two of their arguments).
virtual const Scalar & npsi() const
Factor at the current time step (jtime ).
virtual const Scalar & nn() const
Lapse function N at the current time step (jtime )
virtual const Sym_tensor & hata() const
Conformal representation of the traceless part of the extrinsic curvature: .
virtual const Scalar & psi() const
Conformal factor relating the physical metric to the conformal one: .
virtual Sym_tensor aa() const
Conformal representation of the traceless part of the extrinsic curvature: .
const Scalar & psi4() const
Factor at the current time step (jtime ).
virtual const Metric & tgam() const
Conformal metric Returns the value at the current time step (jtime ).
virtual const Vector & beta() const
shift vector at the current time step (jtime )
virtual const Sym_tensor & hh(Param *par_bc=0x0, Param *par_mat=0x0) const
Deviation of the conformal metric from the flat metric : .
virtual Scalar solve_psi(const Scalar *ener_dens=0x0) const
Solves the elliptic equation for the conformal factor $\Psi$ (Hamiltonian constraint).
virtual Scalar solve_npsi(const Scalar *ener_dens=0x0, const Scalar *trace_stress=0x0) const
Solves the elliptic equation for (maximal slicing condition + Hamiltonian constraint)
virtual Vector solve_beta(int method=6) const
Solves the elliptic equation for the shift vector from (Eq.
Tensor field of valence 1.
Vector poisson(double lambda, int method=6) const
Solves the vector Poisson equation with *this as a source.
Cmp pow(const Cmp &, int)
Power .
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values in the compactified ext...
const Tensor & derive_cov(const Metric &gam) const
Returns the covariant derivative of this with respect to some metric .
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.