% levi's orbit1 test to see how many leves, iterations at each...
%
% Use: test(3,4) displays result of 3 iterations per level and 4 levels


function testoseehowmanylevelsanditerationsatneeded(itperlevel,numberoflevels)
% returns nothing

more off
l=7;     % for example, l=7 does frames 7-1=6 and frame l=7
         % e.g. l=7 is for frames 6 and 7
A = eval(sprintf('loadframe(%d)',l-1));
figure(1);  % keep from running out of memory (limit fig use)
B = eval(sprintf('loadframe(%d)',l));
disp(sprintf('--- estimating translation frame 0 to %d',l))
p_trans=est_trans_fft2(A,B); % now takes B to A
p=est_pchirp2m(A,B,itperlevel,numberoflevels,p_trans,1);
mse=mean2nan(abs(A-pchirp2(p,B)));
s = sprintf('save test_pairwisem p mse');
disp(s)
eval(s)

