Accurate computations with collocation matrices of the Lupaş-type (p,q)-analogue of the Bernstein basis

Abstract

A fast and accurate algorithm to compute the bidiagonal decomposition of collocation matrices of the Lupaş-type (p,q)-analogue of the Bernstein basis is presented. The error analysis of the algorithm and the perturbation theory for the bidiagonal decomposition are also included. Starting from this bidiagonal decomposition, the accurate and efficient solution of several linear algebra problems involving these matrices is addressed: linear system solving, eigenvalue and singular value computation, and computation of the inverse and the Moore-Penrose inverse. The numerical experiments carried out show the good behaviour of the algorithm.