|
0
|
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 |
|