source/unres/src-HCD-5D/matmult.f

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