PROGRAM TEST_MATRIX_MULTIPLY
      IMPLICIT NONE
      INTEGER N
      PARAMETER (N=1024)
      DOUBLE PRECISION A(N,N), B(N,N), C(N,N)
      INTEGER I, J
      
      DO J=1,N
         DO I=1,N
              A(I,J) = 1.d0
              B(I,J) = 1.d0
              C(I,J) = 0.d0
	 END DO
      END DO	      
      
      CALL DGEMM ( 'N', 'N', N, N, N, 1.0d0, A, N, B, N,
     $                   0.0d0, C, N )
     
      END PROGRAM TEST_MATRIX_MULTIPLY