4

2 つの行列をAlgLibで乗算するにはどうすればよいですか

4

3 に答える 3

9

免責事項: 私は AlgLib を使用していません。ドキュメントに書かれているように見えることに従っているだけです。どなたか詳しい方に訂正していただけると幸いです。

とにかく、答えは、次のようにcmatrixgemmor rmatrixgemm(行列が実数か複素数かによって異なります) を使用する必要があるようです。

rmatrixgemm(m,n,k, 1, A,0,0,0, B,0,0,0, 0, C,0,0);

どこ：

  • m, n,kは行列のサイズです ( Ais mby k, Bis kby n, Cis mby n)
  • the1は積に掛ける値です (たまたま AB ではなく 3AB が必要な場合は、3そこに置きます)
  • A,0,0,0およびグループは次のB,0,0,0とおりです: 行列、行オフセット、列オフセット、演算タイプ
    • 演算の種類は、AやBをそのまま使う場合は0、転置を使う場合は1、共役転置を使う場合は2です(もちろん に2を使うことはできませんrmatrixgemm) 。
  • 次0は、結果に 0*C を追加することを示しています (ここに 0 を入力すると、C の初期値は完全に無視されます)。
  • C の後の 2 つ0の s は、行と列のオフセットです

おそらく、そのレベルの一般性は過剰であり、これらのデフォルトを提供するより単純な関数が必要であると考えるかもしれません。私はそれに同意しませんが、私が見る限り、AlgLib にはこれほど単純な関数はありません。独自に記述したい場合があります (rmatrixgemmまたはを呼び出すだけcmatrixgemmです)。

(なぜそんなに一般的なのでしょうか? 効率的な行列乗算を行うには非常に洗練されたコードが必要であり、それは本質的に、より一般的なC=a.f(A).g(B)+b.C操作を行うために必要な非常に洗練されたコードと同じ*matrixgemmであり、より一般的な操作が役立つ場合があるためです。)

役に立つかもしれないいくつかの発言を追加するために編集されました。

  • オフセットは、部分行列で処理できるようにするためのものです。これを実行できることは、一部の数値アルゴリズムで役立ちます。m、n、使用しているサブマトリックスkのサイズです。一般的なケースでは、それらは配列の次元と同じになり、オフセットはゼロになります。
  • rmatrixgemmまたはを呼び出す前に、配列自体が存在し、適切なサイズである必要がありますcmatrixgemm。少なくとも、AそしてB確かにそうです。Cとして渡されるため、エントリrefにある場合、これらの関数がそれを作成する可能性があります。null
  • rmatrixgemmorの署名から考えると、 は参照によって渡されるのにコピーされるかもしれませんが、C cmatrixgemm#のセマンティクスについて完全に混乱していなければ、それらはすべて実質的に (オブジェクト) 参照によって渡されます。ABC
于 2011-04-09T20:39:53.167 に答える
4

ガレスが書いたことを確認するためだけに：

double[,] a = new double[,] {
    {1,2,3},
    {4,5,6}
};
double[,] b = new double[,] {
    {7,8,9,10},
    {11,12,13,14},
    {15,16,17,18}
};
int m = a.GetLength(0);
int n = b.GetLength(1);
int k = a.GetLength(1);
double[,] c = new double[m,n];
alglib.rmatrixgemm(m, n, k, 1, a, 0,0,0, b,0,0,0, 0, ref c, 0,0);
//c = {{74, 80, 86, 92}, {173, 188, 203, 218}}
于 2011-04-09T22:44:58.093 に答える
0

vba では、この関数の複雑なバージョンを使用できました。

Alpha.x = 1: Alpha.y = 0
Beta.x = 0:  Beta.y = 0

Call CMatrixGEMM(4, 1, 4, Alpha, r, 0, 0, 0, x, 0, 0, 0, Beta, RX, 0, 0)

補足として、すべてのalglibモジュールを1つのアクセスデータベースにロードし、関数が必要な現在のデータベースからそのデータベースへの参照を設定することにより、alglib関数のセット全体を任意のアクセスプログラムにロードできます。これにより、作業データベースにとって非常に便利で軽量になります。

于 2013-10-09T20:30:11.660 に答える