28 char 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)
103 assert ((jp>=0) && (jp<N_MET_MAX)) ;
111 int method_poisson)
const {
115 assert ((jp>=0) && (jp<N_MET_MAX)) ;
121 for (
int ic=0; ic<
n_comp; ic++) {
122 assert(
cmp[ic]->check_dzpuis(4)) ;
144 int method_poisson)
const {
148 assert ((jp>=0) && (jp<N_MET_MAX)) ;
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++) {
165 dzp5 = dzp5 || hhh(i).check_dzpuis(5) ;
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" ;
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++){
196 cout <<
diffrel(lapw(i),hhh(i) )(l) <<
" " ;
void set_dependance(const Metric &) const
To be used to describe the fact that the derivatives members have been calculated with met .
Sym_tensor ope_killing(const Metric &gam) const
Computes the Killing operator associated with a given metric.
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...
Class intended to describe valence-2 symmetric tensors.
Metric for tensor calculation.
Tbl diffrel(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (norme version).
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
virtual void dec_dzpuis(int dec=1)
Decreases by dec units the value of dzpuis and changes accordingly the values in the compactified ext...
Scalar ** cmp
Array of size n_comp of pointers onto the components.
int get_place_met(const Metric &) const
Returns the position of the pointer on metre in the array met_depend .
const Scalar & divergence(const Metric &) const
The divergence of this with respect to a Metric .
Sym_tensor_trans * p_transverse[N_MET_MAX]
Array of the transverse part of the tensor with respect to various metrics, transverse meaning diver...
Flat metric for tensor calculation.
Vector poisson(double lambda, int method=6) const
Solves the vector Poisson equation with *this as a source.
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 Map *const mp
Mapping on which the numerical values at the grid points are defined.
virtual void del_deriv() const
Deletes the derived quantities.
const Metric & get_met_div() const
Returns the metric with respect to which the divergence and the trace are defined.
int get_nzone() const
Returns the number of domains.
int get_index_type(int i) const
Gives the type (covariant or contravariant) of the index number i .
const Base_vect * triad
Vectorial basis (triad) with respect to which the tensor components are defined.
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 Tensor_sym & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
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.
Tensor field of valence 1.
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 Tensor & divergence(const Metric &gam) const
Computes the divergence of this with respect to some metric .
Itbl type_indice
1D array of integers (class Itbl ) of size valence containing the type of each index: COV for a cov...
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values in the compactified ext...
Transverse symmetric tensors of rank 2.
const Vector & divergence(const Metric &) const
Returns the divergence of this with respect to a Metric .
int n_comp
Number of stored components, depending on the symmetry.