29char cmp_asymptot_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Cmp/cmp_asymptot.C,v 1.5 2014/10/13 08:52:47 j_novak Exp $" ;
84 for (
int h=0; h<=n0; h++) {
94 for (
int h=0; h<
dzpuis; h++) {
102 for (
int h=
dzpuis; h<=n0; h++) {
108 for (
int l=0; l<nzm1; l++) {
113 for (
int k=0; k<np; k++) {
114 for (
int j=0; j<nt; j++) {
115 vu[h]->
set(nzm1, k, j, 0) = uu(nzm1, k, j, nrm1) ;
124 cout <<
"Term in 1/r^" << h << endl ;
125 cout <<
"-------------" << endl ;
128 double vmin = (*vu[h])(nzm1, 0, 0, 0) ;
131 cout <<
" Values at the point (phi_k, theta_j) : " << endl ;
132 cout.precision(precis) ;
133 cout.setf(ios::showpoint) ;
134 for (
int k=0; k<np; k++) {
135 cout <<
" k=" << k <<
" : " ;
136 for (
int j=0; j<nt; j++) {
137 double xx = (*vu[h])(nzm1, k, j, 0) ;
138 cout <<
" " << setw(precis) << xx ;
140 vmin = ( xx < vmin ) ? xx : vmin ;
141 vmax = ( xx > vmax ) ? xx : vmax ;
146 cout <<
"Minimum value on S^2 : " << vmin << endl ;
147 cout <<
"Maximum value on S^2 : " << vmax << endl ;
148 cout <<
"L^1 norm on S^2 : " << ndec << endl ;
152 for (
int k=0; k<np; k++) {
153 for (
int j=0; j<nt; j++) {
154 double v_inf = (*vu[h])(nzm1, k, j, 0) ;
155 for (
int i=0; i<nr; i++) {
156 uu.
set(nzm1, k, j, i) -= v_inf ;
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
const Map * mp
Reference mapping.
void mult_r_zec()
Multiplication by r in the external compactified domain (ZEC)
int dzpuis
Power of r by which the quantity represented by this must be divided in the external compactified z...
Valeur va
The numerical value of the Cmp
Tbl & set(int l)
Read/write of the value in a given domain.
Valeur ** asymptot(int n, const int flag=0) const
Asymptotic expansion at r = infinity.
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).
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
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).
Base_val base
Bases on which the spectral expansion is performed.