gmNewtonCotesIntegrationTables.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_NEWTON_INT_TABLES_H_
#define _GEMA_NEWTON_INT_TABLES_H_

#include "gmCoreConfig.h"
#include "gmVector.h"

#include <QVector>

class GMC_API_EXPORT GmNewtonIntegrationTables
{
public:

  static const QVector<QPair<double, double>>* lineRuleTable(int rule, bool closed)
  {
    return select(rule, closed, _closedLineTable, _openLineTable);
  }

  static int lineRuleDegree(int rule) { return rule%2 ? rule : rule-1; }

  static const QVector<QPair<GmVector, double>>* triRuleTable(int rule, bool closed)
  {
    return select(rule, closed, _closedTriTable, _openTriTable);
  }

  static int triRuleDegree(int rule) { return rule; }

  static const QVector<QPair<GmVector, double>>* tetRuleTable(int rule, bool closed)
  {
    return select(rule, closed, _closedTetTable, _openTetTable);
  }

  static int tetRuleDegree(int rule) { return rule; }

  static bool initTables();

private:

  template <class T> static const T* select(int rule, bool closed, const T ct[], const T ot[])
  {

    if(rule < 1 || rule > 5)
      return NULL;
    if(closed)
    {
      assert(!ct[rule-1].isEmpty());
      return &ct[rule-1];
    }
    else
    {
      assert(!ot[rule-1].isEmpty());
      return &ot[rule-1];
    }
  }

  static void initLineTables();
  static void initTriTables();
  static void initTetTables();

  static void   simplexRulesPermutations(const char* s, QStringList& permutations);
  static double simplexRulesPointAndWeight(int degree, const char* index, int u, int num, int den, QVector<GmVector>& posList);


  static QVector<QPair<double, double>> _closedLineTable[5];

  static QVector<QPair<double, double>> _openLineTable[5];

  static QVector<QPair<GmVector, double>> _closedTriTable[5];

  static QVector<QPair<GmVector, double>> _openTriTable[5];

  static QVector<QPair<GmVector, double>> _closedTetTable[5];

  static QVector<QPair<GmVector, double>> _openTetTable[5];
};

#endif