Elegant users forum

Hi,

I got a large DA, with a set of sextupole, but found that during the tracking of "find_aperture" all particles are lost.
Then the (x,y)=(+/-1, +/-1) are the default one ? the Area=2 m^2. are they supposed to be zero ?

Thanks

Here are the ouput:




lyyang@debian:~/devel/nsls2/trunk/optim$ sddsprintout -col lat_0000_offm.aper
Printout for SDDS file lat_0000_offm.aper

x y xClipped yClipped
m m m m
------------------------------------------------------------
1.000000e+00 -1.000000e+00 1.000000e+00 -1.000000e+00
1.000000e+00 -1.000000e+00 1.000000e+00 -1.000000e+00
1.000000e+00 -1.000000e+00 1.000000e+00 -1.000000e+00
1.000000e+00 -1.000000e+00 1.000000e+00 -1.000000e+00
1.000000e+00 -1.000000e+00 1.000000e+00 -1.000000e+00
1.000000e+00 -1.000000e+00 -1.000000e+00 -1.000000e+00
-1.000000e+00 -1.000000e+00 -1.000000e+00 -1.000000e+00
-1.000000e+00 -1.000000e+00 -1.000000e+00 -1.000000e+00
-1.000000e+00 -1.000000e+00 -1.000000e+00 -1.000000e+00
-1.000000e+00 -1.000000e+00 -1.000000e+00 -1.000000e+00
-1.000000e+00 -1.000000e+00 -1.000000e+00 -1.000000e+00
lyyang@debian:~/devel/nsls2/trunk/optim$ sddsprintout -par lat_0000_offm.aper
Printout for SDDS file lat_0000_offm.aper

Step = 1 Area (m$a2$n) = 2.000000e+00



tracking step 1
** Starting 11-line aperture search
* Searching line 0
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.000000e-03, 0.000000e+00
particle lost for x=9.990000e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -5.000000e-04, 0.000000e+00
particle lost for x=9.995000e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -2.500000e-04, 0.000000e+00
particle lost for x=9.997500e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.250000e-04, 0.000000e+00
particle lost for x=9.998750e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -6.250000e-05, 0.000000e+00
particle lost for x=9.999375e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -3.125000e-05, 0.000000e+00
particle lost for x=9.999687e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.562500e-05, 0.000000e+00
particle lost for x=9.999844e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -7.812500e-06, 0.000000e+00
particle lost for x=9.999922e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -3.906250e-06, 0.000000e+00
particle lost for x=9.999961e-01, y=-1.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
* Searching line 1
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -9.510565e-04, 4.635255e-05
particle lost for x=9.990489e-01, y=-9.999536e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -4.755283e-04, 2.317627e-05
particle lost for x=9.995245e-01, y=-9.999768e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -2.377641e-04, 1.158814e-05
particle lost for x=9.997622e-01, y=-9.999884e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.188821e-04, 5.794069e-06
particle lost for x=9.998811e-01, y=-9.999942e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -5.944103e-05, 2.897034e-06
particle lost for x=9.999406e-01, y=-9.999971e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -2.972052e-05, 1.448517e-06
particle lost for x=9.999703e-01, y=-9.999986e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.486026e-05, 7.242586e-07
particle lost for x=9.999851e-01, y=-9.999993e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -7.430129e-06, 3.621293e-07
particle lost for x=9.999926e-01, y=-9.999996e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -3.715065e-06, 1.810646e-07
particle lost for x=9.999963e-01, y=-9.999998e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
* Searching line 2
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -8.090170e-04, 8.816779e-05
particle lost for x=9.991910e-01, y=-9.999118e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -4.045085e-04, 4.408389e-05
particle lost for x=9.995955e-01, y=-9.999559e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -2.022542e-04, 2.204195e-05
particle lost for x=9.997977e-01, y=-9.999780e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.011271e-04, 1.102097e-05
particle lost for x=9.998989e-01, y=-9.999890e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -5.056356e-05, 5.510487e-06
particle lost for x=9.999494e-01, y=-9.999945e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -2.528178e-05, 2.755243e-06
particle lost for x=9.999747e-01, y=-9.999972e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.264089e-05, 1.377622e-06
particle lost for x=9.999874e-01, y=-9.999986e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -6.320445e-06, 6.888108e-07
particle lost for x=9.999937e-01, y=-9.999993e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -3.160223e-06, 3.444054e-07
particle lost for x=9.999968e-01, y=-9.999997e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
* Searching line 3
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -5.877853e-04, 1.213525e-04
particle lost for x=9.994122e-01, y=-9.998786e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -2.938926e-04, 6.067627e-05
particle lost for x=9.997061e-01, y=-9.999393e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.469463e-04, 3.033814e-05
particle lost for x=9.998531e-01, y=-9.999697e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -7.347316e-05, 1.516907e-05
particle lost for x=9.999265e-01, y=-9.999848e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -3.673658e-05, 7.584534e-06
particle lost for x=9.999633e-01, y=-9.999924e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.836829e-05, 3.792267e-06
particle lost for x=9.999816e-01, y=-9.999962e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -9.184145e-06, 1.896134e-06
particle lost for x=9.999908e-01, y=-9.999981e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -4.592072e-06, 9.480668e-07
particle lost for x=9.999954e-01, y=-9.999991e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -2.296036e-06, 4.740334e-07
particle lost for x=9.999977e-01, y=-9.999995e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
* Searching line 4
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -3.090170e-04, 1.426585e-04
particle lost for x=9.996910e-01, y=-9.998573e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.545085e-04, 7.132924e-05
particle lost for x=9.998455e-01, y=-9.999287e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -7.725425e-05, 3.566462e-05
particle lost for x=9.999227e-01, y=-9.999643e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -3.862712e-05, 1.783231e-05
particle lost for x=9.999614e-01, y=-9.999822e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.931356e-05, 8.916155e-06
particle lost for x=9.999807e-01, y=-9.999911e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -9.656781e-06, 4.458077e-06
particle lost for x=9.999903e-01, y=-9.999955e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -4.828391e-06, 2.229039e-06
particle lost for x=9.999952e-01, y=-9.999978e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -2.414195e-06, 1.114519e-06
particle lost for x=9.999976e-01, y=-9.999989e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to -1.207098e-06, 5.572597e-07
particle lost for x=9.999988e-01, y=-9.999994e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
* Searching line 5
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 1.500000e-04
particle lost for x=1.000000e+00, y=-9.998500e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 7.500000e-05
particle lost for x=1.000000e+00, y=-9.999250e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 3.750000e-05
particle lost for x=1.000000e+00, y=-9.999625e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 1.875000e-05
particle lost for x=1.000000e+00, y=-9.999813e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 9.375000e-06
particle lost for x=1.000000e+00, y=-9.999906e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 4.687500e-06
particle lost for x=1.000000e+00, y=-9.999953e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 2.343750e-06
particle lost for x=1.000000e+00, y=-9.999977e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 1.171875e-06
particle lost for x=1.000000e+00, y=-9.999988e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
divided search interval to 0.000000e+00, 5.859375e-07
particle lost for x=1.000000e+00, y=-9.999994e-01
Sweep done, particle survived up to x=1.000000e+00, y=-1.000000e+00
* Searching line 6
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 3.090170e-04, 1.426585e-04
particle lost for x=-9.996910e-01, y=-9.998573e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.545085e-04, 7.132924e-05
particle lost for x=-9.998455e-01, y=-9.999287e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 7.725425e-05, 3.566462e-05
particle lost for x=-9.999227e-01, y=-9.999643e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 3.862712e-05, 1.783231e-05
particle lost for x=-9.999614e-01, y=-9.999822e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.931356e-05, 8.916155e-06
particle lost for x=-9.999807e-01, y=-9.999911e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 9.656781e-06, 4.458077e-06
particle lost for x=-9.999903e-01, y=-9.999955e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 4.828391e-06, 2.229039e-06
particle lost for x=-9.999952e-01, y=-9.999978e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 2.414195e-06, 1.114519e-06
particle lost for x=-9.999976e-01, y=-9.999989e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.207098e-06, 5.572597e-07
particle lost for x=-9.999988e-01, y=-9.999994e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
* Searching line 7
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 5.877853e-04, 1.213525e-04
particle lost for x=-9.994122e-01, y=-9.998786e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 2.938926e-04, 6.067627e-05
particle lost for x=-9.997061e-01, y=-9.999393e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.469463e-04, 3.033814e-05
particle lost for x=-9.998531e-01, y=-9.999697e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 7.347316e-05, 1.516907e-05
particle lost for x=-9.999265e-01, y=-9.999848e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 3.673658e-05, 7.584534e-06
particle lost for x=-9.999633e-01, y=-9.999924e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.836829e-05, 3.792267e-06
particle lost for x=-9.999816e-01, y=-9.999962e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 9.184145e-06, 1.896134e-06
particle lost for x=-9.999908e-01, y=-9.999981e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 4.592072e-06, 9.480668e-07
particle lost for x=-9.999954e-01, y=-9.999991e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 2.296036e-06, 4.740334e-07
particle lost for x=-9.999977e-01, y=-9.999995e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
* Searching line 8
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 8.090170e-04, 8.816779e-05
particle lost for x=-9.991910e-01, y=-9.999118e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 4.045085e-04, 4.408389e-05
particle lost for x=-9.995955e-01, y=-9.999559e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 2.022542e-04, 2.204195e-05
particle lost for x=-9.997977e-01, y=-9.999780e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.011271e-04, 1.102097e-05
particle lost for x=-9.998989e-01, y=-9.999890e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 5.056356e-05, 5.510487e-06
particle lost for x=-9.999494e-01, y=-9.999945e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 2.528178e-05, 2.755243e-06
particle lost for x=-9.999747e-01, y=-9.999972e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.264089e-05, 1.377622e-06
particle lost for x=-9.999874e-01, y=-9.999986e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 6.320445e-06, 6.888108e-07
particle lost for x=-9.999937e-01, y=-9.999993e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 3.160223e-06, 3.444054e-07
particle lost for x=-9.999968e-01, y=-9.999997e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
* Searching line 9
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 9.510565e-04, 4.635255e-05
particle lost for x=-9.990489e-01, y=-9.999536e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 4.755283e-04, 2.317627e-05
particle lost for x=-9.995245e-01, y=-9.999768e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 2.377641e-04, 1.158814e-05
particle lost for x=-9.997622e-01, y=-9.999884e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.188821e-04, 5.794069e-06
particle lost for x=-9.998811e-01, y=-9.999942e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 5.944103e-05, 2.897034e-06
particle lost for x=-9.999406e-01, y=-9.999971e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 2.972052e-05, 1.448517e-06
particle lost for x=-9.999703e-01, y=-9.999986e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.486026e-05, 7.242586e-07
particle lost for x=-9.999851e-01, y=-9.999993e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 7.430129e-06, 3.621293e-07
particle lost for x=-9.999926e-01, y=-9.999996e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 3.715065e-06, 1.810646e-07
particle lost for x=-9.999963e-01, y=-9.999998e-01
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
* Searching line 10
particle lost for x=0.000000e+00, y=0.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.000000e-03, 0.000000e+00
particle lost for x=-9.990000e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 5.000000e-04, 0.000000e+00
particle lost for x=-9.995000e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 2.500000e-04, 0.000000e+00
particle lost for x=-9.997500e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.250000e-04, 0.000000e+00
particle lost for x=-9.998750e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 6.250000e-05, 0.000000e+00
particle lost for x=-9.999375e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 3.125000e-05, 0.000000e+00
particle lost for x=-9.999687e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 1.562500e-05, 0.000000e+00
particle lost for x=-9.999844e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 7.812500e-06, 0.000000e+00
particle lost for x=-9.999922e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
divided search interval to 3.906250e-06, 0.000000e+00
particle lost for x=-9.999961e-01, y=-1.000000e+00
Sweep done, particle survived up to x=-1.000000e+00, y=-1.000000e+00
Finished all tracking steps.
Finished dynamic aperture search.
End of input data encountered.
statistics: ET: 00:00:01 CP: 1.20 BIO:0 DIO:0 PF:0 MEM:2365
=====================================================================================
Thanks for using elegant. Please cite the following reference in your publications:
M. Borland, "elegant: A Flexible SDDS-Compliant Code for Accelerator Simulation,"
Advanced Photon Source LS-287, September 2000.
If you use a modified version, please indicate this in all publications.
=====================================================================================

I think this is a bug. I never prepared the algorithm to deal with all particles being lost.

If you can post the input files (or email them to me), it would make it easier for me to understand what's happening.

--Michael

Hi Michael,

Thanks for looking at this problem.

Here are the lattice file and kickmap.


Lingyun

Lingyun,

It is indeed a bug. The attached file contains the fix, which you can use if you build from source.

Otherwise, a workaround is to postprocess the file, which can be done by adding this command at the end of the input file
&subprocess command = "
sddsprocess %s.aper -redefine=parameter,Area,\"Area 2 == ? 0 : Area $ \"
" &end

--Michael

Michael,

Thanks for this quick fix, I didn't build from source code, so the post process would be a best solution for me now.

Lingyun