100 assert (ci.
get_etat() != ETATNONDEF) ;
107 assert( ci.
get_etat() == ETATQCQ ) ;
112 resu = (ci.
va).dsdx() ;
114 (resu.
va).base = (ci.
va).
dsdx().base ;
128 assert (uu.
get_etat() != ETATNONDEF) ;
135 assert( uu.
get_etat() == ETATQCQ ) ;
166 tmp_ced = tmp_ced.mult_x() ;
167 tmp_ced = tmp_ced /
xsr ;
169 tmp_ced.annule(0, nz-2) ;
170 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
192 assert (ci.
get_etat() != ETATNONDEF) ;
199 assert( ci.
get_etat() == ETATQCQ ) ;
204 resu = (ci.
va).dsdx() *
dxdr ;
206 (resu.
va).base = (ci.
va).dsdx().base ;
220 assert (uu.
get_etat() != ETATNONDEF) ;
227 assert( uu.
get_etat() == ETATQCQ ) ;
255 tmp_ced = tmp_ced *
dxdr ;
256 tmp_ced.set_base(sauve_base) ;
258 tmp_ced = tmp_ced.mult_x() ;
259 tmp_ced = tmp_ced /
xsr ;
261 tmp_ced.annule(0, nz-2) ;
262 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
279 assert (uu.
get_etat() != ETATNONDEF) ;
286 assert( uu.
get_etat() == ETATQCQ ) ;
314 tmp_ced = tmp_ced *
dxdr ;
315 tmp_ced.set_base(sauve_base) ;
317 tmp_ced = tmp_ced.mult_x() ;
318 tmp_ced = tmp_ced /
xsr ;
320 tmp_ced.annule(0, nz-2) ;
321 tmp_ced.
set(nzm1) -= dzp * uuva(nzm1) ;
342 assert (ci.
get_etat() != ETATNONDEF) ;
350 assert( ci.
get_etat() == ETATQCQ ) ;
359 srdfdt = srdfdt.
dsdt() ;
360 srdfdt = srdfdt.
sx() ;
364 srdfdt = srdfdt *
xsr ;
366 srdfdt.
base = sauve_base ;
373 adfdx = adfdx.
dsdx() ;
375 sauve_base = adfdx.
base ;
377 adfdx.
base = sauve_base ;
382 resu = srdfdt - adfdx ;
396 assert (uu.
get_etat() != ETATNONDEF) ;
404 assert( uu.
get_etat() == ETATQCQ ) ;
416 srdfdt = srdfdt.
dsdt() ;
418 srdfdt = srdfdt.
sx() ;
422 srdfdt = srdfdt *
xsr ;
424 srdfdt.
base = sauve_base ;
430 adfdx = adfdx.
dsdx() ;
432 sauve_base = adfdx.
base ;
434 adfdx.
base = sauve_base ;
441 resu = srdfdt - adfdx ;
456 Valeur tmp = srdfdt - adfdx ;
462 Valeur tmp_ced = - adfdx ;
466 tmp_ced = tmp_ced.
mult_x() ;
468 tmp_ced = tmp_ced /
xsr ;
470 tmp_ced = tmp_ced + uuva.
dsdt() ;
474 resu = tmp + tmp_ced ;
490 assert (ci.
get_etat() != ETATNONDEF) ;
508 srstdfdp = srstdfdp.
dsdp() ;
509 srstdfdp = srstdfdp.
ssint() ;
510 srstdfdp = srstdfdp.
sx() ;
514 srstdfdp = srstdfdp *
xsr ;
516 srstdfdp.
base = sauve_base ;
522 bdfdx = bdfdx.
dsdx() ;
524 sauve_base = bdfdx.
base ;
526 bdfdx.
base = sauve_base ;
531 resu = srstdfdp - bdfdx ;
545 assert (uu.
get_etat() != ETATNONDEF) ;
553 assert( uu.
get_etat() == ETATQCQ ) ;
566 srstdfdp = srstdfdp.
dsdp() ;
567 srstdfdp = srstdfdp.
ssint() ;
568 srstdfdp = srstdfdp.
sx() ;
572 srstdfdp = srstdfdp *
xsr ;
574 srstdfdp.
base = sauve_base ;
580 bdfdx = bdfdx.
dsdx() ;
582 sauve_base = bdfdx.
base ;
584 bdfdx.
base = sauve_base ;
591 resu = srstdfdp - bdfdx ;
605 Valeur tmp = srstdfdp - bdfdx ;
610 Valeur tmp_ced = - bdfdx ;
613 tmp_ced = tmp_ced.
mult_x() ;
615 tmp_ced = tmp_ced /
xsr ;
617 tmp_ced = tmp_ced + uuva.
dsdp().
ssint() ;
621 resu = tmp + tmp_ced ;
634 assert (ci.
get_etat() != ETATNONDEF) ;
644 assert( ci.
get_etat() == ETATQCQ ) ;
667 resu = dfdt - adfdx ;
679 assert (ci.
get_etat() != ETATNONDEF) ;
687 assert( ci.
get_etat() == ETATQCQ ) ;
711 resu = stdfdp - adfdx ;
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 ***.
const Cmp & dsdx() const
Returns of *this , where .
void coef() const
Computes the coeffcients of *this.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
virtual void dsdt(const Scalar &uu, Scalar &resu) const
Computes of a Scalar.
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.
virtual void srstdsdp(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
const Valeur & sx() const
Returns (r -sampling = RARE ) \ Id (r sampling = FIN ) \ (r -sampling = UNSURR ) ...
virtual void dsdxi(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
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 )
Coord srstdrdp
in the nucleus and in the non-compactified shells; \ in the compactified outer domain.
void set_dzpuis(int)
Modifies the dzpuis flag.
virtual void stdsdp(const Scalar &uu, Scalar &resu) const
Computes of a Scalar.
const Valeur & mult_x() const
Returns (r -sampling = RARE ) \ Id (r sampling = FIN ) \ (r -sampling = UNSURR ) ...
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Coord stdrdp
in the nucleus and in the non-compactified shells; \ in the compactified external domain (CED)...
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.
virtual void dsdr(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
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...
Coord drdt
in the nucleus and in the non-compactified shells; \ in the compactified external domain (CED)...
virtual void dsdradial(const Scalar &uu, Scalar &resu) const
Computes of a Scalar if the description is affine and if it is logarithmic.
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.
bool check_dzpuis(int dzi) const
Returns false if the last domain is compactified and *this is not zero in this domain and dzpuis is n...
virtual void srdsdt(const Cmp &ci, Cmp &resu) const
Computes of a Cmp.
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.
Coord srdrdt
in the nucleus and in the non-compactified shells; \ in the compactified outer domain.
const Map & get_mp() const
Returns the mapping.
Valeur va
The numerical value of the Cmp.
const Valeur & get_spectral_va() const
Returns va (read only version)
const Base_val & get_base() const
Return the bases for spectral expansions (member base )