Subject | F bench2.sad on iMac late 2012 (3.4 GHz i7). |
Article No | F 1095 |
Date | F 2012/12/29(Sat) 06:34:02 |
Contributor | F K. Oide |
*** Welcome to SAD Ver.1.0.10.6.5.2k64 built at 2012-12-29 06:27:48 +0900 *** *** Today: 06:32:25 Saturday 12/29/2012 *** cpu time= 5.4680E-03(sec) dt= 5.468(msec) free area:: 1794 OFF LOG ECHO;READ 77 ; 23 cpu time= 5.5600E-03(sec) dt= 5.560(msec) free area:: 1794 cpu time= 5.5920E-03(sec) dt= 5.592(msec) free area:: 1793 *** SADScript Initialization: /Users/oide/SAD/oldsad/Packages/init.n *** RFSW RADCOD RAD FLUC INTRA POL COD DAPER EMIOU CMPLO FOURI SMEAR F F F T F F T F F F F T Design orbit length = 76.180924396860036 NPARA=1; EMFLG='M'; FLENG='ON'; FSUPP='OFF'; FINIT='OFF'; cell cal Matched. ( 0.000 ) DP = 0.02000 DP0 = 0.00000 ExponentOfResidual = 2.0 OffMomentumWeight = 1.000 $$$ f AX ####### # 1.554E-15 $$$ f BX ####### # 6.981218 $$$ f NX ####### # 1.250000 $$$ f AY ####### # -4.44E-15 $$$ f BY ####### # 27.858289 $$$ f NY ####### # 1.250000 $$$ f LENG ####### # 76.180924 vquad=Element['VALUE','Q*'];nquad=Element['NAME','Q*']; ! setquad[name_,val_]:=Module[{},Element['VALUE',name]=val]; ! MapThread[Element['VALUE',#1 ]=#2 &,{nquad,vquad}] ! Do[Element['VALUE',nquad[ii]]=vquad[ii],{ii,Length[nquad]}] fit SD1.2 nx .25 ny .25 fit SF1.2 SF1.4 nx 0.5 ny 0.5 fit B.2 ex .183 fit ex -0.065 fit leng 76.1809 nx 1.25 ny 1.25 fit qf1.2 qf3.2 bx 10 fit qf1.2 qf2.2 bx 10 fit qd1 qd2.2 by 10 fit qd1 qd3.2 by 10 b 2*Pi/112 b l .86 l7 .52 free q* l0 l1 l2 l3 l4 l5 l6 fix qe2 CONVERGENCE=1e-28; IF EMFLG=='H' fix * free qe2 qd2 qd3 qf2 qf3 fit rej leng nx ny;fit sd1.2 rej * ;fit qf1.2 qf2.2 rej bx;fit qf1.2 qf3.2 rej bx; fit qd1 qd2.2 rej by fit qd1 qd3.2 rej by;fit qd1 qd2.2 by 10; fit ex -.30 fit b.2 ex 0.35 ENDIF IF EMFLG=='L' fix * free qe2 qd2 qd3 qf2 qf3 fit rej leng nx ny;fit sd1.2 rej * ;fit qf1.2 qf2.2 rej bx;fit qf1.2 qf3.2 rej bx; fit qd1 qd2.2 rej by fit qd1 qd3.2 rej by;fit qd2.2 qd3.2 by 10; fit ex 0.11 fit b.2 ex 0.05 ENDIF ! go org nmark 0 0 0 360/112*2 0 0 ; calc; Residual = 6.5569E-12 DP = 0.02000 DP0 = 0.00000 ExponentOfResidual = 2.0 OffMomentumWeight = 1.000 QD1.1/QD2.2 BY 27.85829 1 27.858289 QD1.1/QD3.f BY 27.85829 1 27.858289 QF1.2/QF2.2 BX 26.11497 1 26.114975 QF1.2/QF3.2 BX 26.11497 1 26.114975 SD1.2 NX .25 1 .250000 SD1.2 NY .25 1 .250000 SF1.2/SF1.4 NX .5 1 .500000 SF1.2/SF1.4 NY .5 1 .500000 B.2 EX .183 1 .183000 $$$ AX ####### # 2.109E-15 $$$ BX ####### # 6.981218 $$$ NX 1.25 1 1.250000 $$$ AY ####### # -2.66E-15 $$$ BY ####### # 27.858289 $$$ NY 1.25 1 1.250000 $$$ EX -.065 1 -.065000 $$$ LENG 76.1809 1 76.180924 xxx=LINE['GEO','***'][1,1]/2/Tan[Pi*2/112*2]-LINE['GEO','QF3.2'][1,2]; yyy=LINE['GEO','***'][1,1]/2/Sin[Pi*2/112*2]; delr=(xxx+yyy)/2-yyy Out[33]:= .2332992167380894 Clear[xxx];Clear[yyy]; stop FFS end: CPU time = 0.002 sec RFSW RADCOD RAD FLUC INTRA POL COD DAPER EMIOU CMPLO FOURI SMEAR F F F T F F T F F F F T Design orbit length = 3016.3033401153457 CONVERGENCE=1e-25; fit rfmark2 bx @ by @ ax 0 ay 0 fit rfmark2.2 bx @ by @ ax 0 ay 0 fit sfeo.2 rfmark2.2 bxm 80 bym 80 fit nx nx0 ny ny0 free qr* qv* ! go sdsq:=Apply[Plus,Element["K2","SD*"]^2]; coda="DP=0.024 fitp 24 go DP=0.02 fitp 20 go DP=0.016 fitp 16 go DP=0.02 fitp 20 go"; cell MatchingAmplitude:={{-0.026,xa,xa/3},{0.026,xa,xa/3},{0.022,xa*1.2,0},{-0.022,xa*1.2,0},{-0.024,-xa*1.1,0},{0.024,-xa*1.1,0}};xa=8; fix * rej total; ! fit rfmark2 bx @ 20 by @ 20 ax 0 20 ay 0 20 fit rfmark2.2 bx @ 20 by @ 20 ax 0 20 ay 0 20 fit nx nx0 20 ny ny0 20 ax 0 20 ay 0 20 bx .33 20 by by0 20 free s{fdc}*; nodifcod Out[68]:= nodifcod ExponentOfResidual=0.5; DP=0.02 Out[69]:= .02 gym="p=12;rep 4 DP=p*0.001 fitp p go p=p+4 until; DP=0.02 fitp 20 go" Out[70]:= "p=12;rep 4 DP=p*0.001 fitp p go p=p+4 until; DP=0.02 fitp 20 go" sf* min 0 sf* max 2.5 sd* max 0 sd* min -6 !!!! Function Test !!!!! FunctionTest=Class[{},{},{t0,a,b,c,d,e,f,g,C,a1,x,a2,b2,x2, ndim=400}, Constructor[]:=( C=Literal[Class][{},{},{aa},]; ); Fact[xx_]:=If[xx<=1,1,xx*Fact[xx-1]]; Fun[xx_,z__,xx_]:=xx+z; Test[]:=( t0=TimeUsed[]; SeedRandom[35]; Print[Timing[Do[c=C[];c@aa=1,{ndim}];"Instance"]]; a1=a=Table[Random[],{ndim},{ndim}]; b=Table[Join[Random[200],Table["b",{200}]],{ndim}]; Print[Timing[Do[D[Sin[x2]*x2^3,x2],{800}];"D"]]; Print[Timing[Do[ Fit[a[[,{1,2}]],a2*x2+b2,x2,{a2,0},{b2,0}],{300}];"Fit"]]; Print[Timing[Do[Sort[Flatten[a]],{2}];"Sort"]]; Print[Timing[Do[Plot[{Sin[x2],Cos[x2]},{x2,-Pi,Pi}, DisplayFunction->Identity],{6}];"Plot"]]; Print[Timing[Do[ListPlot[a[[,{1,2}]], DisplayFunction->Identity],{400}];"ListPlot"]]; Print[Timing[Do[Spline[a[[,{1,2}]]],{700}];"Spline"]]; Print[Timing[Do[Join[a,b],{2200}];"Join"]]; Print[Timing[Do[Append[a,b],{8000}];"Append"]]; Print[Timing[Do[ToExpression["a=1+Pi"],{60000}];"ToExpression"]]; Print[Timing[Do[ToString/@[a,{-1}],{1}];"ToString"]]; Print[Timing[Do[Fun[1,2,3,1],{40000}];"Fun"]]; Print[Timing[Do[b/.{x_String:>x//"c"},{2}];"Replace"]]; Print[Timing[Do[Do[b[[,i]],{i,ndim}],{10}];"ThreadPart"]]; Print[Timing[ Do[Do[a1[[i,j]]=a1[[j,i]],{i,ndim},{j,i}],{2}];"Part"]]; Print[Timing[Do[Fact/@[a,{-1}],{1}];"Map"]]; Print[Timing[Do[Thread[a],{30}];"Thread Real List"]]; Print[Timing[Do[Thread[b],{15}];"Thread Non-Real List"]]; Print[Timing[Do[Fact[200],{500}];"Fact"]]; Print[Timing[Do[Flatten[b],{16}];"Flatten"]]; Print[" \n******* Function Test CPU used = ",tf=TimeUsed[]-t0," sec\n"]); ]; ft=FunctionTest[]; ft@Test[]; {.13316699489951134,"Instance"} {.07584898173809052,"D"} {.06892898678779602,"Fit"} {.04306900501251221,"Sort"} {.06711101531982422,"Plot"} {.07419297099113464,"ListPlot"} {.051777005195617676,"Spline"} {.01518702507019043,"Join"} {.02437502145767212,"Append"} {.04338604211807251,"ToExpression"} {.17225104570388794,"ToString"} {.013862967491149902,"Fun"} {.08092701435089111,"Replace"} {.03201699256896973,"ThreadPart"} {.1032480001449585,"Part"} {.035827040672302246,"Map"} {.006170988082885742,"Thread Real List"} {.016003966331481934,"Thread Non-Real List"} {.08439195156097412,"Fact"} {.021400928497314453,"Flatten"} ******* Function Test CPU used = 1.1767779476940632 sec ft=.; MemoryCheck[1] Out[72]:= {8702089,10485760,592,0} !!!! Optics Calculation Test !!!!!!! t0=TimeUsed[] Out[73]:= 1.2128390073776245 Do[FFS["CALC"],{15}]; Print[" \n******* Optics Calculation Test: CPU time used = ",toc=TimeUsed[]-t0," sec\n"]; ******* Optics Calculation Test: CPU time used = 1.1532069444656372 sec !!! Tracking Test !!!!!!! t0=TimeUsed[] Out[89]:= 2.366075038909912 RFSW RADCOD e=Emittance[SaveEMIT->True] Out[92]:= {Stable->1,Tunes->{-.47942117581376276,.08030829333156965,-.0139431323245126},EnergyLossU0->841899.6811929208, RfVoltageVc->6669375.398487202,EquilibriumPosition->.011851105149355131,MomentumCompaction->.00012620782330407236, OrbitDilation->0,BucketHeight->.038987664907056,HarmonicNumber->5120.067344058354, OrbitAtExit->{-1.4012452055814659e-08,1.5048702833196652e-07,-1.1493899025236514e-09,-2.953640703863867e-07, -3.824391941056511e-14,-6.053749691366422e-07}, DampingRate->{-.00012030204052793035,-.00012026427661264297,-.00024044716498004433}, Emittances->{1.8589831294798087e-08,9.872011117942012e-13,2.369377526735907e-06},MomentumSpread->.0007388954736169904, BunchLength->.0032066480372587585,TuneShiftByRadiation->{-1.4152151073185874e-07,9.935168750626089e-07,-7.903717944824638e-09}} DynamicApertureSurvey[{{0,50},{0,50},Table[i,{i,-16,16,4}]},60,Output->6]; Maximum number of particles =126 Range Xmin: 0.000 Xmax: 50.000 (Ymin: 0.000 Ymax: 50.000) Zmin: -16.000 Zmax: 16.000 Display: 10 turns/character NZ 0----!----1----!----2----!----3----!----4----!----5 -16.00 20 *************66666660000101010000000000000000000000 -12.00 23 ****************66666660060060000000000000000000000 -8.00 28 *********************666666600000006001000000000000 -4.00 30 ***********************6666666335653000000000000000 0.00 35 ****************************66666662111111112361111 4.00 28 *********************666666646636633323332211121111 8.00 26 *******************66666663333343322222222222222111 12.00 23 ****************66666663335333332221100000000000000 16.00 21 **************6666666333332000000000000000000000000 NZ 0----!----1----!----2----!----3----!----4----!----5 Score: 234 Print[" \n******* Tracking Test CPU used = ",tt=TimeUsed[]-t0," sec\n"]; ******* Tracking Test CPU used = 1.2453789710998535 sec !!!! Matching Test !!!!!!! t0=TimeUsed[] Out[93]:= 3.6114840507507324 CONVERGENCE=2e-4; dp0=0.020; dp=dp0; FFS["fitp dp/0.001;go",6]; Matched. ( 8.3268E-04) DP = 0.02000 DP0 = 0.00000 ExponentOfResidual = 0.5 OffMomentumWeight = 1.000 DP -.020 -.018 -.016 -.014 -.012 -.010 -.008 -.006 .000 .006 .008 .010 .012 .014 .016 .018 .020 Res. .027 .027 .025 .024 .022 .020 .018 .016 .000 .018 .021 .022 .024 .025 .026 .027 .028 RFMARK2.2 AX 0 20 -.074 -.052 -.031 -.012 .002 .010 .012 .010 -3E-7 .007 .006 .002 -.006 -.016 -.025 -.032 -.033 RFMARK2.2 BX 4.362 20 4.285 4.271 4.276 4.299 4.332 4.366 4.393 4.405 4.362 4.373 4.401 4.425 4.435 4.424 4.387 4.323 4.235 RFMARK2.2 AY 0 20 -.091 -.109 -.120 -.124 -.120 -.110 -.094 -.074 2.E-8 .054 .062 .063 .059 .050 .038 .023 .009 RFMARK2.2 BY 25.103 20 23.53 23.31 23.29 23.38 23.52 23.68 23.87 24.10 25.10 26.32 26.61 26.76 26.75 26.55 26.15 25.56 24.79 $$$ fAX 0 20 -.010 -.015 -.019 -.021 -.020 -.016 -.010 -.004 3.E-7 -.011 -.013 -.014 -.016 -.018 -.022 -.026 -.031 $$$ fBX .33 20 .341 .336 .330 .325 .322 .319 .319 .321 .330 .335 .336 .338 .342 .348 .355 .364 .374 $$$ fNX 46.52 20 46.54 46.53 46.53 46.53 46.52 46.52 46.52 46.52 46.52 46.52 46.52 46.53 46.53 46.53 46.54 46.54 46.55 $$$ fAY 0 20 .032 .033 .030 .023 .015 .007 -1E-4 -.005 -3E-8 .031 .044 .056 .066 .073 .075 .072 .061 $$$ fBY .008 20 .008 .008 .008 .008 .007 .007 .008 .008 .008 .009 .009 .009 .009 .009 .009 .009 .009 $$$ fNY 46.08 20 46.10 46.09 46.09 46.09 46.08 46.08 46.08 46.08 46.08 46.08 46.08 46.08 46.08 46.08 46.07 46.07 46.07 $$$ fLENG ##### # 3016.3033 Print[" \n******* Matching Test: CPU time used = ",tm=TimeUsed[]-t0," sec\n"]; ******* Matching Test: CPU time used = .07752084732055664 sec Print["\n *********** summary of the bench mark \n", "\tFunction Test \t",tf,"\n\tOptics Test\t",toc, "\n\tTracking Test\t",tt, "\n\tMatching Test\t",tm, "\n\tOverall Score (smaller is faster): ",tf/42+(toc+tt+tm)/96];
*********** summary of the bench mark Function Test 1.1767779476940632 Optics Test1.1532069444656372 Tracking Test1.2453789710998535 Matching Test.07752084732055664 Overall Score (smaller is faster): .05381130134420735 Abort
|