基于坐标分割的聚集型代数多重网格预条件研究-计算机应用与软件.PDF

第３５卷第７期　　　 计算机应用与软件 Ｖｏｌ３５Ｎｏ．７ ２０１８年７月　　 ＣｏｍｐｕｔｅｒＡｐｐｌｉｃａｔｉｏｎｓａｎｄＳｏｆｔｗａｒｅ Ｊｕｌ．２０１８ 基于坐标分割的聚集型代数多重网格预条件研究 吴建平　银福康　彭　军　杨锦辉 （国防科学技术大学海洋科学与工程研究院　湖南长沙４１００７３） 摘　要　　针对基于坐标分割的聚集型代数多重网格预条件，给出了三种进行坐标分割的方法，即正方分割、最 小界面分割与逐步单向分割，并对其进行了高效实现。正方分割以每个子图接近于正方体或正方形的方式进行 分割。最小界面分割遍历所有可能的分割，并以每个子图表面积或周长之和最短的方式进行实际分割。逐步单 向分割以分割数的素因子分解为基础，并按素因子从大到小的顺序，每次沿不同坐标数最大的方向进行分割，直 到所有素因子遍历完为止。之后对从模型偏微分方程离散得到的稀疏线性方程组，通过Ｖ型、Ｗ型与Ｋ型等多 种循环，从多重网格预条件共轭斜量法的效率上，对这三种分割算法进行了实验对比分析。结果表明，逐步单向 分割更适合于Ｊａｃｏｂｉ光滑、Ｋ循环与强各向异性等情形。最小界面分割算法更适合于ＧａｕｓｓＳｅｉｄｅｌ光滑、系数矩 阵具有较多非零元素等情形。 关键词　　稀疏线性方程组　聚集型代数多重网格　预条件　共轭斜量法　坐标分割 中图分类号　ＴＰ３０１　　　　文献标识码　Ａ　　　　ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１０００３８６ｘ．２０１８．０７．０４９ ＲＥＳＥＡＲＣＨＯＦＴＨＥＣＯＯＲＤＩＮＡＴＥＳＰＡＲＴＩＴＩＯＮＩＮＧＢＡＳＥＤＡＧＧＲＥＧＡＴＩＯＮＴＹＰＥ ＡＬＧＥＢＲＡＩＣＭＵＬＴＩＧＲＩＤＰＲＥＣＯＮＤＩＴＩＯＮＳ ＷｕＪｉａｎｐｉｎｇ　ＹｉｎＦｕｋａｎｇ　ＰｅｎｇＪｕｎ　ＹａｎｇＪｉｎｈｕｉ （ＡｃａｄｅｍｙｏｆＯｃｅａｎＳｃｉｅｎｃｅａｎｄＥｎｇｉｎｅｅｒｉｎｇ，ＮａｔｉｏｎａｌＵｎｉｖｅｒｓｉｔｙｏｆＤｅｆｅｎｓｅＴｅｃｈｎｏｌｏｇｙ，Ｃｈａｎｇｓｈａ４１００７３，Ｈｕｎａｎ，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ　　Ｉｎｔｈｉｓｐａｐｅｒ，ｗｅｐｒｏｐｏｓｅｄｔｈｒｅｅｍｅｔｈｏｄｓｆｏｒｃｏｏｒｄｉｎａｔｅｐａｒｔｉｔｉｏｎｉｎｇｂａｓｅｄｏｎａｇｇｒｅｇａｔｅｔｙｐｅａｌｇｅｂｒａｉｃ ｍｕｌｔｉｇｒｉｄｐｒｅｃｏｎｄｉｔｉｏｎｓ，ｉｅ，ｓｑｕａｒｅｐａｒｔｉｔｉｏｎｉｎｇ，ｍｉｎｉｍｕｍｉｎｔｅｒｆａｃｅｐａｒｔｉｔｉｏｎｉｎｇ，ａｎｄｇｒａｄｕａｌｕｎｉｄｉｒｅｃｔｉｏｎａｌｐａｒｔｉｔｉｏｎｉｎｇ， ａｎｄｉｍｐｌｅｍｅｎｔｔｈｅｍｅｆｆｉｃｉｅｎｔｌｙ．Ｓｑｕａｒｅｐａｒｔｉｔｉｏｎｉｎｇｗａｓｐｅｒｆｏｒｍｅｄｉｎｓｕｃｈａｗａｙｔｈａｔｅａｃｈｓｕｂｇｒａｐｈｗａｓｃｌｏｓｅｔｏａｃｕｂｅ ｏｒａｓｑｕａｒｅ．Ｔｈｅｍｉｎｉｍｕｍｉｎｔｅｒｆａｃｅｐａｒｔｉｔｉｏｎｔｒａｖｅｒｓｅｄａｌｌｐｏｓｓｉｂｌｅｐａｒｔｉｔｉｏｎｓａｎｄｐｅｒｆｏｒｍｅｄｔｈｅａｃｔｕａｌｓｅｇｍｅｎｔａｔｉｏｎｉｎ ｔｈｅｓｈｏｒｔｅｓｔｗａｙｏｆｔｈｅｓｕｍｏｆｔｈｅｓｕｒｆａｃｅａｒｅａｓｏｒｐｅｒｉｍｅｔｅｒｓｏｆｅａｃｈｓｕｂｇｒａｐｈ．Ｓｔｅｐｗｉｓｅｏｎｅｗａｙｓｅｇｍｅｎｔａｔｉｏｎｗａｓ ｂａｓｅｄｏｎｔｈｅｐｒｉｍｅｆａｃｔｏｒｄｅｃｏｍｐｏｓｉｔｉｏｎｏｆｔｈｅｓｅｇｍｅｎｔａｔｉｏｎｎｕｍｂｅｒ，ａｎｄｔｈｅｓｅｇｍｅｎｔａｔｉｏｎｗａｓｐｅｒｆｏｒｍｅｄａｌｏｎｇｔｈｅ ｄｉｒｅｃｔｉｏｎｗｉｔｈｔｈｅｌａｒｇｅｓｔｎｕｍｂｅｒｏｆｄｉｆｆｅｒｅｎｔｃｏｏｒｄｉｎａｔｅｓｉｎｏｒｄｅｒｏｆｔｈｅｐｒｉｍｅｆａｃｔｏｒｓｆｒｏｍｔｈｅｌａｒｇｅｓｔｔｏｔｈｅｓｍａｌｌｅｓｔ， ｕｎｔｉｌａｌｌｐｒｉｍｅｆａｃｔｏｒｓｈａｄｂｅｅｎｔｒａｖｅｒｓｅｄ．Ａｆｔｅｒｗａｒｄｓ，ｆｏｒｔｈｅｓｐａｒｓｅｌｉｎｅａｒｅｑｕａｔｉｏｎｓｄｉｓｃｒｅｔｉｚｅｄｆｒｏｍｔｈｅｐａｒｔｉａｌ ｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｏｆｔｈｅｍｏｄｅｌ，ｔｈｒｏｕｇｈＶｔｙｐｅ，Ｗｔｙｐｅ，ａｎｄＫｔｙｐｅｃｙｃｌｅｓ，ｆｒｏｍｔｈｅ