改进的SSOR-PCG迭代法在接触问题研究中应用

SSOR-PCG方法对于大型对称正定问题具有很高的求解效率,但采用求解静动力接触问题的L agrange乘子法导致结构刚度矩阵对应L agrange乘子的对角元为零,不满足传统的SSOR-PCG方法的应用条件.为此通过建立联系L agrange乘子的罚函数矩阵,提出了SSOR-PCG罚函数方法,并通过主动自由度和被动自由度的关系,提出了SSOR-PCG变量替换法.数值例题证明SSOR-PCG变量替换法具有良好的精度和效率.

基于集群的EBE-PCG算法在物探中的应用

采用有限元法进行电法勘探时,会产生大型稀疏线性方程组,如何提高方程组的求解效率成为物探研究的关键。提出了在Beowulf集群环境下,采用粗粒度EBE-PCG算法处理物探问题。在集群环境下,该算法具有机器间相互通信少,易于并行实现等优点。实验结果表明,采用EBE-PCG算法获得了良好的并行效果。

关于H-矩阵的H-预处理子（英文）

设A为一实对称正定的严格对角占优矩阵.设A=D-B为A的Jacobi分裂.为了求解线性方程组Ax=b,在新提出的预处理子的基础上,我们采用预处理共轭梯度方法(PCG)来求解该问题.新提出的预处理子Pv=D+νvv^T,其中v=|B|e,e=(1,...,1)~T,ν=v^TBv/||v||_2~4,且ν使||cvv^T-B||_F达到极小.我们得到了预处理矩阵P_v^(-1)A特征值的上下界,它的界比JIN提出的预处理子的界简单紧凑.数值结果表明我们的预处理子的有效性.