多自由度块行压缩存储技术及大型稀疏方程组的求解.pdfVIP

2∞9 年 3 )费 第 31 .4.第 2 期 辑4率也4率计算就术 1 文章蜡寄: 1001…1749(2 9)02-0108….Q5 ∞ 多自由度块行压缩存储技术及大型稀疏 方程细的求解 张继锋，汤井凹，王烨，肖晓 (中南大学信息物搜工程学挠，湖南长 41∞83) 摘 要:在地球物嘿电磁勘探领域有限元数值模拟中，最后都舍得到一个大型稀疏的提系数域 性方程姐，受计算叔，内存空间的限制，据规根据有?在走，树皮娃阵的稀疏性时其进行压缩存错。由 于电磁场有限元，计算的自由皮大都在三个以上，因而提出了适合多自由皮的块按行路缩稀就存 储方棠，并遏地存储格式的转换，把块按行压缩方式转换成流行的，大型稀腻知阵的行压缩存储 格式，以便于求解。用求掰大型稀疏方程姐的险y!，梆子空尚方法中的稳定双共辄梯皮(Bicg- s阔的方法，收钱边皮快，精度高，而且稳定性好，结合ilu 预处理拭术，可以大大提高点解大型稀 疏方粗细的效华。 关键词:压辅存储:制度报阵;有限元;稀疏挫降;稳定双典辄棉皮 中剧分提号: 0242.21 文献棉识码:A 包，都是基于该存储方束。目前，关于该方法的存 储模式，巳有许部文献陈边。作者在本文中，从有 。 前宵 限克制度矩阵的组装原理出发，分绍了如何捏一个 关于大型稀疏姐阵的存储方式，已有许多学者 有限元刚度知阵，组装成为一个按付压缩存储的一 维数组，然后提出了多自由度的块按行压捕存储模 作了大攘的研究，目前已有的存储格式就达十多种 以上L1410 上个教铠七十年R拐，出现了以二维等 式，并对该模式租单身自度行压缩存锵模式之间的 带宽存储和一维变带宽存储的方案，当时许多有限 格式转换作了研究，并蝙制了相应的程序。大型稀 冗软件都采用了这二种存储方法。但是，由于工维 疏线性方程组的求解也般有限冗计算的一个关髓 问趣，许多学者作了研究协1210 作者通过读处理 帮带宽存储量均半带宽的大小辑部相关，它有可能 因为周部带宽过大，而使整体刚度短阵的存储量大 技术降低了短阵的条件嫩，然后结合求解大型稀疏 大增加;而一维变带宽存储鼠然会比等带宽节省内 线性方程组的Krylov 子空间方法，稳定~共辄棉 度法Bic伊始，应用到地球物理电磁法数值摸摸 存空i碍，但街存储了许多零元素。丽且从某种意义 上说，这;二种存储方案的存储簸，取决于单光和节 中。 点的编号模式。近t年来，稀疏存储方案逐渐取代 了传统的存储方案，该存储模式只需要始黯度矩醉 1 按行恶结银装整体酣皮矩阵 的非等元素分配存储空间，矩阵的运算金部在索引 形式下进行，大大节省了所需内存空间。所以，许 在实际的数值模拟中，常常形成…个大型稀疏