KADATH
spheric_time.hpp
1 /*
2  Copyright 2017 Philippe Grandclement
3 
4  This file is part of Kadath.
5 
6  Kadath is free software: you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation, either version 3 of the License, or
9  (at your option) any later version.
10 
11  Kadath is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with Kadath. If not, see <http://www.gnu.org/licenses/>.
18 */
19 
20 #ifndef __SPHERIC_TIME_HPP_
21 #define __SPHERIC_TIME_HPP_
22 
23 #include "space.hpp"
24 
25 namespace Kadath {
26 
45 class Domain_spheric_time_nucleus : public Domain {
46 
47  private:
48  double alpha ;
49  double tmin;
50  double tmax ;
51 
52  public:
62  Domain_spheric_time_nucleus (int num, int ttype, double tmmin, double tmmax, double radius, const Dim_array& nbr) ;
63  Domain_spheric_time_nucleus (const Domain_spheric_time_nucleus& so) ;
69  Domain_spheric_time_nucleus (int num, FILE* ff) ;
70 
71  virtual ~Domain_spheric_time_nucleus() ;
72  virtual void save (FILE*) const ;
73 
74  private:
75  virtual void do_absol () const ;
76 
77  private:
78  virtual void set_cheb_base(Base_spectral&) const ;
79  virtual void set_legendre_base(Base_spectral&) const ;
80  virtual void do_coloc () ;
81  virtual void do_radius () const ;
82 
83  public:
84  virtual bool is_in(const Point&xx, double prec=1e-13) const ;
85  virtual const Point absol_to_num(const Point&) const;
86  virtual void do_der_abs_from_der_var(const Val_domain *const *const der_var, Val_domain **const der_abs) const ;
87  virtual Base_spectral mult (const Base_spectral&, const Base_spectral&) const ;
91  double get_alpha () const {return alpha ;} ;
92 
93  public:
94  virtual Val_domain div_x (const Val_domain&) const ;
95  virtual Val_domain mult_r (const Val_domain&) const ;
96  virtual Val_domain div_r (const Val_domain&) const ;
97  virtual Val_domain laplacian (const Val_domain&, int) const ;
98  virtual Val_domain der_r (const Val_domain&) const ;
99  virtual Val_domain srdr (const Val_domain&) const ;
100  virtual Val_domain ddt (const Val_domain&) const ;
101  virtual Val_domain dt (const Val_domain&) const ;
102 
103  virtual double val_boundary (int, const Val_domain&, const Index&) const ;
104  virtual void find_other_dom (int, int, int&, int&) const ;
105  virtual Val_domain der_normal (const Val_domain&, int) const ;
106 
107  virtual int nbr_unknowns (const Tensor&, int) const ;
114  int nbr_unknowns_val_domain (const Val_domain& so) const ;
115  virtual Array<int> nbr_conditions_array (const Tensor&, int, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
123  int nbr_conditions_val_domain_array (const Val_domain& so, const Array<int>& order) const ;
124  virtual Array<int> nbr_conditions_boundary_array (const Tensor&, int, int, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
134  int nbr_conditions_val_domain_boundary_array (const Val_domain& eq, int bound, const Array<int>& order) const ;
135  virtual void export_tau_array (const Tensor&, int, const Array<int>&, Array<double>&, int&, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
145  void export_tau_val_domain_array (const Val_domain& eq, const Array<int>& order, Array<double>& res, int& pos_res, int ncond) const ;
146  virtual void export_tau_boundary_array (const Tensor&, int, int, const Array<int>&, Array<double>&, int&, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
157  void export_tau_val_domain_boundary_array (const Val_domain& eq, int bound, const Array<int>& order, Array<double>& res, int& pos_res, int ncond) const ;
158  virtual void affecte_tau (Tensor&, int, const Array<double>&, int&) const ;
166  void affecte_tau_val_domain (Val_domain& so, const Array<double>& cf, int& pos_cf) const ;
167  virtual void affecte_tau_one_coef (Tensor&, int, int, int&) const ;
175  void affecte_tau_one_coef_val_domain (Val_domain& so, int cc, int& pos_cf) const ;
176 
177 public:
178  virtual ostream& print (ostream& o) const ;
179 } ;
180 
199 class Domain_spheric_time_shell : public Domain {
200 
201  private:
202  double alpha ;
203  double beta ;
204  double tmin ;
205  double tmax ;
206 
207  public:
218  Domain_spheric_time_shell (int num, int ttype, double tmmin, double tmmax , double r_int, double r_ext, const Dim_array& nbr) ;
219  Domain_spheric_time_shell (const Domain_spheric_time_shell& so) ;
225  Domain_spheric_time_shell (int num, FILE* ff) ;
226 
227  virtual ~Domain_spheric_time_shell() ;
228  virtual void save (FILE*) const ;
229 
230  private:
231  virtual void do_absol () const ;
232 
233  private :
234  virtual void set_cheb_base(Base_spectral&) const ;
235  virtual void set_legendre_base(Base_spectral&) const ;
236  virtual void do_coloc () ;
237  virtual void do_radius () const ;
238 
239  public:
240  virtual bool is_in(const Point& xx, double prec=1e-13) const ;
241  virtual const Point absol_to_num(const Point&) const;
242  virtual void do_der_abs_from_der_var(const Val_domain *const *const der_var, Val_domain **const der_abs) const ;
243 
244  virtual Base_spectral mult (const Base_spectral&, const Base_spectral&) const ;
248  double get_alpha () const {return alpha ;} ;
249 
250  public:
251  virtual Val_domain mult_r (const Val_domain&) const ;
252  virtual Val_domain div_r (const Val_domain&) const ;
253  virtual Val_domain laplacian (const Val_domain&, int) const ;
254  virtual Val_domain der_r (const Val_domain&) const ;
255  virtual Val_domain ddt (const Val_domain&) const ;
256  virtual Val_domain dt (const Val_domain&) const ;
257 
258  virtual void find_other_dom (int, int, int&, int&) const ;
259  virtual Val_domain der_normal (const Val_domain&, int) const ;
260  virtual double val_boundary (int, const Val_domain&, const Index&) const ;
261 
262  virtual int nbr_unknowns (const Tensor&, int) const ;
269  int nbr_unknowns_val_domain (const Val_domain& so) const ;
270  virtual Array<int> nbr_conditions_array (const Tensor&, int, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
278  int nbr_conditions_val_domain_array (const Val_domain& so, const Array<int>& order) const ;
279  virtual Array<int> nbr_conditions_boundary_array (const Tensor&, int, int, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
289  int nbr_conditions_val_domain_boundary_array (const Val_domain& eq, int bound, const Array<int>& order) const ;
290  virtual void export_tau_array (const Tensor&, int, const Array<int>&, Array<double>&, int&, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
300  void export_tau_val_domain_array (const Val_domain& eq, const Array<int>& order, Array<double>& res, int& pos_res, int ncond) const ;
301  virtual void export_tau_boundary_array (const Tensor&, int, int, const Array<int>&, Array<double>&, int&, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
312  void export_tau_val_domain_boundary_array (const Val_domain& eq, int bound, const Array<int>& order, Array<double>& res, int& pos_res, int ncond) const ;
313  virtual void affecte_tau (Tensor&, int, const Array<double>&, int&) const ;
321  void affecte_tau_val_domain (Val_domain& so, const Array<double>& cf, int& pos_cf) const ;
322  virtual void affecte_tau_one_coef (Tensor&, int, int, int&) const ;
330  void affecte_tau_one_coef_val_domain (Val_domain& so, int cc, int& pos_cf) const ;
331 
332 public:
333  virtual ostream& print (ostream& o) const ;
334 } ;
335 
354 class Domain_spheric_time_compact : public Domain {
355 
356  private:
357  double alpha ;
358  double tmin ;
359  double tmax ;
360 
361  public:
371  Domain_spheric_time_compact (int num, int ttype, double tmmin, double tmmax , double r_int, const Dim_array& nbr) ;
372  Domain_spheric_time_compact (const Domain_spheric_time_compact& so) ;
378  Domain_spheric_time_compact (int num, FILE* ff) ;
379 
380  virtual ~Domain_spheric_time_compact() ;
381  virtual void save (FILE*) const ;
382 
383  private:
384  virtual void do_absol () const ;
385 
386  private :
387  virtual void set_cheb_base(Base_spectral&) const ;
388  virtual void set_legendre_base(Base_spectral&) const ;
389  virtual void do_coloc () ;
390  virtual void do_radius () const ;
391 
392  public:
393  virtual bool is_in(const Point& xx, double prec=1e-13) const ;
394  virtual const Point absol_to_num(const Point&) const;
395  virtual void do_der_abs_from_der_var(const Val_domain *const *const der_var, Val_domain **const der_abs) const ;
396 
397  virtual Base_spectral mult (const Base_spectral&, const Base_spectral&) const ;
401  double get_alpha () const {return alpha ;} ;
402 
403  public:
404  virtual Val_domain mult_xm1 (const Val_domain&) const ;
405  virtual Val_domain div_xm1 (const Val_domain&) const ;
406  virtual Val_domain mult_r (const Val_domain&) const ;
407  virtual Val_domain div_r (const Val_domain&) const ;
408  virtual Val_domain der_r (const Val_domain&) const ;
409  virtual Val_domain ddt (const Val_domain&) const ;
410  virtual Val_domain dt (const Val_domain&) const ;
411  virtual Val_domain laplacian (const Val_domain&, int) const ;
412 
413  virtual void find_other_dom (int, int, int&, int&) const ;
414  virtual Val_domain der_normal (const Val_domain&, int) const ;
415  virtual double val_boundary (int, const Val_domain&, const Index&) const ;
416 
417  virtual int nbr_unknowns (const Tensor&, int) const ;
424  int nbr_unknowns_val_domain (const Val_domain& so) const ;
425  virtual Array<int> nbr_conditions_array (const Tensor&, int, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
433  int nbr_conditions_val_domain_array (const Val_domain& so, const Array<int>& order) const ;
434  virtual Array<int> nbr_conditions_boundary_array (const Tensor&, int, int, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
444  int nbr_conditions_val_domain_boundary_array (const Val_domain& eq, int bound, const Array<int>& order) const ;
445  virtual void export_tau_array (const Tensor&, int, const Array<int>&, Array<double>&, int&, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
455  void export_tau_val_domain_array (const Val_domain& eq, const Array<int>& order, Array<double>& res, int& pos_res, int ncond) const ;
456  virtual void export_tau_boundary_array (const Tensor&, int, int, const Array<int>&, Array<double>&, int&, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
467  void export_tau_val_domain_boundary_array (const Val_domain& eq, int bound, const Array<int>& order, Array<double>& res, int& pos_res, int ncond) const ;
468  virtual void affecte_tau (Tensor&, int, const Array<double>&, int&) const ;
476  void affecte_tau_val_domain (Val_domain& so, const Array<double>& cf, int& pos_cf) const ;
477  virtual void affecte_tau_one_coef (Tensor&, int, int, int&) const ;
485  void affecte_tau_one_coef_val_domain (Val_domain& so, int cc, int& pos_cf) const ;
486 
487 public:
488  virtual ostream& print (ostream& o) const ;
489 } ;
490 
495 class Space_spheric_time : public Space {
496 
497  protected:
498  double tmin ;
499  double tmax ;
500  bool withcompact ;
501 
502  public:
512  Space_spheric_time (int ttype, const Dim_array& nbr, const Array<double>& bounds, double tmmin, double tmmax, bool wc = false) ;
513  Space_spheric_time (FILE*) ;
514  virtual ~Space_spheric_time() ;
515  virtual void save(FILE*) const ;
516 
520  double get_tmax() const {return tmax ;} ;
521 
525  double get_tmin() const {return tmin ;} ;
526 } ;
527 }
528 #endif
Kadath::Base_spectral
Class for storing the basis of decompositions of a field.
Definition: base_spectral.hpp:69
Kadath::Dim_array
Class for storing the dimensions of an array.
Definition: dim_array.hpp:34
Kadath::Domain_spheric_time_compact
Class for a 2-dimensional compactified spherical domain and a symetry with respect to the plane .
Definition: spheric_time.hpp:354
Kadath::Domain_spheric_time_compact::affecte_tau
virtual void affecte_tau(Tensor &, int, const Array< double > &, int &) const
Affects some coefficients to a Tensor.
Definition: domain_spheric_time_compact_affecte_tau.cpp:43
Kadath::Domain_spheric_time_compact::der_normal
virtual Val_domain der_normal(const Val_domain &, int) const
Normal derivative with respect to a given surface.
Definition: domain_spheric_time_compact.cpp:72
Kadath::Domain_spheric_time_compact::absol_to_num
virtual const Point absol_to_num(const Point &) const
Computes the numerical coordinates from the physical ones.
Definition: domain_spheric_time_compact.cpp:140
Kadath::Domain_spheric_time_compact::div_r
virtual Val_domain div_r(const Val_domain &) const
Division by .
Definition: domain_spheric_time_compact_ope.cpp:59
Kadath::Domain_spheric_time_compact::alpha
double alpha
Relates the numerical to the physical radii.
Definition: spheric_time.hpp:357
Kadath::Domain_spheric_time_compact::do_absol
virtual void do_absol() const
Computes the absolute coordinates.
Definition: domain_spheric_time_compact.cpp:95
Kadath::Domain_spheric_time_compact::set_cheb_base
virtual void set_cheb_base(Base_spectral &) const
Gives the standard base for Chebyshev polynomials.
Definition: domain_spheric_time_compact.cpp:182
Kadath::Domain_spheric_time_compact::Domain_spheric_time_compact
Domain_spheric_time_compact(int num, int ttype, double tmmin, double tmmax, double r_int, const Dim_array &nbr)
Standard constructor :
Definition: domain_spheric_time_compact.cpp:32
Kadath::Domain_spheric_time_compact::div_xm1
virtual Val_domain div_xm1(const Val_domain &) const
Division by .
Definition: domain_spheric_time_compact_ope.cpp:41
Kadath::Domain_spheric_time_compact::affecte_tau_one_coef
virtual void affecte_tau_one_coef(Tensor &, int, int, int &) const
Sets at most one coefficient of a Tensor to 1.
Definition: domain_spheric_time_compact_affecte_tau_one_coef.cpp:55
Kadath::Domain_spheric_time_compact::is_in
virtual bool is_in(const Point &xx, double prec=1e-13) const
Check whether a point lies inside Domain.
Definition: domain_spheric_time_compact.cpp:129
Kadath::Domain_spheric_time_compact::nbr_unknowns
virtual int nbr_unknowns(const Tensor &, int) const
Computes the number of true unknowns of a Tensor, in a given domain.
Definition: domain_spheric_time_compact_nbr_unknowns.cpp:34
Kadath::Domain_spheric_time_compact::nbr_conditions_val_domain_boundary_array
int nbr_conditions_val_domain_boundary_array(const Val_domain &eq, int bound, const Array< int > &order) const
Computes number of discretized equations associated with a given equation on a boundary.
Definition: domain_spheric_time_compact_nbr_conditions_boundary.cpp:28
Kadath::Domain_spheric_time_compact::export_tau_boundary_array
virtual void export_tau_boundary_array(const Tensor &, int, int, const Array< int > &, Array< double > &, int &, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Exports all the residual equations corresponding to one tensorial one on a given boundary It makes us...
Definition: domain_spheric_time_compact_export_tau_boundary.cpp:68
Kadath::Domain_spheric_time_compact::tmin
double tmin
Initial time .
Definition: spheric_time.hpp:358
Kadath::Domain_spheric_time_compact::do_der_abs_from_der_var
virtual void do_der_abs_from_der_var(const Val_domain *const *const der_var, Val_domain **const der_abs) const
Computes the derivative with respect to the absolute Cartesian coordinates from the derivative with r...
Definition: domain_spheric_time_compact.cpp:209
Kadath::Domain_spheric_time_compact::do_coloc
virtual void do_coloc()
Computes the colocation points.
Definition: domain_spheric_time_compact.cpp:152
Kadath::Domain_spheric_time_compact::save
virtual void save(FILE *) const
Saving function.
Definition: domain_spheric_time_compact.cpp:52
Kadath::Domain_spheric_time_compact::tmax
double tmax
Final time .
Definition: spheric_time.hpp:359
Kadath::Domain_spheric_time_compact::get_alpha
double get_alpha() const
Definition: spheric_time.hpp:401
Kadath::Domain_spheric_time_compact::affecte_tau_val_domain
void affecte_tau_val_domain(Val_domain &so, const Array< double > &cf, int &pos_cf) const
Affects some coefficients to a Val_domain.
Definition: domain_spheric_time_compact_affecte_tau.cpp:27
Kadath::Domain_spheric_time_compact::nbr_conditions_array
virtual Array< int > nbr_conditions_array(const Tensor &, int, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Computes number of discretized equations associated with a given tensorial equation in the bulk.
Definition: domain_spheric_time_compact_nbr_conditions.cpp:34
Kadath::Domain_spheric_time_compact::export_tau_array
virtual void export_tau_array(const Tensor &, int, const Array< int > &, Array< double > &, int &, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Exports all the residual equations corresponding to one tensorial one in the bulk.
Definition: domain_spheric_time_compact_export_tau.cpp:48
Kadath::Domain_spheric_time_compact::export_tau_val_domain_boundary_array
void export_tau_val_domain_boundary_array(const Val_domain &eq, int bound, const Array< int > &order, Array< double > &res, int &pos_res, int ncond) const
Exports all the residual equations corresponding to a tensorial one on a given boundary It makes use ...
Definition: domain_spheric_time_compact_export_tau_boundary.cpp:28
Kadath::Domain_spheric_time_compact::find_other_dom
virtual void find_other_dom(int, int, int &, int &) const
Gives the informations corresponding the a touching neighboring domain.
Definition: domain_spheric_time_compact_systems.cpp:25
Kadath::Domain_spheric_time_compact::print
virtual ostream & print(ostream &o) const
Delegate function to virtualize the << operator.
Definition: domain_spheric_time_compact.cpp:62
Kadath::Domain_spheric_time_compact::export_tau_val_domain_array
void export_tau_val_domain_array(const Val_domain &eq, const Array< int > &order, Array< double > &res, int &pos_res, int ncond) const
Exports a residual equation in the bulk.
Definition: domain_spheric_time_compact_export_tau.cpp:28
Kadath::Domain_spheric_time_compact::dt
virtual Val_domain dt(const Val_domain &) const
Compute the derivative with respect to of a scalar field.
Definition: domain_spheric_time_compact_ope.cpp:75
Kadath::Domain_spheric_time_compact::do_radius
virtual void do_radius() const
Computes the generalized radius.
Definition: domain_spheric_time_compact.cpp:114
Kadath::Domain_spheric_time_compact::laplacian
virtual Val_domain laplacian(const Val_domain &, int) const
Computes the ordinary flat Laplacian for a scalar field with an harmonic index m.
Definition: domain_spheric_time_compact_ope.cpp:80
Kadath::Domain_spheric_time_compact::set_legendre_base
virtual void set_legendre_base(Base_spectral &) const
Gives the standard base for Legendre polynomials.
Definition: domain_spheric_time_compact.cpp:196
Kadath::Domain_spheric_time_compact::nbr_conditions_boundary_array
virtual Array< int > nbr_conditions_boundary_array(const Tensor &, int, int, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Computes number of discretized equations associated with a given tensorial equation on a boundary.
Definition: domain_spheric_time_compact_nbr_conditions_boundary.cpp:47
Kadath::Domain_spheric_time_compact::val_boundary
virtual double val_boundary(int, const Val_domain &, const Index &) const
Computes the value of a field at a boundary.
Definition: domain_spheric_time_compact_systems.cpp:38
Kadath::Domain_spheric_time_compact::mult_xm1
virtual Val_domain mult_xm1(const Val_domain &) const
Multiplication by .
Definition: domain_spheric_time_compact_ope.cpp:30
Kadath::Domain_spheric_time_compact::nbr_unknowns_val_domain
int nbr_unknowns_val_domain(const Val_domain &so) const
Computes the number of true unknowns of a Val_domain.
Definition: domain_spheric_time_compact_nbr_unknowns.cpp:28
Kadath::Domain_spheric_time_compact::affecte_tau_one_coef_val_domain
void affecte_tau_one_coef_val_domain(Val_domain &so, int cc, int &pos_cf) const
Sets at most one coefficient of a Val_domain to 1.
Definition: domain_spheric_time_compact_affecte_tau_one_coef.cpp:27
Kadath::Domain_spheric_time_compact::ddt
virtual Val_domain ddt(const Val_domain &) const
Compute the second derivative with respect to of a scalar field.
Definition: domain_spheric_time_compact_ope.cpp:71
Kadath::Domain_spheric_time_compact::der_r
virtual Val_domain der_r(const Val_domain &) const
Compute the radial derivative of a scalar field.
Definition: domain_spheric_time_compact_ope.cpp:66
Kadath::Domain_spheric_time_compact::mult_r
virtual Val_domain mult_r(const Val_domain &) const
Multiplication by .
Definition: domain_spheric_time_compact_ope.cpp:52
Kadath::Domain_spheric_time_compact::mult
virtual Base_spectral mult(const Base_spectral &, const Base_spectral &) const
Method for the multiplication of two Base_spectral.
Definition: domain_spheric_time_compact.cpp:217
Kadath::Domain_spheric_time_compact::nbr_conditions_val_domain_array
int nbr_conditions_val_domain_array(const Val_domain &so, const Array< int > &order) const
Computes number of discretized equations associated with a given tensorial equation in the bulk.
Definition: domain_spheric_time_compact_nbr_conditions.cpp:28
Kadath::Domain_spheric_time_nucleus
Class for a 2-dimensional spherical domain containing the origin and a symetry with respect to the pl...
Definition: spheric_time.hpp:45
Kadath::Domain_spheric_time_nucleus::affecte_tau_val_domain
void affecte_tau_val_domain(Val_domain &so, const Array< double > &cf, int &pos_cf) const
Affects some coefficients to a Val_domain.
Definition: domain_spheric_time_nucleus_affecte_tau.cpp:27
Kadath::Domain_spheric_time_nucleus::affecte_tau_one_coef
virtual void affecte_tau_one_coef(Tensor &, int, int, int &) const
Sets at most one coefficient of a Tensor to 1.
Definition: domain_spheric_time_nucleus_affecte_tau_one_coef.cpp:55
Kadath::Domain_spheric_time_nucleus::print
virtual ostream & print(ostream &o) const
Delegate function to virtualize the << operator.
Definition: domain_spheric_time_nucleus.cpp:63
Kadath::Domain_spheric_time_nucleus::absol_to_num
virtual const Point absol_to_num(const Point &) const
Computes the numerical coordinates from the physical ones.
Definition: domain_spheric_time_nucleus.cpp:136
Kadath::Domain_spheric_time_nucleus::Domain_spheric_time_nucleus
Domain_spheric_time_nucleus(int num, int ttype, double tmmin, double tmmax, double radius, const Dim_array &nbr)
Standard constructor :
Definition: domain_spheric_time_nucleus.cpp:33
Kadath::Domain_spheric_time_nucleus::dt
virtual Val_domain dt(const Val_domain &) const
Compute the derivative with respect to of a scalar field.
Definition: domain_spheric_time_nucleus_ope.cpp:76
Kadath::Domain_spheric_time_nucleus::is_in
virtual bool is_in(const Point &xx, double prec=1e-13) const
Check whether a point lies inside Domain.
Definition: domain_spheric_time_nucleus.cpp:125
Kadath::Domain_spheric_time_nucleus::ddt
virtual Val_domain ddt(const Val_domain &) const
Compute the second derivative with respect to of a scalar field.
Definition: domain_spheric_time_nucleus_ope.cpp:72
Kadath::Domain_spheric_time_nucleus::do_absol
virtual void do_absol() const
Computes the absolute coordinates.
Definition: domain_spheric_time_nucleus.cpp:91
Kadath::Domain_spheric_time_nucleus::set_legendre_base
virtual void set_legendre_base(Base_spectral &) const
Gives the standard base for Legendre polynomials.
Definition: domain_spheric_time_nucleus.cpp:193
Kadath::Domain_spheric_time_nucleus::tmax
double tmax
Final time .
Definition: spheric_time.hpp:50
Kadath::Domain_spheric_time_nucleus::export_tau_val_domain_boundary_array
void export_tau_val_domain_boundary_array(const Val_domain &eq, int bound, const Array< int > &order, Array< double > &res, int &pos_res, int ncond) const
Exports all the residual equations corresponding to a tensorial one on a given boundary It makes use ...
Definition: domain_spheric_time_nucleus_export_tau_boundary.cpp:28
Kadath::Domain_spheric_time_nucleus::nbr_unknowns
virtual int nbr_unknowns(const Tensor &, int) const
Computes the number of true unknowns of a Tensor, in a given domain.
Definition: domain_spheric_time_nucleus_nbr_unknowns.cpp:34
Kadath::Domain_spheric_time_nucleus::div_x
virtual Val_domain div_x(const Val_domain &) const
Division by .
Definition: domain_spheric_time_nucleus_ope.cpp:28
Kadath::Domain_spheric_time_nucleus::set_cheb_base
virtual void set_cheb_base(Base_spectral &) const
Gives the standard base for Chebyshev polynomials.
Definition: domain_spheric_time_nucleus.cpp:179
Kadath::Domain_spheric_time_nucleus::mult
virtual Base_spectral mult(const Base_spectral &, const Base_spectral &) const
Method for the multiplication of two Base_spectral.
Definition: domain_spheric_time_nucleus.cpp:214
Kadath::Domain_spheric_time_nucleus::export_tau_array
virtual void export_tau_array(const Tensor &, int, const Array< int > &, Array< double > &, int &, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Exports all the residual equations corresponding to one tensorial one in the bulk.
Definition: domain_spheric_time_nucleus_export_tau.cpp:48
Kadath::Domain_spheric_time_nucleus::get_alpha
double get_alpha() const
Definition: spheric_time.hpp:91
Kadath::Domain_spheric_time_nucleus::affecte_tau
virtual void affecte_tau(Tensor &, int, const Array< double > &, int &) const
Affects some coefficients to a Tensor.
Definition: domain_spheric_time_nucleus_affecte_tau.cpp:43
Kadath::Domain_spheric_time_nucleus::nbr_conditions_array
virtual Array< int > nbr_conditions_array(const Tensor &, int, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Computes number of discretized equations associated with a given tensorial equation in the bulk.
Definition: domain_spheric_time_nucleus_nbr_conditions.cpp:34
Kadath::Domain_spheric_time_nucleus::find_other_dom
virtual void find_other_dom(int, int, int &, int &) const
Gives the informations corresponding the a touching neighboring domain.
Definition: domain_spheric_time_nucleus_systems.cpp:25
Kadath::Domain_spheric_time_nucleus::srdr
virtual Val_domain srdr(const Val_domain &) const
Compute the of a scalar field .
Definition: domain_spheric_time_nucleus_ope.cpp:68
Kadath::Domain_spheric_time_nucleus::nbr_conditions_boundary_array
virtual Array< int > nbr_conditions_boundary_array(const Tensor &, int, int, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Computes number of discretized equations associated with a given tensorial equation on a boundary.
Definition: domain_spheric_time_nucleus_nbr_conditions_boundary.cpp:44
Kadath::Domain_spheric_time_nucleus::export_tau_boundary_array
virtual void export_tau_boundary_array(const Tensor &, int, int, const Array< int > &, Array< double > &, int &, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Exports all the residual equations corresponding to one tensorial one on a given boundary It makes us...
Definition: domain_spheric_time_nucleus_export_tau_boundary.cpp:60
Kadath::Domain_spheric_time_nucleus::alpha
double alpha
Relates the numerical to the physical radii.
Definition: spheric_time.hpp:48
Kadath::Domain_spheric_time_nucleus::nbr_conditions_val_domain_array
int nbr_conditions_val_domain_array(const Val_domain &so, const Array< int > &order) const
Computes number of discretized equations associated with a given tensorial equation in the bulk.
Definition: domain_spheric_time_nucleus_nbr_conditions.cpp:28
Kadath::Domain_spheric_time_nucleus::do_coloc
virtual void do_coloc()
Computes the colocation points.
Definition: domain_spheric_time_nucleus.cpp:149
Kadath::Domain_spheric_time_nucleus::val_boundary
virtual double val_boundary(int, const Val_domain &, const Index &) const
Computes the value of a field at a boundary.
Definition: domain_spheric_time_nucleus_systems.cpp:38
Kadath::Domain_spheric_time_nucleus::nbr_unknowns_val_domain
int nbr_unknowns_val_domain(const Val_domain &so) const
Computes the number of true unknowns of a Val_domain.
Definition: domain_spheric_time_nucleus_nbr_unknowns.cpp:28
Kadath::Domain_spheric_time_nucleus::do_der_abs_from_der_var
virtual void do_der_abs_from_der_var(const Val_domain *const *const der_var, Val_domain **const der_abs) const
Computes the derivative with respect to the absolute Cartesian coordinates from the derivative with r...
Definition: domain_spheric_time_nucleus.cpp:206
Kadath::Domain_spheric_time_nucleus::affecte_tau_one_coef_val_domain
void affecte_tau_one_coef_val_domain(Val_domain &so, int cc, int &pos_cf) const
Sets at most one coefficient of a Val_domain to 1.
Definition: domain_spheric_time_nucleus_affecte_tau_one_coef.cpp:27
Kadath::Domain_spheric_time_nucleus::save
virtual void save(FILE *) const
Saving function.
Definition: domain_spheric_time_nucleus.cpp:53
Kadath::Domain_spheric_time_nucleus::der_r
virtual Val_domain der_r(const Val_domain &) const
Compute the radial derivative of a scalar field.
Definition: domain_spheric_time_nucleus_ope.cpp:63
Kadath::Domain_spheric_time_nucleus::tmin
double tmin
Initial time .
Definition: spheric_time.hpp:49
Kadath::Domain_spheric_time_nucleus::mult_r
virtual Val_domain mult_r(const Val_domain &) const
Multiplication by .
Definition: domain_spheric_time_nucleus_ope.cpp:39
Kadath::Domain_spheric_time_nucleus::div_r
virtual Val_domain div_r(const Val_domain &) const
Division by .
Definition: domain_spheric_time_nucleus_ope.cpp:51
Kadath::Domain_spheric_time_nucleus::laplacian
virtual Val_domain laplacian(const Val_domain &, int) const
Computes the ordinary flat Laplacian for a scalar field with an harmonic index m.
Definition: domain_spheric_time_nucleus_ope.cpp:81
Kadath::Domain_spheric_time_nucleus::der_normal
virtual Val_domain der_normal(const Val_domain &, int) const
Normal derivative with respect to a given surface.
Definition: domain_spheric_time_nucleus.cpp:73
Kadath::Domain_spheric_time_nucleus::nbr_conditions_val_domain_boundary_array
int nbr_conditions_val_domain_boundary_array(const Val_domain &eq, int bound, const Array< int > &order) const
Computes number of discretized equations associated with a given equation on a boundary.
Definition: domain_spheric_time_nucleus_nbr_conditions_boundary.cpp:28
Kadath::Domain_spheric_time_nucleus::export_tau_val_domain_array
void export_tau_val_domain_array(const Val_domain &eq, const Array< int > &order, Array< double > &res, int &pos_res, int ncond) const
Exports a residual equation in the bulk.
Definition: domain_spheric_time_nucleus_export_tau.cpp:28
Kadath::Domain_spheric_time_nucleus::do_radius
virtual void do_radius() const
Computes the generalized radius.
Definition: domain_spheric_time_nucleus.cpp:110
Kadath::Domain_spheric_time_shell
Class for a 2-dimensional spherical domain bounded between two finite radii and a symetry with respec...
Definition: spheric_time.hpp:199
Kadath::Domain_spheric_time_shell::do_der_abs_from_der_var
virtual void do_der_abs_from_der_var(const Val_domain *const *const der_var, Val_domain **const der_abs) const
Computes the derivative with respect to the absolute Cartesian coordinates from the derivative with r...
Definition: domain_spheric_time_shell.cpp:212
Kadath::Domain_spheric_time_shell::Domain_spheric_time_shell
Domain_spheric_time_shell(int num, int ttype, double tmmin, double tmmax, double r_int, double r_ext, const Dim_array &nbr)
Standard constructor :
Definition: domain_spheric_time_shell.cpp:32
Kadath::Domain_spheric_time_shell::set_legendre_base
virtual void set_legendre_base(Base_spectral &) const
Gives the standard base for Legendre polynomials.
Definition: domain_spheric_time_shell.cpp:199
Kadath::Domain_spheric_time_shell::nbr_conditions_val_domain_array
int nbr_conditions_val_domain_array(const Val_domain &so, const Array< int > &order) const
Computes number of discretized equations associated with a given tensorial equation in the bulk.
Definition: domain_spheric_time_shell_nbr_conditions.cpp:28
Kadath::Domain_spheric_time_shell::beta
double beta
Relates the numerical to the physical radii.
Definition: spheric_time.hpp:203
Kadath::Domain_spheric_time_shell::nbr_unknowns
virtual int nbr_unknowns(const Tensor &, int) const
Computes the number of true unknowns of a Tensor, in a given domain.
Definition: domain_spheric_time_shell_nbr_unknowns.cpp:34
Kadath::Domain_spheric_time_shell::tmax
double tmax
Final time .
Definition: spheric_time.hpp:205
Kadath::Domain_spheric_time_shell::nbr_unknowns_val_domain
int nbr_unknowns_val_domain(const Val_domain &so) const
Computes the number of true unknowns of a Val_domain.
Definition: domain_spheric_time_shell_nbr_unknowns.cpp:28
Kadath::Domain_spheric_time_shell::ddt
virtual Val_domain ddt(const Val_domain &) const
Compute the second derivative with respect to of a scalar field.
Definition: domain_spheric_time_shell_ope.cpp:48
Kadath::Domain_spheric_time_shell::do_radius
virtual void do_radius() const
Computes the generalized radius.
Definition: domain_spheric_time_shell.cpp:115
Kadath::Domain_spheric_time_shell::export_tau_array
virtual void export_tau_array(const Tensor &, int, const Array< int > &, Array< double > &, int &, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Exports all the residual equations corresponding to one tensorial one in the bulk.
Definition: domain_spheric_time_shell_export_tau.cpp:48
Kadath::Domain_spheric_time_shell::der_r
virtual Val_domain der_r(const Val_domain &) const
Compute the radial derivative of a scalar field.
Definition: domain_spheric_time_shell_ope.cpp:44
Kadath::Domain_spheric_time_shell::export_tau_val_domain_array
void export_tau_val_domain_array(const Val_domain &eq, const Array< int > &order, Array< double > &res, int &pos_res, int ncond) const
Exports a residual equation in the bulk.
Definition: domain_spheric_time_shell_export_tau.cpp:28
Kadath::Domain_spheric_time_shell::alpha
double alpha
Relates the numerical to the physical radii.
Definition: spheric_time.hpp:202
Kadath::Domain_spheric_time_shell::nbr_conditions_val_domain_boundary_array
int nbr_conditions_val_domain_boundary_array(const Val_domain &eq, int bound, const Array< int > &order) const
Computes number of discretized equations associated with a given equation on a boundary.
Definition: domain_spheric_time_shell_nbr_conditions_boundary.cpp:28
Kadath::Domain_spheric_time_shell::div_r
virtual Val_domain div_r(const Val_domain &) const
Division by .
Definition: domain_spheric_time_shell_ope.cpp:38
Kadath::Domain_spheric_time_shell::do_coloc
virtual void do_coloc()
Computes the colocation points.
Definition: domain_spheric_time_shell.cpp:155
Kadath::Domain_spheric_time_shell::affecte_tau_val_domain
void affecte_tau_val_domain(Val_domain &so, const Array< double > &cf, int &pos_cf) const
Affects some coefficients to a Val_domain.
Definition: domain_spheric_time_shell_affecte_tau.cpp:27
Kadath::Domain_spheric_time_shell::find_other_dom
virtual void find_other_dom(int, int, int &, int &) const
Gives the informations corresponding the a touching neighboring domain.
Definition: domain_spheric_time_shell_systems.cpp:25
Kadath::Domain_spheric_time_shell::der_normal
virtual Val_domain der_normal(const Val_domain &, int) const
Normal derivative with respect to a given surface.
Definition: domain_spheric_time_shell.cpp:74
Kadath::Domain_spheric_time_shell::export_tau_boundary_array
virtual void export_tau_boundary_array(const Tensor &, int, int, const Array< int > &, Array< double > &, int &, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Exports all the residual equations corresponding to one tensorial one on a given boundary It makes us...
Definition: domain_spheric_time_shell_export_tau_boundary.cpp:68
Kadath::Domain_spheric_time_shell::mult_r
virtual Val_domain mult_r(const Val_domain &) const
Multiplication by .
Definition: domain_spheric_time_shell_ope.cpp:27
Kadath::Domain_spheric_time_shell::laplacian
virtual Val_domain laplacian(const Val_domain &, int) const
Computes the ordinary flat Laplacian for a scalar field with an harmonic index m.
Definition: domain_spheric_time_shell_ope.cpp:56
Kadath::Domain_spheric_time_shell::absol_to_num
virtual const Point absol_to_num(const Point &) const
Computes the numerical coordinates from the physical ones.
Definition: domain_spheric_time_shell.cpp:143
Kadath::Domain_spheric_time_shell::save
virtual void save(FILE *) const
Saving function.
Definition: domain_spheric_time_shell.cpp:53
Kadath::Domain_spheric_time_shell::affecte_tau
virtual void affecte_tau(Tensor &, int, const Array< double > &, int &) const
Affects some coefficients to a Tensor.
Definition: domain_spheric_time_shell_affecte_tau.cpp:43
Kadath::Domain_spheric_time_shell::dt
virtual Val_domain dt(const Val_domain &) const
Compute the derivative with respect to of a scalar field.
Definition: domain_spheric_time_shell_ope.cpp:52
Kadath::Domain_spheric_time_shell::affecte_tau_one_coef_val_domain
void affecte_tau_one_coef_val_domain(Val_domain &so, int cc, int &pos_cf) const
Sets at most one coefficient of a Val_domain to 1.
Definition: domain_spheric_time_shell_affecte_tau_one_coef.cpp:27
Kadath::Domain_spheric_time_shell::val_boundary
virtual double val_boundary(int, const Val_domain &, const Index &) const
Computes the value of a field at a boundary.
Definition: domain_spheric_time_shell_systems.cpp:42
Kadath::Domain_spheric_time_shell::is_in
virtual bool is_in(const Point &xx, double prec=1e-13) const
Check whether a point lies inside Domain.
Definition: domain_spheric_time_shell.cpp:130
Kadath::Domain_spheric_time_shell::affecte_tau_one_coef
virtual void affecte_tau_one_coef(Tensor &, int, int, int &) const
Sets at most one coefficient of a Tensor to 1.
Definition: domain_spheric_time_shell_affecte_tau_one_coef.cpp:55
Kadath::Domain_spheric_time_shell::export_tau_val_domain_boundary_array
void export_tau_val_domain_boundary_array(const Val_domain &eq, int bound, const Array< int > &order, Array< double > &res, int &pos_res, int ncond) const
Exports all the residual equations corresponding to a tensorial one on a given boundary It makes use ...
Definition: domain_spheric_time_shell_export_tau_boundary.cpp:28
Kadath::Domain_spheric_time_shell::do_absol
virtual void do_absol() const
Computes the absolute coordinates.
Definition: domain_spheric_time_shell.cpp:96
Kadath::Domain_spheric_time_shell::get_alpha
double get_alpha() const
Definition: spheric_time.hpp:248
Kadath::Domain_spheric_time_shell::mult
virtual Base_spectral mult(const Base_spectral &, const Base_spectral &) const
Method for the multiplication of two Base_spectral.
Definition: domain_spheric_time_shell.cpp:220
Kadath::Domain_spheric_time_shell::print
virtual ostream & print(ostream &o) const
Delegate function to virtualize the << operator.
Definition: domain_spheric_time_shell.cpp:64
Kadath::Domain_spheric_time_shell::nbr_conditions_boundary_array
virtual Array< int > nbr_conditions_boundary_array(const Tensor &, int, int, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Computes number of discretized equations associated with a given tensorial equation on a boundary.
Definition: domain_spheric_time_shell_nbr_conditions_boundary.cpp:47
Kadath::Domain_spheric_time_shell::tmin
double tmin
Initial time .
Definition: spheric_time.hpp:204
Kadath::Domain_spheric_time_shell::set_cheb_base
virtual void set_cheb_base(Base_spectral &) const
Gives the standard base for Chebyshev polynomials.
Definition: domain_spheric_time_shell.cpp:185
Kadath::Domain_spheric_time_shell::nbr_conditions_array
virtual Array< int > nbr_conditions_array(const Tensor &, int, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Computes number of discretized equations associated with a given tensorial equation in the bulk.
Definition: domain_spheric_time_shell_nbr_conditions.cpp:34
Kadath::Domain
Abstract class that implements the fonctionnalities common to all the type of domains.
Definition: space.hpp:60
Kadath::Domain::radius
Val_domain * radius
The generalized radius.
Definition: space.hpp:78
Kadath::Index
Class that gives the position inside a multi-dimensional Array.
Definition: index.hpp:38
Kadath::Point
The class Point is used to store the coordinates of a point.
Definition: point.hpp:30
Kadath::Space_spheric_time
The Space_oned class fills the space with 2-dimensioanl domains : spherically symmetric and time depe...
Definition: spheric_time.hpp:495
Kadath::Space_spheric_time::get_tmax
double get_tmax() const
Returns the final time.
Definition: spheric_time.hpp:520
Kadath::Space_spheric_time::withcompact
bool withcompact
Indicator of the presence of a compactified Domain.
Definition: spheric_time.hpp:500
Kadath::Space_spheric_time::Space_spheric_time
Space_spheric_time(int ttype, const Dim_array &nbr, const Array< double > &bounds, double tmmin, double tmmax, bool wc=false)
Standard constructor.
Definition: space_spheric_time.cpp:27
Kadath::Space_spheric_time::get_tmin
double get_tmin() const
Returns the initial time.
Definition: spheric_time.hpp:525
Kadath::Space_spheric_time::tmax
double tmax
Final time .
Definition: spheric_time.hpp:499
Kadath::Space_spheric_time::tmin
double tmin
Initial time .
Definition: spheric_time.hpp:498
Kadath::Space_spheric_time::save
virtual void save(FILE *) const
Saving function.
Definition: space_spheric_time.cpp:79
Kadath::Space
The Space class is an ensemble of domains describing the whole space of the computation.
Definition: space.hpp:1362
Kadath::Tensor
Tensor handling.
Definition: tensor.hpp:149
Kadath::Val_domain
Class for storing the basis of decompositions of a field and its values on both the configuration and...
Definition: val_domain.hpp:69