quadrature.h

// Copyright (c) 2020 Chris Richardson
// FEniCS Project
// SPDX-License-Identifier:    MIT
 
#pragma once
 
#include "cell.h"
#include <xtl/xspan.hpp>
#include <utility>
#include <vector>
#include <xtensor/xtensor.hpp>
 
 
namespace basix::quadrature
{
xt::xtensor<double, 2> compute_jacobi_deriv(double a, std::size_t n,
                                            std::size_t nderiv,
                                            const xtl::span<const double>& x);
 
// Computes Gauss-Jacobi quadrature points
std::vector<double> compute_gauss_jacobi_points(double a, int m);
 
std::pair<xt::xarray<double>, std::vector<double>>
compute_gauss_jacobi_rule(double a, int m);
 
std::pair<xt::xarray<double>, std::vector<double>> make_quadrature_line(int m);
 
std::pair<xt::xarray<double>, std::vector<double>>
make_quadrature_triangle_collapsed(std::size_t m);
 
std::pair<xt::xarray<double>, std::vector<double>>
make_quadrature_tetrahedron_collapsed(std::size_t m);
 
std::pair<xt::xarray<double>, std::vector<double>>
make_quadrature(const std::string& rule, cell::type celltype, int m);
 
std::pair<xt::xarray<double>, std::vector<double>> make_gll_line(int m);
 
std::pair<xt::xarray<double>, std::vector<double>> compute_gll_rule(int m);
 
} // namespace basix::quadrature