Fortranで最初の固有値と固有ベクトルを計算する方法

Question

いくつかの例（here）を使用してARPACKを使用しようとしましたが、私の行列を入力する方法も分かりませんでした。このquestionから、PythonとMatlabの実装はARPACKの複雑さを避ける唯一の解決策であるようです。これを行ういくつかの（Fortranベースの）方法がありますか、すべての固有値/固有ベクトルを計算することも避けますか？

Answer 1

行列Aと固有値の数（例えば、最初の固有値のみが必要な場合はnev = 1）を設定する必要があるラッパーを記述しました。このコードは、hereという特定のアプリケーションに適合しています。 私のモジュールのような版のみ行列、機能の数、それらのが取られるべきである固有値との数（「LM」を設定する必要があります - 。最大の大きさ」を

コードを要旨hereとして入手可能である

NB：これは、大規模なスパース行列のために最適化されている間、それは、通常、密行列のために最適化されているため、ARPACKは、LAPACK異なる

ARPACKパッケージは 対称行列、ランチョス の対応する変異体の場合に 暗黙リスタートアーノルディ法（IRAM）またはを使用して、いくつかの固有値を計算し、大規模なスパースまたは構造化行列の 固有ベクトルに対応するように設計されていますアルゴリズム。これは、 の機能を提供するために、SciPy、Mathematica、GNU Octave、MATLABなど多くの一般的な数値計算環境 で使用されています。 （出典：ウィキペディア）