Hostname: page-component-6766d58669-zlvph Total loading time: 0 Render date: 2026-05-19T10:17:15.993Z Has data issue: false hasContentIssue false

Fast algorithms using orthogonal polynomials

Published online by Cambridge University Press:  30 November 2020

Sheehan Olver
Affiliation:
Department of Mathematics, Imperial College, London SW7 2AZ, UK E-mail: s.olver@imperial.ac.uk
Richard Mikaël Slevinsky
Affiliation:
Department of Mathematics, University of Manitoba, Winnipeg, Manitoba, Canada R3T 2M8 E-mail: Richard.Slevinsky@umanitoba.ca
Alex Townsend
Affiliation:
Department of Mathematics, Cornell University, Ithaca, NY 14853, USA E-mail: townsend@cornell.edu

Abstract

We review recent advances in algorithms for quadrature, transforms, differential equations and singular integral equations using orthogonal polynomials. Quadrature based on asymptotics has facilitated optimal complexity quadrature rules, allowing for efficient computation of quadrature rules with millions of nodes. Transforms based on rank structures in change-of-basis operators allow for quasi-optimal complexity, including in multivariate settings such as on triangles and for spherical harmonics. Ordinary and partial differential equations can be solved via sparse linear algebra when set up using orthogonal polynomials as a basis, provided that care is taken with the weights of orthogonality. A similar idea, together with low-rank approximation, gives an efficient method for solving singular integral equations. These techniques can be combined to produce high-performance codes for a wide range of problems that appear in applications.

Information

Type
Research Article
Copyright
© The Author(s), 2020. Published by Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

Article purchase

Temporarily unavailable