【摘 要】
:
随着超级计算机的迅猛发展,高性能计算技术面临着巨大的挑战。由于绝大部分计算机均使用浮点运算标准,在大规模问题中,减弱舍入误差对数值计算产生的影响显得尤为重要。本文以向前向后误差分析理论、动态误差分析理论、无误差变换技术和双倍双精度基本算术运算为基础,提出了几种基于多部分格式的补偿算法,包括商-差算法、Clenshaw-Smith算法、Barrio-Clenshaw-Smith算法、Horner算法
论文部分内容阅读
随着超级计算机的迅猛发展,高性能计算技术面临着巨大的挑战。由于绝大部分计算机均使用浮点运算标准,在大规模问题中,减弱舍入误差对数值计算产生的影响显得尤为重要。本文以向前向后误差分析理论、动态误差分析理论、无误差变换技术和双倍双精度基本算术运算为基础,提出了几种基于多部分格式的补偿算法,包括商-差算法、Clenshaw-Smith算法、Barrio-Clenshaw-Smith算法、Horner算法与Volk-Schumaker算法,并给出了它们的误差分析、数值实验与应用。本文还开发了C语言、Matlab与Scilab平台上的几种基于多部分格式的高精度算法库。本文的主要工作可分为如下三个部分:(1)第一部分(第2章)对商-差算法进行了稳定性分析,并设计了补偿的商-差算法。商-差算法可以应用到亚纯函数求极点、多项式求根、páde逼近和特征值问题等,但是它非常不稳定。本文提出了商-差算法的条件数,利用无误差变换与补偿除法运算设计了补偿商-差算法,并通过向前向后误差分析理论给出了商-差算法及其补偿算法的稳定性分析。数值实验给出了补偿商-差算法及其向前误差界的有效性验证,同时还通过三个应用测试了补偿算法的精确性。(2)第二部分(第3、4章)着重考察了计算多项式序列的问题,包括一维经典正交多项式序列和二维张量积曲面。计算经典正交多项式序列及其k阶导数值需要使用Clenshaw-Smith算法和Barrio-Clenshaw-Smith算法,计算二维幂指数基多项式序列需要双变量Horner算法,计算Bézier张量积曲面需要二维的Volk-Schumaker算法。本文利用无误差变换技术设计了这些算法的补偿算法,并通过向前向后误差分析理论与动态误差分析理论给出了它们的稳定性分析或动态误差界。同时,数值实验给出了补偿算法、向前误差界与动态误差界的有效性验证。另外,本文还将计算经典正交多项式序列及其k阶导数值的各个算法、补偿算法与计算动态误差界的算法整合在了C语言与Matlab平台上,形成了高精度算法库,并通过实际应用中可能存在的问题对算法库进行了精确性、有效性和速度方面的数值测试。(3)第三部分(第5章)设计了混合扩展精度的软件包。本文将双倍双精度、三倍双精度与四倍双精度均定义为扩展精度,结合实数与复数运算,设计了混合扩展精度的算术运算算法。同时,本文在Matlab与Scilab平台上重新定义了扩展精度的类,并对各类算术运算进行了重载,通过初始化设计了两种显示模式与两种精确性模式。数值实验测试了该软件包的精确性与速度。
其他文献
中国空间站将于2022年前后建成,未来空间站的长期在轨运营、复杂科学实验的开展将主要依赖航天员来完成。载人飞行成本高、在轨驻留航天员人数有限,实现航天员工作效率提升,充分利用宝贵的在轨时间开展更多有价值的空间任务,对载人航天活动意义重大。为辅助空间站舱内航天员工作、提升航天员工作效率,本文提出了一款智能的舱内航天员跟随辅助机器人。该机器人能够智能地跟随所服务的航天员并与其进行交互,代替或辅助航天员
强流离子束在离子束驱动快点火、温稠密物质产生以及肿瘤治疗等领域有非常重要的研究价值。本论文采用数值模拟和理论分析的方法研究了超强激光与等离子体薄膜靶相互作用中强流离子束的产生及其在等离子体中的输运过程,主要研究内容如下:一、研究了激光辐射压加速中横向不稳定性的发展过程。当激光作用到调制靶表面时,横向不稳定性迅速激发,质子束密度出现周期性扰动。通过对质子平均面密度的傅里叶分析诊断了横向不稳定性的增长
点阵结构多功能设计是目前飞行器结构设计的重要研究方向,承载/阻尼一体化设计更是航空航天设备结构设计的难题。轻质点阵结构具有高韧性、抗冲击、吸声、电磁波吸收、有效隔热等优异性能,具有十分广泛的应用前景。在航空航天领域引起了许多研究者的关注。各国研究人员对点阵结构进行了广泛的研究。然而,目前的点阵结构优化算法与已有结构结合不足、异质三维点阵研究较少、异质点阵优化拓扑优化方法研究较少,严重制约了点阵结构
随着激光技术的不断发展,超强、超短的激光脉冲与原子、分子、凝聚态相互作用带来了很多新的物理现象。强场物理的研究为人们探测和调控物质的超快动力学过程提供了强有力的技术手段,具有重大的科学意义和应用价值。本文发展了经典轨迹蒙特卡洛方法和德布罗意-波姆力学分析方法,结合数值求解含时薛定谔方程,深入地研究了激光场作用下电子在库伦势和周期势中的超快动力学过程,主要内容由三部分组成:第一部分中,我们研究了太赫
保结构算法是微分方程数值算法的重要研究方向之一,其目的是构造数值积分保持连续系统的相应特征。一切真实的、耗散可忽略不计的物理过程都可以表示成Hamilton系统,它在自然界中有着非常广泛的应用。然而经典力学中研究的大部分系统都不是保守系统,所以很难将这类系统表示为经典的Hamilton力学形式以及最小作用量变分原理形式或者与此等效的Lagrange力学形式,极大地限制了保结构算法在耗散系统中的应用
原子分子在强激光场作用下,可以相干地发射处于极紫外到软X射线之间的高次谐波,也可以相干地发射处于毫米亚毫米波段的太赫兹波。同步探测高次谐波与太赫兹光谱(HATS)是研究原子分子中的电子结构以及强激光场下电子动力学的新型光学方法。由于太赫兹波与高次谐波在能量以及空间尺度上存在的五个量级的巨大差异,同步探测这两种辐射有助于加深对强场下电子运动过程的理解,实现辐射的原位调控。本文首先回顾了强场下的电子运
望远镜不仅是人类探索宇宙奥秘的重要科学工具,也是监测地球轨道、预警飞行器碰撞的守望者。但是大气湍流导致的波前畸变会导致大型地基望远镜的实际分辨率大幅下降,因此世界各大望远镜正竞相发展基于钠导星的自适应光学(AO)技术,以校正大气湍流导致的波前畸变。钠导星作为AO系统的信标源,其亮度是决定AO系统波前探测精度和响应速度的关键因素。尽管当前的钠导星亮度仍然限制着AO系统性能,科学家们已经开始研究下一代
基于格的密码体制作为后量子时代公钥密码体制的主要选择之一,受到了越来越多的关注。格中困难问题求解算法的有效性,在公钥密码体制的安全性分析中起着非常重要的作用。求解格中困难问题,需要设计相应的算法,而对格结构的深刻认识是设计出优异算法的前提。因此,研究格的结构、性质和格中求解困难问题的算法,具有十分重要的意义。本文的第一项工作是研究了计算格中最短向量的约化算法。对比了不同算法对同一个矩阵的不同度量参
燃气轮机具有功率密度高、运行平稳等突出优点,被广泛用做重大军事及民用装备的动力,其安全可靠运行至关重要。压气机、燃烧室、涡轮等气路部件是燃气轮机的主要故障源,对其实现有效的故障在线检测意义重大。与传统检测技术相比,基于静电传感的尾气异常颗粒监测技术具有较强的早期故障检测能力,得到了广泛关注。然而,该技术通常只依靠单个静电传感器,使得灵敏度分布不均匀、有用信息量不足,难以准确重构异常颗粒的数量及电荷
刚刚荣获辽宁省五一劳动奖章的宋秀锋平时很忙,采访从4月份开始联系,最终时间却确定在5月末一个周三的下午。他平均每天至少三台手术,每周二、周三上午还要出专家诊,再加上每天定时查房、会诊、急诊等,他的工作状态当真可以用"连轴转"三个字来形容。