Сумма умножения всех столбцов одной матрицы (n-на-m) на другую матрицу (n-на-n)

Question

Надеюсь, название не требует пояснений. Я надеялся сделать название таким, чтобы другие тоже могли его найти. Я знаю, как выполнить операцию с циклом, но она должна быть быстрее с использованием некоторого умножения матриц, которое мне интересно изучить.

Код с циклом выглядит примерно так

x <- matrix(rexp(300, rate=.1), nrow=20)
y <- matrix(rexp(400, rate=.1), nrow=20)

res <- as.data.frame(matrix(0,ncol = 15, nrow = 20))

for (i in 1:20){

  res <- res + x*y[,i]

}