In this paper we propose a new fast O(n) algorithm for solving Cauchy linear systems of equations. An a priori rounding error analysis, which provides bounds for the forward, backward and residual errors associated with the computed solution, indicates that for the class of totally positive Cauchy matrices the proposed algorithm is forward and backward stable, producing a remarkably high relati...

Among the most important aspects of geophysical data interpretation is the estimation and computation of inverse solution uncertainties. We present a general uncertainty estimation method that allows for the comprehensive search of model posterior space while maintaining computational efficiencies similar to deterministic inverse solutions. Integral to this method is the combination of an effic...

MIUR Center for Science and Application of Advanced Computing Paradigms at the University of Padova, Italy Consorzio Roma Ricerche Roma, Italy International Centre for Mechanical Sciences, Udine, Italy Ente per le Nuove tecnologie, l’Energia e l’Ambiente, Roma, Italy Department of Information Engineering University of Padova, Italy ? IBM Power3@375MHz with 4 GB mem ? HPM Toolkit for performance...