带有双重循环展开的非对称稀疏线性方程组快速直接解法

来源 :第三届全国计算爆炸力学会议 | 被引量 : 0次 | 上传用户:yughg
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文介绍了一种新的直接解法来求解科学与工程计算中生成的大型非对称稀疏线性方程组.该解法从现有的对称解法中演变出来,其分解过程在矩阵的上、下三角阵中对称行进.该解法中的LDU分解算法利用了双重循环展开技术,并且由于其对称行进的求解方式,可以通过修改已有的对称矩阵分解算法的代码来实现,这提供了从对称解法到非对称解法的快捷转换.在数值测试中求解了若干无网格局部Petrov-Galerkin法(MLPG)法生成的矩阵,结果表明本文的方法可以大幅度提高了大型非对称稀疏线性方程组的求解速度.
其他文献
自从60年代开创计算流体力学以来,计算流体力学经历了形成、发展和成熟的几个阶段,计算流体力学已经成为流体力学三个分支:理论流体力学、实验流体力学和计算流体力学中发展
由于凝聚态炸药爆轰机制的极度复杂性导致对爆轰过程的理解至今仍然非常有限。在分子动力学模拟微观研究达到工程应用的水平之前,唯象模型仍然占据主要地位。本文以工程数值模
数值方法是解决数值问题的近似方法,它对于科学的发展是相当重要的,数值软件是对数值方法的一种计算机实现过程.我国在若干领域数值方法研究上处于国际领先地位,但我国的数值
山东德棉成功开发海洋纤维纱线最近,山东德棉股份有限公司与青岛大学阻燃纤维研究所共同研制开发精梳棉/海洋纤维(60/40)混纺40S机织纱和海洋纤维纯纺16S针织纱获得成功。海
期刊
期刊
本文介绍应用于冲击动力学问题中主要的无网格方法.着重阐述这些方法的算法特点及其在冲击动力学中的具体应用现状和最新进展.从形函数的构造方法和微分方程的离散形式两条线
本文对多物质弹塑性流体动力学问题的Euler型有限差分方法进行了研究.针对自行开发的二维及三维爆炸与冲击问题数值仿真程序以及与之相配套的前后处理程序,阐述了其理论基础
维生素K1注射液的严重不良反应/事件的表现如下.全身性损害占49.2%,主要表现为过敏性休克(328例)、过敏样反应、发热、寒战等;呼吸系统损害占26.5%,主要表现为呼吸困难、胸闷、
期刊
本文提出了一种基于PML算法的多尺度模拟方法.在本文的计算工作中,将被研究的区域分成三个部分:MD(MolecularDynamics)区域、MC(MacrosCopic)区域、以及PML区域.粗网格方程用
磁学虽是一门古老的学科,但在当今信息社会中其重要性越来越明显.这是因为信息的记录、存贮、调用和处理都离不开磁学与磁性材料.磁性材料已成为信息时代的主要物质支柱之一.