pond/flow/lsuflow/RectilinearInterpolation.hpp

#ifndef _RectilinearInterpolation_HPP
#define _RectilinearInterpolation_HPP
 
#include <fish/fiber/field/Field.hpp>
#include <eagle/PhysicalSpace.hpp>
#include <memcore/RefPtr.hpp>
#include <fish/fiber/field/DirectProductArrayTypes.hpp>
 
using namespace Fiber;
using namespace Eagle::PhysicalSpace;
 
typedef MemArray<3, tvector>                    VectorArray_t;
typedef CartesianRectilinearArray3D             RectiLinearArray_t;
 
class RectilinearInterpolation
{
        RectiLinearArray_t &rectCoords;
        RefPtr<MemArray<3, double> >u_field,v_field,w_field;
        double timeStepForVectorField;
public:
        RectilinearInterpolation(const RefPtr<RectiLinearArray_t> &Points);
        bool getPoint(const point&pos, const RefPtr<Field> &vectField, tvector&outputField);
        bool find_rectilinear_cell(const point&pos,point &fp_index);
        static void rectilinear_interpolation(const point &fp_index,const RefPtr<Field> &vectField, tvector&outputField);
};
#endif