問題タブ [eigenvalue]

Python 2.6.5 を使用して、大きな行列 (約 1000x1000) の固有値を計算したいと思います。急にできなくなりました。この質問に対処する他のスレッドは見つかりませんでした。

私が走るとき

MATLAB では約 30 秒かかります。Python での同様のテストには 216 秒かかります。RPy を使用して R で実行しても、計算はそれほど高速化されませんでした。Octave でのテストには 93 秒かかりました。速度の違いに少し戸惑います。

私がオンラインで見つけることができるこのような質問の唯一の例は、数年前のthisです。その質問のポスターは、異なる Python ディレクトリ構造を持っているため (間違っている可能性もありますが、これは投稿の年代によるものだと思います)、特派員が投稿した指示に従おうとするほど自信がありませんでした。

パッケージ マネージャーによると、LAPACK がインストールされており、Python の計算に NumPy と SciPy を使用しています。

私はPythonにかなり慣れていないので、ばかげたことをしたかもしれません。さらに情報を提供する必要がある場合はお知らせください。

基本的に、行列の固有値を見つけようとしていますが、約 12 時間かかります。完了すると、すべての固有ベクトルを見つけることができなかったと表示されます (実際にはほとんど見つかりませんでした)。私が実際にできることは、自分のコードを投稿することだけです。誰かが私にいくつかの提案をしてくれることを願っています. 私は mathematica の経験があまりないので、実行時間が遅く、結果が悪いのは、mathematica の能力ではなく、私に関係があるのか​​もしれません。返信してくれた人に感謝します。本当に感謝しています。

私のコードはもう少しありますが、これは実際に速度が低下しているポイントです。m1 と m2 の両方を 0 に設定しても特に問題はありませんが、m1 を定数に設定するとすべてがうまくいきません。

列ベクトル形式の固有値がN個あります。したがって、これらの固有値に対応するN個の固有ベクトルがあり、固有ベクトル行列を形成します。

今、私が取り組んでいる問題では、固有値の列ベクトルを降順で並べ替える必要があります。対応を維持するために、固有ベクトル行列を固有値と同じ順序で並べ替えるにはどうすればよいですか？

サイズの異なる行列の固有値と固有ベクトルを計算しようとしています。私は非常に単純な Fortran90 コードを使用しており、Ubuntu で動作する私のマシンで利用可能な Intel MKL パッケージに含まれる適切な Lapack ライブラリにリンクしてコンパイルしています。メッセージの末尾に「matrix_diag_01.f90」というコードが添付されています。「random」モジュールには、Numerical Recipes の「ran」乱数ジェネレーターが含まれています。コードは次を使用して適切にコンパイルされます

小さい行列が与えられた場合、実行可能ファイルはうまく機能します。ただし、サイズが 3000x3000 の行列の場合、奇妙な動作が発生します。最初にこのエラーが発生します

ただし、SSYEVD の呼び出しには 3 つのパラメーターしかありません。次に、固有ベクトルを返しますが、固有値は返しません。より大きなメモリを搭載した別のマシンでコンパイルして確認しましたが、結果は同じでした。

誰でも助けてもらえますか？

ありがとう！

私は正方行列 A の次元削減を行っています。私の問題は、13000 x 13000 行列 A の eigvalue 分解の計算に問題があること[v d]=eigs(A)です。これはスパース行列であるため、4GB RAM を使用すると「メモリ不足エラー」が発生します。eigsコマンドの実行時にメモリが使い果たされていないため、これは私の PC の問題ではないと確信しています。私がオンラインで見たヘルプは ARPACK に関するものでした。おすすめサイトをチェックしたのですが、ファイルが多くてどれをダウンロードすればいいのかわかりません。また、MATLABでの使い方もわかりませんでした。別のヘルプでは、数値法を使用すると書かれていますが、具体的にどの方法を使用すればよいかわかりません。どんな解決策でも大歓迎です。

これらのエラーの翻訳と修正方法を教えてください。

現在、仕事をするnumpyを使用しています。しかし、私は数千の行/列を持つ行列を扱っており、後でこの数値は数万に達するため、この種の計算をより高速に実行できるパッケージが存在するかどうか疑問に思っていました.

I'm trying to convert a generalized eigenvalue problem into a normal eigenvalue calculation.

I have this code:

[V,D,flag] = eigs(A, T);

Now I convert it into:

Shouldn't I get the same result? From what I understand in the Matlab documentation, the first equation solves:

and the second one solves:

am I missing something?

Thanks!!

Do you know of a good library/way that I can use to solve an eigen system in c#?

My data is 2D/3D and I want to get direction and length of first and second eigen vectors to evaluate how my data is elongated in 2D space.

Thanks

生成された6つの固有値から球のボリュームを簡単に生成できるようにするコードを少し作成しようとしています。

だから私はデータを持っています：

そして、ボリュームに変換したいと思います。すべての値の重みが等しいと想定しています。これはRで可能ですか？

ありがとう