Commit 6090cd30 authored by Simone Liuzzo's avatar Simone Liuzzo

compact and large version of correctors display

parent e1740774
......@@ -3,9 +3,9 @@ tic
if nargin < 5, handles=struct(); end
if isfield(handles,'axes1')
qcheck(qemb.kn,qemb0.kn,handles.axes1);
qcheck(qemb.kn,qemb0.kn,handles.axes1,get(handles.largefigures,'Value'));
else
qcheck(qemb.kn,qemb0.kn);
qcheck(qemb.kn,qemb0.kn,[],get(handles.largefigures,'Value'));
end
toc;
qemb.at=qemat(qemres,qemb,true); % builds the at structure
......@@ -167,9 +167,23 @@ else
end
if isfield(handles,'axes4') % displays corr. strengths
bar(handles.axes4,qemb.cor);
title(handles.axes4,'Corrector strengths');
grid(handles.axes4,'on');
if get(handles.largefigures,'Value')
figure; extfigax = gca;
bar(extfigax,qemb.cor);
legend(extfigax,'initial','corrected');
title(extfigax,'Corrector strengths');
grid(extfigax,'on');
else
Nc = length(qemb.cor);
extracor = mod(Nc,32); % remove extra correctors from injection cells
cor_per_cell_std = squeeze(std(reshape(qemb.cor((1+extracor/2):(Nc-extracor/2),:),32,[],2),1,2));
% cor_per_cell_mean = squeeze(mean(reshape(qemb.cor((1+extracor/2):(Nc-extracor/2),:),32,[],2),2));
bar(handles.axes4,cor_per_cell_std,'EdgeColor','none');
title(handles.axes4,'std Corrector strengths (per cell)');
legend(handles.axes4,'initial','corrected');
grid(handles.axes4,'on');
end
end
% not working, wrong handle. write below code to enable zoom in out in GUI
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment