31char 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 $" ;
66 const Mg3d& mg = *(
mp->get_mg()) ;
76 for (
int h=0; h<=n0; h++) {
86 for (
int h=0; h<
dzpuis; h++) {
88 vu[h]->set_etat_zero() ;
94 for (
int h=
dzpuis; h<=n0; h++) {
98 vu[h]->set_etat_c_qcq() ;
99 vu[h]->c->set_etat_qcq() ;
101 vu[h]->c->t[
l]->set_etat_zero() ;
103 vu[h]->c->t[
nzm1]->set_etat_qcq() ;
105 for (
int k=0;
k<np;
k++) {
106 for (
int j=0;
j<nt;
j++) {
111 vu[h]->set_base(
uu.va.base ) ;
116 cout <<
"Term in 1/r^" << h <<
endl ;
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++) {
144 for (
int k=0;
k<np;
k++) {
145 for (
int j=0;
j<nt;
j++) {
147 for (
int i=0;
i<nr;
i++) {
Time evolution with partial storage (*** under development ***).
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.
Tensor field of valence 0 (or component of a tensorial field).
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...
Values and coefficients of a (real-value) function.
const Map *const mp
Mapping on which the numerical values at the grid points are defined.