e******o
发帖数: 27 | 1 【 以下文字转载自 Computation 讨论区 】
发信人: espresso (espresso), 信区: Computation
标 题: 请问:稀疏矩阵运算有什么好的C++库？
发信站: BBS 未名空间站 (Thu Nov 10 23:45:06 2005), 转信
大概有6000*6000
想求矩阵相乘，SVD之类,有什么好用的库么?谢谢 |
i******e
发帖数: 171 | 2 有时间要求么？
【在 e******o 的大作中提到】
: 【 以下文字转载自 Computation 讨论区 】
: 发信人: espresso (espresso), 信区: Computation
: 标 题: 请问:稀疏矩阵运算有什么好的C++库？
: 发信站: BBS 未名空间站 (Thu Nov 10 23:45:06 2005), 转信
: 大概有6000*6000
: 想求矩阵相乘，SVD之类,有什么好用的库么?谢谢
|
e******o
发帖数: 27 | 3 有多慢?时间应该不要紧
【在 i******e 的大作中提到】
: 有时间要求么？
|
i******e
发帖数: 171 | 4 那就用gsl, gnu scientific library． 几行code就解决了． 但6000的矩阵可能要花
半个小时把所有的eigenvalue, eigenvector算出来．
【在 e******o 的大作中提到】
: 有多慢?时间应该不要紧
|
e******o
发帖数: 27 | 5 gsl没有稀疏矩阵巴，估计gsl用的是malloc,内存分配受得了么?
还是double类型:(
花
【在 i******e 的大作中提到】
: 那就用gsl, gnu scientific library． 几行code就解决了． 但6000的矩阵可能要花
: 半个小时把所有的eigenvalue, eigenvector算出来．
|
i******e
发帖数: 171 | 6 有SVD还有其他分解方法(LU,.....)．但我好像没找到矩阵相乘的函数， 自己写了两个
循环．
花
【在 i******e 的大作中提到】
: 那就用gsl, gnu scientific library． 几行code就解决了． 但6000的矩阵可能要花
: 半个小时把所有的eigenvalue, eigenvector算出来．
|
i******e
发帖数: 171 | 7 内存不是问题．
【在 e******o 的大作中提到】
: gsl没有稀疏矩阵巴，估计gsl用的是malloc,内存分配受得了么?
: 还是double类型:(
:
: 花
|
e******o
发帖数: 27 | 8 恩，不错，试了一个20000*20000没问题。多谢
【在 i******e 的大作中提到】
: 内存不是问题．
|
i******e
发帖数: 171 | 9 昨晚试了试，发现4000X4000的矩阵算eigenvalue,eigenvector花了５个小时(P4 3.0G).
【在 e******o 的大作中提到】
: 恩，不错，试了一个20000*20000没问题。多谢
|
o**o
发帖数: 3964 | 10 Lapack没有稀疏矩阵的么? 4kx4k的double complex矩阵求逆大概在半个小时到一个小时
.
【在 i******e 的大作中提到】
: 昨晚试了试，发现4000X4000的矩阵算eigenvalue,eigenvector花了５个小时(P4 3.0G).
|
|
|
i******e
发帖数: 171 | 11 嗯, 如果本身用c,c++写得话，用GSL就方便．如果时间重要的话，还是要用LAPACK.
时
【在 o**o 的大作中提到】
: Lapack没有稀疏矩阵的么? 4kx4k的double complex矩阵求逆大概在半个小时到一个小时
:
: .
|
o**o
发帖数: 3964 | 12 为什么一定要C++呢，如果不需要改内部的函数。
这个东西做成库就是让别的语言方便调用的
【在 i******e 的大作中提到】
: 嗯, 如果本身用c,c++写得话，用GSL就方便．如果时间重要的话，还是要用LAPACK.
:
: 时
|
i******e
发帖数: 171 | 13 因为题目说了要c++库。 GSL虽然是c,但是基本一样。
【在 o**o 的大作中提到】
: 为什么一定要C++呢，如果不需要改内部的函数。
: 这个东西做成库就是让别的语言方便调用的
|
x***g
发帖数: 52 | 14 sparselib++
【在 e******o 的大作中提到】
: 恩，不错，试了一个20000*20000没问题。多谢
|
N***m
发帖数: 4460 | 15 什么样的矩阵？hermitian的应该很快。
0G).
【在 i******e 的大作中提到】
: 昨晚试了试，发现4000X4000的矩阵算eigenvalue,eigenvector花了５个小时(P4 3.0G).
|
