From 804a239d308bea55722c59b70c95459336230488 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Fri, 10 Apr 2009 19:54:43 +0000
Subject: patch from Moritz Lenz to allow solving transposed problem with
 superlu

---
 test/sparse_solvers.cpp | 8 ++++++++
 1 file changed, 8 insertions(+)

(limited to 'test/sparse_solvers.cpp')

diff --git a/test/sparse_solvers.cpp b/test/sparse_solvers.cpp
index d1090dfed..e1ec1ef35 100644
--- a/test/sparse_solvers.cpp
+++ b/test/sparse_solvers.cpp
@@ -191,6 +191,14 @@ template<typename Scalar> void sparse_solvers(int rows, int cols)
           VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LU: SuperLU");
         }
         // std::cerr << refDet << " == " << slu.determinant() << "\n";
+        if (slu.solve(b, &x, SvTranspose)) {
+          VERIFY(b.isApprox(m2.transpose() * x, test_precision<Scalar>()));
+        }
+
+        if (slu.solve(b, &x, SvAdjoint)) {
+//          VERIFY(b.isApprox(m2.adjoint() * x, test_precision<Scalar>()));
+        }
+
         if (count==0) {
           VERIFY_IS_APPROX(refDet,slu.determinant()); // FIXME det is not very stable for complex
         }
-- 
cgit v1.2.3