polar_periodic.hpp

/*
    Copyright 2017 Philippe Grandclement
 
    This file is part of Kadath.
 
    Kadath is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    Kadath is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.
 
    You should have received a copy of the GNU General Public License
    along with Kadath.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef __POLAR_PERIODIC_HPP_
#define __POLAR_PERIODIC_HPP_
 
#include "space.hpp"
 
#define TO_PI 0
namespace Kadath {
 
class Domain_polar_periodic_nucleus : public Domain {
 
 private:
  double alpha ; 
 
  double ome ; 
  mutable Term_eq* ome_term_eq ; 
 
  int type_time ;
  double maxt ; 
 
 public: 
  Domain_polar_periodic_nucleus (int num, int ttype, double radius, double ome, const Dim_array& nbr) ;
  Domain_polar_periodic_nucleus (const Domain_polar_periodic_nucleus& so) ; 
  Domain_polar_periodic_nucleus (int num, FILE* ff) ;
 
  virtual ~Domain_polar_periodic_nucleus() ;
  virtual void save (FILE*) const ;
 
  private:    
    virtual void do_absol ()  const ;
 
  private:
     virtual void set_cheb_base(Base_spectral&) const ;       
     virtual void set_legendre_base(Base_spectral&) const ;
     virtual void set_anti_cheb_base(Base_spectral&) const ;       
     virtual void set_anti_legendre_base(Base_spectral&) const ;
     virtual void set_cheb_base_with_m(Base_spectral&, int m) const ;       
     virtual void set_legendre_base_with_m(Base_spectral&, int m) const ;
     virtual void set_anti_cheb_base_with_m(Base_spectral&, int m) const ;       
     virtual void set_anti_legendre_base_with_m(Base_spectral&, int m) const ;
     virtual void set_cheb_base_r_spher(Base_spectral&) const ;
     virtual void set_cheb_base_t_spher(Base_spectral&) const ;
     virtual void set_cheb_base_p_spher(Base_spectral&) const ;
     virtual void set_legendre_base_r_spher(Base_spectral&) const ;
     virtual void set_legendre_base_t_spher(Base_spectral&) const ;
     virtual void set_legendre_base_p_spher(Base_spectral&) const ;
 
 
     virtual void do_coloc () ;
     virtual void do_radius () const ;
 
  public:
     virtual bool is_in(const Point&xx, double prec=1e-13) const ;
     virtual const Point absol_to_num(const Point&) const;
     virtual void do_der_abs_from_der_var(const Val_domain *const *const der_var, Val_domain **const der_abs) const ;
     virtual Base_spectral mult (const Base_spectral&, const Base_spectral&) const ;
     double get_ome() const {return ome ;} ;
 
  public:
     virtual Val_domain mult_cos_theta (const Val_domain&) const ;
     virtual Val_domain mult_sin_theta (const Val_domain&) const ;
     virtual Val_domain div_sin_theta (const Val_domain&) const ;
     virtual Val_domain div_x (const Val_domain&) const ;     
     virtual Val_domain mult_r (const Val_domain&) const ;
     virtual Val_domain div_r (const Val_domain&) const ;
     virtual Val_domain der_r (const Val_domain&) const ; 
     virtual Val_domain dt (const Val_domain&) const ;
     virtual Val_domain srdr (const Val_domain&) const ;
     virtual Val_domain dtime (const Val_domain&) const ;
     virtual Val_domain mult_cos_time (const Val_domain&) const ;
     virtual Val_domain mult_sin_time (const Val_domain&) const ;
   
     virtual Term_eq partial_spher (const Term_eq&) const ;
     virtual Term_eq derive_flat_spher (int, char, const Term_eq&, const Metric*) const ;
 
     virtual double val_boundary (int, const Val_domain&, const Index&) const ;
     virtual void find_other_dom (int, int, int&, int&) const ;
     virtual Val_domain der_normal (const Val_domain&, int) const ;
    
 
  //   virtual int nbr_unknowns_from_adapted() const ;
   //  virtual void vars_to_terms() const ;
   //  virtual void affecte_coef (int&, int, bool&) const ;
   //  virtual void xx_to_vars_from_adapted (double, const Array<double>&, int&) const ;
   //  virtual void xx_to_ders_from_adapted (const Array<double>&, int&) const ;
   //  virtual void update_term_eq (Term_eq*) const ;
   //  virtual void update_variable (double, const Scalar&, Scalar&) const ;   
   //  virtual void update_constante (double, const Scalar&, Scalar&) const ;   
   //  virtual void update_mapping(double) ;
 
    // void update() const ;
    
     virtual Term_eq dtime_term_eq (const Term_eq& so) const ;   
     virtual Term_eq ddtime_term_eq (const Term_eq& so) const ;   
 
         virtual int nbr_unknowns (const Tensor&, int) const ;
     int nbr_unknowns_val_domain (const Val_domain& so,  int llim) const ;
 
     virtual Array<int> nbr_conditions (const Tensor&, int, int, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
     int nbr_conditions_val_domain (const Val_domain& so, int llim, int order) const ;
     virtual Array<int> nbr_conditions_boundary (const Tensor&, int, int, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
     int nbr_conditions_val_domain_boundary (const Val_domain& eq) const ;
     virtual void export_tau (const Tensor&, int, int, Array<double>&, int&, const Array<int>&,int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
     void export_tau_val_domain (const Val_domain& eq, int llim, int order, Array<double>& res, int& pos_res, int ncond) const ;
     virtual void export_tau_boundary (const Tensor&, int, int, Array<double>&, int&, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
     void export_tau_val_domain_boundary (const Val_domain& eq, int bound, Array<double>& res, int& pos_res, int ncond) const ;
     virtual void affecte_tau (Tensor&, int, const Array<double>&, int&) const ;
     void affecte_tau_val_domain (Val_domain& so, int llim, const Array<double>& cf, int& pos_cf) const ;
     virtual void affecte_tau_one_coef (Tensor&, int, int, int&) const ;
      void affecte_tau_one_coef_val_domain (Val_domain& so, int llim, int cc, int& pos_cf) const ;
 
 
public:
     virtual ostream& print (ostream& o) const ;
} ;
 
class Domain_polar_periodic_shell : public Domain {
 
 private:
  double alpha ; 
  double beta ; 
 
 
  double ome ; 
  mutable Term_eq* ome_term_eq ; 
 
  int type_time ;
  double maxt ; 
 
 public: 
  Domain_polar_periodic_shell (int num, int ttype, double rin, double rout, double ome, const Dim_array& nbr) ;
  Domain_polar_periodic_shell (const Domain_polar_periodic_shell& so) ; 
  Domain_polar_periodic_shell (int num, FILE* ff) ;
 
  virtual ~Domain_polar_periodic_shell() ;
  virtual void save (FILE*) const ;
 
  private:    
    virtual void do_absol ()  const ;
 
  private:
     virtual void set_cheb_base(Base_spectral&) const ;       
     virtual void set_legendre_base(Base_spectral&) const ;
     virtual void set_anti_cheb_base(Base_spectral&) const ;       
     virtual void set_anti_legendre_base(Base_spectral&) const ;
     virtual void set_cheb_base_with_m(Base_spectral&, int m) const ;       
     virtual void set_legendre_base_with_m(Base_spectral&, int m) const ;
     virtual void set_anti_cheb_base_with_m(Base_spectral&, int m) const ;       
     virtual void set_anti_legendre_base_with_m(Base_spectral&, int m) const ;
     virtual void set_cheb_base_r_spher(Base_spectral&) const ;
     virtual void set_cheb_base_t_spher(Base_spectral&) const ;
     virtual void set_cheb_base_p_spher(Base_spectral&) const ;
     virtual void set_legendre_base_r_spher(Base_spectral&) const ;
     virtual void set_legendre_base_t_spher(Base_spectral&) const ;
     virtual void set_legendre_base_p_spher(Base_spectral&) const ;
 
     virtual void do_coloc () ;
     virtual void do_radius () const ;
 
  public:
     virtual bool is_in(const Point&xx, double prec=1e-13) const ;
     virtual const Point absol_to_num(const Point&) const;
     virtual void do_der_abs_from_der_var(const Val_domain *const *const der_var, Val_domain **const der_abs) const ;
     virtual Base_spectral mult (const Base_spectral&, const Base_spectral&) const ;
 
  public:
     virtual Val_domain mult_cos_theta (const Val_domain&) const ;
     virtual Val_domain mult_sin_theta (const Val_domain&) const ;
     virtual Val_domain div_sin_theta (const Val_domain&) const ;
     virtual Val_domain mult_r (const Val_domain&) const ;
     virtual Val_domain div_r (const Val_domain&) const ;
     virtual Val_domain der_r (const Val_domain&) const ; 
     virtual Val_domain dt (const Val_domain&) const ; 
     virtual Val_domain div_1mrsL (const Val_domain&) const ;
     virtual Val_domain dtime (const Val_domain&) const ;
     virtual Val_domain mult_cos_time (const Val_domain&) const ;
     virtual Val_domain mult_sin_time (const Val_domain&) const ;
   
     virtual Term_eq partial_spher (const Term_eq&) const ;
     virtual Term_eq derive_flat_spher (int, char, const Term_eq&, const Metric*) const ;
 
     virtual double val_boundary (int, const Val_domain&, const Index&) const ;
     virtual void find_other_dom (int, int, int&, int&) const ;
     virtual Val_domain der_normal (const Val_domain&, int) const ;
 
  /*   virtual int nbr_unknowns_from_adapted() const ;
     virtual void vars_to_terms() const ;
     virtual void affecte_coef (int&, int, bool&) const ;
     virtual void xx_to_vars_from_adapted (double, const Array<double>&, int&) const ;
     virtual void xx_to_ders_from_adapted (const Array<double>&, int&) const ;
     virtual void update_term_eq (Term_eq*) const ;
     virtual void update_variable (double, const Scalar&, Scalar&) const ;   
     virtual void update_constante (double, const Scalar&, Scalar&) const ;   
     virtual void update_mapping(double) ;
*/
  //   void update() const ;
    
     virtual Term_eq dtime_term_eq (const Term_eq& so) const ;   
     virtual Term_eq ddtime_term_eq (const Term_eq& so) const ;   
 
 
         virtual int nbr_unknowns (const Tensor&, int) const ;
     int nbr_unknowns_val_domain (const Val_domain& so) const ;
 
     virtual Array<int> nbr_conditions (const Tensor&, int, int, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
     int nbr_conditions_val_domain (const Val_domain& so, int order) const ;
     virtual Array<int> nbr_conditions_boundary (const Tensor&, int, int, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
     int nbr_conditions_val_domain_boundary (const Val_domain& eq) const ;
     virtual void export_tau (const Tensor&, int, int, Array<double>&, int&, const Array<int>&,int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
     void export_tau_val_domain (const Val_domain& eq, int order, Array<double>& res, int& pos_res, int ncond) const ;
     virtual void export_tau_boundary (const Tensor&, int, int, Array<double>&, int&, const Array<int>&, int n_cmp=-1, Array<int>** p_cmp=0x0) const ;
     void export_tau_val_domain_boundary (const Val_domain& eq, int bound, Array<double>& res, int& pos_res, int ncond) const ;
     virtual void affecte_tau (Tensor&, int, const Array<double>&, int&) const ;
     void affecte_tau_val_domain (Val_domain& so, const Array<double>& cf, int& pos_cf) const ;
     virtual void affecte_tau_one_coef (Tensor&, int, int, int&) const ;
      void affecte_tau_one_coef_val_domain (Val_domain& so, int cc, int& pos_cf) const ;
 
 
public:
     virtual ostream& print (ostream& o) const ;
} ;
 
 
class Space_polar_periodic : public Space {
     public:
    Space_polar_periodic (int ttype, double omega, const Dim_array& nbr, const Array<double>& bounds) ;
    Space_polar_periodic (FILE*) ; 
    virtual ~Space_polar_periodic() ;      
    virtual void save(FILE*) const ;
 
    /*virtual int nbr_unknowns_from_variable_domains() const ;
    virtual void affecte_coef_to_variable_domains(int& , int, Array<int>&) const ;
    virtual void xx_to_ders_variable_domains(const Array<double>&, int&) const ;
    virtual void xx_to_vars_variable_domains(System_of_eqs*, const Array<double>&, int&) const ;
*/
    double get_omega() const ;
} ;
 
}
#endif