QuadCellRegion.h

//----------------------------------------------------------------------------
//  XC program; finite element analysis code
//  for structural analysis and design.
//
//  Copyright (C)  Luis Claudio Pérez Tato
//
//  This program derives from OpenSees <http://opensees.berkeley.edu>
//  developed by the  «Pacific earthquake engineering research center».
//
//  Except for the restrictions that may arise from the copyright
//  of the original program (see copyright_opensees.txt)
//  XC 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, either version 3 of the License, or 
//  (at your option) any later version.
//
//  This software 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, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------------
//QuadCellRegion.h

#ifndef QUADCELLREGION_H 
#define QUADCELLREGION_H


#include <iostream>
#include "material/section/repres/geom_section/region/SectRegion.h" 

class Grid2d;

namespace XC {

//
//
//         nDivIJ= 4
// L +---+---+---+---+ K
//   |   |   |   |   |
//   +---+---+---+---+ nDivKL= 2
//   |   |   |   |   |
// I +---+---+---+---+ J
class QuadCellRegion: public SectRegion
  {
    mutable Grid2d *grid; 
  protected:
    int nDivIJ, nDivJK; 
    void free_mem(void) const;
    const Grid2d &alloc(const Grid2d &) const;

  public:
    QuadCellRegion(Material *);
    QuadCellRegion(Material *, int numSubdivIJ, int numSubdivJK);
    QuadCellRegion(const QuadCellRegion &);
    QuadCellRegion &operator=(const QuadCellRegion &);
    virtual ~QuadCellRegion(void);

    void setDiscretization(int numSubdivIJ, int numSubdivJK);

    inline void setNDivIJ(const int &n)
      { nDivIJ= n; }
    inline int getNDivIJ(void) const
      { return nDivIJ; }
    inline void setNDivJK(const int &n)
      { nDivJK= n; }
    inline int getNDivJK(void) const
      { return nDivJK; }
    //inquiring functions
    inline const int &nDivCirc(void) const
      { return nDivIJ; }
    inline int &nDivCirc(void)
      { return nDivIJ; }
    inline const int &nDivRad(void) const
      { return nDivJK; }
    inline int &nDivRad(void)
      { return nDivJK; }
    int getNumCells(void) const;
    void getDiscretization(int &numSubdivIJ, int &numSubdivJK) const;
    virtual const Grid2d &getMesh(void) const=0;
    Vector getCenterOfMassGrid(void) const;
    Vector getVertCoords(const size_t &,const size_t &) const;
    Matrix getCellVertCoords(const size_t &,const size_t &) const;
    

  };
} // end of XC namespace


#endif