共轭梯度法的l~2模收敛性研究

讨论了共轭梯度法(Conjugate Gradient,CG)在l2模意义下的单调性与收敛性.所得结论对结合共轭梯度法来求解大规模线性方程组的各类方法,特别是对研究外推瀑布式多重网格法(EXCMG)按l2模的收敛性有重要作用.     

新外推多网格法研究

提出一种求解大规模有限元方程组的新外推多网格法,即用前两层网格上有限元的新外推值提供第3层细网格上的好初值,然后用几次CG-迭代能得到较准确的有限元解.数值试验表明此方法有很高的效率.     

基于局部加密等级网格的2.5D直流电法有限元模拟

合理截取半圆形计算区域,采取局部加密的?-等级网格,结合对称行索引存贮格式(CSR)及并行稀疏直接求解器PARDISO,提出一种高效、高精度的2.5D直流电阻率法有限元正演方案,并编制相应的Fortran程序,对具有解析解的3个典型地电模型进行计算与分析。研究结果表明:圆形截断边界不仅便于在径向方向上采取?-等级网格剖分,而且能大大简化有限元模拟中单元刚度矩阵的计算;结构化的等级网格避开了通常非结构化网格有限元计算时繁琐的网格剖分及总体刚度阵的集成过程,且能在不增大问题规模的前提下,显著提高2.5D直流电法正演源点附近的模拟精度;Intel MKL的PARDISO求解器能在普通PC机上5 s内求解电法正演有限元离散得到的100万阶稀疏线性方程组,可广泛用于各种地球物理正演问题。