31 char scalar_asymptot_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Tensor/Scalar/scalar_asymptot.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $" ;
76 for (
int h=0; h<=n0; h++) {
86 for (
int h=0; h<
dzpuis; h++) {
94 for (
int h=
dzpuis; h<=n0; h++) {
100 for (
int l=0; l<nzm1; l++) {
105 for (
int k=0; k<np; k++) {
106 for (
int j=0; j<nt; j++) {
116 cout <<
"Term in 1/r^" << h << endl ;
117 cout <<
"-------------" << endl ;
120 double vmin = (*vu[h])(nzm1, 0, 0, 0) ;
123 cout <<
" Values at the point (phi_k, theta_j) : " << endl ;
124 cout.precision(precis) ;
125 cout.setf(ios::showpoint) ;
126 for (
int k=0; k<np; k++) {
127 cout <<
" k=" << k <<
" : " ;
128 for (
int j=0; j<nt; j++) {
129 double xx = (*vu[h])(nzm1, k, j, 0) ;
130 cout <<
" " << setw(precis) << xx ;
132 vmin = ( xx < vmin ) ? xx : vmin ;
133 vmax = ( xx > vmax ) ? xx : vmax ;
138 cout <<
"Minimum value on S^2 : " << vmin << endl ;
139 cout <<
"Maximum value on S^2 : " << vmax << endl ;
140 cout <<
"L^1 norm on S^2 : " << ndec << endl ;
144 for (
int k=0; k<np; k++) {
145 for (
int j=0; j<nt; j++) {
146 double v_inf = (*vu[h])(nzm1, k, j, 0) ;
147 for (
int i=0; i<nr; i++) {
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
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_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial 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.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tensor field of valence 0 (or component of a tensorial field).
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
void mult_r_ced()
Multiplication by r in the compactified external domain (CED), the dzpuis flag is not changed.
Valeur va
The numerical value of the Scalar
Valeur ** asymptot(int n, const int flag=0) const
Asymptotic expansion at r = infinity.
int dzpuis
Power of r by which the quantity represented by this must be divided in the compactified external d...
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Values and coefficients of a (real-value) function.
void set_etat_c_qcq()
Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl c ).
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
void set_base(const Base_val &)
Sets the bases for spectral expansions (member base )
Mtbl * c
Values of the function at the points of the multi-grid
Base_val base
Bases on which the spectral expansion is performed.
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).
const Map *const mp
Mapping on which the numerical values at the grid points are defined.