Fast solving of weighted pairing least-squares systems

  • Authors:
  • Pierre Courrieu

  • Affiliations:
  • LPC, UMR 6146, CNRS-University of Provence, Marseille, France

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2009

Quantified Score

Hi-index 7.29

Visualization

Abstract

This paper presents a generalization of the ''weighted least-squares'' (WLS), named ''weighted pairing least-squares'' (WPLS), which uses a rectangular weight matrix and is suitable for data alignment problems. Two fast solving methods, suitable for solving full rank systems as well as rank deficient systems, are studied. Computational experiments clearly show that the best method, in terms of speed, accuracy, and numerical stability, is based on a special {1, 2, 3}-inverse, whose computation reduces to a very simple generalization of the usual ''Cholesky factorization-backward substitution'' method for solving linear systems.