% this runme script tests out the capabilities of logunrolling
% (logunrolling is logistical logarithmic unrolling of comparagrams).

disp('loading comparagram...');
% load in the comparagram for k = (cbrt(2))^3 = 2 (e.g. K=1)
%load Jsum03.mat
% if you're using matlab, uncompress Jsum03.mat.gz and uncomment above...
load -force Jsum03.txt
% if you're using octave, uncomment above line since octave text more compact

disp('loaded comparagram...');
input('hit any key to display comparagram and logunroll to response function: ');
tvs(log(Jsum03+0.1));
Finv=cls_fast(Jsum03,2);
title('Response function logunrolled from comparagram');
xlabel('Q = log(q)'); ylabel('pixel integer')
axis([min(Finv),max(Finv),0,255]); gset nokey;
plot(Finv,0:255);

axis;

figure;

title('Certainty function recovered from response function');
%xlabel('pixel integer'); ylabel('relative certainty')
%axis([0,255,0,1.05*max(certainty)]); gset nokey;
plotcertaintyfromFinv(Finv);
Jverify=clshow(Finv,2); % verify with slenderized comparagram
input('hit any key to regenerate verifying comparagram from response function :');
tvs(Jverify);