source:
CIVL/mods/dev.civl.abc/examples/fortran/argonne/MXM/ex1c.F
| Last change on this file was aad342c, checked in by , 3 years ago | |
|---|---|
|
|
| File size: 409 bytes | |
| Rev | Line | |
|---|---|---|
| [b7b71e8] | 1 | subroutine mxm(A,N1,B,N2,C,N3) |
| 2 | c | |
| 3 | real a(n1,n2),b(n2,n3),c(n1,n3),s | |
| 4 | ||
| 5 | c$OMP PARALLEL DEFAULT(PRIVATE) SHARED(A,B,C,N1,N2,N3) | |
| 6 | c$OMP DO | |
| 7 | do j=1,n3 | |
| 8 | do i=1,n1 | |
| 9 | s = 0.0 | |
| 10 | do k=1,n2 | |
| 11 | s = s + a(i,k)*b(k,j) | |
| 12 | enddo | |
| 13 | c(i,j) = s | |
| 14 | enddo | |
| 15 | enddo | |
| 16 | c$OMP END DO | |
| 17 | c$OMP END PARALLEL | |
| 18 | ||
| 19 | return | |
| 20 | end |
Note:
See TracBrowser
for help on using the repository browser.
