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 4e372994d1750ada63afb4579ccbbd8274dc881d..6ea1e27454f43ac1c3e7664d48d6cfa9b5a96f13 100644 (file)
--- 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)
       {