#ifndef _GAUSS_H_
#define _GAUSS_H_
#include <petscksp.h>

// Gaussian evaluation points - x-component for gauss5
static const PetscScalar GAUSS_5_X[7] = {
0.33333333333333,
0.47014206410511,
0.47014206410511,
0.05971587178977,
0.10128650732346,
0.10128650732346,
0.79742698535309};

// Gaussian evaluation points - y-component for gauss5
static const PetscScalar GAUSS_5_Y[7] = {
0.33333333333333,
0.47014206410511,
0.05971587178977,
0.47014206410511,
0.10128650732346,
0.79742698535309,
0.10128650732346};

// Gaussian weights for gauss5
static const PetscScalar GAUSS_5_W[7] = {
0.22500000000000,
0.13239415278851,
0.13239415278851,
0.13239415278851,
0.12593918054483,
0.12593918054483,
0.12593918054483};

// number of evaluation points for gauss5
#define GAUSS_5_NG 7
PetscScalar gauss5(PetscScalar (*f)(PetscScalar x, PetscScalar y));

#endif