diff --git a/test/jacobisvd.cpp b/test/jacobisvd.cpp
index f44e66c94..cf8c76055 100644
--- a/test/jacobisvd.cpp
+++ b/test/jacobisvd.cpp
@@ -132,6 +132,11 @@ void jacobisvd_test_all_computation_options(const MatrixType& m)
     jacobisvd_solve<MatrixType, QRPreconditioner>(m, ComputeThinU | ComputeFullV);
     jacobisvd_solve<MatrixType, QRPreconditioner>(m, ComputeThinU | ComputeThinV);
   }
+
+  {
+    JacobiSVD<MatrixType, QRPreconditioner> svd(m, ComputeThinU | ComputeThinV);
+    VERIFY_IS_APPROX(m, svd.matrixU() * svd.singularValues().asDiagonal() * svd.matrixV().transpose());
+  }
 }
 
 template<typename MatrixType>