ISS_VertexField.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_ISS_VertexField_h__)
#define __geom_ISS_VertexField_h__

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

#include "../../defs.h"

#include "../../../ads/array/Array.h"
// Interpolation functions.
#include "../../../numerical/interpolation/simplex.h"

BEGIN_NAMESPACE_GEOM


template<class ISS, 
         typename F = typename ISS::Number>
class ISS_VertexField
{
  //
  // 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:

  const IssType& _iss;
  FieldContainer _fields;
  mutable ads::FixedArray<ISS::M + 1,Vertex> _pos;
  mutable ads::FixedArray<ISS::M + 1,Field> _val;

  //
  // Not implemented.
  //

private:

  ISS_VertexField();

  ISS_VertexField&
  operator=(const ISS_VertexField& x);

public:

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


  template <bool A2>
  ISS_VertexField(const IssType& iss,
                  const ads::Array<1,Field,A2>& fields) :
    _iss(iss),
    _fields(fields),
    _pos(),
    _val() {
    assert(_iss.getVerticesSize() == _fields.size());
  }


  ISS_VertexField(const IssType& iss,
                  const SizeType numVertices, 
                  FieldVoidPointer fields) :
    _iss(iss),
    _fields(numVertices, fields),
    _pos(),
    _val() {
    assert(_iss.getVerticesSize() == numVertices);
  }

  ISS_VertexField(const ISS_VertexField& other) :
    _iss(other._iss),
    _fields(other._fields),
    _pos(),
    _val()
  {}

  ~ISS_VertexField()
  {}

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

  FieldParameterType
  interpolate(int n, const Vertex& x) const;

};

END_NAMESPACE_GEOM

#define __geom_ISS_VertexField_ipp__
#include "ISS_VertexField.ipp"
#undef __geom_ISS_VertexField_ipp__

#endif

---