一种基于HYPRE的高次Lagrange有限元方程的并行代数多层网格法

来源 :湘潭大学 | 被引量 : 0次 | 上传用户:baalzio
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
多层网格法通常可分为几何多层网格(GMG)法和代数多层网格(AMG)法与GMG法相比,AMG法具有更强的普适性和鲁棒性(robustness),它是求解许多大规模科学工程计算问题特别是偏微分方程离散化系统的最为有效的方法之一。HYPRE是目前国际上流行的一种在大规模并行计算机上求解大型稀疏线性方程组的数值软件包其目的是为用户提供高性能并行解法器和预条件子,BoomerAMG是其中最常用的并行代数多层网格解法器之一.本文基于HYPRE平台,针对三维二阶椭圆边值问题在分层基下的高次Lagrange有限元方程讨论其并行AMG求解算法我们的工作主要有1.简要介绍了HYPRE数值软件包重点介绍了几种经典AMG的网格粗化算法,如RS算法,CLJP算法以及一种常用的并行网格粗化算法:Falgout算法,最后介绍了近年发展起来的连续子空间校正算法框架及其收敛性理论2.在四面体网格剖分下,针对一种分层基下的三维二阶椭圆边值问题的高次有限元方程分别设计了生成其总刚度矩阵和总载荷向量的串行算法和一种基于子区域划分的并行算法。对于后者,我们是通过引入与侧面、分划棱及角点的关联矩阵,并进行通信,来获得当前进程上与其它进程关联的刚度矩阵和载荷向量。另外,我们对分层基采用了一种合理的序,不仅为程序设计带来方便而且提高了并行AMG的磨光效率数值实验结果表明,我们所设计的并行算法扩大了刚度矩阵的生成规模,具有较好的可扩展性3.首先针对上述并行生成的高次有限元方程设计并分析了一种基于辅助变分问题的新的AMG(X-AMG)法在拟一致网格剖分下,利用连续子空间校正理论证明了这种新的AMG法的下降率与网格尺寸无关,数值结果也验证了理论的正确性接着为X-AMG法设计了两种并行算法,第一种基于串行刚度矩阵结构,该算法虽然关于迭代次数具有着较好可扩展性但由于其中的并行向量与串行向量的相互转化过于频繁,并且算法过于依赖磨光子的选取从而降低了并行效率因此基于子区域划分的并行矩阵结构,我们给出了一种更加合理的X-AMG并行算法数值实验结果表明,对于求解高次有限元方程这种并行X-AMG算法比BoomerAMG具有更高的效率
其他文献
期刊
本文将Hartmann函数推广到异面情形实现了两个异面直线的光滑拼接,还用两个锥面交线成功光滑拼接了两个异面直线.并以这两种曲线为轴线构造了光滑拼接两个半径相同的轴异面圆
近日,东莞检验检疫局在对一批进境废纸检查时发现,这批进境废纸中竟然有近十种有害生物藏身于中。据悉,此次携带多种危险性生物的废纸是从美国入境,检疫人员在进行现场掏柜检
半无穷区间上二阶边值问题起源于对非线性椭圆微分方程对称径向解以及半直线上中间多漏洞的煤气压力模型的研究[1].现在人们越来越关注半无穷区间边值问题正解的存在性,并取得
公平交换数字签名问题是密码学的一个基本问题,在电子商务中有着广泛的应用,同时签名是解决公平交换的一个新的工具。在一个同时签名方案中,两个实体共同生成两个没有绑定的
Mizar系统是用于证明或计算数学问题的计算机语言系统。它由波兰华沙大学AndrzejTrybulec教授组织的Mizar协会领导,其逻辑框架是基于Jaskowski自然演绎推理的古典逻辑。Mizar
本文基于可加的ψ-模糊偏好结构,我们分别在不可比关系为空集和大偏好关系强完全的条件下,研究了大偏好关系、严格偏好关系、无区别关系的传递性及其传递性指标.  首先,我们简
期刊
本文研究t-模的旋转不变性.首先,我们给出了旋转不变t-模的基本性质,我们讨论旋转不变性与其它t-模的性质之间的关系,主要包括旋转不变性与t-模的左连续性、幂零性、矛盾律、IT
Artin代数表示论的主要目的就是用一个代数的模范畴的性质来刻画这个代数.用模论来研究代数的好处之一是我们可以应用范畴理论和同调代数.从上个世纪40年代开始,同调代数逐渐