下一节:Block Size. 上一级:块Arnoldi方法 上一节:块Arnoldi约化

实用算法

算法7.12列出了一种通用的块隐式重启Arnoldi方法(BIRAM)。本节剩余部分将讨论实现稳健软件所需的一些关键实现问题。特别是,我们将探讨块大小(b)、收敛考虑、缩减、隐式重启以及位移选择和p值选取等问题。

算法 7.12:BIRAM用于NHEP
直到收敛,执行以下步骤:
(1) 将长度为r的块Arnoldi约减扩展p块:
AV_{[r+p]} = V_{[r+p]} H_{[r+p]} + F_{r+p} E_{r+p}^* \text{ 通过算法 7.11}
(2) 确定感兴趣的 k 个Ritz值是否足够接近 A 的值。
(3) 锁定(或消去)满足收敛容忍度的Ritz值。
(4) 隐式重启 p 个移位,并保留长度为 r 的块Arnoldi约减。



小节

Susan Blackford 2000-11-20