critic.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 __CRITIC_HPP_
#define __CRITIC_HPP_
 
#include "space.hpp"
 
namespace Kadath {
 
class Domain_critic_inner : public Domain {
 
  protected:
    mutable Val_domain* p_X ; 
    mutable Val_domain* p_T ; 
    double xlim ; 
 
 public:  
  Domain_critic_inner (int num, int ttype, const Dim_array& nbr, double xl) ;
  Domain_critic_inner (const Domain_critic_inner & so) ; 
  Domain_critic_inner (int num, FILE* ff) ;
 
  virtual ~Domain_critic_inner() ;
  void del_deriv() override;
  virtual void save (FILE*) const ;
   
  private : 
     virtual void set_cheb_base(Base_spectral& so) const ;
     virtual void set_legendre_base(Base_spectral&) const ;
 
     virtual void set_cheb_xodd_base(Base_spectral& so) const ;
     virtual void set_legendre_xodd_base(Base_spectral&) const ;     
     virtual void set_cheb_todd_base(Base_spectral& so) const ;
     virtual void set_legendre_todd_base(Base_spectral&) const ;
     virtual void set_cheb_xodd_todd_base(Base_spectral& so) const ;
     virtual void set_legendre_xodd_todd_base(Base_spectral&) const ;
 
     virtual void do_coloc () ;
     virtual int give_place_var (char*) const ;
     void do_X() const ; 
     void do_T() const ; 
 
   public:
     virtual const Val_domain & get_X() const ;
     virtual const Val_domain & get_T() const ;
     virtual Val_domain div_x (const Val_domain&) const ;     
     
     virtual bool is_in(const Point& xx, double prec=1e-13) const ;
     virtual const Point absol_to_num(const Point& xxx) const ;
     virtual void do_absol () const ;
 
  public:  
     virtual void find_other_dom (int, int, int&, int&) const ;
     virtual Base_spectral mult (const Base_spectral&, const Base_spectral&) 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 int nbr_points_boundary (int, const Base_spectral&) const ;
     virtual void do_which_points_boundary (int, const Base_spectral&, Index**, int) const ;
 
public:
     virtual ostream& print (ostream& o) const ;
} ;
 
 
class Domain_critic_outer : public Domain {
 
  protected:
    mutable Val_domain* p_X ; 
    mutable Val_domain* p_T ; 
    double xlim ; 
 
 public:  
  Domain_critic_outer (int num, int ttype, const Dim_array& nbr, double xl) ;
  Domain_critic_outer (const Domain_critic_outer & so) ; 
  Domain_critic_outer (int num, FILE* ff) ;
 
  virtual ~Domain_critic_outer() ;   
  void del_deriv() override;
  virtual void save (FILE*) const ;
   
  private : 
     virtual void set_cheb_base(Base_spectral& so) const ;
     virtual void set_legendre_base(Base_spectral&) const ;
     virtual void set_cheb_todd_base(Base_spectral& so) const ;
     virtual void set_legendre_todd_base(Base_spectral&) const ;
 
     virtual void do_coloc () ;
     virtual int give_place_var (char*) const ;
     void do_X() const ; 
     void do_T() const ;
 
   public:
     virtual const Val_domain & get_X() const ;
     virtual const Val_domain & get_T() const ;
     
     virtual bool is_in(const Point& xx, double prec=1e-13) const ;
     virtual const Point absol_to_num(const Point& xxx) const ;
     virtual void do_absol () const ;
 
  public:  
     virtual void find_other_dom (int, int, int&, int&) const ;
     virtual Base_spectral mult (const Base_spectral&, const Base_spectral&) 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 int nbr_points_boundary (int, const Base_spectral&) const ;
     virtual void do_which_points_boundary (int, const Base_spectral&, Index**, int) const ;
 
     //friend ostream& operator<< (ostream& o, const Domain_critic_outer& so) ; ///< Display
public:
     virtual ostream& print (ostream& o) const ;
} ;
 
class Space_critic : public Space {
     public:
    Space_critic (int ttype, double xl, const Dim_array& nbr_inner, const Dim_array& nbr_outer) ;
    Space_critic (FILE*) ; 
    virtual ~Space_critic() ;       
    virtual void save(FILE*) const ;
 
    void add_bc_zero (System_of_eqs& syst, const char* eq, int nused=-1, Array<int>** pused=0x0)  ;
    void add_bc_one (System_of_eqs& syst, const char* eq, int nused=-1, Array<int>** pused=0x0)  ;
    void add_eq_zero_mode_center (System_of_eqs& syst, const char* f, int jtarget) ;
         
    virtual Array<int> get_indices_matching_non_std(int, int) const ;
} ;
}
#endif