source/wham/src-restraints-DFA/matmult.f

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