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