Shifted CholeskyQR for sparse matrices

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Guan, Haoran, Fan, Yuwei
Format: Preprint
Veröffentlicht: 2024
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866918163399049216
author Guan, Haoran
Fan, Yuwei
author_facet Guan, Haoran
Fan, Yuwei
contents In this work, we focus on Shifted CholeskyQR (SCholeskyQR) for sparse matrices. We provide a new shifted item $s$ for Shifted CholeskyQR3 (SCholeskyQR3) based on the number of non-zero elements (nnze) and the element with the largest absolute value of the input sparse $X \in \mathbb{R}^{m\times n}$ with $m \ge n$. We do rounding error analysis of SCholeskyQR3 with such an $s$ and show that SCholeskyQR3 is accurate in this case. Therefore, an alternative choice of $s$ can be taken for SCholeskyQR3 with the comparison between our new $s$ and the $s$ shown in the previous work when the input $X$ is sparse, improving the applicability and residual of the algorithm for the ill-conditioned cases. Numerical experiments demonstrate the advantage of SCholeskyQR3 with our alternative choice of $s$ in both applicablity and accuracy over the case with the original $s$, together with the same level of efficiency. This work is also the first to build connections between sparsity and numerical algorithms with detailed rounding error analysis to the best of our knowledge.
format Preprint
id arxiv_https___arxiv_org_abs_2410_06525
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Shifted CholeskyQR for sparse matrices
Guan, Haoran
Fan, Yuwei
Numerical Analysis
65F25, 15A23, 65F50, 65G50
In this work, we focus on Shifted CholeskyQR (SCholeskyQR) for sparse matrices. We provide a new shifted item $s$ for Shifted CholeskyQR3 (SCholeskyQR3) based on the number of non-zero elements (nnze) and the element with the largest absolute value of the input sparse $X \in \mathbb{R}^{m\times n}$ with $m \ge n$. We do rounding error analysis of SCholeskyQR3 with such an $s$ and show that SCholeskyQR3 is accurate in this case. Therefore, an alternative choice of $s$ can be taken for SCholeskyQR3 with the comparison between our new $s$ and the $s$ shown in the previous work when the input $X$ is sparse, improving the applicability and residual of the algorithm for the ill-conditioned cases. Numerical experiments demonstrate the advantage of SCholeskyQR3 with our alternative choice of $s$ in both applicablity and accuracy over the case with the original $s$, together with the same level of efficiency. This work is also the first to build connections between sparsity and numerical algorithms with detailed rounding error analysis to the best of our knowledge.
title Shifted CholeskyQR for sparse matrices
topic Numerical Analysis
65F25, 15A23, 65F50, 65G50
url https://arxiv.org/abs/2410.06525