make cp src-HCD-5D
[unres.git] / source / unres / src-HCD-5D / matmult.f
index e9257cf..748df4d 100644 (file)
@@ -1,11 +1,12 @@
       SUBROUTINE MATMULT(A1,A2,A3)
-      implicit real*8 (a-h,o-z)
+      IMPLICIT NONE
       include 'DIMENSIONS'
-      DIMENSION A1(3,3),A2(3,3),A3(3,3)
-      DIMENSION AI3(3,3)
+      DOUBLE PRECISION A1(3,3),A2(3,3),A3(3,3)
+      DOUBLE PRECISION AI3(3,3),A3IJ
+      integer I,J,K
       DO 1 I=1,3
         DO 2 J=1,3
-          A3IJ=0.0
+          A3IJ=0.0D0
           DO 3 K=1,3
     3       A3IJ=A3IJ+A1(I,K)*A2(K,J)
           AI3(I,J)=A3IJ