
fftw ('planner','measure'); 

Nmax = 5;
for h = 1:Nmax
  for i = 1:Nmax
    for j = 1:Nmax
      for k = 1:Nmax
        fprintf ('row=%d col=%d k=%d h=%d      \r', i, j, k, h);
        x1 = rand (i, j, k, h);
        y = fft (x1);
      %  x2 = ifft (y);
      end
    end
  end
end


