DirectInterpolation.hpp

#include <fish/fiber/field/DirectProductArray.hpp> 
#include <fish/fiber/vector/Interpolate.hpp>
#include <fish/fiber/vector/LinearIpol.hpp>
#include <fish/fiber/baseop/UniGridMapper.hpp>
#include <fish/fiber/baseop/LocalFromWorldPoint.hpp>
#include <fish/fiber/baseop/ExpandBBox.hpp>

#include <fish/fiber/field/Fragment.hpp>
#include <fish/fiber/grid/Chart.hpp>

#include <fish/fiber/baseop/BundleBaseSpaceFunctions.hpp>

#include <eagle/PhysicalSpace.hpp>
//#include <eagle/KDTree.hpp>
#include <eagle/FixedArray.hpp>

#include <memcore/RefPtr.hpp>
#include "FragmentSkeleton.hpp"
#include "PointSearch.hpp"
#include "gridopDllApi.h"

using namespace Fiber;
using namespace Eagle;
using namespace Eagle::PhysicalSpace;


class DirectInterpolate
{
 
  typedef MemArray<3, Eagle::tvector3> VectorArray_t; 
  typedef MemArray<3, point> CoordsArray_t;
  typedef std::list<MultiIndex<3> > IndexList_t; 
  //typedef pair< RefPtr<MemArray<1, point> >, RefPtr<MemArray<1, Eagle::tvector3> > > Pair_t;

  RefPtr<Field> Coords;
  CartesianFragments FragmentProperties;
  //RefPtr<Field> vecField;

  RefPtr<KDTree<3, int> > MyTree;
  //RefPtr<UniGridMapper> uniMap;

public:

  typedef MemArray<1, point> CoordArray_t;
  typedef MemArray<1, Eagle::tvector3> VecArray_t;
  typedef pair<RefPtr<CoordArray_t>,RefPtr<VecArray_t> > pairArray_t;

  DirectInterpolate(const RefPtr<Field>&coordinates, const RefPtr<Grid>&grid);
  ~DirectInterpolate();

  bool getPoint(const point&pos, const RefPtr<Field>&vectField, tvector&outputField);  
  int PointInCell(const point&p, const MultiArray<1, point>&coords);//const point&c1, const point&c2, const point&c3, const point&c4, const double z_up);
  unsigned interpVectorField(const point&p, const MultiArray<3, tvector3>&vecField, const MultiArray<3, point>&coords, tvector&outputField,const RefPtr<UniGridMapper>&uniMap);
  void InterpCell(const point&p, const MultiArray<1, tvector3>&vecField, const MultiArray<1, point>&coords, tvector&outVector);
  double InterpQuad(const point&p, const MultiArray<1, tvector3>&vecField, const MultiArray<1, point>&coords, double v[], int offset, double *zface);
  void RearrangePoints(RefPtr<MemArray<1, point> >&vertexCoord, RefPtr<MemArray<1, tvector3> >&vertexVector, pairArray_t&cell);
   
};