清华大学 大规模科学计算与并行计算9.pdf

第九讲 ：矩阵计算并行算法 陈俊清 清华大学数学科学系 科学计算 引言 密矩阵向量乘并行算法 密矩阵乘并行算法 科学计算 引言 密矩阵向量乘并行算法 密矩阵乘并行算法 科学计算 前言：微分方程 给定一个常微分方程组 ，有两类问题 ： 求解未知量在每个时间点的值 显式格式 隐式格式 求振动频率 特征值 问题求解 科学计算 显式格式 假定常微分方程为 是一个稀疏矩阵 给定初始条件 计算 显式格式比如 向前法 方程近似为： ，稀疏矩阵向量乘积 。 显式格式的优缺点 算法简单 ，只需要计算稀疏矩阵向量乘积 稳定性 问题 ：可能需要步长很小 ，特别地当系统具有刚性时 （具有 特别大的元 ，可能变化很快） 科学计算 隐式格式 假定常微分方程为 是一稀疏矩阵 给定初始条件 计算 隐式格式比如 向后法 方程近似为： ，求解稀疏线性方程组 。 隐式格式的优缺点 可 以用大的步长 ，特别是刚性 问题 算法困难 ：每个时间步都要求解一个稀疏线性方程组 科学计算 特征值 问题 计算振动的模式：求特征值和特征 向量 求方程 具有结构 的解 ，其中 为一个常向量 。 因此 ， 为 的特征值 ， 为对应的特征 向量 。 求解过程涉及稀疏矩阵乘积 ，或者求解稀疏线性方程组 。