外推多网格法——大规模科学计算的新算法

来源 :湖南师范大学 | 被引量 : 0次 | 上传用户:t123
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
快速求解由差分法或有限元法导出的大型线性方程组是大规模科学和工程计算中的重要问题。人们的目标是对N阶方程组用O(N)次乘除运算得到所需精度的解。而多网格法(MG)第一次实现了这个目标,成为求解大规模问题最有效的方法。经典MG的理论方法已经比较成熟,由于使用了3种运算:插值,限制和迭代,求解时需要在多重粗细网格上反复迭代,所以程序较为复杂。1996年德国Bornemann等人提出瀑布式多网格法(CMG),即从粗网到细网的单向计算,只采用了插值与迭代两种运算,程序容易实现,令人注目。   本文研究一类新的瀑布式多重网格法,主要工作和创新点如下。   首次提出了外推瀑布式多网格法(EXCMG)。本论文基于有限元的渐近展开式,导出了新的外推公式.新算法沿用CMG的思想,但将粗网上的线性插值改为新外推与二次插值,为密网提供更好的初值,本质上减小了初始误差,对加速收敛起着关键的作用。新算法对函数和导数都收敛,且都有高精度。本论文在PC机上用EXCMG求解了400万未知数的问题,对函数和导数都达到8-10位精度,仅费时10分钟,进一步证实了这些优势。   首次证明了共轭梯度法(CG)按l2模的有界性和收敛性.对CG法至今只有按能量模的收敛估计,因此这是一个新结果。本文还发现当第i层网格的迭代次数取幂指数形式(mLβL-i)时,基于CG迭代的多重网格法,存在着一个网格层数的“阀值”i0,即当i<i0时,第i层上的误差快速衰减,CG迭代压缩效果非常明显;而当i>i0时,CG主要起磨光作用.此特性对证明EXCMG按离散L2模的收敛性起重要作用。   对正则性较弱的解u∈H3,首次证明了双线性元的渐近展开在节点集上按离散L2模有高精度O(h3),因此EXCMG仍有效.已有的渐近展开式研究主要是对光滑解(如u∈C3),并只得到了在逐点意义下的结果。上述结果拓宽了EXCMG的应用范围。数值试验表明,EXCMG对非光滑解u∈H3甚至H2仍有很高精度。
其他文献
非线性偏微分方程的求解通常是非常困难的.已有某些方法求孤子方程的精确解:例如反散射方法,双线性方法,达布变换方法以及代数几何方法等.在这些方法中,达布变换方法是一种强
学位
教育信息化已经成为教学主流趋势之一,在此背景下不断有新的教育理念、模式应运而生,改变了人们学习和交流的方式,为教育提供了新的可能.翻转课堂是一种新型的教学模式,翻转
随着《高中数学课程标准(实验)》(以下简称《标准》)的颁布实施,国家已经审定通过了五个版本的基于《标准》的高中数学教科书,并于2004年9月开始在各试验区试用。本研究是以人
在新的教学理念下,小组合作学习作为一种新的教学模式被广泛地运用,而要如何才能让小组合作学习模式成为小学英语教学的得力助手呢?本文结合笔者在教学中的实操过程,阐述在英语教
支持向量机是建立在统计学习理论基础之上的一种机器学习方法,它能有效的处理小样本,非线性等学习问题。经典的支持向量机是建立在概率空间上基于实随机样本的,因此,它难以处
物理学、力学和工程技术中的许多问题,都可以归结为求解偏微分方程的初边值问题。但除了极少数的问题可以给出解析解,绝大多数的问题不得不求助于数值方法寻求其近似解。目前,有
所谓微分算子主要研究两个方面的问题,一方面研究微分算子的谱问题,另一方面研究微分算子的逆谱问题。所谓逆谱问题就是由谱数据的信息,尤其是特征值,确定微分算子进而将其重构。
随着国内外星球探测工程的实施和相关技术的快速发展,尤其是我国嫦娥一、二号摄制获得了海量的月球数据需要分析和处理。针对复杂光照和地理环境下月貌撞击坑图像呈现多模态,
图像的局部特征在保留图像重要信息的同时,又有效地减少了图像处理的数据量,极大地提高了运算速度。因此特征提取成为模式识别与计算机视觉等图像处理相关领域的基础性研究内
我国的教育改革进程随着教育体制的改变在不断的深化,而思想品德作为主要的素质教育科目,越来越受到教育领域的重视,如何合理建构小学阶段学生的德育教学,是所有一线教师要在实践