1 #!/usr/bin/env python

     2 # 11.9

     3 import os, sys

     4 sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]

     5 

     6 import pylab as pl

     7 from xdync import xdync

     8 from zpowk import zpowk

     9 from polyfuncs import polsplit3, polmul

    10 

    11 def mv(A, dA, B, dB, C, dC, k, int1):

    12     zk, dzk = zpowk(k)

    13     Bzk, dBzk = polmul(B, dB, zk, dzk)

    14     

    15     Bg, Bb = polsplit3(B)

    16     Bbr = Bb[::-1]

    17 

    18     RHS = pl.convolve(C, pl.convolve(Bg, Bbr))

    19     dRHS = len(RHS) - 1

    20 

    21     Sc, dSc, Rc, dRc = xdync(Bzk, dBzk, A, dA, RHS, dRHS)[:4]

    22     return Sc, dSc, Rc, dRc

    23