integration/quadratureBase.hh

Go to the documentation of this file.

/* 1D quadrature rules
 */

#ifndef quadratureBase_hh
#define quadratureBase_hh

#include "basics/typedefs.hh"

namespace concepts {

  // ************************************************************ C routines **

#ifdef __cplusplus
  extern "C" {
#endif

    void JacobiPol(const double alf, const double bet, const int maxn,
       const double* x, const int m, double* p);

    void JacobiDerivatives(const double alf, const double bet, const int maxn,
         const double* x, const int m, const double* p,
         double* q);

    void JacobiZeros(double* x, int p, double alf, double bet);

    void GaussLobattoAbscWght(double* x, double* w, const uint p,
            const uint j = 0);

    void GaussRadauAbscWght(double* x, double* w, const uint p,
          const uint j = 0);

    void GaussJacobiAbscWght(double* x, double* w, const uint p);

#ifdef __cplusplus
  }
#endif

} // namespace concepts

#endif // quadratureBase_hh