source: PSPA/madxPSPA/tests/test-match-3/test-match-3.ref @ 478

Last change on this file since 478 was 430, checked in by touze, 11 years ago

import madx-5.01.00

File size: 8.1 KB
Line 
1
2  +++++++++++++++++++++++++++++++++++++++++++
3  +    MAD-X 5.00.19  (64 bit, Darwin)      +
4  +    DEBUG Version - use at own risk!     +
5  + Code Modification Date: 2012.12.11      +
6  + Execution Time Stamp: 13.02.13 21.52.09 +
7  +++++++++++++++++++++++++++++++++++++++++++
8TITLE,'Transferline matching with MAD-X';
9
10
11
12option,-info,-warn;
13
14
15
16call file="fv9.opt"; 
17
18// forces and other constants;
19
20
21
22const l.bpm = .3;
23
24bang = .509998807401e-2;
25
26kqf = .872651312e-2;
27
28kqd = -.872777242e-2;
29
30//ksf = .0198492943;
31
32//ksd = -.039621283;
33
34acbv1 = 1.e-4;
35
36dummy := ksf * ksd;
37
38// element definitions;
39
40
41
42mb: rbend, l=14.2, angle:=bang;
43
44mq: quadrupole, l=3.1;
45
46m0:multipole, knl:={0},ksl:={0.0,0.0};
47
48//mm:multipole,lrad=dummy,kn={0,dummy,bang*3,1},ks={1.1,3*acbv1+5.,4,0};
49
50qf.1: mq, k1:=kqf;
51
52qf.2: mq, k1:=kqf;
53
54qf.3: mq, k1:=kqf;
55
56qf.4: mq, k1:=kqf;
57
58qf.5: mq, k1:=kqf;
59
60qd.1: mq, k1:=kqd;
61
62qd.2: mq, k1:=kqd;
63
64qd.3: mq, k1:=kqd;
65
66qd.4: mq, k1:=kqd;
67
68qd.5: mq, k1:=kqd;
69
70bpv: vmonitor, l:=l.bpm;
71
72cbh: hkicker;
73
74cbh.1:cbh,kick:=acbh1;
75
76cbh.2:cbh,kick:=acbh2;
77
78cbh.3:cbh,kick:=acbh3;
79
80cbh.4:cbh,kick:=acbh4;
81
82cbh.5:cbh,kick:=acbh5;
83
84cbv: vkicker;
85
86cbv.1:cbv,kick:=acbv1;
87
88cbv.2:cbv,kick:=acbv2;
89
90cbv.3:cbv,kick:=acbv3;
91
92cbv.4:cbv,kick:=acbv4;
93
94cbv.5:cbv,kick:=acbv5;
95
96mscbh: sextupole, l=1.1, k2:=ksf;
97
98mscbv: sextupole, l=1.1, k2:=ksd;
99
100//bph: hmonitor, l:=l.bpm;
101
102bph: rfcavity,l:=l.bpm,
103
104 volt=.070833333333, harmon=432, shunt=26.166, tfill=90.0;
105
106call file="fv9.seq";
107
108lfive=534.6;
109
110dumzero = 0;
111
112// sequence declaration;
113
114shared ins: sequence, l=5.615,refpos=qd.1;
115
116   bpv.1: bpv, at=5.246000e+01-5.501000e+01,from=qd.1;
117
118   qd.1: qd.1, at=2.8;
119
120   mscbv.1: mscbv, at=5.727500e+01-5.501000e+01,from=qd.1;
121
122endsequence;
123
124
125
126pos = 1.084700e+02-1.550000e+00;
127
128fivecell: sequence, refer=centre, l=lfive;
129
130   seqstart: marker, at = 0;
131
132   qf.1: qf.1, at=dumzero+1.550000e+00, from=seqstart;
133
134   mscbh: mscbh, at=3.815000e+00;
135
136   cbh.1: cbh.1, at=4.365000e+00;
137
138   mb: mb, at=1.262000e+01;
139
140   mb: mb, at=2.828000e+01;
141
142   mb: mb, at=4.394000e+01;
143
144   ins, at=5.501000e+01;
145
146   cbv.1: cbv.1, at=5.782500e+01;
147
148   mb: mb, at=6.608000e+01;
149
150   mb: mb, at=8.174000e+01;
151
152   mb: mb, at=9.740000e+01;
153
154   m0, at=105;
155
156//   mm, at=105;
157
158   bph: bph, at=1.059200e+02, from=seqstart;
159
160   startend: marker, at = pos;
161
162   partstart: marker, at = pos;
163
164   qf.2: qf.2, at=dumzero+1.084700e+02;
165
166   mscbh: mscbh, at=1.107350e+02;
167
168   cbh.2: cbh.2, at=1.112850e+02;
169
170   mb: mb, at=1.195400e+02;
171
172   mb: mb, at=1.352000e+02;
173
174   mb: mb, at=1.508600e+02;
175
176   bpv: bpv, at=1.593800e+02;
177
178   qd.2: qd.2, at=1.619300e+02;
179
180   mscbv: mscbv, at=1.641950e+02;
181
182   cbv.2: cbv.2, at=1.647450e+02;
183
184   mb: mb, at=1.730000e+02;
185
186   mb: mb, at=1.886600e+02;
187
188   mb: mb, at=2.043200e+02;
189
190   bph: bph, at=2.128400e+02;
191
192   partend: marker, at = 213.84;
193
194   qf.3: qf.3, at=2.153900e+02;
195
196   mscbh: mscbh, at=2.176550e+02;
197
198   cbh.3: cbh.3, at=2.182050e+02;
199
200   mb: mb, at=2.264600e+02;
201
202   mb: mb, at=2.421200e+02;
203
204   mb: mb, at=2.577800e+02;
205
206   bpv: bpv, at=2.663000e+02;
207
208   qd.3: qd.3, at=2.688500e+02;
209
210   mscbv: mscbv, at=2.711150e+02;
211
212   cbv.3: cbv.3, at=2.716650e+02;
213
214   mb: mb, at=2.799200e+02;
215
216   mb: mb, at=2.955800e+02;
217
218   mb: mb, at=3.112400e+02;
219
220   bph: bph, at=3.197600e+02;
221
222   qf.4: qf.4, at=3.223100e+02;
223
224   mscbh: mscbh, at=3.245750e+02;
225
226   cbh.4: cbh.4, at=3.251250e+02;
227
228   mb: mb, at=3.333800e+02;
229
230   mb: mb, at=3.490400e+02;
231
232   mb: mb, at=3.647000e+02;
233
234   bpv: bpv, at=3.732200e+02;
235
236   qd.4: qd.4, at=3.757700e+02;
237
238   mscbv: mscbv, at=3.780350e+02;
239
240   cbv.4: cbv.4, at=3.785850e+02;
241
242   mb: mb, at=3.868400e+02;
243
244   mb: mb, at=4.025000e+02;
245
246   mb: mb, at=4.181600e+02;
247
248   bph: bph, at=4.266800e+02;
249
250   qf.5: qf.5, at=4.292300e+02;
251
252   mscbh: mscbh, at=4.314950e+02;
253
254   cbh.5: cbh.5, at=4.320450e+02;
255
256   mb: mb, at=4.403000e+02;
257
258   mb: mb, at=4.559600e+02;
259
260   mb: mb, at=4.716200e+02;
261
262   bpv: bpv, at=4.801400e+02;
263
264   qd.5: qd.5, at=4.826900e+02;
265
266   mscbv: mscbv, at=4.849550e+02;
267
268   cbv.5: cbv.5, at=4.855050e+02;
269
270   mb: mb, at=4.937600e+02;
271
272   mb: mb, at=5.094200e+02;
273
274   mb: mb, at=5.250800e+02;
275
276   bph: bph, at=5.336000e+02;
277
278end_fivecell: marker at=5.346000e+02;
279
280endsequence;
281
282
283
284option,info,warn;
285
286
287
288kqf = 0.0089;
289
290++++++ info: kqf redefined
291kqd = -0.0097;
292
293++++++ info: kqd redefined
294
295
296Beam, particle = proton, sequence=fivecell, energy = 450.0,
297
298          NPART=1.05E11, sige=      4.5e-4 ;
299
300
301
302use,period=fivecell;
303
304
305
306savebeta,label=abc,place=#s;
307
308select,flag=twiss,range=full,column=name,s,betx,bety;
309
310twiss, sequence=fivecell;
311
312enter Twiss module
313 
314iteration:   1 error:   6.134314E-03 deltap:   0.000000E+00
315orbit:   3.446033E-07  4.759082E-09  2.948855E-03 -1.098131E-04 -6.889332E-08  1.519132E-07
316 
317iteration:   2 error:   2.757735E-06 deltap:   0.000000E+00
318orbit:  -2.556164E-05 -3.552748E-07  2.948751E-03 -1.098197E-04 -2.419912E-07 -1.142313E-05
319 
320iteration:   3 error:   2.748594E-10 deltap:   0.000000E+00
321orbit:  -2.556193E-05 -3.552790E-07  2.948751E-03 -1.098197E-04 -2.419989E-07 -1.142338E-05
322
323++++++ table: summ
324
325            length             orbit5               alfa            gammatr
326             534.6    2.419989205e-07    0.0004498974988        47.14582183
327
328                q1                dq1            betxmax              dxmax
329       1.256622886       -1.709613092        186.0593755        2.235741425
330
331             dxrms             xcomax             xcorms                 q2
332       1.711709766     2.55811844e-05    1.957039145e-05        1.468423245
333
334               dq2            betymax              dymax              dyrms
335      -1.988647234        185.2849056      0.05539960915      0.02099810508
336
337            ycomax             ycorms             deltap            synch_1
338    0.008577862743     0.004467121568                  0                  0
339
340           synch_2            synch_3            synch_4            synch_5
341                 0                  0                  0                  0
342
343
344match,sequence=fivecell,beta0=abc;
345
346START MATCHING
347
348number of sequences: 1
349sequence name: fivecell
350number of beta0s: 1
351BETA0 name: abc
352constraint,sequence=fivecell,range=#e,betx=180.0,bety=30.0;
353
354vary,name=kqf,step=1.0e-6;
355
356vary,name=kqd,step=1.0e-6;
357
358lmdif,calls=500,tolerance=1.e-15;
359
360number of variables:    2
361user given constraints: 2
362total constraints:      2
363
364START LMDIF:
365
366Initial Penalty Function =   0.69516383E+02
367
368
369call:       7   Penalty function =   0.68759290E+02
370call:      10   Penalty function =   0.67426945E+02
371call:      13   Penalty function =   0.65332869E+02
372call:      16   Penalty function =   0.61722731E+02
373call:      19   Penalty function =   0.55331866E+02
374call:      22   Penalty function =   0.43829835E+02
375call:      25   Penalty function =   0.33718015E+02
376call:      28   Penalty function =   0.12560706E+02
377call:      31   Penalty function =   0.26435430E+00
378call:      34   Penalty function =   0.73332041E-04
379call:      37   Penalty function =   0.84935053E-09
380call:      40   Penalty function =   0.38717782E-14
381call:      43   Penalty function =   0.56083230E-19
382 ++++++++++ LMDIF ended: converged successfully
383call:      43   Penalty function =   0.56083230E-19
384endmatch;
385
386
387MATCH SUMMARY
388
389Node_Name                  Constraint   Type  Target Value       Final Value        Penalty
390--------------------------------------------------------------------------------------------------
391fivecell$end:1             betx         4     1.80000000E+02     1.80000000E+02     5.60249946E-20
392fivecell$end:1             bety         4     3.00000000E+01     3.00000000E+01     5.82356558E-23
393
394
395Final Penalty Function =   5.60832302e-20
396
397
398
399
400
401Variable                 Final Value  Initial Value Lower Limit  Upper Limit
402--------------------------------------------------------------------------------
403kqf                       8.72106e-03  8.90000e-03 -1.00000e+20  1.00000e+20
404kqd                      -9.19078e-03 -9.70000e-03 -1.00000e+20  1.00000e+20
405
406END MATCH SUMMARY
407
408VARIABLE "TAR" SET TO   5.60832302e-20
409
410
411value, kqf;
412
413kqf =     0.008721063916 ;
414value, kqd;
415
416kqd =    -0.009190778403 ;
417
418  Number of warnings: 0
419
420  ++++++++++++++++++++++++++++++++++++++++++++
421  + MAD-X 5.00.19 (64 bit) finished normally +
422  ++++++++++++++++++++++++++++++++++++++++++++
Note: See TracBrowser for help on using the repository browser.