Rythmos_InterpolationBuffer_decl.hpp

//@HEADER
// ***********************************************************************
//
//                           Rythmos Package
//                 Copyright (2006) Sandia Corporation
//
// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
// license for use of this work by or on behalf of the U.S. Government.
//
// This library is free software; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 2.1 of the
// License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301
// USA
// Questions? Contact Todd S. Coffey (tscoffe@sandia.gov)
//
// ***********************************************************************
//@HEADER
 
#ifndef Rythmos_INTERPOLATION_BUFFER_DECL_H
#define Rythmos_INTERPOLATION_BUFFER_DECL_H
 
#include "Rythmos_InterpolationBufferBase.hpp"
#include "Rythmos_Types.hpp"
#include "Rythmos_DataStore.hpp"
#include "Rythmos_InterpolatorAcceptingObjectBase.hpp"
 
 
 
namespace Rythmos {
 
enum IBPolicy {
  BUFFER_POLICY_INVALID = 0,
  BUFFER_POLICY_STATIC = 1,
  BUFFER_POLICY_KEEP_NEWEST = 2
};
 
 
template<class Scalar>
class InterpolationBuffer :
  virtual public InterpolationBufferBase<Scalar>,
  virtual public InterpolatorAcceptingObjectBase<Scalar>
{
public:
 
  typedef typename Teuchos::ScalarTraits<Scalar>::magnitudeType ScalarMag;
 
 
  RCP<const Thyra::VectorSpaceBase<Scalar> > get_x_space() const;
 
  InterpolationBuffer();
 
  void initialize( const RCP<InterpolatorBase<Scalar> >& interpolator, int storage );
 
 
  void setInterpolator(const RCP<InterpolatorBase<Scalar> >& interpolator);
 
  RCP<InterpolatorBase<Scalar> >
    getNonconstInterpolator();
 
  RCP<const InterpolatorBase<Scalar> >
    getInterpolator() const;
 
  RCP<InterpolatorBase<Scalar> > unSetInterpolator();
 
 
  void setStorage( int storage );
 
  int getStorage() const;
 
  IBPolicy getIBPolicy();
 
  ~InterpolationBuffer() {};
 
  void addPoints(
    const Array<Scalar>& time_vec
    ,const Array<RCP<const Thyra::VectorBase<Scalar> > >& x_vec
    ,const Array<RCP<const Thyra::VectorBase<Scalar> > >& xdot_vec);
 
  void getPoints(
    const Array<Scalar>& time_vec
    ,Array<RCP<const Thyra::VectorBase<Scalar> > >* x_vec
    ,Array<RCP<const Thyra::VectorBase<Scalar> > >* xdot_vec
    ,Array<ScalarMag>* accuracy_vec
    ) const;
 
  TimeRange<Scalar> getTimeRange() const;
 
  void getNodes(Array<Scalar>* time_vec) const;
 
  int getOrder() const;
 
  void removeNodes(Array<Scalar>& time_vec);
 
 
  std::string description() const;
 
  void describe(
    Teuchos::FancyOStream       &out
    ,const Teuchos::EVerbosityLevel      verbLevel
    ) const;
 
 
  void setParameterList(RCP<Teuchos::ParameterList> const& paramList);
 
  RCP<Teuchos::ParameterList> getNonconstParameterList();
 
  RCP<Teuchos::ParameterList> unsetParameterList();
 
  RCP<const Teuchos::ParameterList> getValidParameters() const;
 
private:
 
  RCP<InterpolatorBase<Scalar> > interpolator_;
  int storage_limit_;
  RCP<typename DataStore<Scalar>::DataStoreVector_t> data_vec_;
 
  RCP<Teuchos::ParameterList> paramList_;
 
  IBPolicy policy_;
 
 
  // Private member functions:
  void defaultInitializeAll_();
 
};
 
 
template<class Scalar>
RCP<InterpolationBuffer<Scalar> > interpolationBuffer(
  const RCP<InterpolatorBase<Scalar> >& interpolator = Teuchos::null,
  int storage = 0
  )
{
  RCP<InterpolationBuffer<Scalar> > ib = rcp(new InterpolationBuffer<Scalar>());
  ib->initialize(interpolator, storage);
  return ib;
}
 
 
} // namespace Rythmos
 
 
#endif // Rythmos_INTERPOLATION_BUFFER_DECL_H