102 if ((*
this)(1,1).get_dzpuis() == 2)
109 assert ((jp>=0) && (jp<N_MET_MAX)) ;
117 int method_poisson)
const {
121 assert ((jp>=0) && (jp<N_MET_MAX)) ;
127 for (
int ic=0; ic<
n_comp; ic++) {
128 assert(
cmp[ic]->check_dzpuis(4)) ;
150 int method_poisson)
const {
154 assert ((jp>=0) && (jp<N_MET_MAX)) ;
160 cout <<
"Sym_tensor::longit_pot : the case of a non flat metric" 161 << endl <<
" is not treated yet !" << endl ;
170 for (
int i=1; i<=3; i++) {
171 dzp5 = dzp5 || (hhh(i).get_dzpuis() == 5);
175 if (dynamic_cast<const Map_af*>(
mp) != 0x0)
191 cout <<
"## Sym_tensor::longit_pot : test of Poisson : \n" ;
193 " Max absolute error in each domain on the vector Poisson equation:\n" ;
198 cout <<
" Relative error in each domain on the vector Poisson equation:\n" ;
199 for (
int i=1; i<=3; i++){
200 cout <<
" Comp. " << i <<
" : " ;
201 for (
int l=0; l<nz; l++){
202 cout <<
diffrel(lapw(i),hhh(i) )(l) <<
" " ;
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.
Metric for tensor calculation.
Sym_tensor ope_killing(const Metric &gam) const
Computes the Killing operator associated with a given metric.
int n_comp
Number of stored components, depending on the symmetry.
void set_dependance(const Metric &) const
To be used to describe the fact that the derivatives members have been calculated with met ...
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Flat metric for tensor calculation.
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...
const Base_vect * triad
Vectorial basis (triad) with respect to which the tensor components are defined.
int get_place_met(const Metric &) const
Returns the position of the pointer on metre in the array met_depend .
Tensor field of valence 1.
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 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 & divergence(const Metric &) const
Returns the divergence of this with respect to a Metric .
Itbl type_indice
1D array of integers (class Itbl ) of size valence containing the type of each index: COV for a cova...
virtual void del_deriv() const
Deletes the derived quantities.
int get_index_type(int i) const
Gives the type (covariant or contravariant) of the index number i .
Vector poisson(double lambda, int method=6) const
Solves the vector Poisson equation with *this as a source.
Scalar ** cmp
Array of size n_comp of pointers onto the components.
Vector * p_longit_pot[N_MET_MAX]
Array of the vector potential of the longitudinal part of the tensor with respect to various metrics ...
int get_nzone() const
Returns the number of domains.
const Metric & get_met_div() const
Returns the metric with respect to which the divergence and the trace are defined.
Transverse symmetric tensors of rank 2.
const Tensor & divergence(const Metric &gam) const
Computes the divergence of this with respect to some metric .
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_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
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 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...
const Map *const mp
Mapping on which the numerical values at the grid points are defined.
Class intended to describe valence-2 symmetric tensors.
Sym_tensor_trans * p_transverse[N_MET_MAX]
Array of the transverse part of the tensor with respect to various metrics, transverse meaning diver...