77 assert (lambda != -1) ;
86 assert (
etat != ETATNONDEF) ;
89 if (
etat == ETATZERO) {
94 for (
int i=0; i<3; i++) {
104 for (
int i=0 ; i<3 ; i++)
105 assert ((*
this)(i).check_dzpuis(4)) ;
114 vecteur_div.set_etat_qcq() ;
115 dvect_div.set_etat_qcq() ;
116 souvect_regu.set_etat_qcq() ;
117 souvect_div.set_etat_qcq() ;
123 for (
int i=0 ; i<2 ; i++) {
126 (*this)(i).poisson_regular(k_div, nzet, unsgam1, *par,
128 vecteur_regu.
set(i), vecteur_div.
set(i),
130 souvect_regu.set(i), souvect_div.set(i)) ;
137 (*this)(2).poisson(*par, vecteur.
set(2)) ;
146 assert (source_scal().check_dzpuis(3)) ;
153 Cmp souscal_regu(
mp) ;
154 Cmp souscal_div(
mp) ;
158 dscal_div.set_etat_qcq() ;
160 souscal_regu.std_base_scal() ;
163 source_scal().poisson_regular(k_div, nzet, unsgam1, *par,
165 scalaire_regu.
set(), scalaire_div.
set(),
166 dscal_div, souscal_regu, souscal_div) ;
186 for (
int i=0 ; i<3 ; i++)
187 shift.
set(i) = (lambda+2)/2/(lambda+1) * vecteur(i)
188 - (lambda/2/(lambda+1)) * (dxsi(i) + dp(i)) ;
Itbl type_indice
Array of size valence contening the type of each index, COV for a covariant one and CON for a contrav...
const Map *const mp
Reference mapping.
const Base_vect * triad
Vectorial basis (triad) with respect to which the tensor components are defined.
int get_type_indice(int i) const
Returns the type of the index number i .
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
void dec2_dzpuis()
dzpuis -= 2 ;
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
double & set(int i)
Read/write of a particular element (index i) (1D case)
int etat
Logical state ETATZERO , ETATQCQ or ETATNONDEF.
int get_valence() const
Returns the valence.
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
double poids
For tensor densities: the weight.
virtual Param * donne_para_poisson_vect(Param ¶, int i) const =0
Function intended to be used by Map::poisson_vect and Map::poisson_vect_oohara .
friend Tenseur skxk(const Tenseur &)
Contraction of the last index of (*this) with or , depending on the type of S .
void std_base_scal()
Sets the spectral bases of the Valeur va to the standard ones for a scalar.
Tbl & set(int l)
Read/write of the value in a given domain.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
void dec_dzpuis()
dzpuis -= 1 ;
const Metrique * metric
For tensor densities: the metric defining the conformal factor.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
void set_etat_zero()
Sets the logical state to ETATZERO (zero state).
void poisson_vect_regu(int k_div, int nzet, double unsgam1, double lambda, Param &par, Tenseur &shift, Tenseur &vect, Tenseur &scal) const
Solves the vectorial Poisson equation : .
Tensor handling *** DEPRECATED : use class Tensor instead ***.
const Tenseur & gradient() const
Returns the gradient of *this (Cartesian coordinates)