- 1、本文档共44页,可阅读全部内容。
- 2、本文档付费后,不意味着付费购买了版权,只能用于单位或个人使用,不得用于商业用途(如:【转卖】进行直接盈利和【编辑后售卖】进行间接盈利)。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:18428362892(电话支持时间:9:00-19:00)。
- 稀疏矩阵直接解法(补充内容1) - 喻文健
- 稀疏矩阵直接解法(补充内容1) -
喻文健
from T.A. Davis’s book: Direct Methods for Sparse Linear
Systems , and Y. Saad and J.R. Gilbert’s course slides
\ --“It’s amazing how much numerical linear
algebra is contained in that single character”
大规模稀疏矩阵概述
矩阵填入元与图
稀疏下三角阵的Lsolve算法
稀疏矩阵LU分解与矩阵重排序
稀疏对称正定阵的Cholesky分解
3
稀疏矩阵的应用领域
◦ 结构工程(structural engineering),
计算流体力学(computational fluid dynamics),
油藏数值仿真(reservoir simulation),
电力网络(electrical network),
优化问题(optimation),
Google PageRank,
信息抽取(information retrieval),
有限差分
电路仿真(circuit simulation), 有限元
器件仿真(device simulation) 边界元
4
ARC130 SHERMAN5
非结构化的
PORES3 BP_1000 5
集成电路供电网仿真
◦ DC分析:求解电阻网络
电路中的节点电压
集成电路供电线网的一部分
◦ 一个来自IBM公司的例子
n= 474,524
nnz= 2,020,882
在Matlab中, 占26M内存 ?
“\”求解(矩阵时对称正定矩阵),
t= 15s
“chol”: Cholesky分解系数矩阵
内存不够用! ?
6
稀疏矩阵的运算
矩阵与向量的乘法
7
计算y=Ax (CSR存储格式)
// 第i行的非零元素值
// 第i行非零元的列号
ki[k] x y
文档评论(0)