% Example 4: Fit a spline curve to sample points using chunkerfit, with
% varying tolerance and an option to split between sample points.

mip load --install magland/magland/chunkie;

tic;

rng(0)
n = 20;
tt = sort(2*pi*rand(n,1));
r = chnk.curves.bymode(tt, [2 0.5 0.2 0.7]);

opts = [];
opts.ifclosed = true;
opts.cparams = [];
opts.cparams.eps = 1e-3;
opts.pref = [];
opts.pref.k = 16;
chnkr = chunkerfit(r, opts);

opts.splitatpoints = true;
chnkrsp = chunkerfit(r, opts);

opts.cparams.eps = 1e-9;
opts.splitatpoints = false;
chnkr2 = chunkerfit(r, opts);

opts.splitatpoints = true;
chnkrsp2 = chunkerfit(r, opts);

toc;

figure(7)
clf
tiledlayout(2,2,"TileSpacing","compact");
nexttile; plot(chnkr,'k-x'); title(sprintf("nch = %d, eps=1e-3",chnkr.nch));
hold on; plot(r(1,:),r(2,:),'bd');
nexttile; plot(chnkrsp,'k-x'); title(sprintf("nch = %d, eps=1e-3\n (splitting)",chnkrsp.nch));
hold on; plot(r(1,:),r(2,:),'bd');
nexttile; plot(chnkr2,'k-x'); title(sprintf("nch = %d, eps=1e-9",chnkr2.nch));
hold on; plot(r(1,:),r(2,:),'bd');
nexttile; plot(chnkrsp2,'k-x'); title(sprintf("nch = %d, eps=1e-9\n (splitting)",chnkrsp2.nch));
hold on; plot(r(1,:),r(2,:),'bd');