作为其他工作的一部分，我需要解决相对较大（~1E5x1E5）和稀疏（在几个块中的每个原始元素中~100 个非零元素）厄米特特征系统。通常只需要很少的特征值+向量，但精度很高。目前我正在使用 Arpack（当精度更可取时，Arnoldi 使用 shift-inverse 或在尺寸很重要时使用频谱折叠）。作为一个选项计划使用 TRLan（厚重启 Lanczos）并尝试切比雪夫滤波而不是频谱折叠。

可能，为此目的存在更新的方法？