src/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;

    int getVtx(int x, int y);
    int getVtxX(int vtxID);
    int getVtxY(int vtxID);

    int getXYvtx(vec3 pos);

    float getPosX(int vtxID);
    float getPosY(int vtxID);

    bool isFlipped;

    //TODO for students: improve this
    int getNearestVtx(vec3 pos);

public:
    FlowGeometry();
    ~FlowGeometry();

    vec3* geometryData;


    bool getInterpolationAt(vec3 pos, int* vtxID, float* coef);

    bool readFromFile(char* header, FILE* fp, bool bigEndian);

    vec3 getPos(int vtxID);

    //remember that our grids are curvilinear and only 2D
    int getDimX();
    int getDimY();
    int getDimZ();

    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