FORTRAN可分配数组。这段代码的数学表示是什么？

Question

这段代码中发生了什么？

DOUBLE PRECISION, INTENT(IN) :: x(3,3), y(3,3)
DOUBLE PRECISION :: e(3,3) 

do i=1,3
j=mod(i,3)+1

e(:,i)=x(:,j)-x(:,i)

done

Answer 1

简单地说，这会占用每一列并从它后面的一列中减去它（索引包围数组）。

即。 e的第一列包含x的第二列减去x的第一列，e的第二列包含x的第三列减去e的第二列和第三列包含x的第一列减去第三列。