Randomized conjugate gradient least squares

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Zeng, Yun, Cai, Jian-Feng, Han, Deren, Xie, Jiaxin
Formato: Preprint
Publicado: 2026
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866917529160515584
author Zeng, Yun
Cai, Jian-Feng
Han, Deren
Xie, Jiaxin
author_facet Zeng, Yun
Cai, Jian-Feng
Han, Deren
Xie, Jiaxin
contents We develop a novel randomized conjugate gradient least squares (RCGLS) method for solving least-squares problems, in which iterative sketching is employed at each step to reduce the dimension and hence the computational cost. In particular, we propose a new perspective on the classical CGLS method, where the next descent direction is determined via a constraint correction problem associated with the gradient. Based on this insight, we replace the gradient with a randomized coordinate gradient that naturally satisfies the variance reduction property, leading directly to the proposed RCGLS method. We prove that RCGLS converges linearly in expectation, with a better convergence bound compared to the randomized coordinate descent method. Furthermore, we investigate an implementation of the method that avoids full-dimensional vector operations, which are the major bottleneck of vanilla RCGLS for sparse matrices and render it impractical. We also show how to apply the RCGLS method to solve the ridge regression problem, yielding a lightweight, parallelizable, and accelerated method for such problems. Numerical experiments are provided to confirm our results.
format Preprint
id arxiv_https___arxiv_org_abs_2605_25034
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Randomized conjugate gradient least squares
Zeng, Yun
Cai, Jian-Feng
Han, Deren
Xie, Jiaxin
Numerical Analysis
We develop a novel randomized conjugate gradient least squares (RCGLS) method for solving least-squares problems, in which iterative sketching is employed at each step to reduce the dimension and hence the computational cost. In particular, we propose a new perspective on the classical CGLS method, where the next descent direction is determined via a constraint correction problem associated with the gradient. Based on this insight, we replace the gradient with a randomized coordinate gradient that naturally satisfies the variance reduction property, leading directly to the proposed RCGLS method. We prove that RCGLS converges linearly in expectation, with a better convergence bound compared to the randomized coordinate descent method. Furthermore, we investigate an implementation of the method that avoids full-dimensional vector operations, which are the major bottleneck of vanilla RCGLS for sparse matrices and render it impractical. We also show how to apply the RCGLS method to solve the ridge regression problem, yielding a lightweight, parallelizable, and accelerated method for such problems. Numerical experiments are provided to confirm our results.
title Randomized conjugate gradient least squares
topic Numerical Analysis
url https://arxiv.org/abs/2605.25034