From f93654ae16b261e462ee00c5255072f8dd7d387b Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Mon, 26 Oct 2015 16:00:25 +0100
Subject: bug #1098: fix regression introduced when generalizing some compute()
 methods in changeset 7031a851d45a8526474ac1ac972ad12a48e99f1a .

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

(limited to 'test/eigensolver_generalized_real.cpp')

diff --git a/test/eigensolver_generalized_real.cpp b/test/eigensolver_generalized_real.cpp
index 566a4bdc6..a46a2e50e 100644
--- a/test/eigensolver_generalized_real.cpp
+++ b/test/eigensolver_generalized_real.cpp
@@ -39,6 +39,14 @@ template<typename MatrixType> void generalized_eigensolver_real(const MatrixType
   VectorType realEigenvalues = eig.eigenvalues().real();
   std::sort(realEigenvalues.data(), realEigenvalues.data()+realEigenvalues.size());
   VERIFY_IS_APPROX(realEigenvalues, symmEig.eigenvalues());
+
+  // regression test for bug 1098
+  {
+    GeneralizedSelfAdjointEigenSolver<MatrixType> eig1(a.adjoint() * a,b.adjoint() * b);
+    eig1.compute(a.adjoint() * a,b.adjoint() * b);
+    GeneralizedEigenSolver<MatrixType> eig2(a.adjoint() * a,b.adjoint() * b);
+    eig2.compute(a.adjoint() * a,b.adjoint() * b);
+  }
 }
 
 void test_eigensolver_generalized_real()
-- 
cgit v1.2.3