对于fortran矩阵乘法排列顺序的疑惑

今天的fortran学习中发现一个有趣的现象。

附图：

fortran源码：

 

 运行结果：

 

 利用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)->A(3,3)

这样，我再改一下我的输入数组的顺序：

 

运行结果如图：

 

 这样就能和我们熟悉的矩阵对上了。

 另外，如果想按行输入数据，可以利用隐do循环：

read *, ((A(i,j),j=1,3),i=1,3)

这时输入顺序是按行的：A(1,1)->A(1,2)->A(1,3) => A(2,1)->A(2,2)->A(2,3) => A(3,1)->A(3,2)->A(3,3)。