141 assert (ci.
get_etat() != ETATNONDEF) ;
149 assert( ci.
get_etat() == ETATQCQ ) ;
160 resu = (ci.
va).dsdx() ;
179 tmp_ced.
set(nzm1) -= 2. * ci.
va(nzm1) ;
182 resu = tmp + tmp_ced ;
193 tmp2.
base = (ci.
va).dsdx().base ;
195 tmp2.annule(0, nz-2) ;
196 tmp2.mult_xm1_zec() ;
198 tmp2.
set(nzm1) -= 4*ci.
va(nzm1) ;
200 tmp2.mult_xm1_zec() ;
202 resu = tmp + tmp2 /
xsr ;
208 cerr <<
"Map_af::dsdxi: unexpected value of input dzpuis !\n" 209 <<
" ci.get_dzpuis() = " << ci.
get_dzpuis() << endl ;
217 (resu.
va).set_base( (ci.
va).dsdx().get_base() ) ;
225 assert (uu.
get_etat() != ETATNONDEF) ;
234 assert( uu.
get_etat() == ETATQCQ ) ;
244 resu = uuva.
dsdx() * unity ;
262 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
284 assert (ci.
get_etat() != ETATNONDEF) ;
292 assert( ci.
get_etat() == ETATQCQ ) ;
303 resu = (ci.
va).dsdx() *
dxdr ;
322 tmp_ced.
set(nzm1) -= 2. * ci.
va(nzm1) ;
325 resu = tmp + tmp_ced ;
336 tmp2.
base = (ci.
va).dsdx().base ;
338 tmp2.annule(0, nz-2) ;
339 tmp2.mult_xm1_zec() ;
341 tmp2.
set(nzm1) -= 4*ci.
va(nzm1) ;
343 tmp2.mult_xm1_zec() ;
345 resu = tmp + tmp2 /
xsr ;
351 cerr <<
"Map_af::dsdr: unexpected value of input dzpuis !\n" 352 <<
" ci.get_dzpuis() = " << ci.
get_dzpuis() << endl ;
360 (resu.
va).set_base( (ci.
va).dsdx().get_base() ) ;
368 assert (uu.
get_etat() != ETATNONDEF) ;
376 assert( uu.
get_etat() == ETATQCQ ) ;
404 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
422 assert (uu.
get_etat() != ETATNONDEF) ;
430 assert( uu.
get_etat() == ETATQCQ ) ;
459 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
480 assert (ci.
get_etat() != ETATNONDEF) ;
488 assert( ci.
get_etat() == ETATQCQ ) ;
508 tmp.set_base(sauve_base) ;
529 tmp.set_base(sauve_base) ;
532 resu = tmp + tmp_ced ;
544 tmp_ced.mult_xm1_zec() ;
552 resu = tmp + tmp_ced /
xsr ;
560 cerr <<
"Map_af::srdsdt: unexpected value of input dzpuis !\n" 561 <<
" ci.get_dzpuis() = " << ci.
get_dzpuis() << endl ;
573 assert (uu.
get_etat() != ETATNONDEF) ;
581 assert( uu.
get_etat() == ETATQCQ ) ;
598 tmp.set_base(sauve_base) ;
619 tmp.set_base(sauve_base) ;
622 resu = tmp + tmp_ced ;
638 assert (ci.
get_etat() != ETATNONDEF) ;
646 assert( ci.
get_etat() == ETATQCQ ) ;
667 tmp.set_base(sauve_base) ;
688 tmp.set_base(sauve_base) ;
691 resu = tmp + tmp_ced ;
703 tmp_ced.mult_xm1_zec() ;
711 resu = tmp + tmp_ced /
xsr ;
719 cerr <<
"Map_af::srstdsdp: unexpected value of input dzpuis !\n" 720 <<
" ci.get_dzpuis() = " << ci.
get_dzpuis() << endl ;
734 assert (uu.
get_etat() != ETATNONDEF) ;
742 assert( uu.
get_etat() == ETATQCQ ) ;
762 tmp.set_base(sauve_base) ;
784 tmp.set_base(sauve_base) ;
787 resu = tmp + tmp_ced ;
802 assert (ci.
get_etat() != ETATNONDEF) ;
810 assert( ci.
get_etat() == ETATQCQ ) ;
828 assert (ci.
get_etat() != ETATNONDEF) ;
836 assert( ci.
get_etat() == ETATQCQ ) ;
const Map * get_mp() const
Returns the mapping.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
const Valeur & dsdt() const
Returns of *this.
const Valeur & dsdx() const
Returns of *this.
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
virtual void srstdsdp(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
void coef() const
Computes the coeffcients of *this.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
virtual void dsdradial(const Scalar &, Scalar &) const
Computes of a Scalar.
virtual void srdsdt(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_etat() const
Returns the logical state.
Tensor field of valence 0 (or component of a tensorial field).
const Valeur & dsdp() const
Returns of *this.
void mult_xm1_zec()
Applies the following operator to *this : \ Id (r sampling = RARE, FIN ) \ (r -sampling = UNSURR ) ...
const Valeur & sx() const
Returns (r -sampling = RARE ) \ Id (r sampling = FIN ) \ (r -sampling = UNSURR ) ...
const Scalar & dsdx() const
Returns of *this , where .
Values and coefficients of a (real-value) function.
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
void annule(int l)
Sets the Valeur to zero in a given domain.
void set(const Map *mp, Mtbl *(*construct)(const Map *))
Semi-constructor from a mapping and a method.
void set_base(const Base_val &)
Sets the bases for spectral expansions (member base )
virtual void stdsdp(const Scalar &uu, Scalar &resu) const
Computes of a Scalar.
void set_dzpuis(int)
Modifies the dzpuis flag.
virtual void dsdxi(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
const Valeur & ssint() const
Returns of *this.
Coord dxdr
in the nucleus and in the non-compactified shells; \ in the compactified outer domain.
Base_val base
Bases on which the spectral expansion is performed.
int get_dzpuis() const
Returns dzpuis.
const Valeur & stdsdp() const
Returns of *this.
int get_nzone() const
Returns the number of domains.
Coord xsr
in the nucleus; \ 1/R in the non-compactified shells; \ in the compactified outer domain...
void set_spectral_base(const Base_val &)
Sets the spectral bases of the Valeur va
const Mg3d * mg
Pointer on the multi-grid Mgd3 on which this is defined.
Bases of the spectral expansions.
int get_dzpuis() const
Returns dzpuis.
void set_dzpuis(int)
Set a value to dzpuis.
Valeur & set_spectral_va()
Returns va (read/write version)
int get_type_r(int l) const
Returns the type of sampling in the radial direction in domain no.
const Map & get_mp() const
Returns the mapping.
Valeur va
The numerical value of the Cmp.
virtual void dsdr(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
virtual void dsdt(const Scalar &uu, Scalar &resu) const
Computes of a Scalar.
const Valeur & get_spectral_va() const
Returns va (read only version)
Coord r
r coordinate centered on the grid
const Base_val & get_base() const
Return the bases for spectral expansions (member base )
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).