Commit b71af7bb authored by frey_m's avatar frey_m

Update RingCyclotronMatche results due to improvements on implementation: EV instead of RDM, etc.

parent a406fafd
......@@ -5,6 +5,7 @@ OPTION, PSDUMPEACHTURN = FALSE;
OPTION, REPARTFREQ = 20;
OPTION, CZERO = FALSE;
OPTION, TELL = TRUE;
Option, CLOTUNEONLY=FALSE;
Option, VERSION=10900;
Title,string="OPAL-cycl: test matched distribution PSI Ring";
......@@ -60,7 +61,7 @@ rf4: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
l1: LINE = (ring,rf0,rf1,rf2,rf3,rf4);
Dist1:DISTRIBUTION, TYPE=GAUSSMATCHED, LINE=l1, FMAPFN="ring590_bfld.dat", FMTYPE="CARBONCYCL", MAGSYM= 8,
Dist1:DISTRIBUTION, TYPE=GAUSSMATCHED, LINE=l1, NSTEPS=1440, SECTOR=FALSE,
EX = 1.0e-06,
EY = 1.0e-06,
ET = 1.0e-06;
......
# 4
24996 24997 25429 24578
SDDS1
&description
text="Processor statistics 'RingCyclotronMatched.in' 28/08/201815:24:02",
contents="stat parameters"
&end
&parameter
name=processors,
type=long,
description="Number of Cores used"
&end
&parameter
name=revision,
type=string,
description="git revision of opal"
&end
&parameter
name=flavor,
type=string,
description="OPAL flavor that wrote file"
&end
&column
name=t,
type=double,
units=ns,
description="1 Time"
&end
&column
name=processor-0,
type=long,
units=1,
description="2 Number of particles of processor 0"
&end
&column
name=processor-1,
type=long,
units=1,
description="3 Number of particles of processor 1"
&end
&column
name=processor-2,
type=long,
units=1,
description="4 Number of particles of processor 2"
&end
&column
name=processor-3,
type=long,
units=1,
description="5 Number of particles of processor 3"
&end
&data
mode=ascii,
no_row_counts=1
&end
4
OPAL 2.0.0 git rev. #598d1782c7d3290d015bca8754e2b76c4a48bf1d
opal-cycl
0.000000000000000e+00 26366 23216 23310 27108
dfddffbce2dc55bc24277dace3c0f637 RingCyclotronMatched.lbal
db6c18c994176e2e65d93687c0652625 RingCyclotronMatched.lbal
Ippl> CommMPI: Parent process waiting for children ...
Ippl> CommMPI: Child 2 ready.
Ippl> CommMPI: Child 1 ready.
Ippl> CommMPI: Child 3 ready.
Ippl> CommMPI: Initialization complete.
> ____ _____ ___
> / __ \| __ \ /\ | |
> | | | | |__) / \ | |
> | | | | ___/ /\ \ | |
> | |__| | | / ____ \| |____
> \____/|_| /_/ \_\______|
OPAL>
OPAL> This is OPAL (Object Oriented Parallel Accelerator Library) Version 1.3.99.1
OPAL>
OPAL> (c) PSI, http://amas.web.psi.ch
OPAL>
OPAL> Please send cookies, goodies or other motivations (wine and beer ... )
OPAL> to the OPAL developers opal@lists.psi.ch
OPAL>
OPAL> Time: 08:04:58 date: 07/07/2015
OPAL>
OPAL> Couldn't find startup file "/Users/adelmann/init.opal".
OPAL> Note: this is not mandatory for an OPAL simulation!
OPAL>
OPAL> * Reading input stream "RingCyclotronMatched.in".
OPAL>
OPAL> Current settings of options:
OPAL> OPTION,ECHO=FALSE,INFO=TRUE,TRACE=FALSE,VERIFY=FALSE,WARN=TRUE,
OPAL> SEED=1.23457e+08,TELL=TRUE,PSDUMPFREQ=1,STATDUMPFREQ=10,
OPAL> PSDUMPEACHTURN=FALSE,PSDUMPLOCALFRAME=FALSE,SPTDUMPFREQ=1,
OPAL> REPARTFREQ=20,REBINFREQ=100,SCSOLVEFREQ=1,MTSSUBSTEPS=1,
OPAL> REMOTEPARTDEL=-1,SCAN=FALSE,RHODUMP=FALSE,EBDUMP=FALSE,
OPAL> CSRDUMP=FALSE,AUTOPHASE=0,PPDEBUG=FALSE,SURFDUMPFREQ=-1,
OPAL> NUMBLOCKS=0,RECYCLEBLOCKS=0,NLHS=1,CZERO=FALSE,RNGTYPE="RANDOM",
OPAL> SCHOTTKYCORR=FALSE,SCHOTTKYRENO=-1,ENABLEHDF5=TRUE,ASCIIDUMP=FALSE,
OPAL> BOUNDPDESTROYFQ=10;
OPAL>
OPAL> * ----------------------------------------------------
OPAL> * About to find closed orbit and matched distribution
OPAL> * I= 2.2 (mA) E= 580 (MeV)
OPAL> * EX= 1e-06* EY= 1e-06* ET= 1e-06* FMAPFN ring590_bfld.dat
OPAL> * FMSYM= 8* HN= 6* PHIINIT= 0
OPAL> * ----------------------------------------------------
OPAL> * Converged (Ex, Ey, Ez) = (0.313794, 0.313794, 0.313794) pi mm mrad for E= 580 (MeV)
OPAL> * Sigma-Matrix
OPAL> 3.914 & -0.6486 & 0 & 0 & 0.4542 & 1.362 \\
OPAL> -0.6486 & 0.6396 & 0 & 0 & 0.7265 & -0.2685 \\
OPAL> 0 & 0 & 6.239 & 1.198 & 0 & 0 \\
OPAL> 0 & 0 & 1.198 & 0.3859 & 0 & 0 \\
OPAL> 0.4542 & 0.7265 & 0 & 0 & 2.363 & 0.1752 \\
OPAL> 1.362 & -0.2685 & 0 & 0 & 0.1752 & 0.8944 \\
OPAL> * m before gsl_linalg_cholesky_decomp
OPAL> * r(0, : ) = 1 & -0.4099 & 0 & 0 & 0.1493 & 0.7277 \\
OPAL> * r(1, : ) = -0.4099 & 1 & 0 & 0 & 0.5909 & -0.355 \\
OPAL> * r(2, : ) = 0 & 0 & 1 & 0.7722 & 0 & 0 \\
OPAL> * r(3, : ) = 0 & 0 & 0.7722 & 1 & 0 & 0 \\
OPAL> * r(4, : ) = 0.1493 & 0.5909 & 0 & 0 & 1 & 0.1205 \\
OPAL> * r(5, : ) = 0.7277 & -0.355 & 0 & 0 & 0.1205 & 1 \\
OPAL> * m after gsl_linalg_cholesky_decomp
OPAL> * r(0, : ) = 1 & 0 & 0 & 0 & 0 & 0 \\
OPAL> * r(1, : ) = -0.4099 & 0.9121 & 0 & 0 & 0 & 0 \\
OPAL> * r(2, : ) = 0 & 0 & 1 & 0 & 0 & 0 \\
OPAL> * r(3, : ) = 0 & 0 & 0.7722 & 0.6354 & 0 & 0 \\
OPAL> * r(4, : ) = 0.1493 & 0.715 & 0 & 0 & 0.683 & 0 \\
OPAL> * r(5, : ) = 0.7277 & -0.06211 & 0 & 0 & 0.08229 & 0.6781 \\
OPAL> * Mass of simulation particle= 2544 GeV/c^2
OPAL> * Charge of simulation particle= 4.344e-16 [C]
OPAL> * **********************************************************************************
OPAL> * Selected Tracking Method == CYCLOTRON-T, NEW TRACK
OPAL> * **********************************************************************************
OPAL> * Number of neighbour bunches= 1
OPAL> * DT = 1e-12
OPAL> * MAXSTEPS = 0
OPAL> * Phase space dump frequency = 1
OPAL> * Statistics dump frequency = 10 w.r.t. the time step.
OPAL> * **********************************************************************************
OPAL> * ************* D I S T R I B U T I O N ********************************************
OPAL> *
OPAL> * Distribution type: MATCHEDGAUSS
OPAL> * SIGMAX = 0.001978 [m]
OPAL> * SIGMAY = 0.002498 [m]
OPAL> * SIGMAZ = 0.001537 [m]
OPAL> * SIGMAPX = 4.129e-05 [Beta Gamma]
OPAL> * SIGMAPY = 3.639e-05 [Beta Gamma]
OPAL> * SIGMAPZ = 4.49e-05 [Beta Gamma]
OPAL> * AVRGPZ = 0 [Beta Gamma]
OPAL> * CORRX = -0.4099
OPAL> * CORRY = 0.7722
OPAL> * CORRZ = 0.1205
OPAL> * R61 = 0.7277
OPAL> * R62 = -0.355
OPAL> * R51 = 0.1493
OPAL> * R52 = 0.5909
OPAL> * CUTOFFX = 3 [units of SIGMAX]
OPAL> * CUTOFFY = 3 [units of SIGMAY]
OPAL> * CUTOFFZ = 3 [units of SIGMAZ]
OPAL> * CUTOFFPX = 3 [units of SIGMAPX]
OPAL> * CUTOFFPY = 3 [units of SIGMAPY]
OPAL> * CUTOFFPZ = 3 [units of SIGMAPY]
OPAL> *
OPAL> * Distribution is injected.
OPAL> *
OPAL> * **********************************************************************************
OPAL>
OPAL> * ************* B E A M ************************************************************
OPAL> * BEAM BEAM1
OPAL> * PARTICLE PROTON
OPAL> * CURRENT 0.0022 A
OPAL> * FREQUENCY 50.65 Hz
OPAL> * CHARGE +e * 1
OPAL> * REST MASS 0.9383 GeV
OPAL> * MOMENTUM 1.194
OPAL> * NPART 1e+05
OPAL> * **********************************************************************************
OPAL>
OPAL> * ************* F I E L D S O L V E R **********************************************
OPAL> * FIELDSOLVER FS1
OPAL> * TYPE FFT
OPAL> * N-PROCESSORS 1
OPAL> * MX 32
OPAL> * MY 32
OPAL> * MT 16
OPAL> * BBOXINCR 2
OPAL> * GRRENSF INTEGRATED
OPAL> * XDIM parallel
OPAL> * YDIM parallel
OPAL> * Z(T)DIM parallel
OPAL>
OPAL> * **********************************************************************************
OPAL>
OPAL> * -------------------------- Adding Cyclotron ---------------------------- *
OPAL>
OPAL> * Bunch global starting position:
OPAL> * RINIT = 2130 [mm]
OPAL> * PHIINIT = 0 [deg]
OPAL> * ZINIT = 0 [mm]
OPAL>
OPAL> * Bunch global starting momenta:
OPAL> * Initial gamma = 1.618
OPAL> * Initial beta = 0.7862
OPAL> * Reference total momentum (beta * gamma) = 1272 [MCU]
OPAL> * Reference azimuthal momentum (Pt) = 1272 [MCU]
OPAL> * Reference radial momentum (Pr) = -0.174 [MCU]
OPAL> * Reference axial momentum (Pz) = 0 [MCU]
OPAL>
OPAL> * 8-fold field symmerty
OPAL> * Field map file name = s03av.nar
OPAL> * Type of cyclotron = RING
OPAL> * Radial aperture = 0 ... 1e+04 [mm]
OPAL> * Vertical aperture = -1e+04 ... 1e+04 [mm]
OPAL> * Harmonic number h = 6
OPAL> * ----------------------------------------------
OPAL> * READ IN RING FIELD MAP
OPAL> * (The first data block is useless)
OPAL> * ----------------------------------------------
OPAL> * Minimal radius of measured field map: 1900 [mm]
OPAL> * Stepsize in radial direction: 20 [mm]
OPAL> * Minimal angle of measured field map: 0 [deg.]
OPAL> * Stepsize in azimuth direction: 0.3333 [deg.]
OPAL> * Index in radial direction: 141
OPAL> * Index in azimuthal direction: 135
OPAL> * Accordingly, total grid point along azimuth: 136
OPAL> * Total stored grid point number ( ntetS * nrad ) : 19176
OPAL> * Read-in loop one block per radius
OPAL> * Rescaling of the fields with factor: 1
OPAL> * Field Map read successfully!
OPAL>
OPAL> * --------- RFCavity ------------------------------
OPAL> * Minimal radius of cavity= 1900 [mm]
OPAL> * Maximal radius of cavity= 4500 [mm]
OPAL> * RF frequency (2*pi*f)= 3.182e+08 [rad/s]
OPAL> * RF Field map file name= Cav1.dat
OPAL> * Cavity azimuth position= 35 [deg]
OPAL> * Cavity gap width= 220 [mm]
OPAL> * Cavity Shift distance= 416 [mm]
OPAL> * Initial RF phase (t=0)= 139.4 [deg]
OPAL> * Read cavity voltage profile data
OPAL> * Cavity voltage data read successfully!
OPAL> * --------- RFCavity ------------------------------
OPAL> * Minimal radius of cavity= 1900 [mm]
OPAL> * Maximal radius of cavity= 4500 [mm]
OPAL> * RF frequency (2*pi*f)= 3.182e+08 [rad/s]
OPAL> * RF Field map file name= Cav1.dat
OPAL> * Cavity azimuth position= 125 [deg]
OPAL> * Cavity gap width= 220 [mm]
OPAL> * Cavity Shift distance= 416 [mm]
OPAL> * Initial RF phase (t=0)= 319.4 [deg]
OPAL> * Read cavity voltage profile data
OPAL> * Cavity voltage data read successfully!
OPAL> * --------- RFCavity ------------------------------
OPAL> * Minimal radius of cavity= 1900 [mm]
OPAL> * Maximal radius of cavity= 4500 [mm]
OPAL> * RF frequency (2*pi*f)= 9.547e+08 [rad/s]
OPAL> * RF Field map file name= Cav3.dat
OPAL> * Cavity azimuth position= 170 [deg]
OPAL> * Cavity gap width= 250 [mm]
OPAL> * Cavity Shift distance= 452 [mm]
OPAL> * Initial RF phase (t=0)= 149.4 [deg]
OPAL> * Read cavity voltage profile data
OPAL> * Cavity voltage data read successfully!
OPAL> * --------- RFCavity ------------------------------
OPAL> * Minimal radius of cavity= 1900 [mm]
OPAL> * Maximal radius of cavity= 4500 [mm]
OPAL> * RF frequency (2*pi*f)= 3.182e+08 [rad/s]
OPAL> * RF Field map file name= Cav1.dat
OPAL> * Cavity azimuth position= 215 [deg]
OPAL> * Cavity gap width= 220 [mm]
OPAL> * Cavity Shift distance= 416 [mm]
OPAL> * Initial RF phase (t=0)= 139.4 [deg]
OPAL> * Read cavity voltage profile data
OPAL> * Cavity voltage data read successfully!
OPAL> * --------- RFCavity ------------------------------
OPAL> * Minimal radius of cavity= 1900 [mm]
OPAL> * Maximal radius of cavity= 4500 [mm]
OPAL> * RF frequency (2*pi*f)= 3.182e+08 [rad/s]
OPAL> * RF Field map file name= Cav1.dat
OPAL> * Cavity azimuth position= 305 [deg]
OPAL> * Cavity gap width= 220 [mm]
OPAL> * Cavity Shift distance= 416 [mm]
OPAL> * Initial RF phase (t=0)= 319.4 [deg]
OPAL> * Read cavity voltage profile data
OPAL> * Cavity voltage data read successfully!
OPAL> * -------------------------------------
OPAL> * The selected Beam line elements are :
OPAL> * -> CYCLOTRON
OPAL> * -> CAVITY
OPAL> * -> CAVITY
OPAL> * -> CAVITY
OPAL> * -> CAVITY
OPAL> * -> CAVITY
OPAL> * -------------------------------------
OPAL> * 4th order Runge-Kutta integrator
OPAL> * Before repartition: min local particle number 100000 max local particle number: 100000
OPAL> * After repartition: min local particle number 100000 max local particle number: 100000
OPAL>
OPAL> * *********************** Bunch information in local frame: ************************
OPAL> * ************** B U N C H *********************************************************
OPAL> * NP = 100000
OPAL> * Qtot = 4.34353e-02 [nC] Qi = 4.34353e-07 [nC]
OPAL> * Ekin = 5.80000e+02 [MeV] dEkin = 0.00000e+00 [MeV]
OPAL> * rmax = ( 6.17679e-03 , 7.79303e-03 , 6.44782e-03 ) [m]
OPAL> * rmin = ( -6.17044e-03 , -7.77654e-03 , -6.41690e-03 ) [m]
OPAL> * rms beam size = ( 1.94408e-03 , 2.47775e-03 , 1.51444e-03 ) [m]
OPAL> * rms momenta = ( 4.06990e-05 , 3.59214e-05 , 4.41418e-05 ) [beta gamma]
OPAL> * mean position = ( 3.81344e-15 , -2.54458e-19 , 1.46818e-19 ) [m]
OPAL> * mean momenta = ( 1.72287e-16 , 1.27218e+00 , -8.12833e-18 ) [beta gamma]
OPAL> * rms emittance = ( 5.66347e-08 , 4.42968e-08 , 5.21833e-08 ) (not normalized)
OPAL> * rms correlation = ( -4.13270e-01 , 7.74027e-01 , 1.17579e-01 )
OPAL> * hr = ( 3.98298e-04 , 5.02244e-04 , 8.57648e-04 ) [m]
OPAL> * dh = 2.00000e-02 [m]
OPAL> * t = 0.00000e+00 [s] dT = inf [s]
OPAL> * spos = 1.46818e-19 [m]
OPAL> * **********************************************************************************
OPAL>
OPAL>
OPAL> * Phase space dump 0 (global frame) at integration step 0
OPAL> * T = 0 ns, Live Particles: 100000
OPAL> * E = 580 MeV, beta * gamma = 1.2722
OPAL> * Bunch position: R = 2130 mm, Theta = -0.00013946 Deg, Z = -0.00025179 mm
OPAL> * Local Azimuth = 0.0078358 Deg, Local Elevation = -0.0028944 Deg
OPAL>
OPAL> * *********************** Bunch information in global frame: ***********************
OPAL> * ************** B U N C H *********************************************************
OPAL> * NP = 100000
OPAL> * Qtot = 4.34353e-02 [nC] Qi = 4.34353e-07 [nC]
OPAL> * Ekin = 5.80000e+02 [MeV] dEkin = 3.25616e-02 [MeV]
OPAL> * rmax = ( 6.17679e-03 , 7.79303e-03 , 6.44782e-03 ) [m]
OPAL> * rmin = ( -6.17044e-03 , -7.77654e-03 , -6.41690e-03 ) [m]
OPAL> * rms beam size = ( 1.94408e-03 , 2.47775e-03 , 1.51444e-03 ) [m]
OPAL> * rms momenta = ( 4.06990e-05 , 3.59220e-05 , 4.41418e-05 ) [beta gamma]
OPAL> * mean position = ( 2.12999e+00 , -5.18450e-06 , -2.51787e-07 ) [m]
OPAL> * mean momenta = ( -1.73985e-04 , 1.27218e+00 , -6.42670e-05 ) [beta gamma]
OPAL> * rms emittance = ( 5.66348e-08 , 4.42986e-08 , 5.21833e-08 ) (not normalized)
OPAL> * rms correlation = ( -4.13269e-01 , 7.74013e-01 , 1.17580e-01 )
OPAL> * hr = ( 3.98298e-04 , 5.02244e-04 , 8.57648e-04 ) [m]
OPAL> * dh = 2.00000e-02 [m]
OPAL> * t = 0.00000e+00 [s] dT = inf [s]
OPAL> * spos = -2.51787e-07 [m]
OPAL> * **********************************************************************************
OPAL>
OPAL> * Beginning of this run is at t = 0 [ns]
OPAL> * The time step is set to dt = inf [ns]
OPAL> * Single particle trajectory dump frequency is set to 1
OPAL> * The frequency to solve space charge fields is set to 1
OPAL> * The repartition frequency is set to 20
OPAL>
OPAL> * --------------------------------- Start tracking ------------------------------------ *
OPAL>
OPAL> * ---------------------------- DONE TRACKING PARTICLES -------------------------------- *
OPAL> * Final energy of reference particle = 580 [MeV]
OPAL> * Total phase space dump number(includes the initial distribution) = 1
OPAL> * One can restart simulation from the last dump step (-restart 0)
OPAL> *
OPAL> * Finished during turn 1 (0 turns completed)
OPAL> * Cave: Turn number is not correct for restart mode
OPAL>
OPAL> * *********************** Bunch information in global frame: ***********************
OPAL> * ************** B U N C H *********************************************************
OPAL> * NP = 100000
OPAL> * Qtot = 4.34353e-02 [nC] Qi = 4.34353e-07 [nC]
OPAL> * Ekin = 5.80000e+02 [MeV] dEkin = 3.25616e-02 [MeV]
OPAL> * rmax = ( 6.17679e-03 , 7.79303e-03 , 6.44782e-03 ) [m]
OPAL> * rmin = ( -6.17044e-03 , -7.77654e-03 , -6.41690e-03 ) [m]
OPAL> * rms beam size = ( 1.94408e-03 , 2.47775e-03 , 1.51444e-03 ) [m]
OPAL> * rms momenta = ( 4.06990e-05 , 3.59220e-05 , 4.41418e-05 ) [beta gamma]
OPAL> * mean position = ( 2.12999e+00 , -5.18450e-06 , -2.51787e-07 ) [m]
OPAL> * mean momenta = ( -1.73985e-04 , 1.27218e+00 , -6.42670e-05 ) [beta gamma]
OPAL> * rms emittance = ( 5.66348e-08 , 4.42986e-08 , 5.21833e-08 ) (not normalized)
OPAL> * rms correlation = ( -4.13269e-01 , 7.74013e-01 , 1.17580e-01 )
OPAL> * hr = ( 3.98298e-04 , 5.02244e-04 , 8.57648e-04 ) [m]
OPAL> * dh = 2.00000e-02 [m]
OPAL> * t = 0.00000e+00 [s] dT = inf [s]
OPAL> * spos = -2.51787e-07 [m]
OPAL> * **********************************************************************************
OPAL>
OPAL> * ----------------------------------------------- *
OPAL> * Finalizing i.e. write data and close files :
OPAL> * Finalize cyclotron
OPAL> * ----------------------------------------------- *
OPAL> * End of input stream "RingCyclotronMatched.in".
Timings> -----------------------------------------------------------------
Timings> Timing results for 1 nodes:
Timings> -----------------------------------------------------------------
Timings> mainTimer........... Wall tot = 22.8619, CPU tot = 22.32
Timings>
Timings> Binaryrepart........ Wall max = 0.009949, CPU max = 0.01
Timings> Wall avg = 0.009949, CPU avg = 0.01
Timings> Wall min = 0.009949, CPU min = 0.01
Timings>
Timings> Boundingbox......... Wall max = 0.019488, CPU max = 0.02
Timings> Wall avg = 0.019488, CPU avg = 0.02
Timings> Wall min = 0.019488, CPU min = 0.02
Timings>
Timings> ComputePotential.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> CreatDistr.......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Dump................ Wall max = 0.089778, CPU max = 0.04
Timings> Wall avg = 0.089778, CPU avg = 0.04
Timings> Wall min = 0.089778, CPU min = 0.04
Timings>
Timings> Fast inside test.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Frametransform...... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> H5PartTimer......... Wall max = 0.082904, CPU max = 0.03
Timings> Wall avg = 0.082904, CPU avg = 0.03
Timings> Wall min = 0.082904, CPU min = 0.03
Timings>
Timings> Histogram........... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Initialize geometry. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Inside test......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Integration......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> LoadDistr........... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Particle Inside..... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Ray tracing......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Secondary emission.. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SelfField total..... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: GreensFTotal.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: IntGreenF1...... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: IntGreenF2...... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: IntGreenF3...... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: MirrorRho1...... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: MirrorRho2...... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: Potential....... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: ShIntGreenF1.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: ShIntGreenF2.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: ShIntGreenF3.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SF: ShIntGreenF4.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Statistics.......... Wall max = 0.008464, CPU max = 0.02
Timings> Wall avg = 0.008464, CPU avg = 0.02
Timings> Wall min = 0.008464, CPU min = 0.02
Timings>
Timings> StatMarkerTimer..... Wall max = 0.004854, CPU max = 0.01
Timings> Wall avg = 0.004854, CPU avg = 0.01
Timings> Wall min = 0.004854, CPU min = 0.01
Timings>
Timings> -----------------------------------------------------------------
Ippl> CommMPI: Started job 1 on host `R2-D2'.
Ippl> CommMPI: Started job 2 on host `R2-D2'.
Ippl> CommMPI: Started job 3 on host `R2-D2'.
{0}> ____ _____ ___
{0}> / __ \| __ \ /\ | |
{0}> | | | | |__) / \ | |
{0}> | | | | ___/ /\ \ | |
{0}> | |__| | | / ____ \| |____
{0}> \____/|_| /_/ \_\______|
OPAL{0}>
OPAL{0}> This is OPAL (Object Oriented Parallel Accelerator Library) Version 2.0.0
OPAL{0}>
OPAL{0}> (c) PSI, http://amas.web.psi.ch
OPAL{0}> The optimiser (former opt-Pilot) is integrated
OPAL{0}>
OPAL{0}> Please send cookies, goodies or other motivations (wine and beer ... )
OPAL{0}> to the OPAL developers opal@lists.psi.ch
OPAL{0}>
OPAL{0}> Time: 15:23:52 date: 28/08/2018
OPAL{0}>
OPAL{0}> Couldn't find startup file "/home/matthias/init.opal".
OPAL{0}> Note: this is not mandatory for an OPAL simulation!
OPAL{0}>
OPAL{0}> * Reading input stream "RingCyclotronMatched.in".
OPAL{0}>
OPAL{0}> Current settings of options:
OPAL{0}> OPTION,ECHO=FALSE,INFO=TRUE,TRACE=FALSE,WARN=TRUE,SEED=1.23457e+08,
OPAL{0}> TELL=TRUE,PSDUMPFREQ=1,STATDUMPFREQ=10,PSDUMPEACHTURN=FALSE,
OPAL{0}> PSDUMPFRAME="GLOBAL",SPTDUMPFREQ=1,REPARTFREQ=20,REBINFREQ=100,
OPAL{0}> SCSOLVEFREQ=1,MTSSUBSTEPS=1,REMOTEPARTDEL=0,RHODUMP=FALSE,
OPAL{0}> EBDUMP=FALSE,CSRDUMP=FALSE,AUTOPHASE=6,PPDEBUG=FALSE,
OPAL{0}> SURFDUMPFREQ=-1,NUMBLOCKS=0,RECYCLEBLOCKS=0,NLHS=1,CZERO=FALSE,
OPAL{0}> RNGTYPE="RANDOM",ENABLEHDF5=TRUE,ASCIIDUMP=FALSE,BOUNDPDESTROYFQ=10,
OPAL{0}> BEAMHALOBOUNDARY=0,CLOTUNEONLY=FALSE,IDEALIZED=FALSE,
OPAL{0}> LOGBENDTRAJECTORY=FALSE,VERSION=10000,MEMORYDUMP=FALSE;
OPAL{0}>
OPAL{0}> * ----------------------------------------------------
OPAL{0}> * About to find closed orbit and matched distribution
OPAL{0}> * I= 2.2 (mA) E= 580 (MeV)
OPAL{0}> * EX= 1e-06 EY= 1e-06 ET= 1e-06
OPAL{0}> * SECTOR: match using all sectors
OPAL{0}> * NSTEPS = 1440
OPAL{0}> * HN= 6 PHIINIT= 0
OPAL{0}> * ----------------------------------------------------
OPAL{0}> * ----------------------------------------------
OPAL{0}> * READ IN RING FIELD MAP
OPAL{0}> * (The first data block is useless)
OPAL{0}> * ----------------------------------------------
OPAL{0}> * Minimal radius of measured field map: 1900 [mm]
OPAL{0}> * Stepsize in radial direction: 20 [mm]
OPAL{0}> * Minimal angle of measured field map: 0 [deg.]
OPAL{0}> * Stepsize in azimuth direction: 0.333333 [deg.]
OPAL{0}> * Index in radial direction: 141
OPAL{0}> * Index in azimuthal direction: 135
OPAL{0}> * Accordingly, total grid point along azimuth: 136
OPAL{0}> * Total stored grid point number ( ntetS * nrad ) : 19176
OPAL{0}> * Read-in loop one block per radius
OPAL{0}> * Rescaling of the fields with factor: 1
OPAL{0}> * Field Map read successfully!
OPAL{0}>
OPAL{0}> * ----------------------------
OPAL{0}> * Closed orbit info:
OPAL{0}> *
OPAL{0}> * average radius: 4.40772 [m]
OPAL{0}> * initial radius: 4.38192 [m]
OPAL{0}> * initial momentum: 0.211084 [Beta Gamma]
OPAL{0}> * frequency error: -0.0136027 [ % ]
OPAL{0}> * horizontal tune: 1.69214
OPAL{0}> * vertical tune: 0.828534
OPAL{0}> * ----------------------------
OPAL{0}>
OPAL{0}> * Converged (Ex, Ey, Ez) = (0.313794, 0.313794, 0.313794) pi mm mrad for E= 580 (MeV)
OPAL{0}> * Sigma-Matrix
OPAL{0}> 2.418 & 0.06913 & 0 & 0 & 0 & 0 \\
OPAL{0}> 0.06913 & 0.5213 & 0 & 0 & 0 & 0.05129 \\
OPAL{0}> -2.506e-16 & 0 & 7.248 & 0 & 0 & 0 \\
OPAL{0}> 1.167e-17 & 0 & 0 & 0.2425 & 0 & 0 \\
OPAL{0}> -0.5205 & 0 & 0 & 0 & 5.674 & 0 \\
OPAL{0}> -0.3361 & 0.05129 & 0 & 0 & 0 & 0.2216 \\
OPAL{0}> * m before gsl_linalg_cholesky_decomp
OPAL{0}> * r(0, : ) = 1 & 0.06158 & 0 & 0 & -0.1405 & -0.4592 \\
OPAL{0}> * r(1, : ) = 0.06158 & 1 & 0 & 0 & -0.4446 & 0.1509 \\
OPAL{0}> * r(2, : ) = 0 & 0 & 1 & -0.6686 & 0 & 0 \\
OPAL{0}> * r(3, : ) = 0 & 0 & -0.6686 & 1 & 0 & 0 \\
OPAL{0}> * r(4, : ) = -0.1405 & -0.4446 & 0 & 0 & 1 & -0.0387 \\
OPAL{0}> * r(5, : ) = -0.4592 & 0.1509 & 0 & 0 & -0.0387 & 1 \\
OPAL{0}> * m after gsl_linalg_cholesky_decomp
OPAL{0}> * r(0, : ) = 1 & 0 & 0 & 0 & 0 & 0 \\
OPAL{0}> * r(1, : ) = 0.06158 & 0.9981 & 0 & 0 & 0 & 0 \\
OPAL{0}> * r(2, : ) = 0 & 0 & 1 & 0 & 0 & 0 \\
OPAL{0}> * r(3, : ) = 0 & 0 & -0.6686 & 0.7436 & 0 & 0 \\
OPAL{0}> * r(4, : ) = -0.1405 & -0.4368 & 0 & 0 & 0.8885 & 0 \\
OPAL{0}> * r(5, : ) = -0.4592 & 0.1795 & 0 & 0 & -0.02795 & 0.8696 \\
OPAL{0}> * Mass of simulation particle= 2544 GeV/c^2
OPAL{0}> * Charge of simulation particle= 4.344e-16 [C]
OPAL{0}> * **********************************************************************************
OPAL{0}> * Selected Tracking Method == CYCLOTRON-T, NEW TRACK
OPAL{0}> * **********************************************************************************
OPAL{0}> * Number of neighbour bunches= 1
OPAL{0}> * DT = 1e-12
OPAL{0}> * MAXSTEPS = 0
OPAL{0}> * Phase space dump frequency = 1
OPAL{0}> * Statistics dump frequency = 10 w.r.t. the time step.
OPAL{0}> * **********************************************************************************
OPAL{0}>
OPAL{0}> * ************* D I S T R I B U T I O N ********************************************
OPAL{0}> *
OPAL{0}> * Number of particles: 100000
OPAL{0}> *
OPAL{0}> * Distribution type: MATCHEDGAUSS
OPAL{0}> * SIGMAX = 0.001555 [m]
OPAL{0}> * SIGMAY = 0.002382 [m]
OPAL{0}> * SIGMAZ = 0.002692 [m]
OPAL{0}> * SIGMAPX = 0.0009185 [Beta Gamma]
OPAL{0}> * SIGMAPY = 0.0009683 [Beta Gamma]
OPAL{0}> * SIGMAPZ = 0.0006265 [Beta Gamma]
OPAL{0}> * AVRGPZ = 0 [Beta Gamma]
OPAL{0}> * CORRX = 0.06158
OPAL{0}> * CORRY = -0.6686
OPAL{0}> * CORRZ = -0.0387
OPAL{0}> * R61 = -0.4592
OPAL{0}> * R62 = 0.1509
OPAL{0}> * R51 = -0.1405
OPAL{0}> * R52 = -0.4446
OPAL{0}> * CUTOFFX = 3 [units of SIGMAX]
OPAL{0}> * CUTOFFY = 3 [units of SIGMAY]
OPAL{0}> * CUTOFFLONG = 3 [units of SIGMAZ]
OPAL{0}> * CUTOFFPX = 3 [units of SIGMAPX]
OPAL{0}> * CUTOFFPY = 3 [units of SIGMAPY]
OPAL{0}> * CUTOFFPZ = 3 [units of SIGMAPY]
OPAL{0}> *
OPAL{0}> * Distribution is injected.
OPAL{0}> *
OPAL{0}> * **********************************************************************************
OPAL{0}>
OPAL{0}> * ************* B E A M ************************************************************
OPAL{0}> * BEAM BEAM1
OPAL{0}> * PARTICLE PROTON
OPAL{0}> * CURRENT 0.0022 A
OPAL{0}> * FREQUENCY 50.65 MHz
OPAL{0}> * CHARGE +e * 1
OPAL{0}> * REST MASS 0.9383 GeV
OPAL{0}> * MOMENTUM 1.194
OPAL{0}> * NPART 1e+05
OPAL{0}> * **********************************************************************************