Error control technique of quadrature-based algorithms for the action of real powers of a Hermitian positive-definite matrix

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Otsuka, Motohiro, Tatsuoka, Fuminori, Sogabe, Tomohiro, Takeda, Kota, Zhang, Shao-Liang
Format: Preprint
Veröffentlicht: 2026
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866917384112046080
author Otsuka, Motohiro
Tatsuoka, Fuminori
Sogabe, Tomohiro
Takeda, Kota
Zhang, Shao-Liang
author_facet Otsuka, Motohiro
Tatsuoka, Fuminori
Sogabe, Tomohiro
Takeda, Kota
Zhang, Shao-Liang
contents This study considers quadrature-based algorithms to compute $A^α\boldsymbol{b}$, the action of a real power of a Hermitian positive-definite matrix $A$ on a vector $ \boldsymbol{b}$. In these algorithms, the computation of an integral representation of $A^α \boldsymbol{b}$ is reduced to solving several tens or hundreds of shifted linear systems. Current approaches usually analyze the quadrature discretization error, but rarely take into account the additional error introduced by solving these shifted linear systems with iterative solvers. Here, we bound this error with the residual of the approximated solution of these linear systems. This allows the derivation of a stopping criterion for iterative solvers to keep the error of $A^α\boldsymbol{b}$ below a prescribed error tolerance. Numerical results demonstrate that the proposed criterion enables the computation of $A^α\boldsymbol{b}$ within prescribed tolerance limits.
format Preprint
id arxiv_https___arxiv_org_abs_2604_03923
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Error control technique of quadrature-based algorithms for the action of real powers of a Hermitian positive-definite matrix
Otsuka, Motohiro
Tatsuoka, Fuminori
Sogabe, Tomohiro
Takeda, Kota
Zhang, Shao-Liang
Numerical Analysis
This study considers quadrature-based algorithms to compute $A^α\boldsymbol{b}$, the action of a real power of a Hermitian positive-definite matrix $A$ on a vector $ \boldsymbol{b}$. In these algorithms, the computation of an integral representation of $A^α \boldsymbol{b}$ is reduced to solving several tens or hundreds of shifted linear systems. Current approaches usually analyze the quadrature discretization error, but rarely take into account the additional error introduced by solving these shifted linear systems with iterative solvers. Here, we bound this error with the residual of the approximated solution of these linear systems. This allows the derivation of a stopping criterion for iterative solvers to keep the error of $A^α\boldsymbol{b}$ below a prescribed error tolerance. Numerical results demonstrate that the proposed criterion enables the computation of $A^α\boldsymbol{b}$ within prescribed tolerance limits.
title Error control technique of quadrature-based algorithms for the action of real powers of a Hermitian positive-definite matrix
topic Numerical Analysis
url https://arxiv.org/abs/2604.03923