|
Last change
on this file since 3567 was
57,
checked in by aslmd, 14 years ago
|
|
mineur LMD_MM_MARS: ajout du GCM ancienne physique, systeme maintenant complet sur SVN (ne manque que la base de donnees d'etats initiaux)
|
|
File size:
623 bytes
|
| Line | |
|---|
| 1 | SUBROUTINE MXVA (SA,IAC,IAR,SB,IB,SC,IC,NRA,NCA) |
|---|
| 2 | IMPLICIT NONE |
|---|
| 3 | |
|---|
| 4 | INTEGER IAC,IAR,IB,IC,NRA,NCA |
|---|
| 5 | REAL SA(*), SB(*), SC(*) |
|---|
| 6 | |
|---|
| 7 | INTEGER I,J |
|---|
| 8 | |
|---|
| 9 | c initialize product |
|---|
| 10 | DO 100 I = 1, NRA |
|---|
| 11 | SC( 1 + (I-1)*IC ) = 0. |
|---|
| 12 | c ( C(i) := 0. ) |
|---|
| 13 | 100 CONTINUE |
|---|
| 14 | c multiply matrix from sa and vector from sb |
|---|
| 15 | DO 220 J = 1, NCA |
|---|
| 16 | DO 210 I = 1, NRA |
|---|
| 17 | SC( 1 + (I-1)*IC ) |
|---|
| 18 | & = SC( 1 + (I-1)*IC ) |
|---|
| 19 | & + SA( 1 + (I-1)*IAC + (J-1)*IAR ) |
|---|
| 20 | & * SB( 1 + (J-1)*IB ) |
|---|
| 21 | c ( C(i) := C(i) + A(i,j)*B(j) ) |
|---|
| 22 | 210 CONTINUE |
|---|
| 23 | 220 CONTINUE |
|---|
| 24 | RETURN |
|---|
| 25 | END |
|---|
Note: See
TracBrowser
for help on using the repository browser.
