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