main
| Line | |
|---|
| 1 | PROGRAM MXMDRIVER
|
|---|
| 2 |
|
|---|
| 3 | INTEGER N
|
|---|
| 4 | PARAMETER (N=4)
|
|---|
| 5 | real A(N,N),B(N,N),C(N,N)
|
|---|
| 6 | INTEGER I,J
|
|---|
| 7 |
|
|---|
| 8 | DO J=1,N
|
|---|
| 9 | DO I=1,N
|
|---|
| 10 | A(I,J) = 1.0*(I+J)
|
|---|
| 11 | B(I,J) = 1.0*(I-J)
|
|---|
| 12 | ENDDO
|
|---|
| 13 | ENDDO
|
|---|
| 14 |
|
|---|
| 15 | c DO I=1,N
|
|---|
| 16 | c DO J=1,N
|
|---|
| 17 | c C(I,J) = 0.0
|
|---|
| 18 | c DO K=1,N
|
|---|
| 19 | c C(I,J) = C(I,J) + A(I,K)*B(K,J)
|
|---|
| 20 | c ENDDO
|
|---|
| 21 | c ENDDO
|
|---|
| 22 | c ENDDO
|
|---|
| 23 |
|
|---|
| 24 | call mxm44_0(A,N,B,N,C,N)
|
|---|
| 25 |
|
|---|
| 26 | DO J=1,N
|
|---|
| 27 | DO I=1,N
|
|---|
| 28 | PRINT *,'C(',I,',',J,') = ',C(I,J)
|
|---|
| 29 | ENDDO
|
|---|
| 30 | ENDDO
|
|---|
| 31 |
|
|---|
| 32 | END
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.
