From 2da3a5bf4bcaa83229e450e3faacad96076a1216 Mon Sep 17 00:00:00 2001
From: Sidafa Conde <sconde@sandia.gov>
Date: Thu, 18 Oct 2018 14:27:32 -0600
Subject: [PATCH] Piro: add getter to access nonlinear solver

---
 packages/piro/src/Piro_RythmosSolver.hpp     |  2 ++
 packages/piro/src/Piro_RythmosSolver_Def.hpp | 14 ++++++++++++++
 2 files changed, 16 insertions(+)

diff --git a/packages/piro/src/Piro_RythmosSolver.hpp b/packages/piro/src/Piro_RythmosSolver.hpp
index f37115a50d94..f43050ee359b 100644
--- 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
index 38f3c5c970e1..4b9cbedb1844 100644
--- 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> >