Format of p_central_mev(run_setup) affects calculation of response matrices

Moderators: michael_borland, soliday

Post Reply
Meykopff
Posts: 3
Joined: 19 Oct 2017, 01:56
Location: DESY, Hamburg

Format of p_central_mev(run_setup) affects calculation of response matrices

Post by Meykopff » 19 Oct 2017, 03:30

Dear all,
at the European XFEL we provide a server based online model of our machine. For the calculation of the linear optics we use elegant. Currently i found a problem with the calculation of the response matrices. If i use 4.99999 MeV as start parameter the calculation works. But if i use 5.0 the calculation fails. The calculation of the transport matrices or other data is not affected. I'm using elegant 33.0.3 on Linux. See the attachment for input files.
Output of *.vrm (last line):
BPMD.479.B2D 4.564451674999980e+02 1.033521988232392e+01 9.490213451379422e+00 8.951785317392810e+00 2.747407033117273e+00 -7.168767277097474e+01 -5.683774308147104e+01 -5.527164157547934e+01 -5.699850022820646e+01 ...

BPMD.479.B2D 4.564451674999980e+02 -5.111274173440000e+17 -7.737148194800000e+16 1.447186591320000e+17 5.050675530920000e+17 9.606133432540000e+17 -7.686501400580000e+17 1.660610899320000e+17 9.603252500800000e+16 ...
Attachments
B2D_b.txt
(695 Bytes) Downloaded 139 times
B2D_a.txt
(700 Bytes) Downloaded 126 times
B2D.lte
(113.06 KiB) Downloaded 166 times

michael_borland
Posts: 1927
Joined: 19 May 2008, 09:33
Location: Argonne National Laboratory
Contact:

Re: Format of p_central_mev(run_setup) affects calculation of response matrices

Post by michael_borland » 23 Oct 2017, 09:29

Sorry for the delayed reply; I was away from work for a while.

There seem to be several causes for this problem. First, you shouldn't have "fixed_length=1" for a transport line or linac. More importantly, the correction matrix is ill-conditioned, probably because there are several correctors between some BPMs. You can see the problem by removing the method = "one-to-one" setting. The condition number is >1e17, whereas <~1e3 is recommended. I suggest using

Code: Select all

&correct
        mode = trajectory
        keep_largest_SVs[0] = 2*50
        corrector_tweek[0] = 2*1e-5
        correction_fraction[0] = 2*0.1,
&end
The fact that it worked for some values of p_central_mev appears to be pure luck.

--Michael

Post Reply