Tempus  Version of the Day
Time Integration
Tempus_InterpolatorLagrange_decl.hpp
Go to the documentation of this file.
1 // @HEADER
2 // ****************************************************************************
3 // Tempus: Copyright (2017) Sandia Corporation
4 //
5 // Distributed under BSD 3-clause license (See accompanying file Copyright.txt)
6 // ****************************************************************************
7 // @HEADER
8 
9 #ifndef Tempus_InterpolatorLagrange_decl_hpp
10 #define Tempus_InterpolatorLagrange_decl_hpp
11 
12 #include "Tempus_config.hpp"
13 #include "Tempus_Interpolator.hpp"
14 
15 namespace Tempus {
16 
20 template<class Scalar>
21 class InterpolatorLagrange : virtual public Interpolator<Scalar>
22 {
23 public:
24 
27 
30 
32 
33 
35  void setNodes(
36  const Teuchos::RCP<const std::vector<Teuchos::RCP<SolutionState<Scalar> > > > & nodes) { nodes_ = nodes; }
37 
39  void interpolate(const Scalar& t, SolutionState<Scalar>* state_out) const;
40 
42  int order() const { return order_; }
43 
45 
47 
48  std::string description() const { return "Tempus::InterpolatorLagrange"; }
49  void describe(Teuchos::FancyOStream &out,
50  const Teuchos::EVerbosityLevel /* verbLevel */) const
51  {
52  out.setOutputToRootOnly(0);
53  out << description() << "::describe" << std::endl;
54  }
56 
58 
59  void setParameterList(const Teuchos::RCP<Teuchos::ParameterList>& paramList);
60  Teuchos::RCP<Teuchos::ParameterList> getNonconstParameterList();
61  Teuchos::RCP<Teuchos::ParameterList> unsetParameterList();
62  Teuchos::RCP<const Teuchos::ParameterList> getValidParameters() const;
64 
65 private:
66 
67  void lagrange(const int p, const Scalar& t,
68  SolutionState<Scalar>* state_out) const;
69 
70  Teuchos::RCP<const std::vector<Teuchos::RCP<SolutionState<Scalar> > > > nodes_;
71  Teuchos::RCP<Teuchos::ParameterList> pl_;
72  int order_;
73 
74 };
75 
76 // Nonmember constructor
77 template<class Scalar>
78 Teuchos::RCP<InterpolatorLagrange<Scalar> > lagrangeInterpolator()
79 {
80  return Teuchos::rcp(new InterpolatorLagrange<Scalar>);
81 }
82 
83 } // namespace Tempus
84 
85 #endif // Tempus_InterpolatorLagrange_decl_hpp
int order() const
Return the order of the interpolation.
Teuchos::RCP< const std::vector< Teuchos::RCP< SolutionState< Scalar > > > > nodes_
void interpolate(const Scalar &t, SolutionState< Scalar > *state_out) const
Perform an interpolation.
Teuchos::RCP< Teuchos::ParameterList > pl_
Teuchos::RCP< Teuchos::ParameterList > getNonconstParameterList()
void lagrange(const int p, const Scalar &t, SolutionState< Scalar > *state_out) const
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel) const
Base strategy class for interpolation functionality.
Teuchos::RCP< Teuchos::ParameterList > unsetParameterList()
void setNodes(const Teuchos::RCP< const std::vector< Teuchos::RCP< SolutionState< Scalar > > > > &nodes)
Store pointer to interpolation nodes.
Teuchos::RCP< InterpolatorLagrange< Scalar > > lagrangeInterpolator()
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
Concrete implemenation of Interpolator that does simple lagrange interpolation.
Solution state for integrators and steppers. SolutionState contains the metadata for solutions and th...
void setParameterList(const Teuchos::RCP< Teuchos::ParameterList > &paramList)