LAPACK dgetrf输入和输出应该转置吗?

时间:2012-08-24 22:16:53

标签: c lapack blas atlas

我通过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]

如果我在调用dgetrf之前转换A,然后再转置输出,我会得到正确的结果。

我正在使用的电话是:clapack_dgetrf(CblasRowMajor, 3, 3, A, 3, ipiv); ipiv当然是一个大小为3的整数数组,A是一个双打数组,大小为9:{{1} }。

我已经尝试将订单参数切换为[4,9,2,3,5,7,8,1,6],以防万一我疯了,而且也没有按预期工作。

这是应该的吗?

0 个答案:

没有答案