xfemMeshCopier.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera 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_XFEM_MESH_COPIER_H_
#define _GEMA_XFEM_MESH_COPIER_H_
 
class XfemMesh;
class XfemFracture;
class XfemElement;
class XfemSolver;
 
class GmCell;
class GmElementMesh;
class GmLogCategory;
 
class GmValueAccessor;
class GmCellAccessor;
class GmGaussAccessor;
 
#include <QVector>
#include <gmVector.h>
#include <gmMatrix.h>
#include <xfemEnrichedElementData.h>
 
class XfemMeshCopier
{
public:
        XfemMeshCopier(const XfemMesh* srcMesh, GmElementMesh* newMesh, XfemSolver* solver, const GmLogCategory& logger);
 
  virtual ~XfemMeshCopier();
 
  bool copy(bool splitFractures);
  bool copy(const QStringList& nodeAttributes, const QStringList& elementAttributes, 
            bool createMissingAttributes, bool splitFractures);
 
private:
  bool copyAux(const QStringList& nodeAttributes, const QStringList& elementAttributes, 
               bool createMissingAttributes, bool splitFractures);
 
  bool setupAccessors(const QStringList& nodeAttributes, const QStringList& elementAttributes, bool createMissingAttributes);
  bool setupNodeAccessors   (const QStringList& nodeAttributes,    bool createMissingAttributes);
  bool setupElementAccessors(const QStringList& elementAttributes, bool createMissingAttributes);
  bool setupEnrichedAccessors();
 
  bool setupPropertyMap(QList<int>& psIndex);
 
  struct SubElementNodeData
  {
    int          _elem;    
    int          _subElem; 
    QVector<int> _nodes;   
 
    SubElementNodeData(int e, int se) : _elem(e), _subElem(se) {}
  };
 
  struct NewNodesData
  {
    int      _elem;         
    GmVector _naturalCoord; 
 
    NewNodesData(int e, const GmVector& coord) : _elem(e), _naturalCoord(coord) {}
  };
 
  bool buildSubElementNodeLists(QList<SubElementNodeData>& nodeDataList, QList<NewNodesData>& newNodesList);
 
  bool setNodeData(const QList<NewNodesData>& newNodesList, bool updateGhostNodes);
  bool setElementData(const QList<SubElementNodeData>& nodeDataList, const QList<int>& psIndex, const int* typeIndexMap);
 
  typedef QPair<int, int> SubElementId;
 
  typedef QPair<const XfemFracture*, int> FractureEdgeId;
 
  struct SplitEdgeData
  {
    int _srcNode;         
    int _inElementIndex;  
    int _outElementIndex; 
    int _fractureIndex;   
 
    SplitEdgeData(int srcNode) { _srcNode = srcNode; _inElementIndex = _outElementIndex = -1;  _fractureIndex = -1; }
  };
 
  struct SplitNodeData
  {
    bool                  _borderNode;  
    QList<FractureEdgeId> _fractures;   
    QList<SubElementId>   _elements;    
    QList<SplitEdgeData>  _edgeList;    
 
    SplitNodeData() { _borderNode = false; }
  };
 
  bool splitMesh(const int* typeIndexMap);
 
  void fillSplitInfoMaps(const int* typeIndexMap, QMap<int, SplitNodeData>& nodeInfo,
                         QMap<SubElementId, int>& elementTranslationMap) const;
 
  void fillEdgeList(int node, const QList<SubElementId>& elements, const QList<FractureEdgeId>& fractures,
                    const QMap<SubElementId, int>& elementTranslationMap, QList<SplitEdgeData>& edges) const;
  int  sortEdgeList(int node, QList<SplitEdgeData>& edges, int* lastEdge) const;
 
  void replaceElementNode(int node, int newNode, const QList<int>& elements);
 
  void addEdgeInfo(int node, int srcNode, bool incoming, int elemIndex, const QList<FractureEdgeId>& fractures, QList<SplitEdgeData>& edges) const;
 
  void addNodeInfo(QMap<int, SplitNodeData>& nodeInfo, QMap<SubElementId, int>& elementTranslationMap,
                   const XfemFracture* f, int fractureSegment, const XfemElement* e, int node, int nodeEdge) const;
 
  const XfemElement* segmentInfo(const XfemFracture* f, int segment, int* n1, int* n2, int* e1 = NULL, int* e2 = NULL) const;
 
  void updateNodeAttributes(int srcNode, int newNode, const FractureEdgeId& fracture, 
                            const SubElementId* nodeElement, const SubElementId* otherElement);
 
  void interpolateSplittedNodeValue(GmValueAccessor* ac, QList<GmValueAccessor*> eacList, int srcNode, int newNode, const FractureEdgeId& fracture,
                                    const SubElementId* nodeElement, const SubElementId* otherElement, GmVector& result) const;
 
  //void interpolateSplittedNodeValue(GmValueAccessor* ac, GmValueAccessor* eAc, int srcNode, int newNode, const FractureEdgeId& fracture,
         // const SubElementId* nodeElement, const SubElementId* otherElement, GmVector& result) const;
 
  void interpolateInternalNodeValue(const XfemElement* e, GmValueAccessor* ac, QList<GmValueAccessor*> eacList,
                                    const GmVector& naturalCoord, GmVector& result) const;
 
  //void interpolateInternalNodeValue(const XfemElement* e, GmValueAccessor* ac, GmValueAccessor* eAc,
         // const GmVector& naturalCoord, GmVector& result) const;
 
  void polygonCentroid(GmMatrix polygonCoords, GmVector& Xcg) const;
 
  void fillNodalLevelSetMatrix(GmMatrix X, GmMatrix Xg, GmMatrix& H) const;
 
  void fillGaussLevelSetMatrix(GmMatrix X, GmMatrix Xg, GmMatrix& H) const;
 
  void fillInteractingCraksNodalLevelSets(const XfemElement* e, XfemEnrichedElementData* exData, int nEnrich, int nCracks, QList<double>crackIndex, QList<GmMatrix>& nodalLevelSet)const;
 
  void fillInteractingCraksGaussLevelSets(const XfemElement* e, XfemEnrichedElementData* exData, int nEnrich, int nCracks, QList<double>crackIndex, QList<GmMatrix>& gaussLevelSet)const;
 
  const XfemMesh* _srcMesh; 
  //XfemMesh* _srcMesh; //!< The source mesh storing the elements to be copied
  GmElementMesh*  _newMesh; 
  XfemSolver*     _solver;  
  const GmLogCategory&  _logger;  
 
  QList<GmValueAccessor*> _nodeAcList;     
  QList<QList<GmValueAccessor*>> _enrichedAcList; 
  //QList<GmValueAccessor*> _enrichedAcList; //!< The set of corresponding enriched node accessors for each entrie in _nodeAcList (or NULL)
  QList<GmValueAccessor*> _newNodeAcList;  
  QList<GmCellAccessor*>  _cellAcList;     
  QList<GmCellAccessor*>  _newCellAcList;  
  QList<GmGaussAccessor*> _gaussAcList;    
  QList<GmGaussAccessor*> _newGaussAcList; 
 
  GmValueAccessor* _coordAc;    
  GmValueAccessor* _newCoordAc; 
};
 
#endif