4 #include "utilitaires.h"
9 #include "param_elliptic.h"
12 #include "excision_surf.h"
28 for (
int ii=1; ii<=3; ii++)
29 for(
int jj=1; jj<=3; jj++)
38 const Scalar& psi4 = ppsi*ppsi*ppsi*ppsi;
42 const Scalar qq = nn*ppsi*ppsi;
80 sym_tmp = nn * (8.* tdln_psi_u * tdln_psi_u)
81 + 4.*( tdln_psi_u * tdnn_u + tdnn_u * tdln_psi_u )
82 - 0.3333333333333333 *
83 ( nn * ( 8.*
contract(dln_psi, 0, tdln_psi_u, 0) )
84 + 8.*
contract(dln_psi, 0, tdnn_u, 0) ) *tgam_uu ;
96 sym_tmp -= 0.3333333333333333 * tmp *tgam_uu ;
98 ss -= sym_tmp / (psi4*ppsi*ppsi) ;
100 for (
int i=1; i<=3; i++) {
101 for (
int j=1; j<=i; j++) {
103 for (
int k=1; k<=3; k++) {
104 for (
int l=1; l<=3; l++) {
105 tmp += tgam_dd(k,l) * aa(i,k) * aa(j,l) ;
108 sym_tmp.
set(i,j) = tmp ;
112 tmp = psi4 * strain_tens.
trace(ff) ;
114 ss += (2.*nn) * ( sym_tmp );
115 Sym_tensor ss2 =2.*nn*( qpig*(psi4*strain_tens - 0.33333333333333 * tmp * tgam_uu));
140 source_hh += 0.6666666666666666* div_beta * l_beta - sym_tmp ;
144 dNdt = -source_hh(1,1)/2.*aa(1,1);
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.
const Tensor & derive_cov(const Metric &gam) const
Returns the covariant derivative of this with respect to some metric .
Cmp log(const Cmp &)
Neperian logarithm.
Scalar * p_get_BC_lapse_4
Source of Dirichlet condtion on , based on einstein equations (conservation of isotropic gauge)
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
Metric gamij
The 3-d metric on the slice.
Scalar derive_lie(const Vector &v) const
Computes the derivative of this along a vector field v.
const Scalar & divergence(const Metric &) const
The divergence of this with respect to a Metric .
Tensor up_down(const Metric &gam) const
Computes a new tensor by raising or lowering all the indices of *this .
Flat metric for tensor calculation.
Tensor field of valence 0 (or component of a tensorial field).
const Map & get_mp() const
Returns the mapping.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
Sym_tensor ope_killing_conf(const Metric &gam) const
Computes the conformal Killing operator associated with a given metric.
Valeur & set_spectral_va()
Returns va (read/write version)
Standard units of space, time and mass.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
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)
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
Spherical orthonormal vectorial bases (triads).
Vector shift
The Shift 3-vector on the slice.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Sym_tensor derive_lie(const Vector &v) const
Computes the Lie derivative of this with respect to some vector field v.
Tensor field of valence 1.
Tensor trace(int ind1, int ind2) const
Trace on two different type indices.
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
double delta_t
The time step for evolution in parabolic drivers.
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 void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
void annule_hard()
Sets the Scalar to zero in a hard way.
Scalar conf_fact
The value of the conformal factor on the 3-slice.
Sym_tensor Kij
The 3-d extrinsic curvature on the slice.
Scalar lapse
The lapse defined on the 3 slice.
void ylm()
Computes the coefficients of *this.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .