Skip to main content
CenXiv.org
此网站处于试运行阶段,支持我们!
我们衷心感谢所有贡献者的支持。
贡献
赞助
cenxiv logo > math > arXiv:2504.09891

帮助 | 高级搜索

数学 > 数值分析

arXiv:2504.09891 (math)
[提交于 2025年4月14日 (v1) ,最后修订 2025年4月16日 (此版本, v2)]

标题: 适用于任意奇异系统和最小二乘问题的NR-SSOR右预处理RRGMRES

标题: NR-SSOR right preconditioned RRGMRES for arbitrary singular systems and least squares problems

Authors:Kouta Sugihara, Ken Hayami
摘要: GMRES 被认为可以确定$ A x = b $的最小二乘解,其中$ A \in R^{n \times n} $在任意$ b \in R^n $和初始迭代$ x_0 \in R^n $下不会出现中断,当且仅当$ A $是范围对称的,即 $ R(A^T) = R(A) $, where $ A $ may be singular and $ b $ may not be in the range space $ R(A) $ of $ A $. In this paper, we propose applying the Range Restricted GMRES (RRGMRES) to $ A C A^T z = b $, where $ C \in R^{n \times n} $ is symmetric positive definite. 这确定了在任意(奇异)矩阵$ A \in R^{n \times n} $和$ b, x_0 \in R^n $的情况下,$ A x = b $的最小二乘解$ x = C A^T z $不会中断,并且与将 GMRES、RRGMRES 和 MINRES-QLP 应用于$ A x = b $的不一致问题相比更加稳定和准确,当$ b \notin R(A) $时。 特别是,我们建议将NR-SSOR作为内部迭代右预条件器应用,这在最小二乘问题 $ \min_{x \in R^n} \| b - A x\|_2 $ 对于 $ A \in R^{m \times n} $ 和任意 $ b \in R^m $ 时也表现出高效的性能。数值实验验证了所提出方法的有效性。
摘要: GMRES is known to determine a least squares solution of $ A x = b $ where $ A \in R^{n \times n} $ without breakdown for arbitrary $ b \in R^n $, and initial iterate $ x_0 \in R^n $ if and only if $ A $ is range-symmetric, i.e. $ R(A^T) = R(A) $, where $ A $ may be singular and $ b $ may not be in the range space $ R(A) $ of $ A $. In this paper, we propose applying the Range Restricted GMRES (RRGMRES) to $ A C A^T z = b $, where $ C \in R^{n \times n} $ is symmetric positive definite. This determines a least squares solution $ x = C A^T z $ of $ A x = b $ without breakdown for arbitrary (singular) matrix $ A \in R^{n \times n} $ and $ b, x_0 \in R^n $, and is much more stable and accurate compared to GMRES, RRGMRES and MINRES-QLP applied to $ A x = b $ for inconsistent problems when $ b \notin R(A) $. In particular, we propose applying the NR-SSOR as the inner iteration right preconditioner, which also works efficiently for least squares problems $ \min_{x \in R^n} \| b - A x\|_2 $ for $ A \in R^{m \times n} $ and arbitrary $ b \in R^m $. Numerical experiments demonstrate the validity of the proposed method.
评论: arXiv管理员注释:与arXiv:2310.16442存在文本重叠
主题: 数值分析 (math.NA)
MSC 类: 65F10, 65F08, 65F20, 65F50
引用方式: arXiv:2504.09891 [math.NA]
  (或者 arXiv:2504.09891v2 [math.NA] 对于此版本)
  https://doi.org/10.48550/arXiv.2504.09891
通过 DataCite 发表的 arXiv DOI

提交历史

来自: Ken Hayami [查看电子邮件]
[v1] 星期一, 2025 年 4 月 14 日 05:34:02 UTC (306 KB)
[v2] 星期三, 2025 年 4 月 16 日 05:44:33 UTC (224 KB)
全文链接:

获取论文:

    查看标题为《》的 PDF
  • 查看中文 PDF
  • 查看 PDF
  • HTML(实验性)
  • TeX 源代码
  • 其他格式
许可图标 查看许可
当前浏览上下文:
cs
< 上一篇   |   下一篇 >
新的 | 最近的 | 2025-04
切换浏览方式为:
cs.NA
math
math.NA

参考文献与引用

  • NASA ADS
  • 谷歌学术搜索
  • 语义学者
a 导出 BibTeX 引用 加载中...

BibTeX 格式的引用

×
数据由提供:

收藏

BibSonomy logo Reddit logo

文献和引用工具

文献资源探索 (什么是资源探索?)
连接的论文 (什么是连接的论文?)
Litmaps (什么是 Litmaps?)
scite 智能引用 (什么是智能引用?)

与本文相关的代码,数据和媒体

alphaXiv (什么是 alphaXiv?)
CatalyzeX 代码查找器 (什么是 CatalyzeX?)
DagsHub (什么是 DagsHub?)
Gotit.pub (什么是 GotitPub?)
Hugging Face (什么是 Huggingface?)
带有代码的论文 (什么是带有代码的论文?)
ScienceCast (什么是 ScienceCast?)

演示

复制 (什么是复制?)
Hugging Face Spaces (什么是 Spaces?)
TXYZ.AI (什么是 TXYZ.AI?)

推荐器和搜索工具

影响之花 (什么是影响之花?)
核心推荐器 (什么是核心?)
IArxiv 推荐器 (什么是 IArxiv?)
  • 作者
  • 地点
  • 机构
  • 主题

arXivLabs:与社区合作伙伴的实验项目

arXivLabs 是一个框架,允许合作伙伴直接在我们的网站上开发和分享新的 arXiv 特性。

与 arXivLabs 合作的个人和组织都接受了我们的价值观,即开放、社区、卓越和用户数据隐私。arXiv 承诺这些价值观,并且只与遵守这些价值观的合作伙伴合作。

有一个为 arXiv 社区增加价值的项目想法吗? 了解更多关于 arXivLabs 的信息.

这篇论文的哪些作者是支持者? | 禁用 MathJax (什么是 MathJax?)
  • 关于
  • 帮助
  • contact arXivClick here to contact arXiv 联系
  • 订阅 arXiv 邮件列表点击这里订阅 订阅
  • 版权
  • 隐私政策
  • 网络无障碍帮助
  • arXiv 运营状态
    通过...获取状态通知 email 或者 slack

京ICP备2025123034号