Merge pull request #3668 from sconde/rythmos-get-solver

Piro: add getter to access nonlinear solver
trilinos · Oct 29, 2018 · a68b12e · a68b12e
2 parents 09eeed6 + 2da3a5b
commit a68b12e
Show file tree

Hide file tree

Showing 2 changed files with 16 additions and 0 deletions.
diff --git a/packages/piro/src/Piro_RythmosSolver.hpp b/packages/piro/src/Piro_RythmosSolver.hpp
@@ -119,6 +119,8 @@ class RythmosSolver
 
   Teuchos::RCP<const Rythmos::IntegratorBase<Scalar> > getRythmosIntegrator() const;
 
+  Teuchos::RCP<const Thyra::NonlinearSolverBase<Scalar> > getTimeStepSolver() const;
+
   /** \name Overridden from Thyra::ModelEvaluatorBase. */
   //@{
   /** \brief . */

diff --git a/packages/piro/src/Piro_RythmosSolver_Def.hpp b/packages/piro/src/Piro_RythmosSolver_Def.hpp
@@ -530,6 +530,20 @@ Piro::RythmosSolver<Scalar>::getRythmosIntegrator() const
   return fwdStateIntegrator;
 }
 
+
+#ifdef ALBANY_BUILD
+template <typename Scalar, typename LocalOrdinal, typename GlobalOrdinal, typename Node>
+Teuchos::RCP<const Thyra::NonlinearSolverBase<Scalar> >
+Piro::RythmosSolver<Scalar, LocalOrdinal, GlobalOrdinal, Node>::getTimeStepSolver() const
+#else
+template <typename Scalar>
+Teuchos::RCP<const Thyra::NonlinearSolverBase<Scalar> >
+Piro::RythmosSolver<Scalar>::getTimeStepSolver() const
+#endif
+{
+  return fwdTimeStepSolver;
+}
+
 #ifdef ALBANY_BUILD
 template <typename Scalar, typename LocalOrdinal, typename GlobalOrdinal, typename Node>
 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >