diff -r 000000000000 -r 0efde00f9229 minv/mv.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/minv/mv.py Fri May 27 14:24:59 2011 +0530 @@ -0,0 +1,20 @@ +#!/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 +