矩阵乘以行

Question

有没有办法更快地计算矩阵c？

a=matrix(runif(10),2,5)
b=matrix(runif(15),3,5)

c=matrix(,nrow(a)*nrow(b),5)
k=0
for(i in 1:nrow(a)){
  for(j in 1:nrow(b)){
    k=k+1
    c[k,]=a[i,]*b[j,]
  }
}

Answer 1

这是我的版本：

c1 =  a[ rep(1:nrow(a), each = nrow(b)), ] * 
      b[ rep(1:nrow(b), times = nrow(a)), ];
all.equal(c, c1);

> TRUE