Fortran语言——矩阵乘法排列顺序问题

相信很多小伙伴在fortran学习的过程中都会发现一些有趣的现象。

如下列一段代码，对两个矩阵做乘法：

integer :: A(2,3)=(/1,2,3,4,5,6/),B(3,2)=(/1,2,3,4,5,6/)    !定义A,B两个数组
integer :: C(2,2)      !定义数组C
C=matmul(A,B)      !C返回A,B两个矩阵的乘积
print *, C
print *, '*********'
print *, C(1,:)           !打印C第一行
print *, C(2,:)           !打印C第二行
pause
end

 运行结果如下：

 为了验证，用matlab计算正确的结果：

那么问题来了，fortran中直接打印C数组是横着排列的，排列顺序是C(1,1)，C(1,2)，C(2,1)，C(2,2)，而利用print *，C(1,:)和print *，C(2,:)将C排列成MATLAB中的结果，发现结果并不一样，这是为什么呢？

print *, '***********************************'

后来才知道怎么回事儿：

fortran编译器在读数时存储数据是按列来的，即读入顺序为A(1,1)->A(2,1)->A(3,1) => A(1,2)->A(2,2)->A(3,2) => A(1,3)->A(2,3)-&