为什么矩阵的索引式列主序的

[VeritasFutureKF]

如题

[zhaosiyuan1098]

我对这里也比较感兴趣，个人理解是：

• cublas为了适配Fortran，默认使用了列优先存储 。在test.cu的实际调用cublasSgemm(err, CUBLAS_OP_N, CUBLAS_OP_N, m, n, k, &alpha, dA, m, dB, k, &beta, dC_ref, m);中，根据api可推断确实使用的是列优先。
• 当把 C / C++ 中行优先保存的二维数组 A 以一维形式输入 cublas 时，会被 cublas 理解为列优先存储。
• 为了保证和cublas的一致性，且避免行列的混乱，统一使用了列主序，从kernel1.cu中的#define A(i,j) A[(i) + (j)*lda] #define B(i,j) B[(i) + (j)*ldb] #define C(i,j) C[(i) + (j)*ldc]中也可推断出kernel中的矩阵是列优先的。

不知道理解的对不对，烦请作者@yzhaiustc解答 ，谢谢。