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