A library of functions for polynomial bases used in spectral element methods using the quadrature rules from
FastGaussQuadrature.jl for Float64 and root finding
via the Newton algorithm for other scalar types (such as BigFloat). The algorithms for interpolation and
differentiation use barycentric weights as described in the book "Implementing Spectral Methods for PDEs"
by David Kopriva. If SymPy.jl/SymEngine.jl
is loaded, symbolic computations using SymPy.Sym/SymEngine.Basic are supported.
A brief tutorial is given as notebook.