9 #ifndef Tempus_IntegratorObserverBasic_impl_hpp
10 #define Tempus_IntegratorObserverBasic_impl_hpp
12 #include "Tempus_Stepper.hpp"
16 template <
class Scalar>
21 template <
class Scalar>
26 template <
class Scalar>
30 std::time_t begin = std::time(
nullptr);
34 *out <<
"\nTempus - IntegratorBasic\n"
35 << std::asctime(std::localtime(&begin)) <<
"\n"
36 <<
" Stepper = " << integrator.
getStepper()->description() <<
"\n"
37 <<
" Simulation Time Range ["
40 <<
" Simulation Index Range ["
43 <<
"===================================================================="
45 <<
" Step Time dt Abs Error Rel Error Order nFail "
50 template <
class Scalar>
56 template <
class Scalar>
62 template <
class Scalar>
68 template <
class Scalar>
74 template <
class Scalar>
80 template <
class Scalar>
87 if ((cs->getOutputScreen() ==
true) || (cs->getOutput() ==
true) ||
96 *out << std::scientific << std::setw(6) << std::setprecision(3)
97 << cs->getIndex() << std::setw(11) << std::setprecision(3)
98 << cs->getTime() << std::setw(11) << std::setprecision(3)
99 << cs->getTimeStep() << std::setw(11) << std::setprecision(3)
100 << cs->getErrorAbs() << std::setw(11) << std::setprecision(3)
101 << cs->getErrorRel() << std::fixed << std::setw(7)
102 << std::setprecision(1) << cs->getOrder() << std::scientific
103 << std::setw(7) << std::setprecision(3) << cs->getNFailures()
104 << std::setw(11) << std::setprecision(3) << steppertime << std::endl;
108 template <
class Scalar>
112 std::string exitStatus;
117 exitStatus =
"Time integration FAILURE!";
120 exitStatus =
"Time integration complete.";
122 std::time_t end = std::time(
nullptr);
127 *out <<
"===================================================================="
129 <<
" Total runtime = " << runtime <<
" sec = " << runtime / 60.0
131 << std::asctime(std::localtime(&end)) <<
"\nNumber of Accepted Steps = "
133 <<
"\nNumber of Failures = "
137 ->getNRunningFailures()
139 << exitStatus <<
"\n"
144 #endif // Tempus_IntegratorObserverBasic_impl_hpp
virtual Teuchos::RCP< Teuchos::Time > getStepperTimer() const =0
virtual Teuchos::RCP< Teuchos::Time > getIntegratorTimer() const =0
Returns the IntegratorTimer_ for this Integrator.
virtual void observeEndIntegrator(const Integrator< Scalar > &integrator) override
Observe the end of the time integrator.
virtual void observeNextTimeStep(const Integrator< Scalar > &integrator) override
Observe after the next time step size is selected.
virtual Teuchos::RCP< const TimeStepControl< Scalar > > getTimeStepControl() const =0
Returns the TimeStepControl for this Integrator.
IntegratorObserverBasic()
Constructor.
virtual void observeStartIntegrator(const Integrator< Scalar > &integrator) override
Observe the beginning of the time integrator.
virtual Tempus::Status getStatus() const =0
Get the Status.
virtual void observeAfterCheckTimeStep(const Integrator< Scalar > &integrator) override
basic_FancyOStream & setOutputToRootOnly(const int rootRank)
Thyra Base interface for time integrators. Time integrators are designed to advance the solution from...
virtual void observeStartTimeStep(const Integrator< Scalar > &integrator) override
Observe the beginning of the time step loop.
virtual void observeAfterTakeStep(const Integrator< Scalar > &integrator) override
Observe after Stepper takes step.
virtual void observeEndTimeStep(const Integrator< Scalar > &integrator) override
Observe the end of the time step loop.
double totalElapsedTime(bool readCurrentTime=false) const
virtual void observeBeforeTakeStep(const Integrator< Scalar > &integrator) override
Observe before Stepper takes step.
virtual RCP< FancyOStream > getOStream() const
virtual Teuchos::RCP< Stepper< Scalar > > getStepper() const =0
Get the stepper.
virtual ~IntegratorObserverBasic()
Destructor.
virtual Teuchos::RCP< const SolutionHistory< Scalar > > getSolutionHistory() const =0
Returns the SolutionHistory for this Integrator.