gmDiscontinuitySet.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2023 by Carlos Augusto Teixeira Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/

#ifndef _GEMA_DISCONTINUITY_SET_H_
#define _GEMA_DISCONTINUITY_SET_H_

#include "gmLog.h"

#include <QVector>
#include <unit.h>

class  GmDiscontinuity;
class  GmDiscontinuityAccessor;
struct GmSegmentCellIntersection;
class  GmCompact2DSegmentCellIntersection;
class  GmSegmentCell2DIntersectionCoordAccessor;
class  GmCell;
class  GmCellMesh;
class  GmValueSet;
class  GmPropertySet;
class  GmValueInfo;
class  GmSimulationData;
class  GmSpatialIndex;
class  LuaTable;


class GMC_API_EXPORT GmDiscontinuitySet
{
public:
  enum DiscontinuityType
  {
    GM_DISC_POLYLINE,    
    GM_DISC_TRISURFACE,  
    GM_DISC_EDGELIST,    
    GM_DISC_FACELIST,    
    // -- If changes are made to this enum, please remember to adjust 
    //    discontinuityTypeToStr()
  };

  typedef QVarLengthArray<GmCompact2DSegmentCellIntersection*, 1> SegmentCell2DIntersectionList;

  ~GmDiscontinuitySet();

  QString id() const { return _id; }

  QString description() const { return _description; }

  DiscontinuityType type() const { return _type; }

  GmCellMesh* mesh() const { return _mesh; }

  GmSpatialIndex* spatialIndex() const { return _index; }

  bool autoAddElements() const { return _addElements; }

  bool extraDofsInAddedElements() const { return _addedInterfElemHaveExtraDofs; }

  double snapTolerance() const { return _snapTol; }

  bool noInternalNodes() const { return _noInternalNodes; }             

  bool internalIntersections() const { return _internalIntersections; }

  bool discardOriginalGeometry() const { return _discardOriginalGeometry; }

  int numDiscontinuities() const { return _data.size(); }

  const QVector<GmDiscontinuity*>& discontinuities() const { return _data; }

  int discontinuityIndex(QString id) const;

  GmDiscontinuity* discontinuity(int index) { return _data[index]; }

  const GmDiscontinuity* discontinuity(int index) const { return _data[index]; }

  QStringList              attributeIds() const;
  GmValueInfo*             attributeInfo     (QString id) const;
  GmDiscontinuityAccessor* attributeAccessor (QString id, Unit desiredUnit, const GmLogCategory& logger) const;
  GmDiscontinuityAccessor* attributeAccessor (QString id, int snum, bool locked, Unit desiredUnit, const GmLogCategory& logger) const;
  int                      numAttributeStates(QString id) const;

  int numPropertySets() const { return _propertySets.size(); }

  const QVector<GmPropertySet*>& propertySets() const { return _propertySets; }

  int propertySetIndex(QString id) const;

  GmValueInfo*             propertyInfo    (QString id) const;
  GmDiscontinuityAccessor* propertyAccessor(QString id, Unit desiredUnit, const GmLogCategory& logger) const;

  bool buildIntersectionMaps(const GmLogCategory& logger);

  const SegmentCell2DIntersectionList&      cell2DPolylineIntersections(const GmCell* c) const;
  GmSegmentCell2DIntersectionCoordAccessor* cell2DPolylineCoordAccessor(Unit desiredUnit, const GmLogCategory& logger) const;

  static GmDiscontinuitySet* loadFromTable(LuaTable& tab, GmSimulationData* data, const GmLogCategory& logger);

  void printParameters(const GmLogCategory& logger, GmLogLevel level) const;
  void print          (const GmLogCategory& logger, GmLogLevel level) const;

  size_t printMemoryStatistics(const GmLogCategory& logger, GmLogLevel level) const;
  size_t usedMemory() const;

  static QString discontinuityTypeToStr(DiscontinuityType type);

private:
  GmDiscontinuitySet(QString id, QString description, DiscontinuityType type, GmCellMesh* mesh,
                     GmSpatialIndex* index, bool addElements, bool extraDofs, double snapTol,
                     bool noInternalNodes, bool internalIntersections, bool discardOriginalGeometry);

  bool add2DSegmentCellIntersectionsToMesh(const GmLogCategory& logger);
  bool break2DPolylinesAtIntersections(const GmLogCategory& logger);

  bool loadGeometry(GmDiscontinuity* d, LuaTable& tab, DiscontinuityType type, const GmCellMesh* mesh, const GmLogCategory& logger);
  template <class T, class V> T* loadIntoVector(LuaTable& tab, V validator, int n, QString discontinuityId, const char* msg, const GmLogCategory& logger);

  size_t usedElementMapMemory() const;

  Q_DISABLE_COPY(GmDiscontinuitySet);

  QString           _id;          
  QString           _description; 
  DiscontinuityType _type;        
  GmCellMesh*       _mesh;        
  GmSpatialIndex*   _index;       
  bool              _addElements; 
  bool              _addedInterfElemHaveExtraDofs; 
  double            _snapTol;                      
  bool              _noInternalNodes;              
  bool              _internalIntersections;        
  bool              _discardOriginalGeometry;      

  QVector<GmValueSet*>      _attributes;    
  QVector<GmPropertySet*>   _propertySets;  
  QVector<GmDiscontinuity*> _data;          
  QMap<QString, int>        _attributesMap; 
  QMap<QString, int>        _dataMap;       

  QMap<int, SegmentCell2DIntersectionList> _elementPolylineIntersectionsMap;
};


#endif // _GEMA_DISCONTINUITY_SET_H_