Accurate and fast computations with Green matrices✩

This paper provides a linear time complexity method to obtain the bidiagonal decomposition of Green matrices with high relative accuracy. In addition, when the Green matrix is nonsingular and totally positive, this bidiagonal decomposition can be used to compute the eigenvalues, the inverse and the solution of some linear system of equations with high relative accuracy. A numerical example illustrates the advantages of this method. & COPY; 2023 The Author(s). Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).