changeset 0 | 0efde00f9229 |
-1:000000000000 | 0:0efde00f9229 |
---|---|
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 |