aboutsummaryrefslogtreecommitdiffhomepage
path: root/test/svd_common.h
diff options
context:
space:
mode:
Diffstat (limited to 'test/svd_common.h')
-rw-r--r--test/svd_common.h6
1 files changed, 5 insertions, 1 deletions
diff --git a/test/svd_common.h b/test/svd_common.h
index e300336c6..d8611b541 100644
--- a/test/svd_common.h
+++ b/test/svd_common.h
@@ -33,6 +33,7 @@ void svd_check_full(const MatrixType& m, const SvdType& svd)
};
typedef typename MatrixType::Scalar Scalar;
+ typedef typename MatrixType::RealScalar RealScalar;
typedef Matrix<Scalar, RowsAtCompileTime, RowsAtCompileTime> MatrixUType;
typedef Matrix<Scalar, ColsAtCompileTime, ColsAtCompileTime> MatrixVType;
@@ -40,7 +41,10 @@ void svd_check_full(const MatrixType& m, const SvdType& svd)
sigma.diagonal() = svd.singularValues().template cast<Scalar>();
MatrixUType u = svd.matrixU();
MatrixVType v = svd.matrixV();
- VERIFY_IS_APPROX(m, u * sigma * v.adjoint());
+ RealScalar scaling = m.cwiseAbs().maxCoeff();
+ if(scaling<=(std::numeric_limits<RealScalar>::min)())
+ scaling = RealScalar(1);
+ VERIFY_IS_APPROX(m/scaling, u * (sigma/scaling) * v.adjoint());
VERIFY_IS_UNITARY(u);
VERIFY_IS_UNITARY(v);
}