nDiki : 疎行列

PDL と 疎行列

計算しようとしている行列が 2000 x 2000 次元を超えているので、PDLで素直に計算させるとしんどい。 PDLはネイティブに疎行列をサポートしていないので、独自に実装する必要があるので面倒だな(一応 CPAN に PDL::Sparse というのもあるにはある)。

外部で計算させようかとも思ったが、GNU Octave もネイティブサポートなし。

やはりPDL::PPを使って自前でPDL上に実装するしかないか?

Meschach

疎行列の計算に Meschach を使ってみる。 sp_lanczos2 で固有値は計算できたが、これだと固有ベクトルを返さない。

固有値と固有ベクトルの組を得るにはどうすれば。