各位大侠：谁有求大型对称稀疏矩阵特征值的算法（c 语言）？矩阵的储存形式为三列式（行，列，元素）。

矩阵如
1.5 0.0 0.0 ...... 0.3 0.0
0.0 3.2 0.0 ...... 0.0 0.5
0.0 0.0 3.8 ...... 0.6 0.4
.   .   .  .      .   .
.   .   .   .     .   .
.   .   .    .    .   .
0.3 0.0 0.6   ......4.4 0.0
0.0 0.5 0.4   ......0.0 2.2
维数可能是10000
存储形式是：
  1     1   1.5
  1  9999   0.3
  2     2   3.2
  2 10000   0.5
等等
如何求这个矩阵的特征值呢？
本人急用，多谢了！