A <- array(1, dim = c(10, 3, 6))
B <- array(1:90, dim = c(3, 10, 3))

for (i in 2:6) {
  foreach (j = 1:3, .combine=cbind) %dopar% {
     colSums(t(A[,,i-1]) * B[,,j])
  }
}