// QUADRATURE.H
// Gaussian quadrature rules for numerically integrate 'smooth' functions over
// simplexes, triangle for 2d and tetrahedra for 3d.

#ifndef QUADRATURE_H
#define QUADRATURE_H

#define one_third 0.333333333333333
#define two_third 0.666666666666667
#define one_sixth 0.166666666666667
#include "../dgtd-performance.hpp"

void GetFormula(int, int, fp_t *, fp_t *, fp_t *, fp_t *);
void GetFormula2dPtr(int, fp_t (*&)[4]);
void GetFormula3d(int, int, fp_t *, fp_t *, fp_t *, fp_t *, fp_t *);
void GetFormula3dPtr(int, fp_t (*&)[5]);
void GetGauss(int, int, int, fp_t *, fp_t *, fp_t *, fp_t *);
void GetDuffy(int, int, int, int, fp_t *, fp_t *, fp_t *, fp_t *);

#endif