79 for (
int h=0; h<=n0; h++) {
89 for (
int h=0; h<
dzpuis; h++) {
97 for (
int h=
dzpuis; h<=n0; h++) {
103 for (
int l=0; l<nzm1; l++) {
108 for (
int k=0; k<np; k++) {
109 for (
int j=0; j<nt; j++) {
119 cout <<
"Term in 1/r^" << h << endl ;
120 cout <<
"-------------" << endl ;
123 double vmin = (*vu[h])(nzm1, 0, 0, 0) ;
126 cout <<
" Values at the point (phi_k, theta_j) : " << endl ;
127 cout.precision(precis) ;
128 cout.setf(ios::showpoint) ;
129 for (
int k=0; k<np; k++) {
130 cout <<
" k=" << k <<
" : " ;
131 for (
int j=0; j<nt; j++) {
132 double xx = (*vu[h])(nzm1, k, j, 0) ;
133 cout <<
" " << setw(precis) << xx ;
135 vmin = ( xx < vmin ) ? xx : vmin ;
136 vmax = ( xx > vmax ) ? xx : vmax ;
141 cout <<
"Minimum value on S^2 : " << vmin << endl ;
142 cout <<
"Maximum value on S^2 : " << vmax << endl ;
143 cout <<
"L^1 norm on S^2 : " << ndec << endl ;
147 for (
int k=0; k<np; k++) {
148 for (
int j=0; j<nt; j++) {
149 double v_inf = (*vu[h])(nzm1, k, j, 0) ;
150 for (
int i=0; i<nr; i++) {
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.
Tensor field of valence 0 (or component of a tensorial field).
Values and coefficients of a (real-value) function.
int dzpuis
Power of r by which the quantity represented by this must be divided in the compactified external do...
void set_base(const Base_val &)
Sets the bases for spectral expansions (member base )
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
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.
Base_val base
Bases on which the spectral expansion is performed.
void mult_r_ced()
Multiplication by r in the compactified external domain (CED), the dzpuis flag is not changed...
Mtbl * c
Values of the function at the points of the multi-grid.
int get_nzone() const
Returns the number of domains.
Valeur va
The numerical value of the Scalar.
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.
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
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.
const Map *const mp
Mapping on which the numerical values at the grid points are defined.
Tbl ** t
Array (size nzone ) of pointers on the Tbl 's.
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).