Fortran中的矩陣生成

Question

問題：我需要生成一個M * N矩陣，其中M表示行，N表示列。 這個矩陣的元素應該像（i-1.0）* N + j。 此代碼不起作用，我不明白爲什麼。 我使用的是Compaq Visual Fortran 6.

我除了這個代碼必須生成一個矩陣。生成後的基質應印像

1.0 2.0 3.0

4.0 5.0 6.0

爲3×2的矩陣。

現在，這段代碼根本行不通，顯示了一串Windows執行錯誤。

INTEGER M, N, I, J 
COMMON /BLK1/ MATA (256, 256) 
READ (7, 100) M 
READ (7, 100) N 
100 FORMAT (I3) 
CALL MATGEN (M, N) 
DO 90003 I = 1, M 
DO 90004 J = 1, N 
WRITE (8, 101) MATA(I, J) 
101 FORMAT (F8.1, ' ', \) 
90004 CONTINUE 
WRITE (8, 102) 
102 FORMAT (/, /) 
90003 CONTINUE 
END 



SUBROUTINE MATGEN (M, N) 
REAL DUM 
INTEGER I, J 
COMMON /BLK1/ MATA (256, 256) 
DUM = 1.0 
DO 90001 I = 1, M 
DO 90002 J = 1, N 
MATA (I, J) = DUM + 1.0 
90002 CONTINUE 
90001 CONTINUE 
RETURN 
END 

COMMON BLK1 
REAL MATA (256, 256) 
END 

Answer 1

有幾個問題與您的代碼：

• 該方案並沒有真正啓動（在本文下面的代碼是PROGRAM TEST）
• MATA（隱含）的整數。我添加了聲明。未來請使用IMPLICIT NONE！
• 有在不真正屬於那裏結束一些奇怪的代碼...
• 作爲francescalus指出，分配不增加，所以MATA充滿2。
• 我的版本的gfortran不接受反斜槓來抑制換行符。我不得不使用$。詳細信息請參見here。
• Agentp表示你沒有打開第七單元的輸入。我改爲從STDIN改爲讀取。
• 縮進不會傷害！

 PROGRAM TEST 
     IMPLICIT NONE 
     INTEGER M, N, I, J 
     REAL MATA 
     COMMON /BLK1/ MATA(256, 256) 
     WRITE(*,*) 'Please enter M:' 
     READ (*, 100) M 
     WRITE(*,*) 'Please enter N:' 
     READ (*, 100) N 
100 FORMAT (I3) 

     CALL MATGEN (M, N) 
     DO 90003 I = 1, M 
     DO 90004 J = 1, N 
      WRITE (8, 101) MATA(I, J) 

101 FORMAT (F8.1, ' ', $) 
90004 CONTINUE 
     WRITE (8, 102) 
102 FORMAT (/, /) 
90003 CONTINUE 
     END 

     SUBROUTINE MATGEN (M, N) 
     IMPLICIT NONE 
     INTEGER M, N 
     REAL DUM 
     INTEGER I, J 
     REAL MATA 
     COMMON /BLK1/ MATA(256, 256) 
     DUM = 0.0 
     DO 90001 I = 1, M 
     DO 90002 J = 1, N 
      DUM = DUM + 1. 
      MATA (I, J) = DUM 
90002 CONTINUE 
90001 CONTINUE 
     RETURN 
     END 

這導致：

$ ./a.out 
Please enter M: 
2 
Please enter N: 
3 

$ cat fort.8 
    1.0  2.0  3.0 


    4.0  5.0  6.0