9#ifndef Tempus_Interpolator_hpp
10#define Tempus_Interpolator_hpp
13#include "Teuchos_VerboseObject.hpp"
14#include "Teuchos_Describable.hpp"
15#include "Teuchos_ParameterListAcceptorDefaultBase.hpp"
16#include "Teuchos_RCP.hpp"
17#include "Teuchos_ScalarTraits.hpp"
20#include "Tempus_config.hpp"
21#include "Tempus_SolutionState.hpp"
31 :
virtual public Teuchos::Describable,
32 virtual public Teuchos::ParameterListAcceptor,
33 virtual public Teuchos::VerboseObject<Interpolator<Scalar> >
56 const Teuchos::RCP<
const std::vector<Teuchos::RCP<SolutionState<Scalar> > > > & nodes) = 0;
61 SolutionState<Scalar>* state_out)
const = 0;
73 SolutionState<Scalar>* state_out)
75 interpolator.interpolate(t, state_out);
81 const Teuchos::RCP<
const std::vector<Teuchos::RCP<SolutionState<Scalar> > > >& nodes,
83 SolutionState<Scalar>* state_out)
85 interpolator.setNodes(nodes);
86 interpolator.interpolate(t, state_out);
90template <
typename Scalar>
94 typedef Teuchos::ScalarTraits<Scalar> ST;
95 typedef typename ST::magnitudeType mag_type;
97 const mag_type tol = ST::magnitude(scale)*ST::eps()*mag_type(1000.0);
98 const mag_type diff = ST::magnitude(a-b);
100 return diff <= tol || diff <= ST::sfmin();
Base strategy class for interpolation functionality.
virtual void setNodes(const Teuchos::RCP< const std::vector< Teuchos::RCP< SolutionState< Scalar > > > > &nodes)=0
Store pointer to interpolation nodes.
virtual void interpolate(const Scalar &t, SolutionState< Scalar > *state_out) const =0
Perform an interpolation.
virtual int order() const =0
Return the order of the interpolation.
void interpolate(const Interpolator< Scalar > &interpolator, const Scalar &t, SolutionState< Scalar > *state_out)
Nonmember functions.
bool floating_compare_equals(const Scalar &a, const Scalar &b, const Scalar &scale)
Helper function for comparing times.