gmLobattoIntegrationTables.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_LOBATTO_INT_TABLES_H_
#define _GEMA_LOBATTO_INT_TABLES_H_

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

#include <QVector>

class GMC_API_EXPORT GmLobattoIntegrationTables
{
public:
  static const QVector<QPair<double, double>>* lineRuleTable(int rule)
  {
    if(rule < 1 || rule > 11)
      return NULL;
    assert(!_lineTable[rule-1].isEmpty());
    return &_lineTable[rule-1];
  }

  static int lineRuleDegree(int rule) { return 2 * rule - 3; }

  static bool initTables();

private:
  static void initLineTable();

  static double jacobi(double eps, int p, double alpha, double beta);
  static double dJacobi(double eps, int order, double alpha, double beta);
  static double rootChevyshev(int indexRoot, int orderPol);
  static void   rootJacobi(int order, double alpha, double beta, GmVector& zeros);
  static void   quadPoints(int Q, GmMatrix & Quad);

  static QVector<QPair<double, double>> _lineTable[11];
};

#endif