numerics/gauss_quadrature.h

// -*- c++ -*-

// +--------------------------------------------------------------------+
// | This file is part of MathTL - the Mathematical Template Library    |
// |                                                                    |
// | Copyright (c) 2002-2009                                            |
// | Thorsten Raasch, Manuel Werner                                     |
// +--------------------------------------------------------------------+

#ifndef _MATHTL_GAUSS_QUADRATURE_H
#define _MATHTL_GAUSS_QUADRATURE_H

#include <utils/array1d.h>
#include <numerics/quadrature.h>
#include <numerics/ortho_poly.h>

namespace MathTL
{
  class GaussLegendreRule
    : public QuadratureRule<1>
  {
  public:
    GaussLegendreRule(const unsigned int N);
  };

  class GaussRule
    : public QuadratureRule<1>
  {
  public:
    GaussRule(const OrthogonalPolynomial& P,
              const double a, const double b,
              const unsigned int N);

    GaussRule(const Array1D<double>& moments,
              const double a, const double b,
              const unsigned int N);

    GaussRule(const Array1D<double>& moments,
              const OrthogonalPolynomial& T,
              const double a, const double b,
              const unsigned int N);

  private:
    void init(const OrthogonalPolynomial& P,
              const double a, const double b,
              const unsigned int N);
  };
}

// include implementation of inline functions
#include <numerics/gauss_quadrature.cpp>

#endif