I’m running Julia 0.5 on Linux, and I found this weird thing, which I think could be a bug:
x = randn(10); # input value x2 = copy(x); # copy of input value y = zeros(10); # output A = plan_idct(x); # plan an inverse discrete cosine transform A_mul_B!(y,A,x); # apply transform in place in y (1) norm(x-x2) # here x should not be modified >2.7... # but it's not!!
In line (1) if change
plan_ifft (in the latter 2 cases you need complex vector)
this issue does not happen.
The problem is solved when using
y .= A*x. But my feeling is that something is wrong with
Please let me know if you experience the same issue or if this was already fixed in julia 0.6!