求助请问如何用dsaupd和dseupd求解大型稀疏矩阵的

目前需要求解超大型对称矩阵的几个较低的本征值及对应的本征态，若用dsyev子程序就要处理很多0项太慢了。目前想用arpack中的dsaupd和dseupd求解，但是查到的参数设置都不同，搞不懂代表着什么含义。
目前矩阵向量乘法子程序我已写出，希望有大佬能举例演示一下。