Optimisation on an R-matrix element

[astecpete]

Hi Michael,

Hope you had a good xmas and new year!

I think i *may* have found another bug. I want to optimise for R56 = 0 at a particular marker. When i try i get no errors but the optimiser reports that the constraint is satisfied e.g.
Terms of equation:
AR1-A2-MARK#1.R56 0 0.000000001 sene: 0.0000000E+0

But looking in the .matr file at that point i see that R56 is around 3E-03.
Or am i daft?

Pete

[astecpete]

Nope - i was daft! I just forgot to regenerate the .matr file. D'oh

[simone.dimitri]

Hi all,
I am working with transfer matrix elements. Let me focus on R56 for a symmetric, 4 dipoles chicane (no active quads inside).
I have a chicane called LH followed by another one called BC1.
R56 is calculated correctly for each chicane alone, while numbers are different for the last chicane when I put LH switched off and then BC1, on a same beam line. Please look at tables below.
I guess it depends on the total path length that is different in the two cases but I am not sure. Am I missing something?
Thanks,
S.


BC1 alone
s ElementName R56 T566
m m/1 m
----------------------------------------------------------
0.000000e+000 _BEG_ 0.000000e+000 0.000000e+000
5.000000e-001 B_BC1.1 6.822414e-004 -6.263285e-016
8.000000e-001 DR_BC1.1A 6.822414e-004 1.235277e-003
9.000000e-001 Q_BC1.1 6.822414e-004 1.647036e-003
3.000000e+000 DR_BC1.1B 6.822414e-004 1.029397e-002
3.500000e+000 B_BC1.2 -2.332700e-002 3.526781e-002
3.600000e+000 DR_BC1.2A -2.332700e-002 3.526781e-002
3.693000e+000 SCRPH_BC1.1 -2.332700e-002 3.526781e-002
3.813000e+000 DR_BC1.2B -2.332700e-002 3.526781e-002
3.913000e+000 DR_BC1.2C -2.332700e-002 3.526781e-002
4.113000e+000 BPM_BC1.1 -2.332700e-002 3.526781e-002
4.295000e+000 DR_BC1.2D -2.332700e-002 3.526781e-002
4.375000e+000 MSCR_BC1.1 -2.332700e-002 3.526781e-002
4.500000e+000 DR_BC1.2E -2.332700e-002 3.526781e-002
5.000000e+000 B_BC1.3 -4.733625e-002 6.023657e-002
7.100000e+000 DR_BC1.3A -4.733625e-002 6.888351e-002
7.200000e+000 Q_BC1.2 -4.733625e-002 6.929527e-002
7.500000e+000 DR_BC1.3B -4.733625e-002 7.053054e-002
8.000000e+000 B_BC1.4 -4.665401e-002 7.055605e-002


LH off followed by BC1

3.601240e+001 DR_L1.SIS15 0.000000e+000 0.000000e+000
3.631440e+001 DR_L1.14 0.000000e+000 0.000000e+000
3.681440e+001 B_BC1.1 2.645912e-004 3.567194e-009
3.711440e+001 DR_BC1.1A 2.645912e-004 1.858006e-004
3.721440e+001 Q_BC1.1 2.645912e-004 2.477330e-004
3.931440e+001 DR_BC1.1B 2.645912e-004 1.548312e-003
3.981440e+001 B_BC1.2 -9.046826e-003 5.304482e-003
3.991440e+001 DR_BC1.2A -9.046826e-003 5.304482e-003
4.000740e+001 SCRPH_BC1.1 -9.046826e-003 5.304482e-003
4.012740e+001 DR_BC1.2B -9.046826e-003 5.304482e-003
4.022740e+001 DR_BC1.2C -9.046826e-003 5.304482e-003
4.042740e+001 BPM_BC1.1 -9.046826e-003 5.304482e-003
4.060940e+001 DR_BC1.2D -9.046826e-003 5.304482e-003
4.068940e+001 MSCR_BC1.1 -9.046826e-003 5.304482e-003
4.081440e+001 DR_BC1.2E -9.046826e-003 5.304482e-003
4.131440e+001 B_BC1.3 -1.835824e-002 9.059890e-003
4.341440e+001 DR_BC1.3A -1.835824e-002 1.036047e-002
4.351440e+001 Q_BC1.2 -1.835824e-002 1.042240e-002
4.381440e+001 DR_BC1.3B -1.835824e-002 1.060820e-002
4.431440e+001 B_BC1.4 -1.809365e-002 1.061204e-002

[michael_borland]

Simone,

Can you post the lattice file? It will be easier to diagnose that way.

Do you have accelerating structures between the chicanes? If so, note that R56 for the full system refers to the initial energy.

--Michael

[simone.dimitri]

Yes Michael, I have accelerating sections.
I have got it. Now it is OK.
Thank you,
S.