English

A MATLAB package computing simultaneous Gaussian quadrature rules for Multiple Orthogonal Polynomials

Numerical Analysis 2024-07-09 v1 Numerical Analysis

Abstract

The aim of this paper is to describe a Matlab package for computing the simultaneous Gaussian quadrature rules associated with a variety of multiple orthogonal polynomials. Multiple orthogonal polynomials can be considered as a generalization of classical orthogonal polynomials, satisfying orthogonality constraints with respect to rr different measures, with r1r \ge 1. Moreover, they satisfy (r+2)(r+2)--term recurrence relations. In this manuscript, without loss of generality, rr is considered equal to 22. The so-called simultaneous Gaussian quadrature rules associated with multiple orthogonal polynomials can be computed by solving a banded lower Hessenberg eigenvalue problem. Unfortunately, computing the eigendecomposition of such a matrix turns out to be strongly ill-conditioned and the \texttt{Matlab} function \texttt{balance.m} does not improve the condition of the eigenvalue problem. Therefore, most procedures for computing simultaneous Gaussian quadrature rules are implemented with variable precision arithmetic. Here, we propose a \texttt{Matlab} package that allows to reliably compute the simultaneous Gaussian quadrature rules in floating point arithmetic. It makes use of a variant of a new balancing procedure, recently developed by the authors of the present manuscript, that drastically reduces the condition of the Hessenberg eigenvalue problem.

Keywords

Cite

@article{arxiv.2406.11269,
  title  = {A MATLAB package computing simultaneous Gaussian quadrature rules for Multiple Orthogonal Polynomials},
  author = {Teresa Laudadio and Nicola Mastronardi and Walter Van Assche and Paul Van Dooren},
  journal= {arXiv preprint arXiv:2406.11269},
  year   = {2024}
}

Comments

32 pages, 1 figure, 6 tables. Includes an appendix with matlab codes

R2 v1 2026-06-28T17:08:14.687Z