数值计算概述
项目 | 说明 | 备注 |
---|---|---|
高精度数值积分 | DE积分方法 | 计算库伦势能矩阵: 隐式计算矩阵和向量乘法 |
迭代求解超大规模稀疏矩阵 | 动能算符对应矩阵可以分解成两个稀疏矩阵和,采用CGS+预处理器方法迭代求解 | 动能算符主矩阵是个超级三对角阵,可表示为Kronecker张量形式,采用LU分解,或是直接高斯消元都可求解,用作CGS中的预处理矩阵M |
ARPACK/PARPACK迭代求解特征值 | 用迭代法求解最大特征值 | 无需矩阵显示形式,提供矩阵x向量操作即可,支持MPI并行模式 |
库伦矩阵x向量 | 利用Jacobi多项式递推关系和DE积分隐式计算 | 无需显示计算 |
收敛性分析 | 对渐增的截断特征值结果进行外推插值,获得高精度结果 | 寻找序列渐进规律,用E算法外推 |
混合精度计算策略 | 用低精度(REAL(4)) 计算近似特征向量,用高精度(REAL(8)/REAL(16))矫正 | 可明显加速计算效率 |
GPU加速策略 | 用GPU计算库伦矩阵向量乘法 | Jacobi多项式递推关系、块状矩阵乘法 |