我正在尝试使用 ARPACK 对稀疏 Hermitian 矩阵进行对角化。

我的程序执行没有错误，但返回不正确的特征值。我试图追溯问题。

这项工作由两个功能完成：

znaupd以及我定义av的用于执行矩阵向量乘法的函数。

znaupd将向量作为输入并将结果向量传递给av. ARPACK 选择了一些起始默认向量。

查看我的程序的输出，似乎默认向量基本上被选为零向量，并且每次都znaupd将相同的零向量传递给它。av