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