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