bin_bh.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 __BIN_BH_HPP_
#define __BIN_BH_HPP_
 
#include "space.hpp"
#include "spheric.hpp"
#include "homothetic.hpp"
#include "bispheric.hpp"
namespace Kadath  {
 
class Space_bin_bh : public Space {
  
  protected:
    double a_minus ; 
    double a_plus ; 
  
     public:  
 
    Space_bin_bh(int ttype, double dist, double rbh1, double rbh2, double rext, int nr) ;
    
    Space_bin_bh (FILE*) ; 
     
        virtual ~Space_bin_bh() ; 
    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 ;
    virtual Array<int> get_indices_matching_non_std(int dom, int bound) const ;
    
    void add_eq (System_of_eqs& syst, const char* eq, const char* rac, const char* rac_der, int nused=-1, Array<int>** pused=0x0)  ; 
    void add_bc_sphere_one (System_of_eqs& syst, const char* eq, int nused=-1, Array<int>** pused=0x0)  ;
    void add_bc_sphere_two (System_of_eqs& syst, const char* eq, int nused=-1, Array<int>** pused=0x0)  ;
    void add_bc_outer (System_of_eqs& syst, const char* eq, int nused=-1, Array<int>** pused=0x0)  ;
     
    void add_eq_int_inf (System_of_eqs& syst, const char* eq) ;
    void add_eq_int_sphere_one (System_of_eqs& syst, const char* eq) ;
    void add_eq_int_sphere_two (System_of_eqs& syst, const char* eq) ;
    void add_eq_zero_mode_inf (System_of_eqs& syst, const char* f, int jtarget, int ktarget) ;
} ;
}
#endif