Redian新闻
>
请问:稀疏矩阵运算有什么好的C++库?
avatar
请问:稀疏矩阵运算有什么好的C++库?# Programming - 葵花宝典
e*o
1
【 以下文字转载自 Computation 讨论区 】
发信人: espresso (espresso), 信区: Computation
标 题: 请问:稀疏矩阵运算有什么好的C++库?
发信站: BBS 未名空间站 (Thu Nov 10 23:45:06 2005), 转信
大概有6000*6000
想求矩阵相乘,SVD之类,有什么好用的库么?谢谢
avatar
i*e
2
有时间要求么?

【在 e******o 的大作中提到】
: 【 以下文字转载自 Computation 讨论区 】
: 发信人: espresso (espresso), 信区: Computation
: 标 题: 请问:稀疏矩阵运算有什么好的C++库?
: 发信站: BBS 未名空间站 (Thu Nov 10 23:45:06 2005), 转信
: 大概有6000*6000
: 想求矩阵相乘,SVD之类,有什么好用的库么?谢谢

avatar
e*o
3
有多慢?时间应该不要紧

【在 i******e 的大作中提到】
: 有时间要求么?
avatar
i*e
4
那就用gsl, gnu scientific library. 几行code就解决了. 但6000的矩阵可能要花
半个小时把所有的eigenvalue, eigenvector算出来.

【在 e******o 的大作中提到】
: 有多慢?时间应该不要紧
avatar
e*o
5
gsl没有稀疏矩阵巴,估计gsl用的是malloc,内存分配受得了么?
还是double类型:(



【在 i******e 的大作中提到】
: 那就用gsl, gnu scientific library. 几行code就解决了. 但6000的矩阵可能要花
: 半个小时把所有的eigenvalue, eigenvector算出来.

avatar
i*e
6
有SVD还有其他分解方法(LU,.....).但我好像没找到矩阵相乘的函数, 自己写了两个
循环.



【在 i******e 的大作中提到】
: 那就用gsl, gnu scientific library. 几行code就解决了. 但6000的矩阵可能要花
: 半个小时把所有的eigenvalue, eigenvector算出来.

avatar
i*e
7
内存不是问题.

【在 e******o 的大作中提到】
: gsl没有稀疏矩阵巴,估计gsl用的是malloc,内存分配受得了么?
: 还是double类型:(
:
: 花

avatar
e*o
8
恩,不错,试了一个20000*20000没问题。多谢

【在 i******e 的大作中提到】
: 内存不是问题.
avatar
i*e
9
昨晚试了试,发现4000X4000的矩阵算eigenvalue,eigenvector花了5个小时(P4 3.0G).

【在 e******o 的大作中提到】
: 恩,不错,试了一个20000*20000没问题。多谢
avatar
o*o
10
Lapack没有稀疏矩阵的么? 4kx4k的double complex矩阵求逆大概在半个小时到一个小时

.

【在 i******e 的大作中提到】
: 昨晚试了试,发现4000X4000的矩阵算eigenvalue,eigenvector花了5个小时(P4 3.0G).
avatar
i*e
11
嗯, 如果本身用c,c++写得话,用GSL就方便.如果时间重要的话,还是要用LAPACK.



【在 o**o 的大作中提到】
: Lapack没有稀疏矩阵的么? 4kx4k的double complex矩阵求逆大概在半个小时到一个小时
:
: .

avatar
o*o
12
为什么一定要C++呢,如果不需要改内部的函数。
这个东西做成库就是让别的语言方便调用的

【在 i******e 的大作中提到】
: 嗯, 如果本身用c,c++写得话,用GSL就方便.如果时间重要的话,还是要用LAPACK.
:
: 时

avatar
i*e
13
因为题目说了要c++库。 GSL虽然是c,但是基本一样。

【在 o**o 的大作中提到】
: 为什么一定要C++呢,如果不需要改内部的函数。
: 这个东西做成库就是让别的语言方便调用的

avatar
x*g
14
sparselib++

【在 e******o 的大作中提到】
: 恩,不错,试了一个20000*20000没问题。多谢
avatar
N*m
15
什么样的矩阵?hermitian的应该很快。

0G).

【在 i******e 的大作中提到】
: 昨晚试了试,发现4000X4000的矩阵算eigenvalue,eigenvector花了5个小时(P4 3.0G).
相关阅读
logo
联系我们隐私协议©2024 redian.news
Redian新闻
Redian.news刊载任何文章,不代表同意其说法或描述,仅为提供更多信息,也不构成任何建议。文章信息的合法性及真实性由其作者负责,与Redian.news及其运营公司无关。欢迎投稿,如发现稿件侵权,或作者不愿在本网发表文章,请版权拥有者通知本网处理。