xfemIntegrationRule.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_INTRULE_H_
#define _GEMA_XFEM_INTRULE_H_
 
#include <gmCellType.h>
#include <gmIntegrationRule.h>
#include <gmMatrix.h>
 
#include <QVector>
 
class GmShape;
class XfemMesh;
 
class XfemIntegrationRule: public GmIntegrationRule
{
public:
 
  XfemIntegrationRule(int ncoord, int npoints = 0);
 
  virtual ~XfemIntegrationRule();
 
  virtual bool isValid() const { return numPoints() > 0; }
 
  // See comments on the base class
  virtual int numPoints() const { return _ipList.count(); }
 
  virtual void integrationPoint(int index, GmVector& naturalCoord, double* weight) const;
 
  virtual int subElement(int index) const { return _subElemList.at(index); }
 
  // See comments on the base class
  virtual int numNaturalCoord() const { return _ncoord; }
 
  virtual int degree() const { return 0; }
 
  // Returns the integration rule type. Always returns 0. Can we do better?
  virtual GmIntegrationRuleType ruleType() const { return (GmIntegrationRuleType)0; }
 
  virtual QString ruleName() const { return "xfem"; }
 
  // See comments on the base class
  virtual int cacheKey() const { return GM_USER_IR_KEY + _internalId; /* TODO: Is this really a good implementation? Do we have another way? */ }
 
  // See comments on the base class
  virtual int rulePar(int num) const { Q_UNUSED(num);  return -1; }
 
  // See comments on the base class
  virtual int ruleParNumPoints(int num) const { Q_UNUSED(num); return -1; }
 
  void clearPoints();
  void setNumPoints(int npoints);
 
  void setIntegrationPoint    (int index, const GmVector& naturalCoord, double weight, int subElement);
  void addIntegrationPoint    (const GmVector& naturalCoord, double weight, int subElement);
  void removeIntegrationPoints(int firstIndex, int num, int seOffset);
 
  void addIntegrationPoints(GmCellType type, GmCellType seType, const GmMatrix& seNaturalCoords,
                            const GmIntegrationRule* seIr, int subElement);
  void setIntegrationPoints(int firstIndex, GmCellType type, GmCellType seType, 
                            const GmMatrix& seNaturalCoords, const GmIntegrationRule* seIr, int subElement);
 
private:
  int nextCacheKey() const { static int globalXfemIrId = 0; return globalXfemIrId++; }
 
  void fillSubelementIntegrationPoint(GmCellType type, GmCellType seType, const GmShape* seShape,
                                      const GmMatrix& seNaturalCoords, const GmIntegrationRule* seIr, 
                                      int ipIndex, QPair<GmVector, double>& ipData);
 
  int                               _internalId;  
  int                               _ncoord;      
  QVector<QPair<GmVector, double> > _ipList;      
  QVector<int>                      _subElemList; 
};
 
 
#endif