近似三对角Toeplitz方程组的并行化实现

来源 :吉林农业·C版 | 被引量 : 0次 | 上传用户:zyr2007
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  摘要:近似三对角Toeplitz矩阵在时间序列分析、线性逼近、图论的特征值分析等多个领域经常出现。本文对近似三角Toeplitz 矩阵求解算法的串、并行算法进行了分析,并在流处理器平台上实现了算法的并行化,实验结果说明并行算法在流处理器平台上获得较好的加速。
  关键词:Toeplitz;方程组;并行化;流处理器
  中图分类号:TP301 文献标识码:A 文章编号:1674-0432(2013)-06-0280-2
  0 引言
  線性方程组的求解是许多科学和工程计算的基本问题之一。在核物理、流体力学、石油地震、数值天气预报等大规模数据处理中都会遇到三对角系统的求解问题,其中很重要的一类是以近似三角Toeplitz矩阵为系数矩阵的线性方程组的求解。由于在理论和实际应用上的重要性,近20年来三对角方程组的并行算法研究十分活跃,对于这种三对角矩阵算法的研究也取得了一系列的成果,本文在前人研究的基础上[1],在流处理器平台上实现了算法的并行化。
  1 近似三对角Toeplitz方程组串行求解的一般思想
  我们选用的实验平台式是Imagine流处理器的时钟精确模拟器,并用流编程模型来编程。Imagine处理器[4][5][6]有8个运算簇,可以开发流程序的数据级并行性。软件管理的片上存储系统可以很好的捕获程序的数据局部性。我们使用的流编程模型包括StreamC/KernelC两级编程语言。流编程语言的主要思想是把相关数据的集合组成一个记录,再把有序的记录构成流。应用程序则被程序员分解成一些对大量流数据进行操作的计算函数(kernel),对输入流中数据进行计算,并把计算结果作为下一个计算函数的输入。
  3 并行结果
  在流处理器平台上,我们获得了图1中的数据,可以看出随着输入集规模的不断增加,流处理器也能获得不错的加速比。一方面流处理器中的大量的计算单元可以完成众多的计算操作,另一方面流处理器中软件管理的片上存储系统很好的开发了计算函数间的数据局部性,这两点使得并行算法在流处理器获得了不错的加速效果。
  4 结论
  本文首先讨论了求解近似三对角Toeplitz方程组的串行算法,然后在其分布式算法的基础上,描述了在流处理器上并行化算法的步骤,最后给出了并行化结果,说明了并行算法在流处理器平台上获得了较高的并行效率。
  参考文献
  [1] 单润红,高峰.近似三对角Toeplitz 方程组的快速分布式并行算法.
  [2] 面向结构的并行算法——设计与分析.
  [3] 赵自春,李晓梅1 Toeplitz 循环三对角方程组的并行解法1 国防系统分析与软件,1989,(2):50-55.
  [4] Scott Rixner1,William J. Dally. A Bandwidth-Efficient Architecture for Media Processing. Micro-31 (1998).
  [5] Ujval J. Kapasi, William J. Dally, Scott Rixner. The Imagine Stream Processor.
  [6] Jung Ho Ahn, William J. Dally, Brucek Khailany. Evaluating the Imagine Stream Architecture.
  作者简介:张传龙(1982-),男,安徽淮南人,安徽理工大学计算机学院2010级工程在读硕士研究生,研究方向:计算机科学与技术。
其他文献
【正】我国玉米收获机械化水平相对较低,是玉米生产,甚至整个粮食生产机械化的薄弱环节。我国玉米机械收获蕴藏着巨大的农机制造业和服务业市场,但是,当前市场上的玉米收获机
对中信实业银行向零售银行业务转型的猜测近半年来一直持续不断。近日,流传在坊间的传言终于得到证实,中信实业银行将用最直观的方式宣布自己的变身:去掉“实业”两字,改名为“中
安全生产是农电部门各项工作的基础,是电力企业一个永恒不变的主题。文章主要通过分析安全生产的重要性,提出了农电部门安全生产的重点措施。
2004年10月,刚回到四川长虹的赵勇以一场“虹色十月”揭开了自己大刀阔斧整顿旧长虹的序幕。随后在产业布局、公司组织以及高层人事方面,赵勇均展示出了自己的快刀风范。一年时
学生经常产生一些似是而非的错误,如: 例1 求函数y=x+(x<sup>2</sup>-3x+2)<sup>1/2</sup>的值域。 错解 由y-x=x<sup>2</sup>-3x+2)<sup>1/2</sup> 可得 (y-x)<sup>2</sup>=x
《中学数学月刊》文[1]研究了2n位(不妨称为第Ⅰ类)双色完全平方数的构造,本文给出2n+1位(第Ⅱ类)双色完全平方数的定义及其构造。 定义 若2n+l位自然数a<sub>2n+1</sub>a<sub>2
《中学数学月刊》1997年第5期《一类新三角不等式》一文绘出了下面几个三角不等式.
高中课本《平面解析几何》中的定比分点坐标公式,深刻揭示了在直角坐标系中直线上相关点之间的内在联系,是平面解析几何理论的重要基石之一,它不仅为我们更深层次地研究平几
继全球范围内率先发布10款基于AMD Opteron双核芯处理器的服务器新品后,曙光近日又宣布正式启动“双核升级计划”,开始布局中国双核服务器计算。此举不仅标志着曙光正式掀开了