quadrature.h

// Copyright (c) 2020 Chris Richardson
// FEniCS Project
// SPDX-License-Identifier:    MIT
 
#pragma once
 
#include "cell.h"
#include "polyset.h"
#include <array>
#include <concepts>
#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,
};

template <std::floating_point T>
std::array<std::vector<T>, 2> gauss_jacobi_rule(T a, int m);

template <std::floating_point T>
std::array<std::vector<T>, 2> make_quadrature(const quadrature::type rule,
                                              cell::type celltype,
                                              polyset::type polytype, int m);

quadrature::type get_default_rule(cell::type celltype, int m);

template <std::floating_point T>
std::vector<T> get_gll_points(int m);

template <std::floating_point T>
std::vector<T> get_gl_points(int m);
 
} // namespace basix::quadrature