|
Last change
on this file since 937 was
38,
checked in by emillour, 14 years ago
|
|
Ajout du modè Martien (mon LMDZ.MARS.BETA, du 28/01/2011) dans le rértoire mars, pour pouvoir suivre plus facilement les modifs.
EM
|
|
File size:
729 bytes
|
| Line | |
|---|
| 1 | SUBROUTINE MXM(A,NAR,B,NAC,C,NBC) |
|---|
| 2 | c |
|---|
| 3 | c |
|---|
| 4 | c Multiplication of matrix A by matrix B ,and result in matrix C |
|---|
| 5 | c |
|---|
| 6 | c |
|---|
| 7 | c NAR : number of rows for matrix A and matrix C |
|---|
| 8 | c NAC : number of columns for matrix A and number of rows |
|---|
| 9 | c for matrix B |
|---|
| 10 | c NBC : number of columns for matrix B and matrix C |
|---|
| 11 | c |
|---|
| 12 | c |
|---|
| 13 | IMPLICIT NONE |
|---|
| 14 | c |
|---|
| 15 | INTEGER NAR,NAC,NBC, I,J,K |
|---|
| 16 | REAL A(NAR,NAC), B(NAC,NBC), C(NAR,NBC) |
|---|
| 17 | c |
|---|
| 18 | DO 10 J = 1,NBC |
|---|
| 19 | DO 10 I = 1,NAR |
|---|
| 20 | C(I,J)=0. |
|---|
| 21 | 10 CONTINUE |
|---|
| 22 | c |
|---|
| 23 | DO 50 J = 1,NBC |
|---|
| 24 | DO 40 I = 1,NAR |
|---|
| 25 | |
|---|
| 26 | DO 20 K = 1,NAC |
|---|
| 27 | C(I,J) = C(I,J) + A(I,K) * B(K,J) |
|---|
| 28 | 20 CONTINUE |
|---|
| 29 | |
|---|
| 30 | 40 CONTINUE |
|---|
| 31 | 50 CONTINUE |
|---|
| 32 | |
|---|
| 33 | RETURN |
|---|
| 34 | END |
|---|
Note: See
TracBrowser
for help on using the repository browser.
