Tempus_InterpolatorLagrange_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_InterpolatorLagrange_decl_hpp
#define Tempus_InterpolatorLagrange_decl_hpp

#include "Tempus_config.hpp"
#include "Tempus_Interpolator.hpp"

namespace Tempus {

template<class Scalar>
class InterpolatorLagrange : virtual public Interpolator<Scalar>
{
public:

  InterpolatorLagrange() : order_(0) {}

  ~InterpolatorLagrange() {}



  void setNodes(
    const Teuchos::RCP<const std::vector<Teuchos::RCP<SolutionState<Scalar> > > > & nodes) { nodes_ = nodes; }

  void interpolate(const Scalar& t, SolutionState<Scalar>* state_out) const;

  int order() const { return order_; }



  std::string description() const { return "Tempus::InterpolatorLagrange"; }
  void describe(Teuchos::FancyOStream &out,
                const Teuchos::EVerbosityLevel /* verbLevel */) const
  { out << description() << "::describe" << std::endl; }


  void setParameterList(const Teuchos::RCP<Teuchos::ParameterList>& paramList);
  Teuchos::RCP<Teuchos::ParameterList> getNonconstParameterList();
  Teuchos::RCP<Teuchos::ParameterList> unsetParameterList();
  Teuchos::RCP<const Teuchos::ParameterList> getValidParameters() const;

private:

  void lagrange(const int p, const Scalar& t,
                SolutionState<Scalar>* state_out) const;

  Teuchos::RCP<const std::vector<Teuchos::RCP<SolutionState<Scalar> > > > nodes_;
  Teuchos::RCP<Teuchos::ParameterList> pl_;
  int order_;

};

// Nonmember constructor
template<class Scalar>
Teuchos::RCP<InterpolatorLagrange<Scalar> > lagrangeInterpolator()
{
  return Teuchos::rcp(new InterpolatorLagrange<Scalar>);
}

} // namespace Tempus

#endif // Tempus_InterpolatorLagrange_decl_hpp