grid.h

Go to the documentation of this file.

#ifndef grid_h
#define grid_h

#include "floatmatrix.h"
#include "heap.h"

namespace oofem {
class Grid
{
public:
    Grid(int N, int M);
    ~Grid();

    int giveSizeHorizontal() { return n; }
    int giveSizeVertical() { return m; }
    int giveSize(int dir);

    void unFreeze();

    void setMethod(int o, double i, int c) { order = o; initDiag = i; centDiff = c; }

    //void setPrescribedField(double *PF);
    void setPrescribedFieldValue(int i, int j, double val) { F->at(i, j) = val; }
    FloatMatrix *givePrescribedField() { return F; }

    void setZeroValues(FloatMatrix *gridCoords);
    void setSolutionValueAt(int i, int j, double value);

    double giveSolutionValueAt(int i, int j);
    void printSolutionAsData();


private:
    int n, m;

    int order;
    double initDiag;
    int centDiff;

    bool solutionAvailable;

    FloatMatrix *T;

    FloatMatrix *F;

    bool *Frozen;

    Heap *narrowBand;

    void fastMarch(int &eFlag);

    double calcTime(int i, int j, double Fij, int ord, int &eFlag);

    int ij2ind(int i, int j, int m) { return ( ( i - 1 ) + ( j - 1 ) * m ); }
    int ind2i(int ind, int m) { return 1 + ind - ( ind / m ) * m; }
    int ind2j(int ind, int m) { return ( ( ind ) / m ) + 1; }
    bool isInDomain(int i, int j, int m, int n) { return ( ( i >= 1 ) && ( j >= 1 ) && ( i <= m ) && ( j <= n ) ); }
};
} // end namespace oofem
#endif // grid_h