dgetrfを介して、行優先の行列で LAPACK から呼び出していclapack.hます。
私はマトリックスを持っていA = [4,9,2; 3,5,7; 8,1,6]ます。
dgetrf を呼び出すと、結果は[9.0, 0.222222, 0.444444; 5.0, 5.888888, 0.132075; 1.0, 5.777777, 6.792453]. ただし、結果は.[8.0, 1.0, 6.0; 0.5, 8.5, -1.0; 0.375, 0.544118, 5.294118]
を呼び出す前に A をdgetrf転置し、出力をもう一度転置すると、正しい結果が得られます。
私が使用している呼び出し is:clapack_dgetrf(CblasRowMajor, 3, 3, A, 3, ipiv); ipivはもちろん、サイズ 3 の整数配列であり、 A はサイズ 9: の double の配列です[4,9,2,3,5,7,8,1,6]。
私が狂った場合に備えて、 order 引数を に切り替えてみましたがCblasColMajor、それも期待どおりに機能しません。
これはあるべき姿ですか？