Tempus_Stepper_ErrorNorm_decl.hpp

Go to the documentation of this file.

// @HEADER
// ****************************************************************************
//                Tempus: Copyright (2017) Sandia Corporation
//
// Distributed under BSD 3-clause license (See accompanying file Copyright.txt)
// ****************************************************************************
// @HEADER
 
#ifndef Tempus_Stepper_ErrorNorm_decl_hpp
#define Tempus_Stepper_ErrorNorm_decl_hpp
 
#include "Tempus_config.hpp"
 
#include "Teuchos_RCPDecl.hpp"
#include "Thyra_VectorBase.hpp"
#include "Thyra_VectorSpaceFactoryBase.hpp"
namespace Tempus {

template<class Scalar>
class Stepper_ErrorNorm
{
public:

  Stepper_ErrorNorm();

  Stepper_ErrorNorm(const Scalar relTol, const Scalar absTol);

  ~Stepper_ErrorNorm() {};

  Scalar computeWRMSNorm(
    const Teuchos::RCP<const Thyra::VectorBase<Scalar>> &x,
    const Teuchos::RCP<const Thyra::VectorBase<Scalar>> &xNext,
    const Teuchos::RCP<const Thyra::VectorBase<Scalar>> &err);

  Scalar errorNorm(const Teuchos::RCP<const Thyra::VectorBase<Scalar>> &x);
 
  void setRelativeTolerance(const Scalar relTol) { relTol_ = relTol; }
  void setAbsoluteTolerance(const Scalar absTol) { abssTol_ = absTol; }
 
 
protected:
 
  Scalar relTol_;
  Scalar abssTol_;
  Teuchos::RCP<Thyra::VectorBase<Scalar>> u_;
  Teuchos::RCP<Thyra::VectorBase<Scalar>> uNext_;
  Teuchos::RCP<Thyra::VectorBase<Scalar>> errorWeightVector_;
  Teuchos::RCP<Thyra::VectorBase<Scalar>> scratchVector_;
 
};
 
} // namespace Tempus
#endif //Tempus_Stepper_ErrorNorm_decl_hpp