]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Set default problem type for SLEPc - this should be the GNHEP.
authoryoung <young@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 3 Oct 2013 13:38:48 +0000 (13:38 +0000)
committeryoung <young@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 3 Oct 2013 13:38:48 +0000 (13:38 +0000)
git-svn-id: https://svn.dealii.org/trunk@31084 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/lac/slepc_solver.cc

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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.