Random sketching to enhance the numerical stability of block orthogonalization algorithms for s-step GMRES

IF 1.1 4区 计算机科学 Q2 COMPUTER SCIENCE, THEORY & METHODS
Parallel Computing Pub Date : 2026-06-01 Epub Date: 2026-05-22 DOI:10.1016/j.parco.2026.103197
Ichitaro Yamazaki , Andrew J. Higgins , Erik G. Boman , Daniel B. Szyld
{"title":"Random sketching to enhance the numerical stability of block orthogonalization algorithms for s-step GMRES","authors":"Ichitaro Yamazaki ,&nbsp;Andrew J. Higgins ,&nbsp;Erik G. Boman ,&nbsp;Daniel B. Szyld","doi":"10.1016/j.parco.2026.103197","DOIUrl":null,"url":null,"abstract":"<div><div>We integrate random sketching techniques into block orthogonalization schemes needed for <span><math><mi>s</mi></math></span>-step GMRES. The resulting, one-stage and two-stage, block orthogonalization schemes generate the basis vectors whose overall orthogonality error is bounded by machine precision as long as each of the corresponding block vectors are numerically full rank. We implement these randomized block orthogonalization schemes using standard distributed-memory linear algebra kernels for <span><math><mi>s</mi></math></span>-step GMRES available in the Trilinos software packages. Our performance results on the Perlmutter supercomputer (with four NVIDIA A100 GPUs per node) demonstrate that these randomized techniques can enhance the numerical stability of the orthogonalization and overall solver, without a significant increase in the execution time.</div></div>","PeriodicalId":54642,"journal":{"name":"Parallel Computing","volume":"128 ","pages":"Article 103197"},"PeriodicalIF":1.1000,"publicationDate":"2026-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Parallel Computing","FirstCategoryId":"94","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0167819126000153","RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"2026/5/22 0:00:00","PubModel":"Epub","JCR":"Q2","JCRName":"COMPUTER SCIENCE, THEORY & METHODS","Score":null,"Total":0}
引用次数: 0

Abstract

We integrate random sketching techniques into block orthogonalization schemes needed for s-step GMRES. The resulting, one-stage and two-stage, block orthogonalization schemes generate the basis vectors whose overall orthogonality error is bounded by machine precision as long as each of the corresponding block vectors are numerically full rank. We implement these randomized block orthogonalization schemes using standard distributed-memory linear algebra kernels for s-step GMRES available in the Trilinos software packages. Our performance results on the Perlmutter supercomputer (with four NVIDIA A100 GPUs per node) demonstrate that these randomized techniques can enhance the numerical stability of the orthogonalization and overall solver, without a significant increase in the execution time.
随机绘制提高s步GMRES分块正交化算法数值稳定性
我们将随机素描技术整合到s步GMRES所需的块正交化方案中。所得到的一阶段和两阶段块正交化方案产生的基向量,只要每个相应的块向量在数值上是满秩的,其总体正交性误差就以机器精度为界。我们使用Trilinos软件包中提供的s步GMRES的标准分布式内存线性代数内核来实现这些随机块正交化方案。我们在Perlmutter超级计算机(每个节点有四个NVIDIA A100 gpu)上的性能结果表明,这些随机化技术可以增强正交化和整体求解器的数值稳定性,而不会显著增加执行时间。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
Parallel Computing
Parallel Computing 工程技术-计算机:理论方法
CiteScore
3.50
自引率
7.10%
发文量
49
审稿时长
4.5 months
期刊介绍: Parallel Computing is an international journal presenting the practical use of parallel computer systems, including high performance architecture, system software, programming systems and tools, and applications. Within this context the journal covers all aspects of high-end parallel computing from single homogeneous or heterogenous computing nodes to large-scale multi-node systems. Parallel Computing features original research work and review articles as well as novel or illustrative accounts of application experience with (and techniques for) the use of parallel computers. We also welcome studies reproducing prior publications that either confirm or disprove prior published results. Particular technical areas of interest include, but are not limited to: -System software for parallel computer systems including programming languages (new languages as well as compilation techniques), operating systems (including middleware), and resource management (scheduling and load-balancing). -Enabling software including debuggers, performance tools, and system and numeric libraries. -General hardware (architecture) concepts, new technologies enabling the realization of such new concepts, and details of commercially available systems -Software engineering and productivity as it relates to parallel computing -Applications (including scientific computing, deep learning, machine learning) or tool case studies demonstrating novel ways to achieve parallelism -Performance measurement results on state-of-the-art systems -Approaches to effectively utilize large-scale parallel computing including new algorithms or algorithm analysis with demonstrated relevance to real applications using existing or next generation parallel computer architectures. -Parallel I/O systems both hardware and software -Networking technology for support of high-speed computing demonstrating the impact of high-speed computation on parallel applications
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术官方微信
小红书