Face.h

// -*-c++-*-
//----------------------------------------------------------------------------
//  XC program; finite element analysis code
//  for structural analysis and design.
//
//  Copyright (C)  Luis C. 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 "preprocessor/multi_block_topology/entities/1d/CmbEdge.h"

class Polygon3d;
class Ref2d3d;
namespace XC {
class Body;
class PolygonalFace;

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

    void insert_body(Body *b);
    Node *getNode(const size_t &i);
    const Node *getNode(const size_t &i) const;
  public:
    typedef std::deque<PolygonalFace *> dq_holes;
    typedef dq_holes::iterator hole_iterator;
    typedef dq_holes::const_iterator hole_const_iterator;
  protected:
    bool hole; 
    size_t ndivj; 
    dq_holes holes; 
    int get_index_opposite_side(const int &) const;
    virtual const Edge *get_opposite_side(const Edge *l) const;
    Edge *get_opposite_side(const Edge *l);
    size_t calc_ndiv_opposite_sides(const Edge *, const size_t &) const;
    void set_ndiv_opposite_sides(const size_t &, const size_t &);
    void set_ndiv_i(const size_t &);
    void set_ndiv_j(const size_t &);
    
    std::vector<XC::Node *> create_interior_nodes(const std::vector<Pos3d> &);
    int create_elements_from_quads(const std::deque<std::vector<int> > &);
    void create_line_nodes(void);
    int create_gmsh_loop(void) const;
    std::vector<int> create_gmsh_loops_for_holes(void) const;
    std::vector<int> create_gmsh_loops(void) const;
    int create_gmsh_surface(void) const;
  public:
    Face(void);
    Face(Preprocessor *m,const size_t &ndivI= 4, const size_t &ndivJ= 4);
    Face(const std::string &name,Preprocessor *m,const size_t &ndivI= 4, const size_t &ndivJ= 4);
    virtual bool operator==(const Face &) const;
    inline virtual unsigned short int GetDimension(void) const
      { return 2; }

    // Topology
    void update_topology(void);
    const Edge *getOppositeEdge(const Edge &) const;
    Edge *getOppositeEdge(const Edge &);
    std::deque<std::pair<const Edge *, const Edge *> > getOppositeEdges(void) const;

    // Surface orientation.
    Vector3d getKVector(void) const;
    Matrix getLocalAxes(void) const;
    Ref2d3d getRef(void) const;
    
    // Number of divisions
    void setNDiv(const size_t &);
    inline size_t NDivI(void) const
      { return ndiv; }
    virtual void setNDivI(const size_t &);
    inline size_t NDivJ(void) const
      { return ndivj; }
    virtual void setNDivJ(const size_t &);
    virtual void ConciliaNDivIJ(void);
    virtual void SetElemSizeI(const double &sz);
    virtual void SetElemSizeJ(const double &sz);
    virtual void SetElemSizeIJ(const double &,const double &);
    virtual void SetElemSize(const double &sz, bool mustBeEven= true);

    bool isHole(void) const;
    void setHole(const bool &);
    void addHole(PolygonalFace *);
    hole_iterator findHole(PolygonalFace *);
    hole_const_iterator findHole(PolygonalFace *) const;
    const PolygonalFace *findHolePtr(PolygonalFace *) const;
    PolygonalFace *findHolePtr(PolygonalFace *);
    const std::deque<PolygonalFace *> &getHoles(void) const;
    boost::python::list getPyHoles(void) const;
    std::deque<Side *> findSides(const Pos3d &);

    size_t getNumberOfVertices(void) const
      { return getNumberOfEdges(); }
    virtual const Pnt *getVertex(const size_t &i) const;
    virtual Pnt *getVertex(const size_t &i);
    Pnt *findVertex(const Pos3d &);
    Pos3d getCentroid(void) const;
    Plane getPlane(void) const;
    bool isFlat(const double &tol= 1e-6) const;
    bool isWarped(const double &tol= 1e-6) const;
    Polyline3d getContour(void) const;
    Polygon3d getPolygon(void) const;
    double getSquaredDistanceTo(const Pos3d &pt) const;    
    double getArea(void) const;
    const std::set<const Body *> &getConnectedBodies(void) const
      { return bodies_surf; }
    size_t CommonEdge(const Face &other) const;
    int SenseOfEdge(const Edge *l,const Face &other) const;
    bool isConnectedTo(const Body &b) const;
    std::deque<const Edge *> getNDivErrors(const size_t &) const;
    bool checkNDivs(const size_t &) const;
    std::deque<const Edge *> getNDivErrors(void) const;
    virtual bool checkNDivs(void) const;

    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 *> &);

    //Surface orientation.
    bool clockwise(const Pos3d &) const;
    bool counterclockwise(const Pos3d &) const;
    std::string orientation(const Pos3d &) const;
    virtual void reverse(void);
    void revolve(const int &);
    void setKOrientation(const Vector3d &);
    int setIOrientation(const Vector3d &);

    int getVtkCellType(void) const;
    
    boost::python::dict getPyDict(void) const;
    void setPyDict(const boost::python::dict &);
  };

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

} //end of XC namespace.

#endif