111 #include "time_slice.h" 113 #include "graphique.h" 122 const Scalar* p_trace_stress)
const {
128 if (p_ener_dens != 0x0) ener_dens = *(p_ener_dens) ;
131 Scalar trace_stress(map) ;
132 if (p_trace_stress != 0x0) trace_stress = *(p_trace_stress) ;
146 - 0.5 *
contract( dhh, 0, 1, dtgam, 0, 2 ).trace(
tgam()) ;
150 source_npsi +=
npsi()*( 0.5*qpig*(ener_dens + 2.*trace_stress )/(
psi()*
psi() )
151 + 0.875*aa_quad/(
psi4()*
psi4() ) + 0.125*tildeR ) ;
163 "Absolute error in the resolution of the equation for N") ;
181 if (p_ener_dens != 0x0) ener_dens = *(p_ener_dens) ;
197 - 0.5 *
contract( dhh, 0, 1, dtgam, 0, 2 ).trace(
tgam()) ;
201 source_psi -= 0.5*qpig*ener_dens/
psi()
202 + 0.125*( aa_quad*
pow(
psi(), -7) - tildeR*
psi() ) ;
214 "Absolute error in the resolution of the equation for Psi") ;
243 Vector beta_new = source_beta.
poisson(0.3333333333333333,
ff, method) ;
250 maxabs(test_beta - source_beta,
251 "Absolute error in the resolution of the equation for beta") ;
virtual const Vector & beta() const
shift vector at the current time step (jtime )
virtual const Scalar & psi() const
Conformal factor relating the physical metric to the conformal one: .
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.
const Scalar & psi4() const
Factor at the current time step (jtime ).
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
virtual const Metric & tgam() const
Conformal metric Returns the value at the current time step (jtime ).
Standard units of space, time and mass.
const Scalar & laplacian(int ced_mult_r=4) const
Returns the Laplacian of *this.
Tensor field of valence 0 (or component of a tensorial field).
Scalar poisson() const
Solves the scalar Poisson equation with *this as a source.
Base class for coordinate mappings.
Tensor up_down(const Metric &gam) const
Computes a new tensor by raising or lowering all the indices of *this .
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
Tensor field of valence 1.
const Tensor_sym & derive_cov(const Metric &gam) const
Returns the covariant derivative of this with respect to some metric .
virtual Scalar solve_psi(const Scalar *ener_dens=0x0) const
Solves the elliptic equation for the conformal factor $$ (Hamiltonian constraint).
const Vector & divergence(const Metric &) const
Returns the divergence of this with respect to a Metric .
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...
Vector poisson(double lambda, int method=6) const
Solves the vector Poisson equation with *this as a source.
const Scalar & divergence(const Metric &) const
The divergence of this with respect to a Metric .
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 const Scalar & npsi() const
Factor at the current time step (jtime ).
Cmp pow(const Cmp &, int)
Power .
virtual Sym_tensor aa() const
Conformal representation of the traceless part of the extrinsic curvature: .
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
const Tensor & derive_cov(const Metric &gam) const
Returns the covariant derivative of this with respect to some metric .
virtual const Sym_tensor & hata() const
Conformal representation of the traceless part of the extrinsic curvature: .
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values in the compactified ext...
virtual const Sym_tensor & cov() const
Read-only access to the covariant representation.
const Tensor & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
const Metric_flat & ff
Pointer on the flat metric with respect to which the conformal decomposition is performed.
virtual const Sym_tensor & hh(Param *par_bc=0x0, Param *par_mat=0x0) const
Deviation of the conformal metric from the flat metric : .
Symmetric tensors (with respect to two of their arguments).
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
const Map & get_mp() const
Returns the mapping.
virtual Vector solve_beta(int method=6) const
Solves the elliptic equation for the shift vector from (Eq.
Class intended to describe valence-2 symmetric tensors.
virtual const Scalar & nn() const
Lapse function N at the current time step (jtime )