Level.h

#ifndef PBAT_SIM_VBD_MULTIGRID_LEVEL_H
#define PBAT_SIM_VBD_MULTIGRID_LEVEL_H
 
#include "HyperReduction.h"
#include "pbat/Aliases.h"
#include "pbat/sim/vbd/Data.h"
#include "pbat/sim/vbd/Mesh.h"
 
namespace pbat {
namespace sim {
namespace vbd {
namespace multigrid {
 
struct Level
{
    Level(Data const& data, VolumeMesh mesh);
    void Prolong(Data& data) const;
    void Smooth(Scalar dt, Index iters, Data& data);
    void Reduce(Data const& data);

    VolumeMesh mesh;         
    MatrixX u;               
    IndexVectorX colors;     
    IndexVectorX Pptr, Padj; 
 
    using BoolVector = Eigen::Vector<bool, Eigen::Dynamic>;

    IndexMatrixX ecVE; 
    MatrixX
        NecVE; 
    IndexMatrixX ilocalE;      
    IndexVectorX GEptr, GEadj; 

    IndexVectorX ecK; 
    MatrixX NecK;     
    IndexVectorX GKptr, GKadj, GKilocal; 

    BoolVector bIsDirichletVertex;

    IndexVectorX GCFptr, GCFadj, GCFrank, GCFparent; 
    // HyperReduction HR; ///< Energy hyper reduction scheme
};
 
} // namespace multigrid
} // namespace vbd
} // namespace sim
} // namespace pbat
 
#endif // PBAT_SIM_VBD_MULTIGRID_LEVEL_H