nDiki : 固有ベクトル

固有値 行列 PDL 固有値分解 ベクトル 疎行列 Perl Octave usr ランチ インストール 線形空間 Meschach バグ 単体テスト CPAN CD-ROM Red Hat Linux 開発 Perl v5.8.2 GD ビルド

■ [ Perl ] PDL::Math::eigens

各固有ベクトルは返り値の行列に「行ベクトル」として格納されているのか。 列ベクトルだと思ってコーディングしてあったので、残りの全ての計算が全部嘘だったという事が判明。 アウチ。

修正自体は、1個所 xchg(0, 1) するだけ。

• [ Perl ] PDL 2.4.1 から eigens が固有ベクトルを... (2004-03-05)
• SQLite とか DbUnit とか (2005-05-23)
• [ Perl ] PDL::PP で C extension を書く (2004-02-19)
• PDL と 疎行列 (2004-06-14)
• 昼食難民になるところだったラゾーナ川崎プラザ (2008-08-03)

■ [ Perl ] PDL 2.4.1 から eigens が固有ベクトルを列ベクトルで返すように

PDLを使う開発中のプログラムを「速いマシンで使いたい」という要望から、オフィスにある Red Hat Linux 8.0 BOX へ昨日インストール作業を行った。

あいかわらず rpm を探してきてインストールするのが面倒。GDを使うのに freetype の開発パッケージを入れようとしたら、CD-ROM に入っている rpm セットだと freetype とその devel のバージョンが違っていてはいらないし。

結局がしがし /usr/local にインストール。 Perl v5.8.2 自体もそのプログラム専用に別途インストールしてCPANモジュールでがしがし依存するものをインストール。

でようやく環境が整って開発中のプログラムをビルドしたら単体テストでこける。 今日数値をチェックしたら、PDL の eigens 関数が今までと転置した固有ベクトル行列を返すようになっていた。 コレ、以前、列ベクトルだと思ったら行ベクトルだったってやつ。それがまた逆になったらしい。

PDLのドキュメントを見ると確かに、2.4.1 でバグフィックスとして列ベクトルに修正されたらしい。 しかし古いPDLもサポートせねばならない。ということでPDLのバージョンを見て 2.4.0 以下のみ xchg かけるように修正して、問題収束。

• CPAN のPerl モジュールを RPM パッケージ化する (2003-04-22)
• 野良パッケージと依存 Perl モジュールのインストールセット をCPAN... (2006-02-11)
• PAR を ActivePerl 5.6.1 build 638 に (2004-07-20)
• SQLite とか DbUnit とか (2005-05-23)
• 今日のさえずり - 京都の小学校のコンピュータ室にいったら、Squeak が (2008-03-06)

■ Meschach

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

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

• [ Perl ] PDL 2.4.1 から eigens が固有ベクトルを... (2004-03-05)
• PDL と 疎行列 (2004-06-14)
• [ Perl ] PDL::Math::eigens (2003-12-04)