|  | SUBROUTINE SMXM(A,N,B,M,C,L) | 
|  | ** | 
|  | **    C = A * B | 
|  | **    A ET B MATRICES A(N,M)  B(M,L) ==> C(N,L) | 
|  | ** | 
|  | *>A      PREMIERE MATRICE | 
|  | *>N      PREMIERE DIMENSION DE A ET DE C | 
|  | *>B      DEUXIEME MATRICE | 
|  | *>M      DEUXIEME DIMENSION DE A ET PERMIERE DE B | 
|  | *<C      MATRICE PRODUIT DE A ET DE B | 
|  | *>L      DEUXIEME DIMENSION DE B ET DE C | 
|  | *A R. SANCHEZ ( EARLY WINTER 1987 ) | 
|  | *V M.F. ROBEAU | 
|  | *M AM BAUDRON - AVRIL 94 | 
|  | *: ERREUR DANS L'APPEL A L'UTILITAIRE SGEMM | 
|  | *: APPEL A L'UTILITAIRE SGEMM DE LA LIBRAIRIE BLAS SUR HP | 
|  | *M AM BAUDRON - NOVEMBRE 1991 | 
|  | *: ERREUR ( SOMME SUR LES TERMES PAS FAITE ) | 
|  | *: APPEL A L'UTILITAIRE SGEMM DE LA LIBRAIRIE BLAS SUR RISC | 
|  | *M AM BAUDRON - MAI 1993 | 
|  | *: CHANGEMENT DES %IF LOCAL SUN MIPS SUITE A INTRODUCTION VERSION IBM | 
|  | DIMENSION A(N,M),B(M,L),C(N,L) | 
|  | DO 20 I=1,N | 
|  | DO 20 J=1,L | 
|  | R=0. | 
|  | DO 10 K=1,M | 
|  | R=R+A(I,K)*B(K,J) | 
|  | 10    CONTINUE | 
|  | C(I,J)=R | 
|  | 20 CONTINUE | 
|  | RETURN | 
|  | END |