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.setOutputToRootOnly(0);
    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