dirname = dirOuter;
recMat = reshape(hsolCVU,[seti.nROI, seti.nROI]);
figure(101); imagesc(real(recMat)); title('reconstruction (real)'); colorbar;
axis xy;
savePngFigSimple(101,dirname,'rec_real',iPda);
figure(102); imagesc(imag(recMat)); title('reconstruction (imag)'); colorbar;
axis xy;
savePngFigSimple(102,dirname,'rec_imag',iPda);
v1 = real(recMat(round(seti.nROI/2),:));
qExaMat = reshape(seti.qROIexact,[seti.nROI,seti.nROI]);
v2 = real(qExaMat(round(seti.nROI/2),:));
figure(103); plot(1:length(v1),v1,1:length(v2),v2);legend('rec','exact');
savePngFigSimple(103,dirname,'rec_schnitt',iPda);
if seti.invNo ~= 1
figure(104); plot(1:iPda,errInPda(1:iPda)); title('error in pda');
savePngFigSimple(104,dirname,'err',iPda);
figure(105);
xa = 1:iPda;
[ax,h1,h2] = plotyy([xa',xa',xa'],[minf.fs(xa),minf.fg(xa),minf.fp(xa)],xa',minf.fd(xa));
legend([h1;h2],'fs','fg','fp','fd');
clear ps xa;
axis(ax,'square');
xlabel('pda (inner) iterations');
title('parts of min functional');
savePngFigSimple(105,dirname,'minf',iPda);
end
figure(106);
y1 = seti.R(wx+vsp);
y1mat = reshape(y1,[seti.nROI, seti.nROI]);
y1 = y1mat(round(seti.nROI/2),:);
y2 = seti.shrkRe(seti.R(wx+vsp),tau*seti.alpha*seti.dVinv);
y2mat = reshape(y2,[seti.nROI, seti.nROI]);
y2 = y2mat(round(seti.nROI/2),:);
y3 = seti.R(xnpRVD);
y3mat = reshape(y3,[seti.nROI, seti.nROI]);
y3 = y3mat(round(seti.nROI/2),:);
xa = 1:length(y1mat);
plot(xa,y1,xa,y2,xa,y3);
legend('R(wx+vsp) (nROI/2,:)','shrkRe(R(wx+vsp)) (nROI/2,:)','proxGsp (nROI/2,:)');
savePngFigSimple(106,dirname,'proxGsp',iPda);
figure(107);
y1 = seti.R(seti.R(wygVal));
y1mat = reshape(y1,[seti.nROI, seti.nROI]);
y1 = y1mat(round(seti.nROI/2),:);
y2 = seti.R(seti.R(ygnp));
y2mat = reshape(y2,[seti.nROI, seti.nROI]);
y2 = y2mat(round(seti.nROI/2),:);
xa = 1:length(y1mat);
plot(xa,y1,xa,y2);
legend('wygVal (round(nROI/2),:)','proxFgPlus (round(nROI/2),:)');
savePngFigSimple(107,dirname,'wygVal_proxFgPlus',iPda);
v1 = seti.R(seti.Kd(xmnRVD,JA,JB));
v2 = seti.R(seti.R(seti.Kg(xmnRVD)));
figure(108); imagesc(v1); colorbar;
axis xy;
savePngFigSimple(108,dirname,'R_Kd',iPda);
figure(109); imagesc(squeeze(v2)); colorbar;
axis xy;
savePngFigSimple(109,dirname,'R_Kg',iPda);
v1 = seti.R(KdAdjyd);
v2 = seti.R(KgAdjyg);
figure(110); imagesc(real(reshape(v1,[seti.nROI, seti.nROI]))); colorbar;
axis xy;
savePngFigSimple(110,dirname,'R_KdAdj',iPda);
figure(111); imagesc(squeeze(real(reshape(v2,[seti.nROI, seti.nROI])))); colorbar;
axis xy;
savePngFigSimple(111,dirname,'R_KgAdj',iPda);