Commit 84e9c137 authored by Simone Liuzzo's avatar Simone Liuzzo

recompute dispersion and tune derivatives at each iteration

parent 7fed62d8
......@@ -58,10 +58,10 @@ dresp=qemderivAnalytic(qemb.at,dpp,@setk,0.0001,qemres.qpidx(okfit),qemres.bpmid
% % % lines for tests
%drespA=qemderivAnalytic(qemb.at,dpp,@setk,0.0001,qemres.qpidx(okfit),qemres.bpmidx,shidx,svidx,disparg{:});
%drespN= qemderiv(qemb.at,dpp,@setk,0.0001,qemres.qpidx(okfit),qemres.bpmidx,shidx,svidx,disparg{:});
%save('dispderivativecompare','drespA','drespN');
%figure; plot(drespA(:));hold on; plot(drespN(:)); plot(drespA(:)-drespN(:)); xlabel('RM derivative element'); ylabel('RM derivative [m/rad/K_{quad}]'); legend('analytic','qempanel','difference')
% drespA=qemderivAnalytic(qemb.at,dpp,@setk,0.0001,qemres.qpidx(okfit),qemres.bpmidx,shidx,svidx,disparg{:});
% drespN= qemderiv(qemb.at,dpp,@setk,0.0001,qemres.qpidx(okfit),qemres.bpmidx,shidx,svidx,disparg{:});
% save('dispderivativecompare','drespA','drespN');
% figure; plot(drespA(:));hold on; plot(drespN(:)); plot(drespA(:)-drespN(:)); xlabel('RM derivative element'); ylabel('RM derivative [m/rad/K_{quad}]'); legend('analytic','qempanel','difference')
[rh,rh2v,rv2h,rv,frh,frv]=qemdecode(qemb.at,qemres.ct,qemres,...
qemres.resph,semres.respv,semres.resph,qemres.respv,...
......
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