SimpleEquationInterface.hpp

#ifndef EQUATION_SIMPLEEQUATIONINTERFACE_H
#define EQUATION_SIMPLEEQUATIONINTERFACE_H


#include <fem/FemKernel.hpp>


namespace imaging
{
  template<class fem_types_t>
  class SimpleEquationInterface
  {
  public:
    typedef fem_types_t fem_types;
  
    typedef ublas::fixed_matrix<float_t, fem_types::data_dimension, fem_types::data_dimension> matrix_coefficient_t;

    enum boundary_data_types {
      NO_BOUNDARY_DATA ,
      IMPLICIT_NEUMANN_DATA ,
      NEUMANN_DATA ,
      DIRICHLET_DATA ,
      MIXED_DATA  
    };

    static const bool a_active = true;
    
    static const bool b_active = true;
    
    static const bool c_active = true;
    
    static const bool f_active = true;
    
    static const bool g_active = true;
    
    static const size_t boundary_data_type;  
    
    void stiffness_matrix(std::size_t integrator_node,
                          const FemKernel<fem_types> & kernel,
                          matrix_coefficient_t & A,
                          ublas::fixed_vector<float_t, fem_types::data_dimension> & a,
                          ublas::fixed_vector<float_t, fem_types::data_dimension> & b,
                          float_t & c) const;
      
    
    void force_vector(std::size_t integrator_node,
                      const FemKernel<fem_types> & kernel,
                      float_t & f,
                      ublas::fixed_vector<float_t, fem_types::data_dimension> & g) const; 
    
    void stiffness_matrix_at_boundary(std::size_t integrator_node,
                                      const FemKernel<fem_types> & kernel,
                                      float_t & h) const {}
    
    
    void force_vector_at_boundary(std::size_t integrator_node,
                                  const FemKernel<fem_types> & kernel,
                                  float_t & v) const {}
        
    bool sanity_check_stiffness_matrix(const FemKernel<fem_types> & kernel, std::string & error_message) const { return true; }  
        
    bool sanity_check_force_vector(const FemKernel<fem_types> & kernel, std::string & error_message) const { return true; }      
  };

}

#endif

---