35 complex a(m*k,np),b(k*n,np)
38 complex cdotu,one,zero
39 parameter(one = 1e0, zero = 0e0)
40 external cdotu,cgemv,cgemm
49 c(1,i) = cdotu(k,a(1,i),1,b(1,i),1)
53 call cgemv(
"T",k,n,one,b(1,i),k,a(1,i),1,zero,c(1,i),1)
59 call cgemv(
"N",m,k,one,a(1,i),m,b(1,i),1,zero,c(1,i),1)
63 call cgemm(
"N",
"N",m,n,k,
64 + one,a(1,i),m,b(1,i),k,zero,c(1,i),m)