一种求解特征值问题的广义共轭梯度算法

摘要

本文基于阻尼块反幂法与子空间投影算法设计了一种求解特征值问题的广义共轭梯度算法,同时也实现了相应的计算软件包.然后对算法和计算过程进行一系列的优化来提高算法的稳定性、计算效率和并行可扩展性,使得本文的算法适合在并行计算环境下求解大规模稀疏矩阵的特征值.所形成的软件包不依赖于矩阵和向量的具体结构,可以应用于任意的矩阵向量结构.针对几种典型矩阵的测试结果表明,本文的算法和软件包不但具有良好的数值稳定性和可扩展性,同时相比于SLEPc软件包中的LOBPCG (locally optimal block preconditioned conjugate gradient)和Jacobi-Davidson解法器有2至6倍的效率提升.软件包的网址是https://github.com/pase2017/GCGE-1.0.