70 "Map_radial::reevaluate : the mapping mp_prev does not belong" 72 cout <<
" to the class Map_radial !" << endl ;
80 assert(uu.
get_mp() ==
this) ;
82 assert(uu.
get_etat() != ETATNONDEF) ;
83 assert(mp_prev->
mg ==
mg) ;
101 (uu.
va).set_etat_c_qcq() ;
105 Mtbl& va_c = *((uu.
va).c) ;
121 double precis = 1.e-15 ;
130 for (
int l=0; l<nzet; l++) {
138 double* ptx = (va_c.
t[l])->t ;
140 double* pr = (rc.
t[l])->t ;
144 for (
int k=0; k<np; k++) {
145 for (
int j=0; j<nt; j++) {
146 for (
int i=0; i<nr; i++) {
153 mp_prev->
val_lx_jk(*pr, j, k, par, l0, xi0) ;
183 if (mp_prev == 0x0) {
185 "Map_radial::reevaluate : the mapping mp_prev does not belong" 187 cout <<
" to the class Map_radial !" << endl ;
195 assert(uu.
get_mp() == *
this) ;
197 assert(uu.
get_etat() != ETATNONDEF) ;
198 assert(mp_prev->
mg ==
mg) ;
236 double precis = 1.e-15 ;
245 for (
int l=0; l<nzet; l++) {
253 double* ptx = (va_c.
t[l])->t ;
255 double* pr = (rc.
t[l])->t ;
259 for (
int k=0; k<np; k++) {
260 for (
int j=0; j<nt; j++) {
261 for (
int i=0; i<nr; i++) {
268 mp_prev->
val_lx_jk(*pr, j, k, par, l0, xi0) ;
const Map * get_mp() const
Returns the mapping.
Mtbl_cf * c_cf
Coefficients of the spectral expansion of the function.
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
void add_int(const int &n, int position=0)
Adds the address of a new int to the list.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
void coef() const
Computes the coeffcients of *this.
void annule(int l)
Sets the Cmp to zero in a given domain.
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
int get_etat() const
Returns the logical state.
Tensor field of valence 0 (or component of a tensorial field).
Base class for coordinate mappings.
Mtbl * c
The coordinate values at each grid point.
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
double val_point_jk(int l, double x, int j, int k) const
Computes the value of the field represented by *this at an arbitrary point in , but collocation point...
virtual void reevaluate(const Map *mp_prev, int nzet, Cmp &uu) const
Recomputes the values of a Cmp at the collocation points after a change in the mapping.
int get_dzpuis() const
Returns dzpuis.
Mtbl * c
Values of the function at the points of the multi-grid.
Base class for pure radial mappings.
int get_nzone() const
Returns the number of domains.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
void fait() const
Computes, at each point of the grid, the value of the coordinate or mapping derivative represented by...
virtual void val_lx_jk(double rr, int j, int k, const Param &par, int &l, double &xi) const =0
Computes the domain index l and the value of corresponding to a point of arbitrary r but collocation...
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
const Mg3d * mg
Pointer on the multi-grid Mgd3 on which this is defined.
int get_dzpuis() const
Returns dzpuis.
void add_double(const double &x, int position=0)
Adds the the address of a new double to the list.
Coefficients storage for the multi-domain spectral method.
void set_etat_c_qcq()
Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl c )...
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Valeur & set_spectral_va()
Returns va (read/write version)
Tbl ** t
Array (size nzone ) of pointers on the Tbl 's.
const Map & get_mp() const
Returns the mapping.
Valeur va
The numerical value of the Cmp.
void add_int_mod(int &n, int position=0)
Adds the address of a new modifiable int to the list.
Coord r
r coordinate centered on the grid