算法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约减。