Ich habe zwei Matrizen in R, die ich multiplizieren möchte:

a = matrix(rnorm(20*10000, mean=0, sd=5), 20, 10000)
b = matrix(rnorm(20*10000, mean=0, sd=5), 20, 10000)
t(a)%*%b

Gegeben, dass die Dimension dieser Matrixmultiplikation viel Zeit in Anspruch nimmt, gibt es eine spezielle Möglichkeit, Berechnungen zu beschleunigen? Und gibt es irgendwelche eingebauten Funktionen in R, um solche Multiplikationen schneller zu machen?