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