Question

fortran数组顺序问题

初学fortran，做了个数组相乘的算法
program diagonale
integer a(2,2),b(2,2),c(2,2)
print*,'输入矩阵a'
read*,a
print*,'输入矩阵b'
read*,b
c=0
do i=1,2
do k=1,2
do j=1,2
c(i,k)=c(i,k)+a(i,j)*b(j,k)
enddo
enddo
enddo
print*,c
end
发现这样算出来是错的,然后尝试发现比如一个3x3的数组123456789
按理说a(2,3)应该是6的，但是用fortran输出居然是8，i,j的含义颠倒了，查了书，自己没记错啊。求指教

Answer 1

Fortran 数组是列优先的。即先存储前面的维度。

一个3*3二维数组，在内存中是：
1 2 3 4 5 6 7 8 9
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)

因而，a(2,3)=8