source:
CIVL/mods/dev.civl.abc/examples/fortran/simple/mxm_naive.f
| Last change on this file was aad342c, checked in by , 3 years ago | |
|---|---|
|
|
| File size: 518 bytes | |
| Rev | Line | |
|---|---|---|
| [aad342c] | 1 | PROGRAM MXM_NAIVE |
| 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 | DO I=1,N | |
| 15 | DO J=1,N | |
| 16 | C(I,J) = 0.0 | |
| 17 | DO K=1,N | |
| 18 | C(I,J) = C(I,J) + A(I,K)*B(K,J) | |
| 19 | ENDDO | |
| 20 | ENDDO | |
| 21 | ENDDO | |
| 22 | ||
| 23 | DO J=1,N | |
| 24 | DO I=1,N | |
| 25 | PRINT *,"C(",I,",",J,") = ",C(I,J) | |
| 26 | ENDDO | |
| 27 | ENDDO | |
| 28 | END |
Note:
See TracBrowser
for help on using the repository browser.
