gmTriIntegrationRule.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_TRI_INT_RULE_H_
#define _GEMA_TRI_INT_RULE_H_
 
#include "gmIntegrationRule.h"
#include "gmBorderIntegrationRule.h"
#include "gmBarIntegrationRule.h" // Needed by edge rules
 
#include "gmGaussIntegrationTables.h"
#include "gmNewtonCotesIntegrationTables.h"
 
 //--------------------------------------------------------
 // Element rules
 //--------------------------------------------------------
 
 
class GMC_API_EXPORT GmTriGaussIntegrationRule : public GmBarycentricSingleTableIntegrationRule<3>
{
public:
  GmTriGaussIntegrationRule(int rule)
    : GmBarycentricSingleTableIntegrationRule<3>(GM_GAUSS_RULE_TYPE, GM_TRI_GIR_KEY, rule,
                                                 GmGaussIntegrationTables::triRuleTable(rule), 
                                                 GmGaussIntegrationTables::triRuleDegree(rule))
  {
  }
 
  GmTriGaussIntegrationRule(GmIntegrationRuleType irType, int rule1, int rule2, int rule3)
    : GmTriGaussIntegrationRule(rule1)
  {
    assert(irType == GM_GAUSS_RULE_TYPE); Q_UNUSED(irType); Q_UNUSED(rule2); Q_UNUSED(rule3);
  }
};
 
class GMC_API_EXPORT GmTriNewtonIntegrationRule : public GmBarycentricSingleTableIntegrationRule<3>
{
public:
  GmTriNewtonIntegrationRule(int rule, bool closed)
    : GmBarycentricSingleTableIntegrationRule<3>(closed ? GM_CLOSED_NEWTON_RULE_TYPE : GM_OPEN_NEWTON_RULE_TYPE, 
                                                 GM_TRI_NIR_KEY, rule,
                                                 GmNewtonIntegrationTables::triRuleTable(rule, closed), 
                                                 GmNewtonIntegrationTables::triRuleDegree(rule))
  {
  }
 
  GmTriNewtonIntegrationRule(GmIntegrationRuleType irType, int rule1, int rule2, int rule3)
    : GmTriNewtonIntegrationRule(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 GmTriGaussEdgeIntegrationRule : public GmProxyBorderIntegrationRule
{
public:
  GmTriGaussEdgeIntegrationRule(GmCellType type, int rule)
    : GmProxyBorderIntegrationRule(type, 3, new GmLineGaussIntegrationRule(rule), false) {}
 
  GmTriGaussEdgeIntegrationRule(GmCellType type, GmIntegrationRuleType irType, int rule)
    : GmTriGaussEdgeIntegrationRule(type, rule) { Q_UNUSED(irType); }
};
 
 
#endif