87 for (
int h=0; h<=n0; h++) {
97 for (
int h=0; h<
dzpuis; h++) {
105 for (
int h=
dzpuis; h<=n0; h++) {
111 for (
int l=0; l<nzm1; l++) {
116 for (
int k=0; k<np; k++) {
117 for (
int j=0; j<nt; j++) {
118 vu[h]->
set(nzm1, k, j, 0) = uu(nzm1, k, j, nrm1) ;
127 cout <<
"Term in 1/r^" << h << endl ;
128 cout <<
"-------------" << endl ;
131 double vmin = (*vu[h])(nzm1, 0, 0, 0) ;
134 cout <<
" Values at the point (phi_k, theta_j) : " << endl ;
135 cout.precision(precis) ;
136 cout.setf(ios::showpoint) ;
137 for (
int k=0; k<np; k++) {
138 cout <<
" k=" << k <<
" : " ;
139 for (
int j=0; j<nt; j++) {
140 double xx = (*vu[h])(nzm1, k, j, 0) ;
141 cout <<
" " << setw(precis) << xx ;
143 vmin = ( xx < vmin ) ? xx : vmin ;
144 vmax = ( xx > vmax ) ? xx : vmax ;
149 cout <<
"Minimum value on S^2 : " << vmin << endl ;
150 cout <<
"Maximum value on S^2 : " << vmax << endl ;
151 cout <<
"L^1 norm on S^2 : " << ndec << endl ;
155 for (
int k=0; k<np; k++) {
156 for (
int j=0; j<nt; j++) {
157 double v_inf = (*vu[h])(nzm1, k, j, 0) ;
158 for (
int i=0; i<nr; i++) {
159 uu.
set(nzm1, k, j, i) -= v_inf ;
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
void mult_r_zec()
Multiplication by r in the external compactified domain (ZEC)
Values and coefficients of a (real-value) function.
void set_base(const Base_val &)
Sets the bases for spectral expansions (member base )
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
int dzpuis
Power of r by which the quantity represented by this must be divided in the external compactified zo...
Base_val base
Bases on which the spectral expansion is performed.
Mtbl * c
Values of the function at the points of the multi-grid.
int get_nzone() const
Returns the number of domains.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Tbl & set(int l)
Read/write of the value in a given domain.
const Map * mp
Reference mapping.
Valeur ** asymptot(int n, const int flag=0) const
Asymptotic expansion at r = infinity.
void set_etat_c_qcq()
Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl c )...
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_type_r(int l) const
Returns the type of sampling in the radial direction in domain no.
Tbl ** t
Array (size nzone ) of pointers on the Tbl 's.
Valeur va
The numerical value of the Cmp.
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).