-- -*- M2-comint -*- {* hash: -898418298 *}

i1 : A = QQ[x,y,z]

o1 = A

o1 : PolynomialRing

i2 : M = matrix(A, {{1,2,3},{4,5,6},{7,8,9}})

o2 = | 1 2 3 |
     | 4 5 6 |
     | 7 8 9 |

             3       3
o2 : Matrix A  <--- A

i3 : Hom(M,A^2)

o3 = | 1 0 4 0 7 0 |
     | 0 1 0 4 0 7 |
     | 2 0 5 0 8 0 |
     | 0 2 0 5 0 8 |
     | 3 0 6 0 9 0 |
     | 0 3 0 6 0 9 |

             6       6
o3 : Matrix A  <--- A

i4 : Hom(A^2,M)

o4 = | 1 2 3 0 0 0 |
     | 4 5 6 0 0 0 |
     | 7 8 9 0 0 0 |
     | 0 0 0 1 2 3 |
     | 0 0 0 4 5 6 |
     | 0 0 0 7 8 9 |

             6       6
o4 : Matrix A  <--- A

i5 : contraHom = (M, s) -> (
         (n,m) := (numgens target M, numgens source M);
         R := mutableMatrix(ring M, s*n, s*m);
         for b from 0 to m-1 do
           for a from 0 to s-1 do
             for c from 0 to n-1 do
               R_(a*n+c,a*m+b) = M_(b,c);
         matrix R
         )

o5 = contraHom

o5 : FunctionClosure

i6 : contraHom(M,2)

o6 = | 1 4 7 0 0 0 |
     | 2 5 8 0 0 0 |
     | 3 6 9 0 0 0 |
     | 0 0 0 1 4 7 |
     | 0 0 0 2 5 8 |
     | 0 0 0 3 6 9 |

             6       6
o6 : Matrix A  <--- A

i7 :