28char map_af_deriv_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Map/map_af_deriv.C,v 1.13 2014/10/13 08:53:02 j_novak Exp $" ;
138 assert (ci.
get_etat() != ETATNONDEF) ;
146 assert( ci.
get_etat() == ETATQCQ ) ;
157 resu = (ci.
va).dsdx() ;
176 tmp_ced.
set(nzm1) -= 2. * ci.
va(nzm1) ;
179 resu = tmp + tmp_ced ;
190 tmp2.
base = (ci.
va).dsdx().base ;
195 tmp2.
set(nzm1) -= 4*ci.
va(nzm1) ;
199 resu = tmp + tmp2 /
xsr ;
205 cerr <<
"Map_af::dsdxi: unexpected value of input dzpuis !\n"
206 <<
" ci.get_dzpuis() = " << ci.
get_dzpuis() << endl ;
214 (resu.
va).set_base( (ci.
va).dsdx().get_base() ) ;
222 assert (uu.
get_etat() != ETATNONDEF) ;
231 assert( uu.
get_etat() == ETATQCQ ) ;
241 resu = uuva.
dsdx() * unity ;
259 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
281 assert (ci.
get_etat() != ETATNONDEF) ;
289 assert( ci.
get_etat() == ETATQCQ ) ;
300 resu = (ci.
va).dsdx() *
dxdr ;
319 tmp_ced.
set(nzm1) -= 2. * ci.
va(nzm1) ;
322 resu = tmp + tmp_ced ;
333 tmp2.
base = (ci.
va).dsdx().base ;
338 tmp2.
set(nzm1) -= 4*ci.
va(nzm1) ;
342 resu = tmp + tmp2 /
xsr ;
348 cerr <<
"Map_af::dsdr: unexpected value of input dzpuis !\n"
349 <<
" ci.get_dzpuis() = " << ci.
get_dzpuis() << endl ;
357 (resu.
va).set_base( (ci.
va).dsdx().get_base() ) ;
365 assert (uu.
get_etat() != ETATNONDEF) ;
373 assert( uu.
get_etat() == ETATQCQ ) ;
401 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
419 assert (uu.
get_etat() != ETATNONDEF) ;
427 assert( uu.
get_etat() == ETATQCQ ) ;
456 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
477 assert (ci.
get_etat() != ETATNONDEF) ;
485 assert( ci.
get_etat() == ETATQCQ ) ;
529 resu = tmp + tmp_ced ;
549 resu = tmp + tmp_ced /
xsr ;
557 cerr <<
"Map_af::srdsdt: unexpected value of input dzpuis !\n"
558 <<
" ci.get_dzpuis() = " << ci.
get_dzpuis() << endl ;
570 assert (uu.
get_etat() != ETATNONDEF) ;
578 assert( uu.
get_etat() == ETATQCQ ) ;
619 resu = tmp + tmp_ced ;
635 assert (ci.
get_etat() != ETATNONDEF) ;
643 assert( ci.
get_etat() == ETATQCQ ) ;
688 resu = tmp + tmp_ced ;
708 resu = tmp + tmp_ced /
xsr ;
716 cerr <<
"Map_af::srstdsdp: unexpected value of input dzpuis !\n"
717 <<
" ci.get_dzpuis() = " << ci.
get_dzpuis() << endl ;
731 assert (uu.
get_etat() != ETATNONDEF) ;
739 assert( uu.
get_etat() == ETATQCQ ) ;
784 resu = tmp + tmp_ced ;
799 assert (ci.
get_etat() != ETATNONDEF) ;
807 assert( ci.
get_etat() == ETATQCQ ) ;
825 assert (ci.
get_etat() != ETATNONDEF) ;
833 assert( ci.
get_etat() == ETATQCQ ) ;
Bases of the spectral expansions.
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
int get_etat() const
Returns the logical state.
Valeur va
The numerical value of the Cmp
int get_dzpuis() const
Returns dzpuis.
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
void set_dzpuis(int)
Set a value to dzpuis.
const Map * get_mp() const
Returns the mapping.
virtual void srstdsdp(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
virtual void dsdradial(const Scalar &, Scalar &) const
Computes of a Scalar.
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.
virtual void stdsdp(const Scalar &uu, Scalar &resu) const
Computes of a Scalar.
virtual void srdsdt(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
virtual void dsdxi(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
Coord xsr
in the nucleus; \ 1/R in the non-compactified shells; \ in the compactified outer domain.
Coord dxdr
in the nucleus and in the non-compactified shells; \ in the compactified outer domain.
Coord r
r coordinate centered on the grid
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Mg3d * mg
Pointer on the multi-grid Mgd3 on which this is defined
int get_nzone() const
Returns the number of domains.
int get_type_r(int l) const
Returns the type of sampling in the radial direction in domain no.
Tensor field of valence 0 (or component of a tensorial field).
int get_dzpuis() const
Returns dzpuis.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Valeur & set_spectral_va()
Returns va (read/write version)
const Valeur & get_spectral_va() const
Returns va (read only version)
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
void set_dzpuis(int)
Modifies the dzpuis flag.
void set_spectral_base(const Base_val &)
Sets the spectral bases of the Valeur va
Values and coefficients of a (real-value) function.
const Valeur & dsdp() const
Returns of *this.
const Valeur & sx() const
Returns (r -sampling = RARE ) \ Id (r sampling = FIN ) \ (r -sampling = UNSURR )
void mult_xm1_zec()
Applies the following operator to *this : \ Id (r sampling = RARE, FIN ) \ (r -sampling = UNSURR )
const Base_val & get_base() const
Return the bases for spectral expansions (member base )
const Valeur & stdsdp() const
Returns of *this.
void set_base(const Base_val &)
Sets the bases for spectral expansions (member base )
const Valeur & dsdt() const
Returns of *this.
void annule(int l)
Sets the Valeur to zero in a given domain.
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).
const Valeur & dsdx() const
Returns of *this.
const Valeur & ssint() const
Returns of *this.
void coef() const
Computes the coeffcients of *this.
Base_val base
Bases on which the spectral expansion is performed.
const Map & get_mp() const
Returns the mapping.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.