if len(x_r) == 0:
return
- #print ("outer butterfly")
+ print ("outer butterfly", mode, SVSHAPE.skip, "submode", SVSHAPE.submode2)
# I-DCT, reference (read/write) the in-place data in *reverse-bit-order*
ri = list(range(n))
if jle == 0b111: # all loops end
break
- print("transform2 result", vec)
+ print("inverse_transform2 result", vec)
return vec
# j schedule
SVSHAPE0 = SVSHAPE()
SVSHAPE0.lims = [xdim, 0b0000010, 0]
- SVSHAPE0.submode2 = 0b100
+ SVSHAPE0.submode2 = 0b11
SVSHAPE0.mode = 0b11
SVSHAPE0.skip = 0b00
SVSHAPE0.offset = 0 # experiment with different offset, here
SVSHAPE1 = SVSHAPE()
SVSHAPE1.lims = [xdim, 0b0000010, 0]
SVSHAPE1.mode = 0b11
- SVSHAPE1.submode2 = 0b100
+ SVSHAPE1.submode2 = 0b11
SVSHAPE1.skip = 0b01
SVSHAPE1.offset = 0 # experiment with different offset, here
SVSHAPE1.invxyz = [1,0,1] # inversion if desired