相信很多小伙伴在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) !定义数组CC=matmul(A,B) !C返回A,B两个矩阵的乘积print *, Cprint *, '*********'print *, C(1,:) !打印C第一行print *, C(2,:) !打印C第二行pauseend运行结果如下:
为了验证,用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)
这样,再改一下原输入数组的顺序:
!integer :: A(2,3)=(/1,2,3,4,5,6/),B(3,2)=(/1,2,3,4,5,6/) integer :: A(2,3)=(/1,4,2,5,3,6/),B(3,2)=(/1,3,5,2,4,6/) !定义A,B两个数组integer :: C(2,2) !定义数组CC=matmul(A,B) !C返回A,B两个矩阵的乘积print *, Cprint *, '*********'print *, C(1,:) !打印C第一行print *, C(2,:) !打印C第二行pauseend运行结果如下:
这样就能和我们熟悉的矩阵对上了。
另外,如果想按行输入数据,可以利用隐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)。