quadrature.h

// Copyright (c) 2020 Chris Richardson
// FEniCS Project
// SPDX-License-Identifier:    MIT
 
#pragma once
 
#include "cell.h"
#include <utility>
#include <vector>
#include <xtensor/xtensor.hpp>
 
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::pair<xt::xarray<double>, std::vector<double>>
make_quadrature(const quadrature::type rule, cell::type celltype, int m);
 
std::pair<xt::xarray<double>, std::vector<double>>
make_quadrature(cell::type celltype, int m);
 
quadrature::type get_default_rule(cell::type celltype, int m);
 
xt::xtensor<double, 1> get_gll_points(int m);
 
xt::xtensor<double, 1> get_gl_points(int m);
 
} // namespace basix::quadrature