quadrature.h

// Copyright (c) 2020 Chris Richardson
// FEniCS Project
// SPDX-License-Identifier:    MIT
 
#pragma once
 
#include "cell.h"
#include <array>
#include <vector>
 
namespace basix::quadrature
{
 
enum class type
{
  Default = 0,
  gauss_jacobi = 1,
  gll = 2,
  xiao_gimbutas = 3,
  zienkiewicz_taylor = 20,
  keast = 21,
  strang_fix = 22,
};
 
std::array<std::vector<double>, 2> make_quadrature(const quadrature::type rule,
                                                   cell::type celltype, int m);
 
std::array<std::vector<double>, 2> make_quadrature(cell::type celltype, int m);
 
quadrature::type get_default_rule(cell::type celltype, int m);
 
std::vector<double> get_gll_points(int m);
 
std::vector<double> get_gl_points(int m);
 
} // namespace basix::quadrature