gmTetIntegrationRule.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_TET_INT_RULE_H_
#define _GEMA_TET_INT_RULE_H_

#include "gmIntegrationRule.h"
#include "gmBorderIntegrationRule.h"
#include "gmBarIntegrationRule.h"  // Needed by edge rules
#include "gmTriIntegrationRule.h"  // Needed by face rules

#include "gmGaussIntegrationTables.h"
#include "gmNewtonCotesIntegrationTables.h"

//--------------------------------------------------------
// Element rules
//--------------------------------------------------------


class GMC_API_EXPORT GmTetGaussIntegrationRule : public GmBarycentricSingleTableIntegrationRule<4>
{
public:
  GmTetGaussIntegrationRule(int rule)
    : GmBarycentricSingleTableIntegrationRule<4>(GM_GAUSS_RULE_TYPE, GM_TET_GIR_KEY, rule,
                                                 GmGaussIntegrationTables::tetRuleTable(rule), 
                                                 GmGaussIntegrationTables::tetRuleDegree(rule))
  {
  }

  GmTetGaussIntegrationRule(GmIntegrationRuleType irType, int rule1, int rule2, int rule3)
    : GmTetGaussIntegrationRule(rule1)
  {
    assert(irType == GM_GAUSS_RULE_TYPE); Q_UNUSED(irType); Q_UNUSED(rule2); Q_UNUSED(rule3);
  }
};

class GMC_API_EXPORT GmTetNewtonIntegrationRule : public GmBarycentricSingleTableIntegrationRule<4>
{
public:
  GmTetNewtonIntegrationRule(int rule, bool closed)
    : GmBarycentricSingleTableIntegrationRule<4>(closed ? GM_CLOSED_NEWTON_RULE_TYPE : GM_OPEN_NEWTON_RULE_TYPE, 
                                                 GM_TET_NIR_KEY, rule,
                                                 GmNewtonIntegrationTables::tetRuleTable(rule, closed), 
                                                 GmNewtonIntegrationTables::tetRuleDegree(rule))
  {
  }

  GmTetNewtonIntegrationRule(GmIntegrationRuleType irType, int rule1, int rule2, int rule3)
    : GmTetNewtonIntegrationRule(rule1, irType == GM_CLOSED_NEWTON_RULE_TYPE)
  {
    assert(irType == GM_CLOSED_NEWTON_RULE_TYPE || irType == GM_OPEN_NEWTON_RULE_TYPE); Q_UNUSED(rule2); Q_UNUSED(rule3);
  }
};

//--------------------------------------------------------
// Border rules
//--------------------------------------------------------


class GMC_API_EXPORT GmTetGaussFaceIntegrationRule : public GmProxyBorderIntegrationRule
{
public:
  GmTetGaussFaceIntegrationRule(GmCellType type, int rule)
    : GmProxyBorderIntegrationRule(type, 4, new GmTriGaussIntegrationRule(rule), true) {}

  GmTetGaussFaceIntegrationRule(GmCellType type, GmCellFaceType faceType, GmIntegrationRuleType irType, int rule1, int rule2)
    : GmTetGaussFaceIntegrationRule(type, rule1) { Q_UNUSED(faceType); Q_UNUSED(irType); Q_UNUSED(rule2); }
};

class GMC_API_EXPORT GmTetGaussEdgeIntegrationRule : public GmProxyBorderIntegrationRule
{
public:
  GmTetGaussEdgeIntegrationRule(GmCellType type, int rule)
    : GmProxyBorderIntegrationRule(type, 4, new GmLineGaussIntegrationRule(rule), false) {}

  GmTetGaussEdgeIntegrationRule(GmCellType type, GmIntegrationRuleType irType, int rule)
    : GmTetGaussEdgeIntegrationRule(type, rule) { Q_UNUSED(irType); }
};


#endif