source/cluster/wham/src/matmult.f

   1       SUBROUTINE MATMULT(A1,A2,A3)
   2       include 'DIMENSIONS'
   3       DIMENSION A1(3,3),A2(3,3),A3(3,3)
   4       DIMENSION AI3(3,3)
   5       DO 1 I=1,3
   6         DO 2 J=1,3
   7           A3IJ=0.0
   8           DO 3 K=1,3
   9     3       A3IJ=A3IJ+A1(I,K)*A2(K,J)
  10           AI3(I,J)=A3IJ
  11     2   CONTINUE
  12     1 CONTINUE
  13       DO 4 I=1,3
  14       DO 4 J=1,3
  15     4   A3(I,J)=AI3(I,J)
  16       RETURN
  17       END