PerSimplexFunction.hpp

Go to the documentation of this file.

#pragma once
 
#include <wmtk/Primitive.hpp>
 
#include <wmtk/attribute/MeshAttributeHandle.hpp>
 
#include <Eigen/Core>
 
namespace wmtk::function {
 
class PerSimplexFunction
{
public:
    PerSimplexFunction(
        const Mesh& mesh,
        const PrimitiveType primitive_type,
        const attribute::MeshAttributeHandle& variable_attribute_handle);
    virtual ~PerSimplexFunction() {}
 
    virtual double get_value(const simplex::Simplex& domain_simplex) const = 0;
    virtual Eigen::VectorXd get_gradient(
        const simplex::Simplex& domain_simplex,
        const simplex::Simplex& variable_simplex) const
    {
        throw std::runtime_error("Gradient not implemented");
    }
    virtual Eigen::MatrixXd get_hessian(
        const simplex::Simplex& domain_simplex,
        const simplex::Simplex& variable_simplex) const
    {
        throw std::runtime_error("Hessian not implemented");
    }
 
    inline const Mesh& mesh() const { return m_mesh; }
    inline const attribute::MeshAttributeHandle& attribute_handle() const
    {
        assert(m_handle.is_valid());
        return m_handle;
    }
 
    int64_t embedded_dimension() const;
 
private:
    attribute::MeshAttributeHandle m_handle;
    const Mesh& m_mesh;
 
protected:
    const PrimitiveType m_primitive_type;
};
} // namespace wmtk::function