fourD_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 __fourD_PERIODIC_HPP_
#define __fourD_PERIODIC_HPP_
 
#include "space.hpp"
 
#define TO_PI 0
namespace Kadath {
 
class Domain_fourD_periodic_nucleus : public Domain {
 
 private:
  double alpha ; 
 
  double ome ; 
 
  int type_time ;
  double maxt ; 
 
 public: 
  Domain_fourD_periodic_nucleus (int num, int ttype, double radius, double ome, const Dim_array& nbr) ;
  Domain_fourD_periodic_nucleus (const Domain_fourD_periodic_nucleus& so) ; 
  Domain_fourD_periodic_nucleus (int num, FILE* ff) ;
 
  virtual ~Domain_fourD_periodic_nucleus() ;
  virtual void save (FILE*) const ;
 
  private:    
    virtual void do_absol ()  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 ;} ;
 
     virtual Val_domain mult_cos_phi (const Val_domain&) const ;
     virtual Val_domain mult_sin_phi (const Val_domain&) const ;
     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 ;
 
 
     Val_domain translate (const Val_domain& so) const ;
     
     virtual Tensor change_basis_spher_to_cart (int dd, const Tensor&) const ;
public:
     virtual ostream& print (ostream& o) const ;
} ;
 
 
class Domain_fourD_periodic_shell : public Domain {
 
 private:
  double alpha ; 
  double beta ; 
 
 
  double ome ; 
 
  int type_time ;
  double maxt ; 
 
 public: 
  Domain_fourD_periodic_shell (int num, int ttype, double rin, double rout, double ome, const Dim_array& nbr) ;
  Domain_fourD_periodic_shell (const Domain_fourD_periodic_shell& so) ; 
  Domain_fourD_periodic_shell (int num, FILE* ff) ;
 
  virtual ~Domain_fourD_periodic_shell() ;
  virtual void save (FILE*) const ;
 
  private:    
     virtual void do_absol ()  const ;
     virtual void do_coloc () ;
     virtual void do_radius () const ;
     virtual void do_cart ()  const ;
 
     virtual void set_cheb_base_x_cart(Base_spectral& so) const ;
     virtual void set_cheb_base_y_cart(Base_spectral& so) const ;
     virtual void set_legendre_base_x_cart(Base_spectral& so) const ;
     virtual void set_legendre_base_y_cart(Base_spectral& so) 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 Val_domain mult_cos_phi (const Val_domain&) const ;
     virtual Val_domain mult_sin_phi (const Val_domain&) const ;
     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_xm1 (const Val_domain&) const ;    
 
     Val_domain translate (const Val_domain& so) const ;
     virtual Tensor change_basis_spher_to_cart (int dd, const Tensor&) const ;
 
     Array<double> integ_outer_boundary (const Val_domain& so) const ;
public:
     virtual ostream& print (ostream& o) const ;
} ;
 
class Space_fourD_periodic : public Space {
     public:
    Space_fourD_periodic (int tt, double omega, const Dim_array& nbr, const Array<double>& bounds) ;
    
    Space_fourD_periodic (FILE*) ; 
    virtual ~Space_fourD_periodic() ;      
    virtual void save(FILE*) const ;
    double get_omega() const ;
} ;
 
}
#endif