Skip to Main content Skip to Navigation
Reports

Efficient computation of condition estimates for linear least squares problems

Abstract : Linear least squares (LLS) is a classical linear algebra problem in scientific computing, arising for instance in many parameter estimation problems. In addition to computing efficiently LLS solutions, an important issue is to assess the numerical quality of the computed solution. The notion of conditioning provides a theoretical framework that can be used to measure the numerical sensitivity of a problem solution to perturbations in its data. We recall some results for least squares conditioning and we derive a statistical estimate for the conditioning of an LLS solution. We present numerical experiments to compare exact values and statistical estimates. We also propose performance results using new routines on top of the multicore-GPU library MAGMA. This set of routines is based on an efficient computation of the variance-covariance matrix for which, to our knowledge, there is no implementation in current public domain libraries LAPACK and ScaLAPACK.
Document type :
Reports
Complete list of metadata

Cited literature [23 references]  Display  Hide  Download

https://hal.inria.fr/hal-00731136
Contributor : Marc Baboulin <>
Submitted on : Thursday, September 13, 2012 - 10:04:46 PM
Last modification on : Thursday, July 8, 2021 - 3:47:42 AM
Long-term archiving on: : Friday, December 16, 2016 - 1:02:46 PM

File

RR-8065.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : hal-00731136, version 2

Citation

Marc Baboulin, Serge Gratton, Rémi Lacroix, Alan J. Laub. Efficient computation of condition estimates for linear least squares problems. [Research Report] RR-8065, INRIA. 2012, pp.15. ⟨hal-00731136v2⟩

Share

Metrics

Record views

494

Files downloads

1776