28char sym_tensor_decomp_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Tensor/sym_tensor_decomp.C,v 1.14 2014/10/13 08:53:43 j_novak Exp $" ;
96 if ((*
this)(1,1).get_dzpuis() == 2)
154 cout <<
"Sym_tensor::longit_pot : the case of a non flat metric"
155 <<
endl <<
" is not treated yet !" <<
endl ;
164 for (
int i=1;
i<=3;
i++) {
169 if (
dynamic_cast<const Map_af*
>(
mp) != 0x0)
185 cout <<
"## Sym_tensor::longit_pot : test of Poisson : \n" ;
187 " Max absolute error in each domain on the vector Poisson equation:\n" ;
190 int nz =
mp->get_mg()->get_nzone() ;
192 cout <<
" Relative error in each domain on the vector Poisson equation:\n" ;
193 for (
int i=1;
i<=3;
i++){
194 cout <<
" Comp. " <<
i <<
" : " ;
195 for (
int l=0;
l<nz;
l++){
Time evolution with partial storage (*** under development ***).
Flat metric for tensor calculation.
Metric for tensor calculation.
Transverse symmetric tensors of rank 2.
Class intended to describe valence-2 symmetric tensors.
void set_longit_trans(const Vector &v, const Sym_tensor_trans &a)
Assigns the derived members p_longit_pot and p_transverse and updates the components accordingly.
const Sym_tensor_trans & transverse(const Metric &gam, Param *par=0x0, int method_poisson=6) const
Computes the transverse part of the tensor with respect to a given metric, transverse meaning diverg...
const Vector & longit_pot(const Metric &gam, Param *par=0x0, int method_poisson=6) const
Computes the vector potential of longitudinal part of the tensor (see documentation of method transv...
Sym_tensor_trans * p_transverse[N_MET_MAX]
Array of the transverse part of the tensor with respect to various metrics, transverse meaning diver...
virtual void del_deriv() const
Deletes the derived quantities.
Vector * p_longit_pot[N_MET_MAX]
Array of the vector potential of the longitudinal part of the tensor with respect to various metrics ...
const Vector & divergence(const Metric &) const
Returns the divergence of this with respect to a Metric .
Tensor field of valence 1.
Sym_tensor ope_killing(const Metric &gam) const
Computes the Killing operator associated with a given metric.
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 .
Tbl diffrel(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (norme version).
virtual void dec_dzpuis(int dec=1)
Decreases by dec units the value of dzpuis and changes accordingly the values in the compactified ext...
const Map *const mp
Mapping on which the numerical values at the grid points are defined.
int get_place_met(const Metric &) const
Returns the position of the pointer on metre in the array met_depend .
const Tensor & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
void set_dependance(const Metric &) const
To be used to describe the fact that the derivatives members have been calculated with met .
Scalar ** cmp
Array of size n_comp of pointers onto the components.
int n_comp
Number of stored components, depending on the symmetry.
const Tensor_sym & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
Itbl type_indice
1D array of integers (class Itbl ) of size valence containing the type of each index: COV for a cov...
const Base_vect * triad
Vectorial basis (triad) with respect to which the tensor components are defined.
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.