Tempus_InterpolatorFactory.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_InterpolatorFactory_hpp
#define Tempus_InterpolatorFactory_hpp

#include "Teuchos_ParameterList.hpp"
#include "Tempus_config.hpp"
#include "Tempus_InterpolatorLagrange.hpp"

namespace Tempus {

template <class Scalar>
class InterpolatorFactory {
 public:
  static Teuchos::RCP<Interpolator<Scalar> > createInterpolator(
      std::string interpolatorType = "")
  {
    if (interpolatorType == "") interpolatorType = "Lagrange";
    return createInterpolator(interpolatorType, Teuchos::null);
  }

  static Teuchos::RCP<Interpolator<Scalar> > createInterpolator(
      const Teuchos::RCP<Teuchos::ParameterList>& interpolatorPL)
  {
    std::string interpolatorType =
        interpolatorPL->get<std::string>("Interpolator Type", "Lagrange");
    return createInterpolator(interpolatorType, interpolatorPL);
  }

 private:
  static Teuchos::RCP<Interpolator<Scalar> > createInterpolator(
      const std::string& interpolatorType,
      const Teuchos::RCP<Teuchos::ParameterList>& interpolatorPL)
  {
    using Teuchos::rcp;

    Teuchos::RCP<Interpolator<Scalar> > interpolator;
    if (interpolatorType == "Lagrange")
      interpolator = rcp(new InterpolatorLagrange<Scalar>);
    else {
      TEUCHOS_TEST_FOR_EXCEPTION(
          true, std::logic_error,
          "Unknown 'Interpolator Type' = " << interpolatorType);
    }
    interpolator->setParameterList(interpolatorPL);

    return interpolator;
  }
};

}  // namespace Tempus
#endif  // Tempus_InterpolatorFactory_hpp