From: young Date: Thu, 3 Oct 2013 13:38:48 +0000 (+0000) Subject: Set default problem type for SLEPc - this should be the GNHEP. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=3a10add12e7b038fa1520b40e460b17f7cf0d4ed;p=dealii-svn.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) {