% 2d noiselet transform matrix: matrix vector multiplication
%
% (c) by Hannes Nickisch, MPI for Biological Cybernetics, 2010 July 15

function y = mvm(A,x,ctransp)
  x = reshape(x,A.imsz);
  if ctransp
      error('not implemented')
  else
    if isreal(x)
      y = realnoiselet(x);
    else
      y = realnoiselet(real(x)) + 1i*realnoiselet(imag(x));
    end
  end
  y = y/sqrt(A.imsz(1));