clear all
close all
fftw("planner","estimate")
counter = 0;
nrange = 10:10:200;
for n = nrange
  counter = counter + 1;
  A = randn(n,n,n)+1i*randn(n,n,n);
  tic
  for i = 1:50
    fftn(A);
  end
  elapsed(counter) = toc;
end
elapsed_estimate = elapsed;
fftw("planner","measure")
counter = 0;
for n = nrange
  counter = counter + 1;
  A = randn(n,n,n)+1i*randn(n,n,n);
  tic
  for i = 1:50
    fftn(A);
  end
  elapsed(counter) = toc;
end
plot(nrange,elapsed,nrange,elapsed_estimate)
xlabel("n")
ylabel("elapsed time")
legend("planner measure","planner estimate","Location","NorthWest")