gmElementMesh.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_ELEMENT_MESH_H_
#define _GEMA_ELEMENT_MESH_H_
 
#include "gmCell.h"
#include "gmCellMesh.h"
#include "gmIntegrationRule.h"
#include "gmGaussValueSetGroup.h"
 
class GmBorderIntegrationRule;
class GmGaussValueSet;
class GmGaussAccessor;
class GmMeshGaussIndex;
 
class GMC_API_EXPORT GmElementMesh : public GmCellMesh
{
public:
  GmElementMesh(GmSimulationData* simulation, QString id, QString description, const GmLogCategory& logger);
  virtual ~GmElementMesh();
 
  const GmIntegrationRule*       elementIntegrationRule(GmCellType type, int ruleSet) const;
  const GmBorderIntegrationRule* borderIntegrationRule (GmCellType type, int ruleSet, int borderIndex) const;
  const GmBorderIntegrationRule* faceIntegrationRule   (GmCellType type, int ruleSet, int faceIndex) const;
  const GmBorderIntegrationRule* edgeIntegrationRule   (GmCellType type, int ruleSet) const;
 
  int  maxNumIntPoints(int ruleSet) const;
  void addIntegrationRuleSet();
  int  numIntegrationRuleSets() const;
  bool setElementIntegrationRule(GmCellType type, GmIntegrationRule*       rule, int ruleSet);
  bool setFaceIntegrationRule   (GmCellType type, GmBorderIntegrationRule* rule, int ruleSet, int faceType);
  bool setEdgeIntegrationRule   (GmCellType type, GmBorderIntegrationRule* rule, int ruleSet);
 
  GmMeshGaussIndex* gaussIndex    (int ruleSet) const;
  bool              usesGaussIndex(int ruleSet) const;
 
  GmGaussValueSetGroup* gaussValueSetGroup(int ruleSet) const { return _gaussAttributesGroupList[ruleSet]; }
 
  GmGaussValueSet* gaussValueSet(QString id) const
  {
    int rs = _gaussAttributesMap.value(id, -1);
    return (rs == -1) ? NULL : _gaussAttributesGroupList[rs]->gaussValueSet(id);
  }
 
  virtual QStringList      gaussAttributeIds     () const;
  virtual GmValueInfo*     gaussAttributeInfo    (QString id) const;
  virtual GmGaussAccessor* gaussAttributeAccessor(QString id, Unit desiredUnit, const GmLogCategory& logger) const;
  virtual GmGaussAccessor* gaussAttributeAccessor(QString id, int snum, bool locked, Unit desiredUnit, const GmLogCategory& logger) const;
 
  virtual bool addGaussAttributeSet(GmValueInfo* info);
  virtual void removeGaussAttributeSet(QString id);
 
  virtual bool initGaussAttributeSets();
  virtual void beginAddingGaussAttributeSets();
  virtual bool endAddingGaussAttributeSets();
  virtual void clearGaussAttributeSets();
  virtual void clearGaussAttributeSetsData();
 
  virtual bool    saveState(GmValueSet::SaveStateMode mode, const GmLogCategory& logger);
  virtual bool    saveGaussAttributeState   (QString id, GmValueSet::SaveStateMode mode, const GmLogCategory& logger);
  virtual int     numGaussAttributeStates   (QString id) const;
  virtual QString gaussAttributeStateTag    (QString id, int snum) const;
  virtual double  gaussAttributeStateTime   (QString id, int snum) const;
  virtual void    setGaussAttributeStateTag (QString id, int snum, QString tag);
  virtual void    setGaussAttributeStateTime(QString id, int snum, double time);
 
  virtual int hPOrder() const { return 0; }
 
  virtual int hQOrder() const { return 0; }
 
  virtual void printParameters(const GmLogCategory& logger);
  virtual void printValues    (const GmLogCategory& logger, int options = GM_PRINT_DEFAULT);
 
  virtual size_t gaussIndicesMemory() const;
 
protected:
  virtual GmGaussValueSetGroup* createGaussValueSetGroup(GmSimulationData* simulation, GmElementMesh* mesh)
  {
    return new GmGaussValueSetGroup(simulation, mesh);
  }
 
  template <class T> T* getIntegrationRuleWithDCL(QList<QAtomicPointer<T>*>& list, QMutex* mutex,
                                                  T* (*instance)(GmCellType, int, int, int, GmIntegrationRuleType),
                                                  GmCellType type, int ruleSet, int faceType) const;
 
  QVector<GmGaussValueSetGroup*> _gaussAttributesGroupList; 
  QMap<QString, int>             _gaussAttributesMap;       
 
  virtual bool addStateAccessorData    (GmStateDump* state, int groupId);
  virtual bool fillDumpControlMapData  (QVariantMap* map, const GmLogCategory& logger);
  virtual bool dumpControlMapDataLoaded(QVariantMap* map, const GmLogCategory& logger);
 
private:
  QString toRuleStr  (GmIntegrationRule* ir) const;
  bool    fromRuleStr(QString r, int* key, GmIntegrationRuleType* type, int* rule1, int* rule2, int* rule3);
  void    clearRuleSetRules(int rule);
 
  mutable QList<QAtomicPointer<GmIntegrationRule>*>       _elemIntegrationRulesList; 
  mutable QList<QAtomicPointer<GmBorderIntegrationRule>*> _faceIntegrationRulesList; 
  mutable QList<QAtomicPointer<GmBorderIntegrationRule>*> _edgeIntegrationRulesList; 
 
  mutable QMutex _elemIntegrationRulesMutex;  
  mutable QMutex _faceIntegrationRulesMutex;  
  mutable QMutex _edgeIntegrationRulesMutex;  
 
  mutable QVector<QAtomicPointer<GmMeshGaussIndex>> _gaussIndexList;   
  mutable QMutex                                    _gaussIndexMutex;  
};
 
#define GmForeachActiveElement(variable, mesh) GmForeachCellHelper(variable, mesh, GmElementMesh, GmElement, true)
 
#define GmForeachElement(variable, mesh, activeOnly)  GmForeachCellHelper(variable, mesh, GmElementMesh, GmElement, activeOnly)
 
#define GmForeachOrderedElement(variable, mesh, activeOnly, modelOrder) \
GmForeachOrderedCellHelper(variable, mesh, GmElementMesh, GmElement, activeOnly, modelOrder ? (mesh)->pcrObject() : NULL)
 
 
#endif