gmGaussIntegrationTables.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_GAUSS_INT_TABLES_H_
#define _GEMA_GAUSS_INT_TABLES_H_
 
#include "gmCoreConfig.h"
#include "gmVector.h"
 
#include <QVector>
 
class GMC_API_EXPORT GmGaussIntegrationTables
{
public:
 
  static const QVector<QPair<double, double>>* lineRuleTable(int rule)
  {
    if(rule < 1 || rule > 5)
      return NULL;
    assert(!_lineTable[rule-1].isEmpty());    
    return &_lineTable[rule-1];
  }
 
  static int lineRuleDegree(int rule) { return 2 * rule - 1; }
 
  static const QVector<QPair<GmVector, double>>* triRuleTable(int rule)
  {
    if(rule != 1 && rule != 3 && rule != 6)
      return NULL;
    assert(!_triTable[rule-1].isEmpty());
    return &_triTable[rule-1];
  }
 
  static int triRuleDegree(int rule) { return (rule == 1) ? 1 : ((rule == 3) ? 2 : 4); }
 
  static const QVector<QPair<GmVector, double>>* tetRuleTable(int rule)
  {
    if(rule != 1 && rule != 4)
      return NULL;
    assert(!_tetTable[rule-1].isEmpty());
    return &_tetTable[rule-1];
  }
 
  static int tetRuleDegree(int rule) { return (rule == 1) ? 1 : 2; }
 
  static const QVector<QPair<GmVector, double>>* pyraRuleTable(int rule)
  {
    if (rule != 1 && rule != 5 && rule != 13)
      return NULL;
    assert(!_pyraTable[rule - 1].isEmpty());
    return &_pyraTable[rule - 1];
  }
 
  static int pyraRuleDegree(int rule) { return (rule == 1) ? 1 : 2; }
 
  static bool initTables();
 
private:
  static void initLineTable();
  static void initTriTable();
  static void initTetTable();
  static void initPyraTable();
 
  static QVector<QPair<double,   double>> _lineTable[5];
 
  static QVector<QPair<GmVector, double>> _triTable[6];
 
  static QVector<QPair<GmVector, double>> _tetTable[4];
 
  static QVector<QPair<GmVector, double>> _pyraTable[13];
};
 
#endif