1 #!/usr/bin/python

     2 

     3 from polyfuncs import polsize

     4 import scipy as sp

     5 

     6 def transp(Q, degQ):

     7     """ Function to transpose a polynomial matrix. """

     8     rQ, cQ = polsize(Q, degQ)

     9     rP = cQ

    10     cP = rQ

    11     degP = degQ

    12     P = sp.zeros((rP, (degP+1)*cP))

    13     for i in range(degP+1):

    14         P[:, i*cP:(i+1)*cP] = Q[:, i*cQ:(i+1)*cQ].T

    15 

    16     return P, degP