Face.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/>.
//----------------------------------------------------------------------------
//Face.h
//Face entities.

#ifndef FACE_H
#define FACE_H

#include "CmbEdge.h"

namespace XC {
class Body;

class Face: public CmbEdge
  {
    friend class Edge;
    friend class Body;
    std::set<const Body *> cuerpos_sup; 

    void insert_body(Body *b);
    Node *getNode(const size_t &i);
    const Node *getNode(const size_t &i) const;
  protected:
    size_t ndivj; //number of divisions in the j axis.
    virtual const Edge *get_lado_homologo(const Edge *l) const= 0;

  public:
    Face(void);
    Face(Preprocessor *m,const size_t &ndivI= 4, const size_t &ndivJ= 4);
    Face(const std::string &nombre,Preprocessor *m,const size_t &ndivI= 4, const size_t &ndivJ= 4);
    inline virtual unsigned short int GetDimension(void) const
      { return 2; }
    void update_topology(void);
    inline size_t NDivI(void) const
      { return ndiv; }
    virtual void SetNDivI(const size_t &ndi);
    inline size_t NDivJ(void) const
      { return ndivj; }
    virtual void SetNDivJ(const size_t &ndj);
    virtual void ConciliaNDivIJ(void)= 0;
    size_t getNumberOfVertices(void) const
      { return getNumberOfEdges(); }
    virtual const Pnt *getVertex(const size_t &i) const;
    Polyline3d getContour(void) const;
    const std::set<const Body *> &getConnectedBodies(void) const
      { return cuerpos_sup; }
    size_t CommonEdge(const Face &otra) const;
    int SenseOfEdge(const Edge *l,const Face &otra) const;
    bool isConnectedTo(const Body &b) const;
    virtual bool checkNDivs(void) const= 0;

    virtual Node *getNode(const size_t &i1,const size_t &j,const size_t &k);
    virtual const Node *getNode(const size_t &i,const size_t &j,const size_t &k) const;
    virtual Node *getNode(const size_t &i,const size_t &j);
    virtual const Node *getNode(const size_t &i,const size_t &j) const;

    std::set<SetBase *> get_sets(void) const;
    void add_to_sets(std::set<SetBase *> &);

    int getVtkCellType(void) const;
  };

std::set<const Face *> getConnectedSurfaces(const Edge &l);

} //end of XC namespace.

#endif