ISS_Interpolate.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_ISS_Interpolate_h__)
#define __geom_ISS_Interpolate_h__

#if defined(DEBUG_geom) && !defined(DEBUG_ISS_Interpolate)
#define DEBUG_ISS_Interpolate
#endif

#include "ISS_SimplexQuery.h"
#include "ISS_VertexField.h"

BEGIN_NAMESPACE_GEOM


template<class ISS, 
         typename F = typename ISS::Number>
class ISS_Interpolate {
  //
  // Private types.
  //

private:

  typedef ISS IssType;
  typedef typename IssType::SizeType SizeType;

  //
  // Public types.
  //

public:

  typedef typename IssType::Number Number;
  typedef typename IssType::Vertex Vertex;

  //
  // Field types.
  //

  typedef F Field;
  typedef ads::Array<1,Field,false> FieldContainer;

  typedef typename FieldContainer::parameter_type 
  FieldParameterType;

  typedef typename FieldContainer::void_pointer FieldVoidPointer;

  //
  // Data.
  //

private:

  ISS_SimplexQuery<ISS> _simplexQuery;
  ISS_VertexField<ISS,Field> _vertexField;

  //
  // Not implemented.
  //

private:

  ISS_Interpolate();

  ISS_Interpolate&
  operator=(const ISS_Interpolate&);

public:

  //--------------------------------------------------------------------------


  template <bool A2>
  ISS_Interpolate(const IssType& iss,
                  const ads::Array<1,Field,A2>& fields) :
    _simplexQuery(iss),
    _vertexField(iss, fields)
  {}


  ISS_Interpolate(const IssType& iss,
                  const SizeType num_vertices, 
                  FieldVoidPointer fields) :
    _simplexQuery(iss),
    _vertexField(iss, num_vertices, fields)
  {}

  ISS_Interpolate(const ISS_Interpolate& other) :
    _simplexQuery(other._simplexQuery),
    _vertexField(other._vertexField)
  {}

  ~ISS_Interpolate()
  {}

  //--------------------------------------------------------------------------

  FieldParameterType
  operator()(const Vertex& x) const {
    return _vertexField.interpolate
      (_simplexQuery.computeMinimumDistanceIndex(x), x);
  }

};

END_NAMESPACE_GEOM

#endif

---