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);

                vec3 getPos(int vtxID);
                float getPosX(int vtxID);
                float getPosY(int vtxID);    
            
                bool isFlipped;

                //TODO for students: improve this
                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
                inline int getDimX() {return (isFlipped) ? dim[1] : dim[0];}
                inline int getDimY() {return (isFlipped) ? dim[0] : dim[1];}
                inline int getDimZ() {return dim[2];}
            
                float getMinX();
                float getMaxX();
                float getMinY();
                float getMaxY();
            
                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

---