Предположим, у меня есть такая матрица:

a = 
    1    2
    3    4

Я хочу удвоить размер матрицы и создать что-то вроде этого:

aa = 
     1    1    2    2
     1    1    2    2
     3    3    4    4
     3    3    4    4

таким образом, каждый элемент в первой матрице распространяется на четыре элемента в большей матрице.

a(i,j) == aa(2*i-1, 2*j-1)
       == aa(2*i  , 2*j-1)
       == aa(2*i-1, 2*j)
       == aa(2*i  , 2*j)

Есть ли предопределенные функции для этого?

определенно я могу сделать это двумя циклами, но я хочу самый простой и чистый способ!