diff -r 000000000000 -r 0efde00f9229 python/transp.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/python/transp.py Fri May 27 14:24:59 2011 +0530 @@ -0,0 +1,16 @@ +#!/usr/bin/python + +from polyfuncs import polsize +import scipy as sp + +def transp(Q, degQ): + """ Function to transpose a polynomial matrix. """ + rQ, cQ = polsize(Q, degQ) + rP = cQ + cP = rQ + degP = degQ + P = sp.zeros((rP, (degP+1)*cP)) + for i in range(degP+1): + P[..., i*cP:(i+1)*cP] = Q[..., i*cQ:(i+1)*cQ].T + + return P, degP