30char valeur_smooth_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Valeur/valeur_smooth.C,v 1.4 2014/10/13 08:53:51 j_novak Exp $" ;
86 assert(
etat == ETATQCQ) ;
94 pot = *((*this).c_cf) ;
96 Tbl& ccf_nucl = *((pot.
c_cf)->t[nucl]) ;
97 Tbl& ccf_shell = *((pot.
c_cf)->t[nzet]) ;
103 Tbl nucl_kj(np, nt) ;
106 for (
int k=0 ; k<np ; k++) {
107 for (
int j=0 ; j<nt ; j++) {
110 for (
int i=0 ; i<nr ; i++) {
112 tmp += ccf_nucl(k, j, i) ;
115 nucl_kj.
set(k, j) = tmp ;
124 Tbl shell_kj(np, nt) ;
127 for (
int k=0 ; k<np ; k++) {
128 for (
int j=0 ; j<nt ; j++) {
131 for (
int i=0 ; i<nr-1 ; i++) {
133 tmp2 +=
pow(-1., i) * ccf_shell(k, j, i) ;
136 shell_kj.
set(k, j) = tmp2 ;
150 for (
int k=0 ; k<np ; k++) {
151 for (
int j=0 ; j<nt ; j++) {
153 uuva_cf.
set(nzet, k, j, nr-1) = nucl_kj(k, j) - shell_kj(k, j) ;
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Coefficients storage for the multi-domain spectral method.
Tbl & set(int l)
Read/write of the Tbl containing the coefficients in a given domain.
Tbl ** t
Array (size nzone ) of pointers on the Tbl 's which contain the spectral coefficients in each domain.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
double & set(int i)
Read/write of a particular element (index i) (1D case)
Values and coefficients of a (real-value) function.
void set_etat_cf_qcq()
Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl_cf c_c...
const Mg3d * mg
Multi-grid Mgd3 on which this is defined.
void coef_i() const
Computes the physical value of *this.
Mtbl_cf * c_cf
Coefficients of the spectral expansion of the function.
int etat
Logical state (ETATNONDEF , ETATQCQ or ETATZERO ).
void smooth(int nzet, Valeur &uuva) const
Changes the function *this as a smooth one when there exists a discontinuity between the nucleus and ...
Cmp pow(const Cmp &, int)
Power .