ExplicitVectorOperators_fwd.hpp

Go to the documentation of this file.

#ifndef MEDUSA_BITS_OPERATORS_EXPLICITVECTOROPERATORS_FWD_HPP_
#define MEDUSA_BITS_OPERATORS_EXPLICITVECTOROPERATORS_FWD_HPP_
 
#include <medusa/Config.hpp>
#include <Eigen/Core>
#include <medusa/bits/types/traits.hpp>
#include <medusa/bits/approximations/Operators_fwd.hpp>
 
namespace mm {
 
template <class shape_storage_type>
class ExplicitVectorOperators {
  public:
    typedef shape_storage_type shape_storage_t;   
    typedef typename shape_storage_t::vector_t vector_t;  
    typedef typename shape_storage_t::scalar_t scalar_t;  
    enum {  dim = shape_storage_t::dim };
 
  private:
    const shape_storage_t* ss;
 
  public:
    ExplicitVectorOperators() : ss(nullptr) {}
    explicit ExplicitVectorOperators(const shape_storage_t& ss) : ss(&ss) {}
 
    void setShapes(const shape_storage_t& shape_storage_) { ss = &shape_storage_; }
 
    bool hasShapes() const { return ss != nullptr; }
 
    template <class vector_field_t>
    typename vector_type<vector_field_t>::type lap(const vector_field_t& u, int node) const {
        return apply<Lap<dim>, vector_field_t>(u, node, {});
    }
 
    template <typename op_family_t, typename vector_field_t>
    typename vector_type<vector_field_t>::type apply(
            const vector_field_t& u, int node, typename op_family_t::operator_t o) const;
 
 
    template <typename op_family_t, typename vector_field_t>
    typename vector_type<vector_field_t>::type apply(const vector_field_t& u, int node) const {
        return apply<op_family_t, vector_field_t>(u, node, {});
    }
 
 
    template <class vector_field_t>
    Eigen::Matrix<typename vector_type<vector_field_t>::type::scalar_t,
                  vector_type<vector_field_t>::type::dim, dim>
    grad(const vector_field_t& u, int node) const;
 
    template <class vector_field_t>
    typename vector_type<vector_field_t>::type::scalar_t div(
            const vector_field_t& u, int node) const;
 
     template <class vector_field_t>
     typename vector_type<vector_field_t>::type curl(const vector_field_t& u, int node) const;
 
    template <class vector_field_t>
    typename vector_type<vector_field_t>::type graddiv(const vector_field_t& u, int node) const;
    template <class vector_field_t>
    typename vector_type<vector_field_t>::type neumann(const vector_field_t& u, int node,
            const vector_t& normal, typename vector_type<vector_field_t>::type val) const;
 
    template <typename S>
    friend std::ostream& operator<<(std::ostream& os, const ExplicitVectorOperators<S>& op);
};
 
}  // namespace mm
 
#endif  // MEDUSA_BITS_OPERATORS_EXPLICITVECTOROPERATORS_FWD_HPP_