InteractionDiagram.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/>.
//----------------------------------------------------------------------------
//InteractionDiagram.h

#ifndef INTERACTION_DIAGRAM_H
#define INTERACTION_DIAGRAM_H

#include "utility/geom/d2/Trihedron.h"
#include <set>
#include <deque>
#include "ClosedTriangleMesh.h"

class Triang3dMesh;

namespace XC {

class Vector;
class FiberSectionBase;
class InteractionDiagramData;

//
class InteractionDiagram: public ClosedTriangleMesh
  {
  protected:
    typedef std::set<const Trihedron *> set_ptr_trihedrons;

    
    set_ptr_trihedrons quadrant_trihedrons[8];

    void classify_trihedron(const Trihedron &tdro);
    void classify_trihedrons(void);
    void setPositionsMatrix(const Matrix &);
    GeomObj::list_Pos3d get_intersection(const Pos3d &p) const;
  public:
    InteractionDiagram(void);
    InteractionDiagram(const Pos3d &org,const Triang3dMesh &mll);
    InteractionDiagram(const InteractionDiagram &);
    InteractionDiagram &operator=(const InteractionDiagram &);
    virtual InteractionDiagram *getCopy(void) const;

    const Trihedron *findTrihedronPtr(const Pos3d &p) const;
    Pos3d getIntersection(const Pos3d &) const;
    double getCapacityFactor(const Pos3d &) const;
    Vector getCapacityFactor(const GeomObj::list_Pos3d &) const;

    void Print(std::ostream &os) const;
  };

InteractionDiagram calc_interaction_diagram(const FiberSectionBase &,const InteractionDiagramData &);

} // end of XC namespace

#endif