Commit 8445d8ea authored by Simone Liuzzo's avatar Simone Liuzzo

never zero strenghts (1e-12)

parent 3d23d9f7
......@@ -53,6 +53,9 @@ nbpm=length(indbpm);
% dipoles angle and length
K_quad = cellfun(@(a)a.PolynomB(2),r(indquad),'un',1)';
K_quad(K_quad==0) = 1e-12; % deal with K=0 quadrupoles
len_quad = cellfun(@(a)a.Length,r(indquad),'un',1)';
% beta x at bpm and quadrupole
......
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