From: Toby D. Young <tyoung@ippt.pan.pl>
Date: Thu, 3 Oct 2013 13:38:48 +0000 (+0000)
Subject: Set default problem type for SLEPc - this should be the GNHEP.
X-Git-Tag: v8.1.0~650
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=eae9da8784752981828587af9054e9c03f24d357;p=dealii.git

Set default problem type for SLEPc - this should be the GNHEP.

git-svn-id: https://svn.dealii.org/trunk@31084 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/source/lac/slepc_solver.cc b/deal.II/source/lac/slepc_solver.cc
index 4e372994d1..6ea1e27454 100644
--- a/deal.II/source/lac/slepc_solver.cc
+++ b/deal.II/source/lac/slepc_solver.cc
@@ -52,7 +52,7 @@ namespace SLEPcWrappers
     mpi_communicator (mpi_communicator),
     target_eigenvalue (0.),
     set_which (EPS_LARGEST_MAGNITUDE),
-    set_problem (EPS_NHEP),
+    set_problem (EPS_GNHEP),
     opA (NULL),
     opB (NULL),
     initial_vector (NULL),
@@ -111,7 +111,7 @@ namespace SLEPcWrappers
 
   void
   SolverBase::solve (const unsigned int  n_eigenpairs,
-                     unsigned int *n_converged)
+                     unsigned int       *n_converged)
   {
     int ierr;
 
@@ -137,6 +137,10 @@ namespace SLEPcWrappers
         set_solver_type (solver_data->eps);
       }
 
+    // set the problem type
+    ierr = EPSSetProblemType (solver_data->eps, set_problem);
+    AssertThrow (ierr == 0, ExcSLEPcError(ierr));
+
     // set the initial vector(s) if any
     if (initial_vector && initial_vector->size() != 0)
       {