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 | |
| Line | |
|---|---|
| 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.
