y****d
发帖数: 52 | 1 需要求一个矩阵的本征值和本征向量(eigenvalue/eigenvector)，矩阵比较大，1000000*
1000000 (1百万*1百万)，但有一定周期结构，并且绝大多数是0. 而且我只需要本征值最
低的10到100个结果。
各位有没有推荐的子程序？ 最好是Fortran or C写的。
PS，SVD我已经试过了，奇慢。而且没有必要把所有的eigenvalue都接出来。
多谢各位！ |
l*********i
发帖数: 483 | 2 晕,要用好多内存吧?
1000000*
最
【在 y****d 的大作中提到】
: 需要求一个矩阵的本征值和本征向量(eigenvalue/eigenvector)，矩阵比较大，1000000*
: 1000000 (1百万*1百万)，但有一定周期结构，并且绝大多数是0. 而且我只需要本征值最
: 低的10到100个结果。
: 各位有没有推荐的子程序？ 最好是Fortran or C写的。
: PS，SVD我已经试过了，奇慢。而且没有必要把所有的eigenvalue都接出来。
: 多谢各位！
|
t*s
发帖数: 1504 | 3 用那个啥Rayleigh...忘了
1000000*
最
【在 y****d 的大作中提到】
: 需要求一个矩阵的本征值和本征向量(eigenvalue/eigenvector)，矩阵比较大，1000000*
: 1000000 (1百万*1百万)，但有一定周期结构，并且绝大多数是0. 而且我只需要本征值最
: 低的10到100个结果。
: 各位有没有推荐的子程序？ 最好是Fortran or C写的。
: PS，SVD我已经试过了，奇慢。而且没有必要把所有的eigenvalue都接出来。
: 多谢各位！
|
a**a
发帖数: 416 | 4 en, 好像有稀疏矩阵的特征值迭代解法。去netlib上翻翻故纸堆吧。算法是GRMES类的。
【在 t*s 的大作中提到】
: 用那个啥Rayleigh...忘了
:
: 1000000*
: 最
|
