An Object Finite Element Library

Grid.h

Go to the documentation of this file.

/*==============================================================================

                                    O  F  E  L  I

                           Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2008 Rachid Touzani

   This program is free software; you can redistribute it and/or modify it under
   the terms of the GNU General Public License as published by the Free
   Software Foundation; Version 2 of the License.

   This program is distributed in the hope that it will be useful, but WITHOUT
   ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
   FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
   details.

   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================

                 Definition of class Grid for structured grids

  ==============================================================================*/

#ifndef __GRID_H
#define __GRID_H

#include <stdlib.h>
#include <math.h>
#include <string>

#include <iostream>
using std::cout;
using std::ostream;
using std::endl;

#include <iomanip>
using std::setw;

using std::string;

#include "OFELI_Config.h"
#include "Point.h"
#include "Point2D.h"
#include "Node.h"
#include "Element.h"
#include "Side.h"
#include "util.h"

namespace OFELI {

#ifndef DOXYGEN_SHOULD_SKIP_THIS
struct ErrorInGrid {
   void Message(const char *file, size_t line, int code, double p1, double p2);
};
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

class Grid
{
  public:
    
    Grid();

    Grid(double xm, double xM, size_t npx);

    Grid(double xm, double xM, double ym, double yM, size_t npx, size_t npy);

    Grid(double xm, double xM, double ym, double yM, double zm, double zM, 
         size_t npx, size_t npy, size_t npz);

    void setXMin(const Point<double> &x) { _xmin = x; }

    void setXMax(const Point<double> &x) { _xmax = x; }

    void setDomain(double xmin, double xmax);

    void setDomain(double xmin, double xmax, double ymin, double ymax);

    void setDomain(double xmin, double xmax, double ymin, double ymax, double zmin, double zmax);

    const Point<double> &getXMin() const { return _xmin; }

    const Point<double> &getXMax() const { return _xmax; }

    void setNx(size_t nx) { _nx = nx; }

    void setNy(size_t ny) { _ny = ny; }

    void setNz(size_t nz) { _nz = nz; }

    size_t getNx() const { return _nx; }

    size_t getNy() const { return _ny; }

    size_t getNz() const { return _nz; }

    double getHx() const { return _hx; }

    double getHy() const { return _hy; }

    double getHz() const { return _hz; }
    
    Point<double> getCoord(size_t i) const;
    
    Point<double> getCoord(size_t i, size_t j) const;

    Point<double> getCoord(size_t i, size_t j, size_t k) const;

    double getX(size_t i) const;

    double getY(size_t j) const;

    double getZ(size_t k) const;

    Point2D<double> getXY(size_t i, size_t j) const;

    Point<double> getXYZ(size_t i, size_t j, size_t k) const;

    void setCode(int side, int code);

    int getCode(int side) const;

    int getCode(size_t i, size_t j) const;

    int getCode(size_t i, size_t j, size_t k) const;

  private:
    int _cmx, _cMx, _cmy, _cMy, _cmz, _cMz;
    size_t _nx, _ny, _nz;
    double _hx, _hy, _hz;
    Point<double> _xmin, _xmax;
    ErrorInGrid _e;
};

    ostream& operator<<(ostream& s, const Grid &g);

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani