HyperReduction.h

#ifndef PBAT_SIM_VBD_MULTIGRID_HYPERREDUCTION_H
#define PBAT_SIM_VBD_MULTIGRID_HYPERREDUCTION_H
 
#include "pbat/Aliases.h"
 
#include <vector>
 
namespace pbat {
namespace sim {
namespace vbd {
namespace multigrid {
 
struct Hierarchy;
 
struct HyperReduction
{
    static constexpr auto kPolynomialOrder = 1;
    static constexpr auto kDims            = 3;
 
    HyperReduction() = default;
    HyperReduction(Hierarchy const& hierarchy, Index clusterSize = 5);
    void Construct(Hierarchy const& hierarchy, Index clusterSize = 5);
    void AllocateWorkspace(std::size_t nLevels);
    void ConstructHierarchicalClustering(Hierarchy const& hierarchy, Index clusterSize);
    void ComputeClusterQuadratureWeights(Hierarchy const& hierarchy);
    void PrecomputeInversePolynomialMatrices(Hierarchy const& hierarchy);
    void
    ComputeLinearPolynomialErrors(Hierarchy const& hierarchy, Eigen::Ref<MatrixX const> const& u);

    std::vector<IndexVectorX> C; 
    std::vector<IndexVectorX> Cptr;
    std::vector<IndexVectorX> Cadj;
 
    std::vector<MatrixX>
        ApInvC; 
    std::vector<MatrixX> bC;      
    std::vector<VectorX> wC;      
    std::vector<IndexVectorX> eC; 
    std::vector<MatrixX> up;      
    std::vector<VectorX> Ep;      
    Scalar EpMax{1e-6};           
};
 
} // namespace multigrid
} // namespace vbd
} // namespace sim
} // namespace pbat
 
#endif // PBAT_SIM_VBD_MULTIGRID_HYPERREDUCTION_H