VertexCFD_TempusObserver_ErrorNormOutput.hpp

#ifndef VERTEXCFD_TEMPUSOBSERVER_ERRORNORMOUTPUT_HPP
#define VERTEXCFD_TEMPUSOBSERVER_ERRORNORMOUTPUT_HPP
 
#include <observers/VertexCFD_Compute_ErrorNorms.hpp>
 
#include <Tempus_Integrator.hpp>
#include <Tempus_IntegratorObserver.hpp>
 
#include <Teuchos_FancyOStream.hpp>
#include <Teuchos_ParameterList.hpp>
#include <Teuchos_RCP.hpp>
 
#include <limits>
#include <string>
 
namespace VertexCFD
{
namespace TempusObserver
{
//---------------------------------------------------------------------------//
template<class Scalar>
class ErrorNormOutput : virtual public Tempus::IntegratorObserver<Scalar>
{
  public:
    ErrorNormOutput(
        const Teuchos::ParameterList& error_norm_list,
        Teuchos::RCP<ComputeErrorNorms::ErrorNorms<Scalar>> error_norm);
 
    void observeStartIntegrator(
        const Tempus::Integrator<Scalar>& integrator) override;
 
    void
    observeStartTimeStep(const Tempus::Integrator<Scalar>& integrator) override;
 
    void
    observeNextTimeStep(const Tempus::Integrator<Scalar>& integrator) override;
 
    void
    observeBeforeTakeStep(const Tempus::Integrator<Scalar>& integrator) override;
 
    void
    observeAfterTakeStep(const Tempus::Integrator<Scalar>& integrator) override;
 
    void observeAfterCheckTimeStep(
        const Tempus::Integrator<Scalar>& integrator) override;
 
    void
    observeEndTimeStep(const Tempus::Integrator<Scalar>& integrator) override;
 
    void
    observeEndIntegrator(const Tempus::Integrator<Scalar>& integrator) override;
 
  private:
    void
    print_error_norms(const std::string& heading,
                      const std::vector<typename ComputeErrorNorms::ErrorNorms<
                          Scalar>::DofErrorNorm>& error_norm,
                      const double scaling = 1.0);
 
    int _output_freq;
    bool _compute_time_error;
    Teuchos::FancyOStream _ostream;
    Teuchos::RCP<ComputeErrorNorms::ErrorNorms<Scalar>> _error_norm;
    static constexpr int prec = std::numeric_limits<double>::digits10 + 1;
    std::vector<typename ComputeErrorNorms::ErrorNorms<Scalar>::DofErrorNorm>
        _L1_time_error_norms;
    std::vector<typename ComputeErrorNorms::ErrorNorms<Scalar>::DofErrorNorm>
        _L2_time_error_norms;
};
 
//---------------------------------------------------------------------------//
 
} // end namespace TempusObserver
} // end namespace VertexCFD
 
#include "VertexCFD_TempusObserver_ErrorNormOutput_impl.hpp"
 
#endif // end VERTEXCFD_TEMPUSOBSERVER_ERRORNORMOUTPUT_HPP