67 assert( l<mg->get_nzone() ) ;
75 case FIN:
case RARE: {
81 resu = double(1) / (
alpha(l) * xi +
beta(l) ) ;
86 cout <<
"Map_log::val_r: unknown type_r ! " << endl ;
101 cout <<
"Map_log::val_r: unknown type_r ! " << endl ;
109 cout <<
"Map_log::val_r: unknown type_r ! " << endl ;
128 for (
int l=0; l<nz; l++) {
141 cout <<
"Case unknown in Map_log::val_lx" << endl ;
146 if (
mg->
get_type_r(l) == UNSURR) rmax =
double(1)/rmax ;
156 cout.setf(ios::showpoint);
157 cout <<
"Map_log::val_lx: the domain containing r = " << rr <<
158 " has not been found ! " 169 case FIN:
case RARE: {
175 xi = ( double(1)/rr -
beta(lz) ) /
alpha(lz) ;
180 cout <<
"Map_log::val_lx: unknown type_r ! " << endl ;
193 cout <<
"Map_log::val_lx: unknown type_r ! " << endl ;
200 cout <<
"Map_log::val_lx: unknown type_r ! " << endl ;
208 int& lz,
double& xi)
const {
210 val_lx(rr, 0., 0., lz, xi) ;
222 return val_r(l, xi, 0., 0.) ;
231 int& l,
double& xi)
const {
233 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.
Tbl beta
Array (size: mg->nzone ) of the values of in each domain.
Cmp log(const Cmp &)
Neperian logarithm.
Cmp exp(const Cmp &)
Exponential.
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...
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.
int get_nzone() const
Returns the number of domains.
const Mg3d * mg
Pointer on the multi-grid Mgd3 on which this is defined.
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...
int get_type_r(int l) const
Returns the type of sampling in the radial direction in domain no.
virtual double val_r_jk(int l, double xi, int j, int k) const
< Comparison operator