Generation of orthogonal rational functions by procedures for structured matrices

03/08/2021
by   Niel Van Buggenhout, et al.
0

The problem of computing recurrence coefficients of sequences of rational functions orthogonal with respect to a discrete inner product is formulated as an inverse eigenvalue problem for a pencil of Hessenberg matrices. Two procedures are proposed to solve this inverse eigenvalue problem, via the rational Arnoldi iteration and via an updating procedure using unitary similarity transformations. The latter is shown to be numerically stable. This problem and both procedures are generalized by considering biorthogonal rational functions with respect to a bilinear form. This leads to an inverse eigenvalue problem for a pencil of tridiagonal matrices. A tridiagonal pencil implies short recurrence relations for the biorthogonal rational functions, which is more efficient than the orthogonal case. However the procedures solving this problem must rely on nonunitary operations and might not be numerically stable.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/01/2023

Algorithms for Modifying Recurrence Relations of Orthogonal Polynomial and Rational Functions when Changing the Discrete Inner Product

Often, polynomials or rational functions, orthogonal for a particular in...
research
07/12/2023

Fast and stable rational approximation of generalized hypergeometric functions

Rational approximations of generalized hypergeometric functions _pF_q of...
research
09/20/2018

On constructing orthogonal generalized doubly stochastic matrices

A real quadratic matrix is generalized doubly stochastic (g.d.s.) if all...
research
03/06/2021

The Short-term Rational Lanczos Method and Applications

Rational Krylov subspaces have become a reference tool in dimension redu...
research
07/18/2019

Fast approximation of orthogonal matrices and application to PCA

We study the problem of approximating orthogonal matrices so that their ...
research
03/23/2022

Spectral analysis of a mixed method for linear elasticity

The purpose of this paper is to analyze a mixed method for linear elasti...
research
10/19/2020

Twice is enough for dangerous eigenvalues

We analyze the stability of a class of eigensolvers that target interior...

Please sign up or login with your details

Forgot password? Click here to reset