WKS.h

#pragma once

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

namespace Euclid
{
template<typename Mesh>
class WKS
{
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>& wks,
                 unsigned escales = 100,
                 float emin = 0.0f,
                 float emax = -1.0f,
                 float sigma = -1.0f);

private:
    const Mesh* _mesh;
    Mat _phi2;
    Vec _loglambda;
    FT _lambda_max;
    FT _lambda_min;
};

} // namespace Euclid

#include "src/WKS.cpp"