% Testroutine zum Programm
% FFT, Blatt 12

% Daten erzeugen
a = rand(16,1);

% Auf Korrektheit testen
if norm(FFT(a)-fft(a)) < 1e-12
  disp('Ihr Programm arbeitet vermutlich korrekt');
else
  disp('Ihr Programm arbeitet NICHT korrekt');
  return
end

% neue Daten erzeugen
a = rand(2048,1);

% Laufzeit f"ur 100 Aufrufe, Ihr Programm
tic,
for k=1:100
  FFT(a);
end
t_user = toc;

% Laufzeit f"ur 100 Aufrufe, Referenzprogramm
tic,
for k=1:100
  FFT_ref(a);
end
t_ref= toc;

q      = round(t_ref/t_user*100)/100;
t_user = round(t_user*100)/100;
t_ref  = round(t_ref *100)/100;

disp(['  '])
disp(['Laufzeiten:'])
disp(['Ihr Programm:    ' num2str(t_user)])
disp(['Referenzprogramm:' num2str(t_ref)])
disp(['Quotient:        ' num2str(q)])
disp(['  '])
disp(['Laufzeiten unterliegen Zufallsschwankungen, mehrfach testen!'])