HKS.h

#pragma once

#include <Eigen/Core>
#include <Euclid/MeshUtil/MeshDefs.h>

namespace Euclid
{
template<typename Mesh>
class HKS
{
public:
    using FT = FT_t<Mesh>;
    using Vec = Eigen::Matrix<FT, Eigen::Dynamic, 1>;
    using Mat = Eigen::Matrix<FT, Eigen::Dynamic, Eigen::Dynamic>;

public:
    void build(const Mesh& mesh, unsigned k = 300);

    void build(const Mesh& mesh,
               const Vec* eigenvalues,
               const Mat* eigenfunctions);

    template<typename Derived>
    void compute(Eigen::ArrayBase<Derived>& hks,
                 unsigned tscales = 100,
                 float tmin = -1.0f,
                 float tmax = -1.0f);

private:
    const Mesh* _mesh;
    Mat _phi2;     // @f$\phi * \phi@f$.
    Vec _emlambda; // @f$e^{-\lambda}@f$.
    FT _lambda_max;
    FT _lambda_min;
};

} // namespace Euclid

#include "src/HKS.cpp"