迭代法中压缩对角存储的应用框架.pdfVIP

2013 年 1 JI 精操他探计算提术 第 35 卷第 1 期 文章编号: 1001一1749(2013)01--0107--05 迭代法中压缩对角存储的应用框架 邹桂红，梁华庆，尹洪东 (中国石油大学 地球物理与信息工程学院，北京 102249) 摘 要:在许多有限元、有限差分的应用中，待解方程组的系数矩阵是大型稀疏带状阵，方程组 的求解一般使用迭代法。与其它存储方式相比，压缩对角存储由于不存矩阵元素的行列索引，对 内存的使用最为节省。使用压缩对角存储，沿对角线操作可以完成高效的矩阵一向量来。由于 在以往的文献中，没有提及按行、列操作的算法，压缩对角存储的应用范围受到一定的限制。用 行、列、对角版的矩碎一向量来代表普通意义的行、列、对角方向操作模式，通过等价矩阵推导和 伪代码，给出了一个包括按行、列操作算法在内的应用框架。这里运用 C 语言实现不完全Chol­ esky 分解共辄梯度法