a+bt.f

来自「这是一个实用的并行计算源代码」· F 代码 · 共 19 行

      PARAMETER(N=4000, NB=32)
      REAL A(N,N), B(N,N)
*
      DO J = 1, N
         DO I = 1, N
            A(I,J) = 1.D0
            B(I,J) = 1.D0
         ENDDO
      ENDDO
*
      DO J = 1, N
      DO I = 1, N
         A(I,J) = A(I,J) + B(J,I)
      ENDDO
      ENDDO
*
      STOP
      END