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;
53 nn_point.annule_domain(nz -1);
55 nn_point.std_spectral_base();
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 ;
91 sym_tmp.inc_dzpuis() ;
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);
146 dNdt.annule_domain(nz-1);
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
Scalar * p_get_BC_lapse_4
Source of Dirichlet condtion on , based on einstein equations (conservation of isotropic gauge) ...
double delta_t
The time step for evolution in parabolic drivers.
Sym_tensor Kij
The 3-d extrinsic curvature on the slice.
Cmp log(const Cmp &)
Neperian logarithm.
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
Standard units of space, time and mass.
void ylm()
Computes the coefficients of *this.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Flat metric for tensor calculation.
Tensor field of valence 0 (or component of a tensorial field).
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. ...
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
Tensor field of valence 1.
void annule_hard()
Sets the Scalar to zero in a hard way.
Scalar lapse
The lapse defined on the 3 slice.
Sym_tensor ope_killing_conf(const Metric &gam) const
Computes the conformal Killing operator associated with a given 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...
const Scalar & divergence(const Metric &) const
The divergence of this with respect to a Metric .
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 .
Scalar derive_lie(const Vector &v) const
Computes the derivative of this along a vector field v.
Spherical orthonormal vectorial bases (triads).
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Scalar conf_fact
The value of the conformal factor on the 3-slice.
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) ...
Valeur & set_spectral_va()
Returns va (read/write version)
Metric gamij
The 3-d metric on the slice.
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 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.
Vector shift
The Shift 3-vector on the slice.
Class intended to describe valence-2 symmetric tensors.
Tensor trace(int ind1, int ind2) const
Trace on two different type indices.