% Reconstruct images $Id: meshsize_param02.m 4839 2015-03-30 07:44:50Z aadler $
imdl= mk_common_model('c2c2',16);
imdl.RtR_prior= @prior_laplace;
imdl.hyperparameter.value = 0.10;

th= linspace(0,2*pi,100); th(end)= [];
xcirc= 0.5+0.1*cos(th); ycirc= 0.1*sin(th);

subplot(121);
show_fem(img); view(0,60);

subplot(122);
show_fem(inv_solve( imdl, vh, vi(1)));
line(xcirc,ycirc,'Color',[0,0.5,0],'LineWidth',2);
axis image
ylabel(sprintf('No. Elems= %d', size(img.fwd_model.elems,1)));

print_convert meshsize_param02a.png '-density 100'
