LORENE
map_star.h
1 
2 #ifndef __MAP_STAR_H_
3 #define __MAP_STAR_H_
4 
5 namespace Lorene {
6 
7  //------------------------------------//
8  // class Map_star //
9  //------------------------------------//
10 
11 /*
12  * Affine and starlike radial mapping to describe 3D star. \ingroup(map)
13  *
14  * The affine radial mapping is the simplest one between the grid coordinates
15  * \f$(\xi, \theta', \phi')\f$ and the physical coordinates \f$(r, \theta, \phi)\f$.
16  * It comprises a single domain (nucleus, shells to be added in the future)
17  * It is defined by \f$\theta=\theta'\f$, \f$\phi=\phi'\f$ and
18  * \li \f$r=\alpha \xi + \beta\f$, in non-compactified domains,
19  * where \f$\alpha\f$ and \f$\beta\f$ depend upon the angular direction.
20  *
21  *
22  */
23 class Map_star : public Map_radial{
24 
25 
26  // Data :
27  // ----
28  private:
31  Valeur beta ;
32 
33 
34  // Constructors, destructor :
35  // ------------------------
36  public:
46  Map_star(const Mg3d& mgrille, const double* r_limits) ;
56  Map_star(const Mg3d& mgrille, const Tbl& r_limits) ;
57 
58  Map_star(const Map_star& ) ;
59  Map_star(const Mg3d&, const string&) ;
60  Map_star(const Mg3d&, FILE* ) ;
61 
62 
63  virtual ~Map_star() ;
64 
68  virtual const Map_af& mp_angu(int) const ;
69 
70  // Assignment
71  // ----------
72  public:
74  virtual void operator=(const Map_af& ) ;
76  virtual void operator=(const Map_star& ) ;
77 
78  // Memory management
79  // -----------------
80  private:
82  void set_coord() ;
83  //protected:
84  //virtual void reset_coord() ; ///< Resets all the member \c Coords
85 
86  // Extraction of information
87  // -------------------------
88  public:
90  const Valeur& get_alpha() const ;
91  const Valeur& get_beta() const ;
92 
93 
95  void set_alpha(const Tbl& alpha0, int l) ;
96  void set_beta(const Tbl& beta0, int l) ;
97 
98  void set_alpha(const Valeur& alpha0) ;
99  void set_beta(const Valeur& beta0) ;
100 
110  virtual double val_r(int l, double xi, double theta, double pphi) const ;
111 
121  virtual void val_lx(double rr, double theta, double pphi,
122  int& l, double& xi) const ;
123 
133  virtual void val_lx(double rr, double theta, double pphi,
134  const Param& par, int& l, double& xi) const ;
135 
145  virtual double val_r_jk(int l, double xi, int j, int k) const ;
146 
156  virtual void val_lx_jk(double rr, int j, int k, const Param& par,
157  int& l, double& xi) const ;
158 
160  virtual bool operator==(const Map& ) const ;
161 
168  virtual void dsdr (const Scalar& ci, Scalar& resu) const ;
169 
177  virtual void srstdsdp (const Scalar&, Scalar&) const ;
178 
186  virtual void srdsdt (const Scalar&, Scalar&) const ;
187 
192  virtual void dsdt (const Scalar&, Scalar&) const ;
193 
194 
199  virtual void stdsdp (const Scalar&, Scalar&) const ;
200 
206  Scalar interpolate_from_map_af(const Scalar& f_a) const ;
207 
208  // Outputs
209  // -------
210  public:
211  virtual void sauve(FILE* ) const ;
212 
213  private:
214  virtual ostream& operator>>(ostream &) const ;
215 
216  virtual void homothetie (double) ;
217  virtual void resize (int, double) ;
218  virtual void adapt (const Cmp&, const Param&, int) ;
219  virtual void dsdr (const Cmp&, Cmp&) const ;
220  virtual void dsdxi (const Cmp&, Cmp&) const ;
221  virtual void dsdxi (const Scalar&, Scalar&) const ;
222  virtual void dsdradial (const Scalar& uu, Scalar& resu) const ;
223  virtual void srdsdt (const Cmp&, Cmp&) const ;
224  virtual void srstdsdp (const Cmp&, Cmp&) const ;
225  virtual void laplacien (const Scalar&, int, Scalar&) const ;
226  virtual void laplacien (const Cmp&, int, Cmp&) const ;
227  virtual void lapang (const Scalar&, Scalar&) const ;
228  virtual void primr(const Scalar&, Scalar&, bool) const ;
229  virtual Tbl* integrale (const Scalar&) const ;
230  virtual Tbl* integrale (const Cmp&) const ;
231  virtual void poisson (const Cmp&, Param&, Cmp&) const ;
232  virtual void poisson_tau (const Cmp&, Param&, Cmp&) const ;
233  virtual void poisson (const Scalar&, Param&, Scalar&) const ;
234  virtual void poisson_tau (const Scalar&, Param&, Scalar&) const ;
235  virtual void poisson_falloff(const Cmp&, Param&, Cmp&, int) const ;
236  virtual void poisson_ylm(const Cmp&, Param&, Cmp&, int, double*) const ;
237  virtual void poisson_regular (const Cmp&, int, int, double, Param&, Cmp&, Cmp&, Cmp&,
238  Tenseur&, Cmp&, Cmp&) const ;
239  virtual void poisson_angu (const Scalar&, Param&, Scalar&, double=0) const ;
240  virtual void poisson_angu (const Cmp&, Param&, Cmp&, double=0) const ;
241  virtual Param* donne_para_poisson_vect (Param&, int) const ;
242  virtual void poisson_frontiere (const Cmp&, const Valeur&, int, int, Cmp&, double = 0., double = 0.) const ;
243  virtual void poisson_frontiere_double (const Cmp&, const Valeur&, const Valeur&, int, Cmp&) const ;
244  virtual void poisson_interne (const Cmp&, const Valeur&, Param&, Cmp&) const ;
245  virtual void poisson2d (const Cmp&, const Cmp&, Param&, Cmp&) const ;
246  virtual void dalembert (Param&, Scalar&, const Scalar&, const Scalar&, const Scalar&) const ;
247  // Building functions for the Coord's
248  // ----------------------------------
249  friend Mtbl* map_star_fait_r(const Map* ) ;
250  friend Mtbl* map_star_fait_tet(const Map* ) ;
251  friend Mtbl* map_star_fait_phi(const Map* ) ;
252  friend Mtbl* map_star_fait_sint(const Map* ) ;
253  friend Mtbl* map_star_fait_cost(const Map* ) ;
254  friend Mtbl* map_star_fait_sinp(const Map* ) ;
255  friend Mtbl* map_star_fait_cosp(const Map* ) ;
256 
257  friend Mtbl* map_star_fait_x(const Map* ) ;
258  friend Mtbl* map_star_fait_y(const Map* ) ;
259  friend Mtbl* map_star_fait_z(const Map* ) ;
260 
261  friend Mtbl* map_star_fait_xa(const Map* ) ;
262  friend Mtbl* map_star_fait_ya(const Map* ) ;
263  friend Mtbl* map_star_fait_za(const Map* ) ;
264 
265  friend Mtbl* map_star_fait_xsr(const Map* ) ;
266  friend Mtbl* map_star_fait_dxdr(const Map* ) ;
267  friend Mtbl* map_star_fait_drdt(const Map* ) ;
268  friend Mtbl* map_star_fait_stdrdp(const Map* ) ;
269  friend Mtbl* map_star_fait_srdrdt(const Map* ) ;
270  friend Mtbl* map_star_fait_srstdrdp(const Map* ) ;
271  friend Mtbl* map_star_fait_sr2drdt(const Map* ) ;
272  friend Mtbl* map_star_fait_sr2stdrdp(const Map* ) ;
273  friend Mtbl* map_star_fait_d2rdx2(const Map* ) ;
274  friend Mtbl* map_star_fait_lapr_tp(const Map* ) ;
275  friend Mtbl* map_star_fait_d2rdtdx(const Map* ) ;
276  friend Mtbl* map_star_fait_sstd2rdpdx(const Map* ) ;
277  friend Mtbl* map_star_fait_sr2d2rdt2(const Map* ) ;
278 };
279  Mtbl* map_star_fait_r(const Map* ) ;
280  Mtbl* map_star_fait_tet(const Map* ) ;
281  Mtbl* map_star_fait_phi(const Map* ) ;
282  Mtbl* map_star_fait_sint(const Map* ) ;
283  Mtbl* map_star_fait_cost(const Map* ) ;
284  Mtbl* map_star_fait_sinp(const Map* ) ;
285  Mtbl* map_star_fait_cosp(const Map* ) ;
286 
287  Mtbl* map_star_fait_x(const Map* ) ;
288  Mtbl* map_star_fait_y(const Map* ) ;
289  Mtbl* map_star_fait_z(const Map* ) ;
290 
291  Mtbl* map_star_fait_xa(const Map* ) ;
292  Mtbl* map_star_fait_ya(const Map* ) ;
293  Mtbl* map_star_fait_za(const Map* ) ;
294 
295  Mtbl* map_star_fait_xsr(const Map* ) ;
296  Mtbl* map_star_fait_dxdr(const Map* ) ;
297  Mtbl* map_star_fait_drdt(const Map* ) ;
298  Mtbl* map_star_fait_stdrdp(const Map* ) ;
299  Mtbl* map_star_fait_srdrdt(const Map* ) ;
300  Mtbl* map_star_fait_srstdrdp(const Map* ) ;
301  Mtbl* map_star_fait_sr2drdt(const Map* ) ;
302  Mtbl* map_star_fait_sr2stdrdp(const Map* ) ;
303  Mtbl* map_star_fait_d2rdx2(const Map* ) ;
304  Mtbl* map_star_fait_lapr_tp(const Map* ) ;
305  Mtbl* map_star_fait_d2rdtdx(const Map* ) ;
306  Mtbl* map_star_fait_sstd2rdpdx(const Map* ) ;
307  Mtbl* map_star_fait_sr2d2rdt2(const Map* ) ;
308 
309 
310  //------------------------------------//
311  // class Map_eps //
312  //------------------------------------//
313 
314 /*
315  * Affine and starlike radial mapping to describe 3D star. \ingroup(map)
316  *
317  * The affine radial mapping is the simplest one between the grid coordinates
318  * \f$(\xi, \theta', \phi')\f$ and the physical coordinates \f$(r, \theta, \phi)\f$. It comprises a single domain (nucleus, shells to be added in the future)
319  * It is defined by \f$\theta=\theta'\f$, \f$\phi=\phi'\f$ and
320  * \li \f$r=\alpha \xi + \beta\f$, in non-compactified domains,
321  * where \f$\alpha\f$ and \f$\beta\f$ depend upon the angular direction.
322  *
323  *
324  */
325 class Map_eps : public Map_radial{
326 
327 
328  // Data :
329  // ----
330  private:
332  double aa, bb, cc ;
333  Valeur alpha ;
334  Valeur beta ;
335 
336 
337  // Constructors, destructor :
338  // ------------------------
339  public:
349  Map_eps(const Mg3d& mgrille, const double* r_limits) ;
359  Map_eps(const Mg3d& mgrille, double a, double b, double c) ;
360 
361  Map_eps(const Map_eps& ) ;
362  Map_eps(const Mg3d&, const string&) ;
363  Map_eps(const Mg3d&, FILE* ) ;
364 
365 
366  virtual ~Map_eps() ;
367 
371  virtual const Map_af& mp_angu(int) const ;
372 
373  // Assignment
374  // ----------
375  public:
377  virtual void operator=(const Map_af& ) ;
379  virtual void operator=(const Map_eps& ) ;
380 
381  // Memory management
382  // -----------------
383  private:
385  void set_coord() ;
386  //protected:
387  //virtual void reset_coord() ; ///< Resets all the member \c Coords
388 
389  // Extraction of information
390  // -------------------------
391  public:
393  const Valeur& get_alpha() const ;
394  const Valeur& get_beta() const ;
395  double get_aa() const { return aa;} ;
396  double get_bb() const { return bb;} ;
397  double get_cc() const { return cc;} ;
398 
399 
401  void set_alpha(const Tbl& alpha0, int l) ;
402  void set_beta(const Tbl& beta0, int l) ;
403 
404  void set_alpha(const Valeur& alpha0) ;
405  void set_beta(const Valeur& beta0) ;
406 
416  virtual double val_r(int l, double xi, double theta, double pphi) const ;
417 
427  virtual void val_lx(double rr, double theta, double pphi,
428  int& l, double& xi) const ;
429 
439  virtual void val_lx(double rr, double theta, double pphi,
440  const Param& par, int& l, double& xi) const ;
441 
451  virtual double val_r_jk(int l, double xi, int j, int k) const ;
452 
462  virtual void val_lx_jk(double rr, int j, int k, const Param& par,
463  int& l, double& xi) const ;
464 
466  virtual bool operator==(const Map& ) const ;
467 
474  virtual void dsdr (const Scalar& ci, Scalar& resu) const ;
475 
483  virtual void srstdsdp (const Scalar&, Scalar&) const ;
484 
492  virtual void srdsdt (const Scalar&, Scalar&) const ;
493 
498  virtual void dsdt (const Scalar&, Scalar&) const ;
499 
500 
505  virtual void stdsdp (const Scalar&, Scalar&) const ;
506 
507  // Outputs
508  // -------
509  public:
510  virtual void sauve(FILE* ) const ;
511 
512  private:
513  virtual ostream& operator>>(ostream &) const ;
514 
515  virtual void homothetie (double) ;
516  virtual void resize (int, double) ;
517  virtual void adapt (const Cmp&, const Param&, int) ;
518  virtual void dsdr (const Cmp&, Cmp&) const ;
519  virtual void dsdxi (const Cmp&, Cmp&) const ;
520  virtual void dsdxi (const Scalar&, Scalar&) const ;
521  virtual void dsdradial (const Scalar& uu, Scalar& resu) const ;
522  virtual void srdsdt (const Cmp&, Cmp&) const ;
523  virtual void srstdsdp (const Cmp&, Cmp&) const ;
524  virtual void laplacien (const Scalar&, int, Scalar&) const ;
525  virtual void laplacien (const Cmp&, int, Cmp&) const ;
526  virtual void lapang (const Scalar&, Scalar&) const ;
527  virtual void primr(const Scalar&, Scalar&, bool) const ;
528  virtual Tbl* integrale (const Scalar&) const ;
529  virtual Tbl* integrale (const Cmp&) const ;
530  virtual void poisson (const Cmp&, Param&, Cmp&) const ;
531  virtual void poisson_tau (const Cmp&, Param&, Cmp&) const ;
532  virtual void poisson (const Scalar&, Param&, Scalar&) const ;
533  virtual void poisson_tau (const Scalar&, Param&, Scalar&) const ;
534  virtual void poisson_falloff(const Cmp&, Param&, Cmp&, int) const ;
535  virtual void poisson_ylm(const Cmp&, Param&, Cmp&, int, double*) const ;
536  virtual void poisson_regular (const Cmp&, int, int, double, Param&, Cmp&, Cmp&, Cmp&,
537  Tenseur&, Cmp&, Cmp&) const ;
538  virtual void poisson_angu (const Scalar&, Param&, Scalar&, double=0) const ;
539  virtual void poisson_angu (const Cmp&, Param&, Cmp&, double=0) const ;
540  virtual Param* donne_para_poisson_vect (Param&, int) const ;
541  virtual void poisson_frontiere (const Cmp&, const Valeur&, int, int, Cmp&, double = 0., double = 0.) const ;
542  virtual void poisson_frontiere_double (const Cmp&, const Valeur&, const Valeur&, int, Cmp&) const ;
543  virtual void poisson_interne (const Cmp&, const Valeur&, Param&, Cmp&) const ;
544  virtual void poisson2d (const Cmp&, const Cmp&, Param&, Cmp&) const ;
545  virtual void dalembert (Param&, Scalar&, const Scalar&, const Scalar&, const Scalar&) const ;
546  // Building functions for the Coord's
547  // ----------------------------------
548  friend Mtbl* map_eps_fait_r(const Map* ) ;
549  friend Mtbl* map_eps_fait_tet(const Map* ) ;
550  friend Mtbl* map_eps_fait_phi(const Map* ) ;
551  friend Mtbl* map_eps_fait_sint(const Map* ) ;
552  friend Mtbl* map_eps_fait_cost(const Map* ) ;
553  friend Mtbl* map_eps_fait_sinp(const Map* ) ;
554  friend Mtbl* map_eps_fait_cosp(const Map* ) ;
555 
556  friend Mtbl* map_eps_fait_x(const Map* ) ;
557  friend Mtbl* map_eps_fait_y(const Map* ) ;
558  friend Mtbl* map_eps_fait_z(const Map* ) ;
559 
560  friend Mtbl* map_eps_fait_xa(const Map* ) ;
561  friend Mtbl* map_eps_fait_ya(const Map* ) ;
562  friend Mtbl* map_eps_fait_za(const Map* ) ;
563 
564  friend Mtbl* map_eps_fait_xsr(const Map* ) ;
565  friend Mtbl* map_eps_fait_dxdr(const Map* ) ;
566  friend Mtbl* map_eps_fait_drdt(const Map* ) ;
567  friend Mtbl* map_eps_fait_stdrdp(const Map* ) ;
568  friend Mtbl* map_eps_fait_srdrdt(const Map* ) ;
569  friend Mtbl* map_eps_fait_srstdrdp(const Map* ) ;
570  friend Mtbl* map_eps_fait_sr2drdt(const Map* ) ;
571  friend Mtbl* map_eps_fait_sr2stdrdp(const Map* ) ;
572  friend Mtbl* map_eps_fait_d2rdx2(const Map* ) ;
573  friend Mtbl* map_eps_fait_lapr_tp(const Map* ) ;
574  friend Mtbl* map_eps_fait_d2rdtdx(const Map* ) ;
575  friend Mtbl* map_eps_fait_sstd2rdpdx(const Map* ) ;
576  friend Mtbl* map_eps_fait_sr2d2rdt2(const Map* ) ;
577 };
578  Mtbl* map_eps_fait_r(const Map* ) ;
579  Mtbl* map_eps_fait_tet(const Map* ) ;
580  Mtbl* map_eps_fait_phi(const Map* ) ;
581  Mtbl* map_eps_fait_sint(const Map* ) ;
582  Mtbl* map_eps_fait_cost(const Map* ) ;
583  Mtbl* map_eps_fait_sinp(const Map* ) ;
584  Mtbl* map_eps_fait_cosp(const Map* ) ;
585 
586  Mtbl* map_eps_fait_x(const Map* ) ;
587  Mtbl* map_eps_fait_y(const Map* ) ;
588  Mtbl* map_eps_fait_z(const Map* ) ;
589 
590  Mtbl* map_eps_fait_xa(const Map* ) ;
591  Mtbl* map_eps_fait_ya(const Map* ) ;
592  Mtbl* map_eps_fait_za(const Map* ) ;
593 
594  Mtbl* map_eps_fait_xsr(const Map* ) ;
595  Mtbl* map_eps_fait_dxdr(const Map* ) ;
596  Mtbl* map_eps_fait_drdt(const Map* ) ;
597  Mtbl* map_eps_fait_stdrdp(const Map* ) ;
598  Mtbl* map_eps_fait_srdrdt(const Map* ) ;
599  Mtbl* map_eps_fait_srstdrdp(const Map* ) ;
600  Mtbl* map_eps_fait_sr2drdt(const Map* ) ;
601  Mtbl* map_eps_fait_sr2stdrdp(const Map* ) ;
602  Mtbl* map_eps_fait_d2rdx2(const Map* ) ;
603  Mtbl* map_eps_fait_lapr_tp(const Map* ) ;
604  Mtbl* map_eps_fait_d2rdtdx(const Map* ) ;
605  Mtbl* map_eps_fait_sstd2rdpdx(const Map* ) ;
606  Mtbl* map_eps_fait_sr2d2rdt2(const Map* ) ;
607 
608 }
609 #endif
virtual ostream & operator>>(ostream &) const
Operator >>
Definition: map_eps.C:267
virtual void dsdradial(const Scalar &uu, Scalar &resu) const
< Not implemented
virtual Param * donne_para_poisson_vect(Param &, int) const
< Not implemented
virtual const Map_af & mp_angu(int) const
Returns the "angular" mapping for the outside of domain l_zone.
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...
friend Mtbl * map_eps_fait_r(const Map *)
< Not implemented
Definition: map_eps_fait.C:43
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Definition: cmp.h:446
virtual bool operator==(const Map &) const
Comparison operator (egality)
Definition: map_star.C:201
virtual void poisson_tau(const Cmp &, Param &, Cmp &) const
< Not implemented
virtual const Map_af & mp_angu(int) const
Returns the "angular" mapping for the outside of domain l_zone.
Multi-domain array.
Definition: mtbl.h:118
Lorene prototypes.
Definition: app_hor.h:67
virtual void resize(int, double)
< Not implemented
virtual void poisson(const Cmp &, Param &, Cmp &) const
< Not implemented
Tensor field of valence 0 (or component of a tensorial field).
Definition: scalar.h:393
Base class for coordinate mappings.
Definition: map.h:694
virtual void homothetie(double)
Sets a new radial scale.
friend Mtbl * map_star_fait_r(const Map *)
< Not implemented
Definition: map_star_fait.C:43
virtual void operator=(const Map_af &)
Assignment to an affine mapping.
Definition: map_eps.C:161
virtual Tbl * integrale(const Scalar &) const
< Not implemented
virtual void laplacien(const Scalar &, int, Scalar &) const
< Not implemented
virtual void poisson_falloff(const Cmp &, Param &, Cmp &, int) const
< Not implemented
virtual void adapt(const Cmp &, const Param &, int)
< Not implemented
virtual void poisson_regular(const Cmp &, int, int, double, Param &, Cmp &, Cmp &, Cmp &, Tenseur &, Cmp &, Cmp &) const
< Not implemented
Values and coefficients of a (real-value) function.
Definition: valeur.h:297
void set_alpha(const Tbl &alpha0, int l)
Modifies the value of in domain no. l.
Definition: map_eps.C:314
virtual void srdsdt(const Scalar &, Scalar &) const
Computes of a Scalar.
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 void lapang(const Scalar &, Scalar &) const
< Not implemented
virtual void dalembert(Param &, Scalar &, const Scalar &, const Scalar &, const Scalar &) const
< Not implemented
Map_eps(const Mg3d &mgrille, const double *r_limits)
Standard Constructor.
Definition: map_eps.C:52
virtual bool operator==(const Map &) const
Comparison operator (egality)
Definition: map_eps.C:207
virtual void poisson(const Cmp &, Param &, Cmp &) const
< Not implemented
virtual void poisson_frontiere(const Cmp &, const Valeur &, int, int, Cmp &, double=0., double=0.) const
< Not implemented
virtual void operator=(const Map_af &)
Assignment to an affine mapping.
Definition: map_star.C:155
virtual void poisson_frontiere(const Cmp &, const Valeur &, int, int, Cmp &, double=0., double=0.) const
< Not implemented
virtual void poisson_tau(const Cmp &, Param &, Cmp &) const
< Not implemented
virtual void laplacien(const Scalar &, int, Scalar &) const
< Not implemented
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 poisson_interne(const Cmp &, const Valeur &, Param &, Cmp &) const
< Not implemented
virtual Tbl * integrale(const Scalar &) const
< Not implemented
void set_alpha(const Tbl &alpha0, int l)
Modifies the value of in domain no. l.
Definition: map_star.C:307
virtual void adapt(const Cmp &, const Param &, int)
< Not implemented
virtual void poisson_falloff(const Cmp &, Param &, Cmp &, int) const
< Not implemented
virtual void stdsdp(const Scalar &, Scalar &) const
Computes of a Scalar .
const Valeur & get_alpha() const
Returns the reference on the Tbl alpha.
Definition: map_eps.C:243
virtual void homothetie(double)
Sets a new radial scale.
Parameter storage.
Definition: param.h:125
Base class for pure radial mappings.
Definition: map.h:1596
Map_star(const Mg3d &mgrille, const double *r_limits)
Standard Constructor.
Definition: map_star.C:52
virtual void poisson2d(const Cmp &, const Cmp &, Param &, Cmp &) const
< Not implemented
void set_coord()
Assignment of the building functions to the member Coords.
Definition: map_eps.C:171
virtual void poisson_angu(const Scalar &, Param &, Scalar &, double=0) const
< Not implemented
virtual void poisson_ylm(const Cmp &, Param &, Cmp &, int, double *) const
< Not implemented
const Valeur & get_alpha() const
Returns the reference on the Tbl alpha.
Definition: map_star.C:237
virtual void sauve(FILE *) const
Save in a file.
Definition: map_eps.C:254
virtual void poisson_regular(const Cmp &, int, int, double, Param &, Cmp &, Cmp &, Cmp &, Tenseur &, Cmp &, Cmp &) const
< Not implemented
virtual void dsdxi(const Cmp &, Cmp &) const
< Not implemented
virtual ~Map_eps()
Destructor.
Definition: map_eps.C:134
void set_coord()
Assignment of the building functions to the member Coords.
Definition: map_star.C:165
virtual void primr(const Scalar &, Scalar &, bool) const
< Not implemented
virtual void poisson_angu(const Scalar &, Param &, Scalar &, double=0) const
< Not implemented
Multi-domain grid.
Definition: grilles.h:279
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.
virtual void srdsdt(const Scalar &, Scalar &) const
Computes of a Scalar.
Affine radial mapping.
Definition: map.h:2087
virtual void poisson2d(const Cmp &, const Cmp &, Param &, Cmp &) const
< Not implemented
virtual void dsdradial(const Scalar &uu, Scalar &resu) const
< Not implemented
virtual void stdsdp(const Scalar &, Scalar &) const
Computes of a Scalar .
Definition: map_eps_deriv.C:95
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 poisson_ylm(const Cmp &, Param &, Cmp &, int, double *) const
< Not implemented
virtual void poisson_frontiere_double(const Cmp &, const Valeur &, const Valeur &, int, Cmp &) const
< Not implemented
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 void srstdsdp(const Scalar &, Scalar &) const
Computes of a Scalar.
virtual void resize(int, double)
< Not implemented
virtual void srstdsdp(const Scalar &, Scalar &) const
Computes of a Scalar.
virtual void primr(const Scalar &, Scalar &, bool) const
< Not implemented
virtual void poisson_interne(const Cmp &, const Valeur &, Param &, Cmp &) const
< Not implemented
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.
virtual void dsdr(const Scalar &ci, Scalar &resu) const
Computes of a Scalar.
Definition: map_eps_deriv.C:39
Basic array class.
Definition: tbl.h:164
virtual ostream & operator>>(ostream &) const
Operator >>
Definition: map_star.C:261
Map(const Mg3d &)
Constructor from a multi-domain 3D grid.
Definition: map.C:142
virtual void dalembert(Param &, Scalar &, const Scalar &, const Scalar &, const Scalar &) const
< Not implemented
Scalar interpolate_from_map_af(const Scalar &f_a) const
Interpolates from a Scalar defined on a Map_af and returns the new Scalar defined on *this...
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 ~Map_star()
Destructor.
Definition: map_star.C:128
double aa
Array (size: mg->nzone*Nt*Np ) of the values of in each domain.
Definition: map_star.h:332
virtual void poisson_frontiere_double(const Cmp &, const Valeur &, const Valeur &, int, Cmp &) const
< Not implemented
virtual void dsdr(const Scalar &ci, Scalar &resu) const
Computes of a Scalar.
virtual Param * donne_para_poisson_vect(Param &, int) const
< Not implemented
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Definition: tenseur.h:304
virtual void dsdt(const Scalar &, Scalar &) const
Computes of a Scalar .
Valeur alpha
Array (size: mg->nzone*Nt*Np ) of the values of in each domain.
Definition: map_star.h:30
virtual void dsdxi(const Cmp &, Cmp &) const
< Not implemented
virtual void lapang(const Scalar &, Scalar &) const
< Not implemented
virtual void sauve(FILE *) const
Save in a file.
Definition: map_star.C:248
virtual void dsdt(const Scalar &, Scalar &) const
Computes of a Scalar .