35 double precision a(m*k,np),b(k*n,np)
36 double precision c(m*n,np)
38 double precision ddot,one,zero
39 parameter(one = 1d0, zero = 0d0)
40 external ddot,dgemv,dgemm
49 c(1,i) = ddot(k,a(1,i),1,b(1,i),1)
53 call dgemv(
"T",k,n,one,b(1,i),k,a(1,i),1,zero,c(1,i),1)
59 call dgemv(
"N",m,k,one,a(1,i),m,b(1,i),1,zero,c(1,i),1)
63 call dgemm(
"N",
"N",m,n,k,
64 + one,a(1,i),m,b(1,i),k,zero,c(1,i),m)
subroutine dmatm3(m, n, k, np, a, b, c)