LUSolver.h

#pragma once
 
#include <SofaCaribou/Solver/EigenSparseSolver.h>
#include <sofa/helper/OptionsGroup.h>
 
namespace SofaCaribou::solver {
 
template <class EigenSolver>
class LUSolver : public EigenSparseSolver<EigenSolver> {
public:
    SOFA_CLASS(SOFA_TEMPLATE(LUSolver, EigenSolver), SOFA_TEMPLATE(EigenSparseSolver, EigenSolver));
 
    template <typename T>
    using Data = sofa::Data<T>;
 
    LUSolver();
 
    inline bool symmetric() const override {return d_is_symmetric.getValue();}
 
    inline void set_symmetric(bool is_symmetric) override { d_is_symmetric.setValue(is_symmetric); }
 
    void assemble (const sofa::core::MechanicalParams* mparams) override;
 
    // Get the backend name of the class derived from the EigenSolver template parameter
    static std::string BackendName();
private:
    Data<sofa::helper::OptionsGroup> d_backend;
 
    Data<bool> d_is_symmetric;
};
 
} // namespace SofaCaribou::solver