35 real a(m*k,np),b(k*n,np)
39 parameter(one = 1e0, zero = 0e0)
40 external sdot,sgemv,sgemm
49 c(1,i) = sdot(k,a(1,i),1,b(1,i),1)
53 call sgemv(
"T",k,n,one,b(1,i),k,a(1,i),1,zero,c(1,i),1)
59 call sgemv(
"N",m,k,one,a(1,i),m,b(1,i),1,zero,c(1,i),1)
63 call sgemm(
"N",
"N",m,n,k,
64 + one,a(1,i),m,b(1,i),k,zero,c(1,i),m)
subroutine smatm3(m, n, k, np, a, b, c)