论文部分内容阅读
偏微分方程的数值求解是科学与工程计算的重要任务,块三对角线性方程组是偏微分方程离散格式的主要形式之一。本文系统地研究了块三对角线性方程组直接求解方法的可扩展并行算法及其应用,主要工作如下:
1.将三对角线性方程组的RCD(Recursive Doubling)方法,推广到了块三对角线性方程组,得到BRCD(Block Recursive Doubling)方法。该方法将块LU分解的时序性较强的计算转化成矩阵的乘积计算,有利于设计并行算法;同时给出了块LU分解方法、BRCD方法和SBP(Sequential Block Prefix)方法的反向与双向计算公式,在计算中使用不同方向的计算公式可以减少计算量;给出了一个数值稳定的块三对角线性方程组快速BOER(Block Odd-Even Reduction)求解公式,提高了计算精度和速度;
2.对于Larfiba等提出的三对角线性方程组的重叠分割并行方法,一方面将它推广到了块三对角线性方程组,该方法对于强块对角占优块三对角线性方程组的求解有较高的并行效率;另一方面将它与非重叠分割并行方法的缩减方程组结合,提出了块三对角线性方程组的有通信重叠分割并行方法,对于弱块对角占优的块三对角线性方程组,该方法可以用少量通信代替过量的计算,保证了精度和计算效率;
3.对严格块对角占优块三对角线性方程组的两类缩减方程组,提出了近似求解方法和相应的并行算法,该方法的通信次数是可调整的,可以通过它的调整来控制精度;
4.对严格块对角占优的块三对角线性方程组提出了多个可扩展的重叠与非重叠分割并行算法(统称PBDDA算法),其中基于块LU分解的算法统称PBDDA-LU算法、基于BOER的算法统称PBDDA-OE算法。给出了两种通信模式,其特点是:仅使用局部通信;通信复杂度可调整;每次通信都是具有一定间隔的相邻处理器之间的单向数据传递。给出了相对误差与三个可调整的指标k<,0>(通信复杂度指标)、q+?(子问题规模指标)和p(处理器数量即算法并行度指标)的关系。PBDDA-LU算法是对PDD方法的推广,在推广中使用了反向块LU分解,改进了回代方法,比直接使用PDD方法求解块三对角线性方程组减少了将近30%的计算量;
5.提出了一类Topelitz块三对角线性方程组的PBDDA-OE算法公式;使用这个公式进行高维Poisson方程第一边值问题的数值求解;研究了PBDDA-OE算法的性能,提出了保证精度和并行效率的分治策略;
6.利用以上Topelitz块三对角线性方程组的PBDDA-OE算法公式进行高维抛物型方程初始边值问题的数值求解,研究了PBDDA-OE算法的性能,提出了保证精度和并行效率的分治策略;
7.提出了另一类Topelitz块三对角线性方程组的PBDDA-OE算法公式;使用这个公式进行高维双曲型方程初始边值问题的数值求解;给出了控制精度与并行效率的分治策略。