author | Puneeth Chaganti <punchagan@fossee.in> |
Fri, 27 May 2011 14:24:59 +0530 | |
changeset 0 | 0efde00f9229 |
permissions | -rw-r--r-- |
#!/usr/bin/env python # 11.5 import os, sys sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"] import pylab as pl from xdync import xdync from polyfuncs import polmul def mv(A, dA, B, dB, C, dC, k, int1): zk = pl.zeros(k+1) zk[-1] = 1 if int1 >= 1: A, dA = polmul([1, -1], 1, A, dA) Fj, dFj, Ej, dEj = xdync(zk, k, A, dA, C, dC)[:4] Gk, dGk = polmul(Ek, dEk, B, dB) S, dS = Fk, dFk R, dR = Gk, dGk return S, dS, R, dR