优雅代码随硬件发展的演化-北京大学数学科学学院.PDFVIP

How Elegant Code Evolves with Hardware: The Case of Gaussian Elimination 报告人：樊玉伟 Author ： Jack Dongarra Piotr Luszczek 作者简介 Jack Dongarra ： • 田纳西大学计算机科学系的计算机科学杰出教授 • 美国橡树岭国家实验室（ORNL ）杰出研究人员 • 曼彻斯特大学计算机科学与数学学院的Turing Fellow • 美国莱斯大学计算机科学系的副教授 • AAAS，ACM和IEEE的成员 • 美国工程院的院士 • 研究领域：线性代数中的数值算法，并行计算，高级计算 机架构的应用，程序设计方法学，用于并行计算机的工具 • 在ISPACK, LINPACK, the BLAS, LAPACK, ScaLAPACK, Netlib, PVM, MPI, NetSolve, Top500, ATLAS, 和 PAPI上 做出过贡献 Piotr Luszczek： • 田纳西大学诺克斯维尔分校 研究教授 • 研究领域：并行的核外（out-of-core ）库 大型超级计算机安装的标准化评价 • 开发了一个自适应的软件库，能够自动选择最 优的算法来有效地利用现有硬件以及有选择地 处理输入数据 • 对高性能编程语言的设计和实现感兴趣 • 数学上等价的公式在计算上不等价 • “计算上等价”的算法在实现上不等价 真正的适应性Ⅰ • 应用程序层面上 数学模型 算法 软件 • 软件层面上 性能 可移植性 代码可理解性 真正的适应性Ⅱ 真正的适应性应试算法具有很好的可移植性 可迁移性 + 可移植性 = 模块化 + 高效的软件标准 预备知识 • 稀疏矩阵：绝大多数元素为0的矩阵 • 稠密矩阵 • Gauss消元法 • 计算机体系结构 Gauss消元法算法 for k = 1:n-1 find p (k = p = n) s.t. |A(p,k)| = max{|A(i,k)|, i = k:n } A(k, 1:n) -- A(p, 1:n) u(k) = p if A(k,k) =/= 0 A(k+1:n,k) = A(k+1:n,k)/A(k,k) A(k+1:n,k+1:n) = A(k+1:n,k+1:n) - A(k+1:n,k) * A(k,k+1:n) else STOP!!! The matrix is singular! end end 程序实现与运行结果 • MATLAB 内部实现 • MATLAB代码实现 • C++实现 n MATLAB内部 MATLAB代码 C++代码