Saved in:
Bibliographic Details
Main Authors: Suzuki, Kengo, Iwashita, Takeshi
Format: Preprint
Published: 2025
Subjects:
Online Access:https://arxiv.org/abs/2505.20719
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909624147378176
author Suzuki, Kengo
Iwashita, Takeshi
author_facet Suzuki, Kengo
Iwashita, Takeshi
contents Low-precision computing is essential for efficiently utilizing memory bandwidth and computing cores. While many mixed-precision algorithms have been developed for iterative sparse linear solvers, effectively leveraging half-precision (fp16) arithmetic remains challenging. This study introduces a novel nested Krylov approach that integrates the flexible GMRES and Richardson methods in a deeply nested structure, progressively reducing precision from double-precision to fp16 toward the innermost solver. To avoid meaningless computations beyond precision limits, the low-precision inner solvers perform only a few iterations per invocation, while the nested structure ensures their frequent execution. Numerical experiments show that using fp16 in the approach directly enhances solver performance without compromising convergence, achieving speedups of up to 1.65x and 2.42x over double-precision and double-single mixed-precision implementations, respectively. Moreover, the proposed method outperforms or matches other standard Krylov solvers, including restarted GMRES, CG, and BiCGStab methods.
format Preprint
id arxiv_https___arxiv_org_abs_2505_20719
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle A Nested Krylov Method Using Half-Precision Arithmetic
Suzuki, Kengo
Iwashita, Takeshi
Numerical Analysis
G.1.3
Low-precision computing is essential for efficiently utilizing memory bandwidth and computing cores. While many mixed-precision algorithms have been developed for iterative sparse linear solvers, effectively leveraging half-precision (fp16) arithmetic remains challenging. This study introduces a novel nested Krylov approach that integrates the flexible GMRES and Richardson methods in a deeply nested structure, progressively reducing precision from double-precision to fp16 toward the innermost solver. To avoid meaningless computations beyond precision limits, the low-precision inner solvers perform only a few iterations per invocation, while the nested structure ensures their frequent execution. Numerical experiments show that using fp16 in the approach directly enhances solver performance without compromising convergence, achieving speedups of up to 1.65x and 2.42x over double-precision and double-single mixed-precision implementations, respectively. Moreover, the proposed method outperforms or matches other standard Krylov solvers, including restarted GMRES, CG, and BiCGStab methods.
title A Nested Krylov Method Using Half-Precision Arithmetic
topic Numerical Analysis
G.1.3
url https://arxiv.org/abs/2505.20719