28char map_af_radius_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Map/map_af_radius.C,v 1.9 2014/10/13 08:53:03 j_novak Exp $" ;
99 assert( l<mg->get_nzone() ) ;
105 case FIN:
case RARE: {
111 resu = double(1) / (
alpha[l] * xi +
beta[l] ) ;
116 cout <<
"Map_af::val_r: unknown type_r ! " << endl ;
135 for (
int l=0; l<nz; l++) {
141 rmin = double(1)/rmin ;
142 rmax = double(1)/rmax ;
144 if ((rr - rmin >= -1.e-14*fabs(rmin)) && ( rr <= rmax )) {
152 cout.setf(ios::showpoint);
153 cout <<
"Map_af::val_lx: the domain containing r = " << rr <<
154 " has not been found ! "
156 for (
int l=0; l<nz; l++) {
158 if (
mg->
get_type_r(l) == UNSURR) rmin = double(1)/rmin ;
160 cout <<
"domain " << l <<
" : r_min = " << rmin ;
162 if (
mg->
get_type_r(l) == UNSURR) rmax = double(1)/rmax ;
163 cout <<
" : r_max = " << rmax << endl ;
173 case FIN:
case RARE: {
179 xi = ( double(1)/rr -
beta[lz] ) /
alpha[lz] ;
184 cout <<
"Map_af::val_lx: unknown type_r ! " << endl ;
193 int& lz,
double& xi)
const {
195 val_lx(rr, 0., 0., lz, xi) ;
207 return val_r(l, xi, 0., 0.) ;
216 int& l,
double& xi)
const {
218 val_lx(rr, 0., 0., par, l, xi) ;
double * beta
Array (size: mg->nzone ) of the values of in each domain.
double * alpha
Array (size: mg->nzone ) of the values of in each domain.
virtual void val_lx_jk(double rr, int j, int k, const Param &par, int &l, double &xi) const
Computes the domain index l and the value of corresponding to a point of arbitrary r but collocation...
virtual void val_lx(double rr, double theta, double pphi, int &l, double &xi) const
Computes the domain index l and the value of corresponding to a point given by its physical coordina...
virtual double val_r_jk(int l, double xi, int j, int k) const
Returns the value of the radial coordinate r for a given and a given collocation point in in a give...
virtual double val_r(int l, double xi, double theta, double pphi) const
Returns the value of the radial coordinate r for a given in a given domain.
const Mg3d * mg
Pointer on the multi-grid Mgd3 on which this is defined
int get_nzone() const
Returns the number of domains.
int get_type_r(int l) const
Returns the type of sampling in the radial direction in domain no.