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"

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;

  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 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);

  const QList<GmGaussValueSet*>& gaussValueSets() const { return _gaussAttributesData; }

  virtual size_t gaussIndicesMemory() const;

protected:
  virtual GmGaussValueSet* createGaussValueSet(GmValueInfo* info);

  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;

  QMap<QString, int>      _gaussAttributesMap;   
  QList<GmGaussValueSet*> _gaussAttributesData;  

  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 QList<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)

#endif