! define dgemul subroutine

#define dgemul(a,lda,transa,b,ldb,transb,c,ldc,l,m,n) \
dgemm(transa,transb,l,n,m,1.0,a,lda,b,ldb,1.0,c,ldc)

#define DGEMUL(a,lda,transa,b,ldb,transb,c,ldc,l,m,n) \
dgemm(transa,transb,l,n,m,1.0d0,a,lda,b,ldb,0.0d0,c,ldc)