稀疏矩阵压缩 - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Computation - 科学计算

稀疏矩阵压缩

稀疏矩阵压缩# Computation - 科学计算

t*t2002-08-16 07:08

1 楼

这里牛牛多,问个问题.
我有一个阶数很大的BLOCK TRIDIAGONAL矩阵(暂称为D)
每个BLOCK又是TRIDIAGONAL的
按照LAPACK所需的格式我知道怎么压D
问题是等我压完了,猛一抬头看我算法里
实际需要的是 D*D
而这个好象是BLOCK PENTA-DIAGONAL的矩阵
每个BLOCK又是PENTA-DIAGONAL的
哪位有什么好的resource来压这个么?
多谢.

t*t2002-08-16 07:08

2 楼

Problem solved!
I found this package at:
http://www.cs.umn.edu/Research/arpa/SPARSKIT/sparskit.html
SPARSKIT: a basic tool-kit for sparse matrix computations (Version 2)
It is a very nice packge. Using the GMRES from it, I reduced
my computation work to ~1/70 as that of original using dense & direct solver.

【在 t******t 的大作中提到】

: 这里牛牛多,问个问题.
: 我有一个阶数很大的BLOCK TRIDIAGONAL矩阵(暂称为D)
: 每个BLOCK又是TRIDIAGONAL的
: 按照LAPACK所需的格式我知道怎么压D
: 问题是等我压完了,猛一抬头看我算法里
: 实际需要的是 D*D
: 而这个好象是BLOCK PENTA-DIAGONAL的矩阵
: 每个BLOCK又是PENTA-DIAGONAL的
: 哪位有什么好的resource来压这个么?
: 多谢.