From eae9da8784752981828587af9054e9c03f24d357 Mon Sep 17 00:00:00 2001 From: "Toby D. Young" Date: Thu, 3 Oct 2013 13:38:48 +0000 Subject: [PATCH] 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 --- deal.II/source/lac/slepc_solver.cc | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) 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) { -- 2.39.5