28char map_log_radius_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Map/map_log_radius.C,v 1.4 2014/10/13 08:53:06 j_novak Exp $" ;
64 assert( l<mg->get_nzone() ) ;
72 case FIN:
case RARE: {
78 resu = double(1) / (
alpha(l) * xi +
beta(l) ) ;
83 cout <<
"Map_log::val_r: unknown type_r ! " << endl ;
98 cout <<
"Map_log::val_r: unknown type_r ! " << endl ;
106 cout <<
"Map_log::val_r: unknown type_r ! " << endl ;
125 for (
int l=0; l<nz; l++) {
138 cout <<
"Case unknown in Map_log::val_lx" << endl ;
143 if (
mg->
get_type_r(l) == UNSURR) rmax = double(1)/rmax ;
153 cout.setf(ios::showpoint);
154 cout <<
"Map_log::val_lx: the domain containing r = " << rr <<
155 " has not been found ! "
166 case FIN:
case RARE: {
172 xi = ( double(1)/rr -
beta(lz) ) /
alpha(lz) ;
177 cout <<
"Map_log::val_lx: unknown type_r ! " << endl ;
190 cout <<
"Map_log::val_lx: unknown type_r ! " << endl ;
197 cout <<
"Map_log::val_lx: unknown type_r ! " << endl ;
205 int& lz,
double& xi)
const {
207 val_lx(rr, 0., 0., lz, xi) ;
219 return val_r(l, xi, 0., 0.) ;
228 int& l,
double& xi)
const {
230 val_lx(rr, 0., 0., par, l, xi) ;
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.
Itbl type_var
Array (size: mg->nzone ) of the type of variable in each domain.
Tbl alpha
Array (size: mg->nzone ) of the values of in each domain.
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 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 double val_r_jk(int l, double xi, int j, int k) const
< Comparison operator
Tbl beta
Array (size: mg->nzone ) of the values of in each 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.
Cmp exp(const Cmp &)
Exponential.
Cmp log(const Cmp &)
Neperian logarithm.