Estoy luchando por dividir una matriz sin usar bucles. Quiero dividir una matriz de 9x9 en nueve matrices de 3x3 (que no se superpongan). Intenté split()pasar un vector como agrupación de "factores" sin éxito. También intenté usarlo asplit()sin éxito. Creo que no entiendo cómo funciona "MARGIN" asplit(). Puedo obtener los "cubos" de 3x3 usando bucles, pero me pregunto si hay una manera de hacerlo con funciones de utilidad simples.

m <- matrix(1:81, 9, 9)

Quiero almacenar las 9 matrices de 3x3 resultantes en una lista. El resultado esperado para el primer cubo es:

cube[[1]]

     [,1] [,2] [,3]
[1,]   01   10   19
[2,]   02   11   20
[3,]   03   12   21

Lo intenté

# all these split by row
indices <- matrix(1:9, 3, 3)
cubes <- split(m, as.vector(indices))

indices <- array(1:9, dim = c(3, 3, 9))
cubes <- split(m, as.vector(indices))

indices <- array(matrix(1:9, 3, 3), dim = c(3, 3, 9))
cubes <- split(m, as.vector(indices))

# this gives an error
asplit(m, MARGIN = c(3,1))