FlowGeometry.h

#ifndef FLOWGEOMETRY_H
#define FLOWGEOMETRY_H

#include <stdio.h>
#include <iostream>
#include "vec3.h"

class FlowGeometry{
                friend class FlowData;
        private:
                int dim[2]; 
                vec3 boundaryMin;
                vec3 boundaryMax;
                vec3 boundarySize;
            
                vec3* geometryData;

                int getVtx(int x, int y);
                int getVtxX(int vtxID);
                int getVtxY(int vtxID);
            
                int getXYvtx(vec3 pos);


            
                bool isFlipped;

                int getNearestVtx(vec3 pos);
                
        public:
                FlowGeometry();
                ~FlowGeometry();
            

                bool getInterpolationAt(vec3 pos, int* vtxID, float* coef);
                
                bool readFromFile(char* header, FILE* fp, bool bigEndian);
            
                //remember that our grids are curvilinear and only 2D
                int getDimX();
                int getDimY();
                int getDimZ();
            
                float getMinX();
                float getMaxX();
                float getMinY();
                float getMaxY();

                vec3 getPos(int vtxID);
                float getPosX(int vtxID);
                float getPosY(int vtxID);
            
                int getRightNeigh(int vtxID);
                int getTopNeigh(int vtxID);
                int getLeftNeigh(int vtxID);
                int getBottomNeigh(int vtxID);
            
                vec3 normalizeCoords(vec3 pos);
                vec3 unNormalizeCoords(vec3 pos);
};

#endif