oned.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 __ONED_HPP_
#define __ONED_HPP_
 
#include "space.hpp"
#include "val_domain.hpp" 
 
namespace Kadath {
class Domain_oned_ori : public Domain {
 
 private:
  double alpha ; 
  
 public: 
  Domain_oned_ori (int num, int ttype, double radius, const Dim_array& nbr) ;
  Domain_oned_ori (const Domain_oned_ori& so) ; 
  Domain_oned_ori (int num, FILE* ff) ;
 
  virtual ~Domain_oned_ori() ; 
  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_cheb_base_odd(Base_spectral&) const ;       
     virtual void set_legendre_base_odd(Base_spectral&) const ;
     virtual void do_coloc () ;
     virtual int give_place_var (char*) const ;
     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 ;
   
     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 Val_domain der_partial_var (const Val_domain&, int) const ;
   virtual Val_domain div_r (const Val_domain&) 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 ;
 
     virtual double integrale(const Val_domain&) const ; 
 
public:
     virtual ostream& print (ostream& o) const ;    
} ;
 
class Domain_oned_qcq : public Domain {
 
 private:
  double alpha ; 
  double beta ; 
  
 public:  
  Domain_oned_qcq (int num, int ttype, double x_int, double x_ext, const Dim_array& nbr) ;
  Domain_oned_qcq (const Domain_oned_qcq& so) ; 
  Domain_oned_qcq (int num, FILE* ff) ;
 
  virtual ~Domain_oned_qcq() ;
  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_cheb_base_odd(Base_spectral&) const ;       
     virtual void set_legendre_base_odd(Base_spectral&) const ;
     
     virtual void do_coloc () ;
     virtual int give_place_var (char*) const ;
     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 void find_other_dom (int, int, int&, int&) const ;     
     virtual Val_domain der_normal (const Val_domain&, int) const ;
     virtual double val_boundary (int, const Val_domain&, const Index&) const ;
     virtual Val_domain der_partial_var (const Val_domain&, int) 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 ;
 
     virtual Val_domain div_xp1 (const Val_domain&) const ;
     virtual double integrale(const Val_domain&) const ; 
 
public:
     virtual ostream& print (ostream& o) const ;  
} ;
 
class Domain_oned_inf : public Domain {
 
 private:
  double alpha ; 
 
 public:
  Domain_oned_inf (int num, int ttype, double radius, const Dim_array& nbr) ;
  Domain_oned_inf (const Domain_oned_inf& so) ; 
  Domain_oned_inf (int num, FILE* ff) ;
 
  virtual ~Domain_oned_inf() ;
  virtual void save(FILE*) const ;
 
  private:
     virtual void set_cheb_base(Base_spectral&) const ;     
     virtual void set_legendre_base(Base_spectral&) const ;     
     virtual void set_cheb_base_odd(Base_spectral&) const ;       
     virtual void set_legendre_base_odd(Base_spectral&) const ;
     
     virtual void do_coloc()  ;     
 
     virtual void do_absol ()  const ; 
     virtual int give_place_var (char*) const ;
     virtual void do_radius () const ;
 
 public:    
     double get_alpha() const {return alpha ;} ;
   
     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 void set_val_inf (Val_domain& so, double xx) const ;    
 
     virtual void find_other_dom (int, int, int&, int&) const ;     
     virtual Val_domain der_normal (const Val_domain&, int) const ;
     virtual double val_boundary (int, const Val_domain&, const Index&) const ;
     virtual Val_domain der_partial_var (const Val_domain&, int) const ;
     virtual Val_domain mult_xm1 (const Val_domain&) const ;
     virtual Val_domain div_xm1 (const Val_domain&) const ;
     virtual Val_domain mult_x (const Val_domain&) const ;
     virtual Val_domain mult_r (const Val_domain& so) const {return mult_x(so) ;} ;
 
     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 ;
 
     virtual double integrale(const Val_domain&) const ; 
 
public:
     virtual ostream& print (ostream& o) const ;  
} ;
 
 
class Space_oned : public Space {
     public:
    Space_oned (int ttype, const Dim_array& nbr, const Array<double>& bounds) ;
    Space_oned (FILE*) ; 
    virtual ~Space_oned() ;        
    virtual void save(FILE*) const ;
 
    void add_eq_ori (System_of_eqs& syst, const char* eq) ;
} ;
}
#endif