Tempus Version of the Day
Time Integration
Loading...
Searching...
No Matches
Tempus_StepperSubcycling_decl.hpp
Go to the documentation of this file.
1//@HEADER
2// *****************************************************************************
3// Tempus: Time Integration and Sensitivity Analysis Package
4//
5// Copyright 2017 NTESS and the Tempus contributors.
6// SPDX-License-Identifier: BSD-3-Clause
7// *****************************************************************************
8//@HEADER
9
10#ifndef Tempus_StepperSubcycling_decl_hpp
11#define Tempus_StepperSubcycling_decl_hpp
12
13#include "Tempus_config.hpp"
14#include "Tempus_StepperExplicit.hpp"
16#include "Tempus_IntegratorBasic.hpp"
17
18namespace Tempus {
19
52template <class Scalar>
53class StepperSubcycling : virtual public Tempus::Stepper<Scalar> {
54 public:
61
64 const Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >& appModel,
65 const Teuchos::RCP<IntegratorBasic<Scalar> >& integrator, bool useFSAL,
66 std::string ICConsistency, bool ICConsistencyCheck,
67 const Teuchos::RCP<StepperSubcyclingAppAction<Scalar> >&
68 stepperSCAppAction);
69
71
72 virtual void setModel(
73 const Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >& appModel);
74
75 virtual void setNonConstModel(
76 const Teuchos::RCP<Thyra::ModelEvaluator<Scalar> >& appModel);
77
78 virtual Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > getModel() const
79 {
80 return scIntegrator_->getStepper()->getModel();
81 }
82
83 virtual void setAppAction(Teuchos::RCP<StepperSubcyclingAppAction<Scalar> >
84 appAction = Teuchos::null);
85
86 virtual Teuchos::RCP<StepperSubcyclingAppAction<Scalar> > getAppAction() const
87 {
89 }
90
92 virtual void initialize();
93
94 virtual Scalar getInitTimeStep(
95 const Teuchos::RCP<SolutionHistory<Scalar> >& solutionHistory) const;
96
98 virtual void setInitialConditions(
99 const Teuchos::RCP<SolutionHistory<Scalar> >& solutionHistory);
100
102 // thus a no-op for explicit steppers.
103 virtual void setInitialGuess(
104 Teuchos::RCP<const Thyra::VectorBase<Scalar> > initial_guess);
105
106 virtual void setSolver(
107 Teuchos::RCP<Thyra::NonlinearSolverBase<Scalar> > solver = Teuchos::null);
108
109 virtual Teuchos::RCP<Thyra::NonlinearSolverBase<Scalar> > getSolver() const;
110
112 virtual void takeStep(
113 const Teuchos::RCP<SolutionHistory<Scalar> >& solutionHistory);
114
116 virtual Teuchos::RCP<Tempus::StepperState<Scalar> > getDefaultStepperState();
117
118 virtual bool isExplicit() const;
119 virtual bool isImplicit() const;
120 virtual bool isExplicitImplicit() const;
121 virtual bool isOneStepMethod() const;
122 virtual bool isMultiStepMethod() const;
123
124 virtual Scalar getOrder() const;
125 virtual Scalar getOrderMin() const;
126 virtual Scalar getOrderMax() const;
127 virtual OrderODE getOrderODE() const;
129
130 Teuchos::RCP<const Teuchos::ParameterList> getValidParameters() const;
131
133
134 virtual void describe(Teuchos::FancyOStream& out,
135 const Teuchos::EVerbosityLevel verbLevel) const;
137
139
140 virtual void setSubcyclingStepper(Teuchos::RCP<Stepper<Scalar> > stepper);
141 virtual void setSubcyclingMinTimeStep(Scalar MinTimeStep);
142 virtual void setSubcyclingInitTimeStep(Scalar InitTimeStep);
143 virtual void setSubcyclingMaxTimeStep(Scalar MaxTimeStep);
144 virtual void setSubcyclingMaxFailures(int MaxFailures);
145 virtual void setSubcyclingMaxConsecFailures(int MaxConsecFailures);
146 virtual void setSubcyclingScreenOutputIndexInterval(int i);
147 virtual void setSubcyclingScreenOutputIndexList(std::string s);
149 Teuchos::RCP<TimeStepControlStrategy<Scalar> > tscs);
151 Teuchos::RCP<IntegratorObserver<Scalar> > obs);
152 virtual void setSubcyclingPrintDtChanges(bool printDtChanges);
154
156
157 virtual Teuchos::RCP<const Stepper<Scalar> > getSubcyclingStepper() const;
158 virtual Scalar getSubcyclingMinTimeStep() const;
159 virtual Scalar getSubcyclingInitTimeStep() const;
160 virtual Scalar getSubcyclingMaxTimeStep() const;
161 virtual std::string getSubcyclingStepType() const;
162 virtual int getSubcyclingMaxFailures() const;
163 virtual int getSubcyclingMaxConsecFailures() const;
164 virtual int getSubcyclingScreenOutputIndexInterval() const;
165 virtual std::string getSubcyclingScreenOutputIndexList() const;
166 virtual Teuchos::RCP<TimeStepControlStrategy<Scalar> >
168 virtual Teuchos::RCP<IntegratorObserver<Scalar> >
170 virtual bool getSubcyclingPrintDtChanges() const;
172
173 protected:
174 Teuchos::RCP<StepperSubcyclingAppAction<Scalar> > stepperSCAppAction_;
175 Teuchos::RCP<IntegratorBasic<Scalar> > scIntegrator_;
176};
177
179// ------------------------------------------------------------------------
180template <class Scalar>
181Teuchos::RCP<StepperSubcycling<Scalar> > createStepperSubcycling(
182 const Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >& model,
183 Teuchos::RCP<Teuchos::ParameterList> pl);
184
185} // namespace Tempus
186
187#endif // Tempus_StepperSubcycling_decl_hpp
IntegratorObserver class for time integrators.
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
Application Action for StepperSubcycling.
virtual void setSubcyclingStepper(Teuchos::RCP< Stepper< Scalar > > stepper)
virtual Teuchos::RCP< IntegratorObserver< Scalar > > getSubcyclingIntegratorObserver() const
virtual Scalar getInitTimeStep(const Teuchos::RCP< SolutionHistory< Scalar > > &solutionHistory) const
virtual void setSubcyclingPrintDtChanges(bool printDtChanges)
virtual Teuchos::RCP< StepperSubcyclingAppAction< Scalar > > getAppAction() const
Teuchos::RCP< StepperSubcyclingAppAction< Scalar > > stepperSCAppAction_
virtual Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getModel() const
virtual void initialize()
Initialize during construction and after changing input parameters.
virtual Teuchos::RCP< const Stepper< Scalar > > getSubcyclingStepper() const
virtual int getSubcyclingScreenOutputIndexInterval() const
virtual void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
virtual void setAppAction(Teuchos::RCP< StepperSubcyclingAppAction< Scalar > > appAction=Teuchos::null)
virtual void setSolver(Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > solver=Teuchos::null)
Set solver.
virtual void setSubcyclingMaxTimeStep(Scalar MaxTimeStep)
virtual void setModel(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &appModel)
virtual void setNonConstModel(const Teuchos::RCP< Thyra::ModelEvaluator< Scalar > > &appModel)
virtual void setSubcyclingIntegratorObserver(Teuchos::RCP< IntegratorObserver< Scalar > > obs)
virtual void setSubcyclingInitTimeStep(Scalar InitTimeStep)
Teuchos::RCP< IntegratorBasic< Scalar > > scIntegrator_
virtual std::string getSubcyclingStepType() const
virtual void setInitialGuess(Teuchos::RCP< const Thyra::VectorBase< Scalar > > initial_guess)
Pass initial guess to Newton solver (only relevant for implicit solvers)
virtual void setSubcyclingMinTimeStep(Scalar MinTimeStep)
virtual Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > getSolver() const
Get solver.
virtual void takeStep(const Teuchos::RCP< SolutionHistory< Scalar > > &solutionHistory)
Take the specified timestep, dt, and return true if successful.
virtual void setSubcyclingScreenOutputIndexInterval(int i)
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
virtual Teuchos::RCP< TimeStepControlStrategy< Scalar > > getSubcyclingTimeStepControlStrategy() const
virtual std::string getSubcyclingScreenOutputIndexList() const
virtual void setSubcyclingMaxFailures(int MaxFailures)
virtual void setInitialConditions(const Teuchos::RCP< SolutionHistory< Scalar > > &solutionHistory)
Set the initial conditions, make them consistent, and set needed memory.
virtual Teuchos::RCP< Tempus::StepperState< Scalar > > getDefaultStepperState()
Get a default (initial) StepperState.
virtual void setSubcyclingMaxConsecFailures(int MaxConsecFailures)
virtual void setSubcyclingScreenOutputIndexList(std::string s)
virtual void setSubcyclingTimeStepControlStrategy(Teuchos::RCP< TimeStepControlStrategy< Scalar > > tscs)
Thyra Base interface for time steppers.
TimeStepControlStrategy class for TimeStepControl.
Teuchos::RCP< StepperSubcycling< Scalar > > createStepperSubcycling(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &model, Teuchos::RCP< Teuchos::ParameterList > pl)
Nonmember constructor - ModelEvaluator and ParameterList.