diff options
author | Benoit Jacob <jacob.benoit.1@gmail.com> | 2009-08-31 22:09:44 -0400 |
---|---|---|
committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2009-08-31 22:09:44 -0400 |
commit | 29c6b2452dbe82cd49aa701921f2fa5a20017cc0 (patch) | |
tree | 5168dac3cc50ee2dc41442158d4d1a7e014e89c4 /test | |
parent | 5339db61645daa2c829898cbd2c3310b5f80f7a8 (diff) |
simplifications
Diffstat (limited to 'test')
-rw-r--r-- | test/svd.cpp | 6 |
1 files changed, 1 insertions, 5 deletions
diff --git a/test/svd.cpp b/test/svd.cpp index 2ccd94764..e6a32bd3f 100644 --- a/test/svd.cpp +++ b/test/svd.cpp @@ -41,15 +41,11 @@ template<typename MatrixType> void svd(const MatrixType& m) Matrix<Scalar, MatrixType::RowsAtCompileTime, 1>::Random(rows,1); Matrix<Scalar, MatrixType::ColsAtCompileTime, 1> x(cols,1), x2(cols,1); - RealScalar largerEps = test_precision<RealScalar>(); - if (ei_is_same_type<RealScalar,float>::ret) - largerEps = 1e-3f; - { SVD<MatrixType> svd(a); MatrixType sigma = MatrixType::Zero(rows,cols); MatrixType matU = MatrixType::Zero(rows,rows); - sigma.block(0,0,cols,cols) = svd.singularValues().asDiagonal(); + sigma.diagonal() = svd.singularValues(); matU = svd.matrixU(); VERIFY_IS_APPROX(a, matU * sigma * svd.matrixV().transpose()); } |