MockMatrix.hpp

Go to the documentation of this file.

#ifndef IRLBA_MOCK_MATRIX_HPP
#define IRLBA_MOCK_MATRIX_HPP
 
#include "Eigen/Dense"
 
namespace irlba {
 
class MockMatrix {
public:
    MockMatrix(Eigen::MatrixXd x) : my_x(std::move(x)) {}
private:
    Eigen::MatrixXd my_x;
 
public:
 
    Eigen::Index rows() const { return my_x.rows(); }
 
    Eigen::Index cols() const { return my_x.cols(); }
 
public:
    struct Workspace {};
 
    Workspace workspace() const {
        return Workspace();
    }
 
    struct AdjointWorkspace {};
 
    AdjointWorkspace adjoint_workspace() const {
        return AdjointWorkspace();
    }
 
public:
    template<class Right_, class EigenVector_>
    void multiply(const Right_& rhs, [[maybe_unused]] Workspace& work, EigenVector_& out) const {
        out.noalias() = my_x * rhs; 
    }
 
    template<class Right_, class EigenVector_>
    void adjoint_multiply(const Right_& rhs, [[maybe_unused]] AdjointWorkspace& work, EigenVector_& out) const {
        out.noalias() = my_x.adjoint() * rhs;
    }
 
    template<class EigenMatrix_>
    EigenMatrix_ realize() const {
        return EigenMatrix_(my_x);
    }
};
 
}
 
#endif