Commit 177a704b authored by Simone Liuzzo's avatar Simone Liuzzo

potentially different number of h and v steerers

parent 2fe1e7f0
......@@ -38,8 +38,18 @@ v=mode.vnorm;
%w=[ones(orbitrange,1);ones(nbpm,1);1000;1000];
%sk=qemsolvex(dresp(:,qpfit).*(w*v(qpfit)'),resp.*w,bok,4,100);
w=[ones(orbitrange,1);mode.dispweight*ones(nbpm,1);mode.tuneweight*ones(2,1)];
hsets=reshape(1:nhst,mode.nsets,[]);
vsets=reshape(1:nvst,mode.nsets,[]);
try
hsets=reshape(1:nhst,mode.nsets,[]);
vsets=reshape(1:nvst,mode.nsets,[]);
catch
warning('Error dividing in subsets')
disp([ 'NHst/Nsets = ' num2str(nhst) ' / ' num2str(mode.nsets) ' = ' num2str(nhst/mode.nsets) '<-should be integer'])
disp([ 'NVst/Nsets = ' num2str(nvst) ' / ' num2str(mode.nsets) ' = ' num2str(nvst/mode.nsets) '<-should be integer'])
disp('setting nsets = 1');
mode.nsets = 1;
hsets=reshape(1:nhst,mode.nsets,[]);
vsets=reshape(1:nvst,mode.nsets,[]);
end
sk=qemsolvex(dresp.*(w*v),resp.*w,bok,hsets,vsets,mode.neigs);
deltaset=mean(sk,2).*v';
devdev=std(sk,1,2).*v';
......
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