当我有一个具有维度(i,j,k)的数组和一个具有维度(j,q)的矩阵时。我怎样才能将每个(,, k)乘以该矩阵。一个例子更有意义。
A <- array(c(rep(1,20), rep(2,20), rep(3,20)),dim = c(10,2,3))
B <- matrix(c(1:10), nrow = 2)
# multiply each A[,,i]%*%B
C <- array(NA, dim=c(nrow(A), ncol(B), 3))
C[] <- apply(A, 3, function(x) x%*%B)
我可以通过这种方式得到结果,但我正在寻找一种更有效的方法,例如使用ATensor包。我希望有人可以帮我解决这个问题。