test/approximations/JacobiSVDWrapper_test.cpp

#include "medusa/bits/approximations/JacobiSVDWrapper.hpp"
#include <Eigen/SVD>
 
#include "gtest/gtest.h"
 
namespace mm {
 
TEST(Approximations, JacobiSVDWrapper) {
    Eigen::MatrixXd M(7, 9);
    M.setRandom();
 
    JacobiSVDWrapper<double> svd;
    svd.compute(M);
 
    // This code is the save as above, with explicit parameter specification.
    Eigen::JacobiSVD<Eigen::MatrixXd> svd2;
    svd2.compute(M, Eigen::ComputeThinU | Eigen::ComputeThinV);
 
    // Use svd.matrixU(), svd.solve(b) as usual.
 
    EXPECT_TRUE((svd.matrixU() - svd2.matrixU()).isZero(0));
    EXPECT_TRUE((svd.matrixV() - svd2.matrixV()).isZero(0));
}
 
}  // namespace mm