Tempus_UnitTest_Stepper_ErrorNorm.cpp

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

#include "Tempus_UnitTest_Utils.hpp"

#include "Thyra_DefaultSpmdVectorSpace.hpp"

#include "Tempus_NumericalUtils.hpp"
#include "Tempus_Stepper_ErrorNorm.hpp"

#include "../TestModels/DahlquistTestModel.hpp"


namespace Tempus_Unit_Test {

// ************************************************************
// ************************************************************
TEUCHOS_UNIT_TEST(Stepper_ErrorNorm, computeWRMSNorm)
{

  int N = 3;
  Teuchos::RCP<const Thyra::VectorSpaceBase<double> > xSpace =
    Thyra::defaultSpmdVectorSpace<double>(N);
  auto x      = Thyra::createMember(xSpace);
  auto xNext  = Thyra::createMember(xSpace);
  auto eVec   = Thyra::createMember(xSpace);

  auto tol = Tempus::numericalTol<double>();
  auto eNorm = Teuchos::rcp(new Tempus::Stepper_ErrorNorm<double>(tol, 0.0));

  Thyra::assign(x.ptr(),     0.0);
  Thyra::assign(xNext.ptr(), 0.0);
  Thyra::assign(eVec.ptr(),  0.0);
  double rmsNorm = eNorm->computeWRMSNorm(x, xNext, eVec);
  TEST_FLOATING_EQUALITY(rmsNorm, 0.0, 1.0e-12);

  Thyra::assign(x.ptr(),     1.0);
  Thyra::assign(xNext.ptr(), 1.0+tol);
  Thyra::assign(eVec.ptr(),  tol);
  rmsNorm = eNorm->computeWRMSNorm(x, xNext, eVec);
  TEST_FLOATING_EQUALITY(rmsNorm, 1.0/std::sqrt(N), 1.0e-12);
}


// ************************************************************
// ************************************************************
TEUCHOS_UNIT_TEST(Stepper_ErrorNorm, errorNorm)
{

  Teuchos::RCP<const Thyra::VectorSpaceBase<double> > xSpace =
    Thyra::defaultSpmdVectorSpace<double>(3);
  auto x      = Thyra::createMember(xSpace);

  auto tol = Tempus::numericalTol<double>();
  auto eNorm = Teuchos::rcp(new Tempus::Stepper_ErrorNorm<double>(tol, tol));

  Thyra::assign(x.ptr(),     0.0);
  double norm = eNorm->errorNorm(x);
  TEST_FLOATING_EQUALITY(norm, tol, 1.0e-12);

  Thyra::assign(x.ptr(),     1.0);
  norm = eNorm->errorNorm(x);
  TEST_FLOATING_EQUALITY(norm, 1.0/(2.0*tol), 1.0e-12);
}


} // namespace Tempus_Unit_Test