EquationInterface.hpp

#ifndef EQUATION_EQUATIONINTERFACE_H
#define EQUATION_EQUATIONINTERFACE_H


#include <fem/FemKernel.hpp>


namespace imaging
{
  template <class fem_types>
  class EquationInterface
  {
    
  public:
    std::size_t system_size() const;
    
    float_t stiffness_matrix(std::size_t k, std::size_t l,
                             std::size_t i, std::size_t j,
                             std::size_t integrator_node,
                             const FemKernel<fem_types> & kernel) const;
    
    float_t force_vector(std::size_t k,
                         std::size_t i,
                         std::size_t integrator_node,
                         const FemKernel<fem_types> & kernel) const;
  
    float_t stiffness_matrix_at_boundary(std::size_t k, std::size_t l,
                                         std::size_t i, std::size_t j,
                                         std::size_t integrator_node,
                                         const FemKernel<fem_types> & kernel) const;
        
    float_t force_vector_at_boundary(std::size_t k,
                                     std::size_t i,
                                     std::size_t integrator_node,
                                     const FemKernel<fem_types> & kernel) const;
        
    bool sanity_check_stiffnesss_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

---