Tempus Version of the Day
Time Integration
Loading...
Searching...
No Matches
Tempus_StepperSubcyclingModifierXBase.hpp
Go to the documentation of this file.
1// @HEADER
2// ****************************************************************************
3// Tempus: Copyright (2017) Sandia Corporation
4//
5// Distributed under BSD 3-clause license (See accompanying file Copyright.txt)
6// ****************************************************************************
7// @HEADER
8
9#ifndef Tempus_StepperSubcyclingModifierXBase_hpp
10#define Tempus_StepperSubcyclingModifierXBase_hpp
11
12#include "Tempus_config.hpp"
13#include "Tempus_SolutionHistory.hpp"
15
16
17namespace Tempus {
18
37template<class Scalar>
39 : virtual public Tempus::StepperSubcyclingAppAction<Scalar>
40{
41private:
42
43 /* \brief Adaptor execute function
44 *
45 * This is an adaptor function to bridge between the AppAction
46 * interface and the ModifierX interface. It is meant to be private
47 * and non-virtual as deriving from this class should only need to
48 * implement the modify function.
49 *
50 * For the ModifierX interface, this adaptor maps the
51 * StepperSubcyclingAppAction::ACTION_LOCATION to the
52 * StepperSubcyclingModifierX::MODIFIERX_TYPE, and only pass the solution
53 * (\f$x\f$ and/or \f$\dot{x}\f$ and other parameters to the modify
54 * function.
55 */
56 void execute(
57 Teuchos::RCP<SolutionHistory<Scalar> > sh,
58 Teuchos::RCP<StepperSubcycling<Scalar> > stepper,
60 {
61 using Teuchos::RCP;
62
64 RCP<SolutionState<Scalar> > workingState = sh->getWorkingState();
65 const Scalar time = workingState->getTime();
66 const Scalar dt = workingState->getTimeStep();
67 RCP<Thyra::VectorBase<Scalar> > x;
68
69 switch(actLoc) {
71 {
72 modType = X_BEGIN_STEP;
73 x = workingState->getX();
74 break;
75 }
77 {
78 modType = XDOT_END_STEP;
79 if (workingState->getXDot() != Teuchos::null)
80 x = workingState->getXDot();
81 else
82 x = stepper->getStepperXDot();
83 break;
84 }
85 default:
86 TEUCHOS_TEST_FOR_EXCEPTION(true, std::logic_error,
87 "Error - unknown action location.\n");
88 }
89
90 this->modify(x, time, dt, modType);
91 }
92
93public:
94
100
102 virtual void modify(
103 Teuchos::RCP<Thyra::VectorBase<Scalar> > /* x */,
104 const Scalar /* time */, const Scalar /* dt */,
105 const MODIFIER_TYPE modType) = 0;
106
107};
108
109} // namespace Tempus
110
111#endif // Tempus_StepperSubcyclingModifierXBase_hpp
Application Action for StepperSubcycling.
ACTION_LOCATION
Indicates the location of application action (see algorithm).
MODIFIER_TYPE
Indicates the location of application action (see algorithm).
void execute(Teuchos::RCP< SolutionHistory< Scalar > > sh, Teuchos::RCP< StepperSubcycling< Scalar > > stepper, const typename StepperSubcyclingAppAction< Scalar >::ACTION_LOCATION actLoc)
Execute application action for Subcycling Stepper.
virtual void modify(Teuchos::RCP< Thyra::VectorBase< Scalar > >, const Scalar, const Scalar, const MODIFIER_TYPE modType)=0
Modify solution based on the MODIFIER_TYPE.