私は科学的なアプリケーションでいくつかのテストを実行しています。このアプリケーションは、LapackdpotrfLapack関数を使用します。私は線形代数についてはあまり気づいていません。dpotrfへの大きな呼び出しをシミュレートしてから、ランダム行列を生成してdpotrfを実行する必要があります。
このマトリックスには特別なパターンが必要ですか？単純なNxN二重行列を生成できますか？
1 に答える
0
dpotrf.fから、DPOTRF は実対称正定行列 A のコレスキー分解を計算します。
したがって、入力は正定で対称であることが期待されます。
また、dpotrf は入力が配列であることを想定しています。つまり、入力は (N*N)x1 です。私は通常、出力を NxN に整形します。
于 2011-02-22T20:25:36.807 に答える