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