gemmとmatmul

メモ

      real(8):: A(M,K),B(K,N),C(M,N)

      C=matmul(A,B)

と等価になるのが

      call dgemm('N','N',M,N,K,1.0d0,A,M,B,K,1.0d0,C1,M)