Parallel algorithms for downdating the least squares estimator of the regression model

  • Authors:

  • Petko I. Yanev;Erricos J. Kontoghiorghes

  • Affiliations:

  • Institut d'informatique, Université de Neuchítel, CH-2007 Neuchítel, Switzerland and Faculty of Mathematics and Informatics, University of Plovdiv "Paisii Hilendarski", Bulgaria;Department of Public and Business Administration, University of Cyprus, Nicosia, Cyprus and School of Computer Science and Information Systems, Birkbeck College, University of London, UK

  • Venue:
  • Parallel Computing
  • Year:
  • 2008

Quantified Score

Hi-index 0.00

Visualization

Abstract

Computationally efficient parallel algorithms for downdating the least squares estimator of the ordinary linear regression are proposed. The algorithms, which are based on the QR decomposition, are block versions of sequential Givens strategies and efficiently exploit the triangular structure of the data matrices. The first strategy utilizes only part of the orthogonal matrix which is derived from the QR decomposition of the initial data matrix. The rest of the orthogonal matrix is not updated or explicitly computed. A modification of the parallel algorithm, which explicitly computes the whole orthogonal matrix in the downdated QR decomposition, is also considered. An efficient distribution of the matrices over the processors is proposed. Furthermore, the new algorithms do not require any inter-processor communication. The theoretical complexities are derived and experimental results are presented and analyzed. The parallel strategies are scalable and highly efficient for large scale downdating least squares problems. A new parallel block-hyperbolic downdating strategy is developed. The algorithm is rich in BLAS-3 computations, involves negligible duplicated computations and requires insignificant inter-processor communication. It is found to outperform the previous downdating strategies and to be highly efficient for large scale problems. The experimental results confirm the derived theoretical complexities.