Adds SolverOptions structure and SetSolverOptions method.

Author: Edward Palmer

src/problem_operators/problem_operator_base.cpp

@@ -12,16 +12,25 @@ ProblemOperatorBase::ConstructJacobianSolver()
   precond->SetPrintLevel(GetGlobalPrintLevel());
 
   auto solver = std::make_unique<mfem::HypreGMRES>(_problem._comm);
-
-  solver->SetTol(1e-16);
-  solver->SetAbsTol(1e-16);
-  solver->SetMaxIter(1000);
-  solver->SetKDim(10);
-  solver->SetPrintLevel(GetGlobalPrintLevel());
   solver->SetPreconditioner(*precond);
 
   _jacobian_preconditioner = std::move(precond);
   _jacobian_solver = std::move(solver);
+
+  SolverOptions default_options;
+  SetSolverOptions(default_options);
+}
+
+void
+ProblemOperatorBase::SetSolverOptions(SolverOptions & options)
+{
+  auto & solver = static_cast<mfem::HypreGMRES &>(*_jacobian_solver);
+
+  solver.SetTol(options._tolerance);
+  solver.SetAbsTol(options._abs_tolerance);
+  solver.SetMaxIter(options._max_iteration);
+  solver.SetKDim(options._k_dim);
+  solver.SetPrintLevel(options._print_level);
 }
 
 void

src/problem_operators/problem_operator_base.hpp

@@ -26,6 +26,21 @@ class ProblemOperatorBase
     return static_cast<TSolver *>(_jacobian_preconditioner.get());
   }
 
+  /// A structure for setting solver options. Defaults have been set.
+  struct SolverOptions
+  {
+    double _tolerance{1e-16};
+    double _abs_tolerance{1e-16};
+
+    unsigned int _max_iteration{1000};
+
+    int _print_level{GetGlobalPrintLevel()};
+    int _k_dim{10};
+  };
+
+  /// Sets the solver's options. Override in derived classes.
+  virtual void SetSolverOptions(SolverOptions & options);
+
 protected:
   /// Use of protected constructor to only allow construction by derived classes.
   /// All problem operator classes are built on-top of this class and it should not