Commit 1817a7b8 authored by frey_m's avatar frey_m

add AMR test 5

parent 3499ec2a
Option, ECHO=FALSE;
Option, PSDUMPFREQ=4500;
Option, PSDUMPFRAME=BUNCH_MEAN;
Option, SPTDUMPFREQ = 1000000;
Option, ENABLEHDF5=TRUE;
Option, PSDUMPEACHTURN=TRUE;
Option, REPARTFREQ=5;
Option, ECHO=FALSE;
Option, STATDUMPFREQ=1;
Option, CZERO=FALSE;
Option, VERSION=20000;
Option, SCsolveFreq=2;
Option, AMR=TRUE;
Option, AMR_REGRID_FREQ=10;
Option, AMR_YT_DUMP_FREQ=1000000;
Option, MEMORYDUMP=TRUE;
Option, TELL=TRUE;
Title,string="OPAL-cycl: the first turn acceleration in PSI 590MeV Ring";
REAL Edes=.072;
REAL gamma=(Edes+PMASS)/PMASS;
REAL beta=sqrt(1-(1/gamma^2));
REAL gambet=gamma*beta;
REAL P0 = gamma*beta*PMASS;
REAL brho = (PMASS*1.0e9*gambet) / CLIGHT;
//value,{gamma,brho,Edes,beta,gambet};
REAL phi01=139.4281;
REAL phi02=phi01+180.0;
REAL phi04=phi01;
REAL phi05=phi01+180.0;
REAL phi03=phi01+10.0;
REAL volt1st=0.9;
REAL volt3rd=0.9*4.0*0.112;
REAL turns =3;
REAL nstep=360;
REAL frequency=50.650;
REAL frequency3=3.0*frequency;
ring: Cyclotron, TYPE="RING", CYHARMON=6, PHIINIT=0.0,
PRINIT=-0.000174, RINIT=2130.0, SYMMETRY=8.0, RFFREQ=frequency,
FMAPFN="s03av.nar";
rf0: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=35.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi01;
rf1: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=125.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi02;
rf2: RFCavity, VOLT=volt3rd, FMAPFN="Cav3.dat", TYPE="SINGLEGAP",
FREQ=frequency3,RMIN = 1900.0, RMAX = 4500.0, ANGLE=170.0, PDIS = 452.0,
GAPWIDTH = 250.0, PHI0=phi03;
rf3: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=215.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi04;
rf4: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=305.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi05;
l1: Line = (ring,rf0,rf1,rf2,rf3,rf4);
Dist1:DISTRIBUTION, TYPE=gauss,
sigmax = 2.0e-03,
sigmapx = 1.0e-7,
corrx = 0.0,
sigmay = 2.0e-03,
sigmapy = 1.0e-7,
corry = 0.0,
sigmat = 2.0e-03,
sigmapt = 3.394e-4,
corrt=0.0;
Fs1:FIELDSOLVER, FSTYPE=AMR_MG, MX=32, MY=32, MT=32,
PARFFTX=true, PARFFTY=true, PARFFTT=true,
BCFFTX=dirichlet, BCFFTY=dirichlet, BCFFTT=dirichlet,
BBOXINCR=20, AMR_MAXLEVEL=1, AMR_MAXGRIDX=32, AMR_MAXGRIDY=32,
AMR_MAXGRIDZ=32, AMR_BFX=16, AMR_BFY=16, AMR_BFZ=16,
AMR_REFX=2, AMR_REFY=2, AMR_REFZ=2, AMR_DOMAIN_RATIO={1.0, 0.75, 0.75},
AMR_TAGGING=CHARGE_DENSITY, AMR_DENSITY=1.0e-9,
AMR_MG_REBALANCE=TRUE, AMR_MG_NSWEEPS=12, ITSOLVER=SA, AMR_MG_NORM=LINF;
beam1: BEAM, PARTICLE=PROTON, pc=P0, NPART=32*32*32, BCURRENT=2.0E-3, CHARGE=1.0, BFREQ=frequency;
Select, Line=l1;
TRACK,LINE=l1, BEAM=beam1, MAXSTEPS=nstep*turns, STEPSPERTURN=nstep,TIMEINTEGRATOR="RK-4";
run, method = "CYCLOTRON-T", beam=beam1, fieldsolver=Fs1, distribution=Dist1,
MBMODE=FORCE, TURNS=2, MB_BINNING=GAMMA, MB_ETA=0.25;
endtrack;
Stop;
#!/bin/bash
mpirun -np 4 $OPAL_EXE_PATH/opal AmrRingCyclotron.in --info 2 --warn 0 "$@" 2>&1
"Neighboring bunch AMR test where we track two bunches by three turns. It uses the AMR_MG Poisson solver based on Trilinos
packages and one AMR level."
stat "rms_x" last 1E-7 AmrRingCyclotron-bunch-0000.smb # first bunch
stat "rms_y" last 1E-7 AmrRingCyclotron-bunch-0000.smb
stat "rms_z" last 1E-7 AmrRingCyclotron-bunch-0000.smb
stat "emit_x" last 1E-8 AmrRingCyclotron-bunch-0000.smb
stat "emit_y" last 1E-8 AmrRingCyclotron-bunch-0000.smb
stat "emit_z" last 1E-8 AmrRingCyclotron-bunch-0000.smb
stat "halo_x" last 1E-5 AmrRingCyclotron-bunch-0000.smb
stat "halo_y" last 1E-5 AmrRingCyclotron-bunch-0000.smb
stat "halo_z" last 1E-5 AmrRingCyclotron-bunch-0000.smb
stat "rms_x" last 1E-7 AmrRingCyclotron-bunch-0001.smb # second bunch
stat "rms_y" last 1E-7 AmrRingCyclotron-bunch-0001.smb
stat "rms_z" last 1E-7 AmrRingCyclotron-bunch-0001.smb
stat "emit_x" last 1E-8 AmrRingCyclotron-bunch-0001.smb
stat "emit_y" last 1E-8 AmrRingCyclotron-bunch-0001.smb
stat "emit_z" last 1E-8 AmrRingCyclotron-bunch-0001.smb
stat "halo_x" last 1E-5 AmrRingCyclotron-bunch-0001.smb
stat "halo_y" last 1E-5 AmrRingCyclotron-bunch-0001.smb
stat "halo_z" last 1E-5 AmrRingCyclotron-bunch-0001.smb
21
0.0e+00 0.63013 1.08
5.0e-02 0.68626 1.08
1.0e-01 0.73857 1.00
1.5e-01 0.78677 9.20e-1
2.0e-01 0.83060 8.30e-1
2.5e-01 0.86979 7.35e-1
3.0e-01 0.90416 6.36e-1
3.5e-01 0.93349 5.34e-1
4.0e-01 0.95762 4.29e-1
4.5e-01 0.97643 3.21e-1
5.0e-01 0.98981 2.12e-1
5.5e-01 0.99766 1.01e-1
6.0e-01 0.99998 -9.31e-3
6.5e-01 0.99673 -1.20e-1
7.0e-01 0.98793 -2.30e-1
7.5e-01 0.97363 -3.40e-1
8.0e-01 0.95393 -4.47e-1
8.5e-01 0.92891 -5.52e-1
9.0e-01 0.89871 -6.53e-1
9.5e-01 0.86352 -7.51e-1
1.0e+00 0.82354 -7.51e-1
\ No newline at end of file
21
0.0e+00 0.109737 2.91
5.0e-02 0.258532 2.91
1.0e-01 0.401421 2.76
1.5e-01 0.535168 2.55
2.0e-01 0.656550 2.27
2.5e-01 0.762840 1.95
3.0e-01 0.851785 1.58
3.5e-01 0.921153 1.17
4.0e-01 0.969231 7.40e-1
4.5e-01 0.995244 2.92e-1
5.0e-01 0.998514 -1.65e-1
5.5e-01 0.978652 -6.20e-1
6.0e-01 0.936425 -1.05
6.5e-01 0.872880 -1.47
7.0e-01 0.789161 -1.85
7.5e-01 0.687305 -2.19
8.0e-01 0.569821 -2.48
8.5e-01 0.439238 -2.71
9.0e-01 0.298508 -2.88
9.5e-01 0.150978 -2.98
1.0e+00 0.222050 -2.98
\ No newline at end of file
This source diff could not be displayed because it is too large. You can view the blob instead.
ac9f1ba2f0d05627c0f39151864bb5cb AmrRingCyclotron-bunch-0000.smb
This source diff could not be displayed because it is too large. You can view the blob instead.
5425b92b8d7652fd8103b8e98571e1de AmrRingCyclotron-bunch-0001.smb
[R2-D2:18272] mca: base: component_find: unable to open /opt/psi/Compiler/openmpi/1.10.7/gcc/7.3.0/lib/openmpi/mca_btl_openib: librdmacm.so.1: cannot open shared object file: No such file or directory (ignored)
[R2-D2:18270] mca: base: component_find: unable to open /opt/psi/Compiler/openmpi/1.10.7/gcc/7.3.0/lib/openmpi/mca_btl_openib: librdmacm.so.1: cannot open shared object file: No such file or directory (ignored)
[R2-D2:18269] mca: base: component_find: unable to open /opt/psi/Compiler/openmpi/1.10.7/gcc/7.3.0/lib/openmpi/mca_btl_openib: librdmacm.so.1: cannot open shared object file: No such file or directory (ignored)
[R2-D2:18270] mca: base: component_find: unable to open /opt/psi/Compiler/openmpi/1.10.7/gcc/7.3.0/lib/openmpi/mca_mtl_psm: libpsm_infinipath.so.1: cannot open shared object file: No such file or directory (ignored)
[R2-D2:18272] mca: base: component_find: unable to open /opt/psi/Compiler/openmpi/1.10.7/gcc/7.3.0/lib/openmpi/mca_mtl_psm: libpsm_infinipath.so.1: cannot open shared object file: No such file or directory (ignored)
[R2-D2:18271] mca: base: component_find: unable to open /opt/psi/Compiler/openmpi/1.10.7/gcc/7.3.0/lib/openmpi/mca_btl_openib: librdmacm.so.1: cannot open shared object file: No such file or directory (ignored)
[R2-D2:18269] mca: base: component_find: unable to open /opt/psi/Compiler/openmpi/1.10.7/gcc/7.3.0/lib/openmpi/mca_mtl_psm: libpsm_infinipath.so.1: cannot open shared object file: No such file or directory (ignored)
[R2-D2:18271] mca: base: component_find: unable to open /opt/psi/Compiler/openmpi/1.10.7/gcc/7.3.0/lib/openmpi/mca_mtl_psm: libpsm_infinipath.so.1: cannot open shared object file: No such file or directory (ignored)
Ippl> CommMPI: Started job 1 on host `R2-D2'.
Ippl> CommMPI: Started job 3 on host `R2-D2'.
Ippl> CommMPI: Parent process waiting for children ...
Ippl> CommMPI: Child 1 ready.
Ippl> CommMPI: Child 3 ready.
Ippl> CommMPI: Child 2 ready.
Ippl> CommMPI: Initialization complete.
Ippl> CommMPI: Started job 2 on host `R2-D2'.
{0}> ____ _____ ___
{0}> / __ \| __ \ /\ | |
{0}> | | | | |__) / \ | |
{0}> | | | | ___/ /\ \ | |
{0}> | |__| | | / ____ \| |____
{0}> \____/|_| /_/ \_\______|
OPAL{0}>
OPAL{0}> This is OPAL (Object Oriented Parallel Accelerator Library) Version 2.1.0
OPAL{0}> git rev. 88400b8ffdcba529ace97bd6f40d16c55d9e4ca8
OPAL{0}>
OPAL{0}>
OPAL{0}> (c) PSI, http://amas.web.psi.ch
OPAL{0}>
OPAL{0}>
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:22:35 date: 03/04/2020
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 "AmrRingCyclotron.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=4500,STATDUMPFREQ=1,PSDUMPEACHTURN=TRUE,
OPAL{0}> PSDUMPFRAME="BUNCH_MEAN",SPTDUMPFREQ=1e+06,REPARTFREQ=5,
OPAL{0}> REBINFREQ=100,SCSOLVEFREQ=2,MTSSUBSTEPS=1,REMOTEPARTDEL=0,
OPAL{0}> RHODUMP=FALSE,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=20000,AMR=TRUE,
OPAL{0}> AMR_YT_DUMP_FREQ=1e+06,AMR_REGRID_FREQ=10,MEMORYDUMP=TRUE,
OPAL{0}> HALOSHIFT=0,DELPARTFREQ=1,MINBINEMITTED=10,MINSTEPFORREBIN=200;
OPAL{0}>
Ippl{0}> Creating new file for solver information: AmrRingCyclotron.solver
OPAL{0}[2]> * Generation of distribution with seed = 123456789
OPAL{0}[2]> * isn't scalable with number of particles and cores.
OPAL{0}> * m before gsl_linalg_cholesky_decomp
OPAL{0}> * r(0, : ) = 1 & 0 & 0 & 0 & 0 & 0 \\
OPAL{0}> * r(1, : ) = 0 & 1 & 0 & 0 & 0 & 0 \\
OPAL{0}> * r(2, : ) = 0 & 0 & 1 & 0 & 0 & 0 \\
OPAL{0}> * r(3, : ) = 0 & 0 & 0 & 1 & 0 & 0 \\
OPAL{0}> * r(4, : ) = 0 & 0 & 0 & 0 & 1 & 0 \\
OPAL{0}> * r(5, : ) = 0 & 0 & 0 & 0 & 0 & 1 \\
OPAL{0}> * m after gsl_linalg_cholesky_decomp
OPAL{0}> * r(0, : ) = 1 & 0 & 0 & 0 & 0 & 0 \\
OPAL{0}> * r(1, : ) = 0 & 1 & 0 & 0 & 0 & 0 \\
OPAL{0}> * r(2, : ) = 0 & 0 & 1 & 0 & 0 & 0 \\
OPAL{0}> * r(3, : ) = 0 & 0 & 0 & 1 & 0 & 0 \\
OPAL{0}> * r(4, : ) = 0 & 0 & 0 & 0 & 1 & 0 \\
OPAL{0}> * r(5, : ) = 0 & 0 & 0 & 0 & 0 & 1 \\
OPAL{0}> * Mass of simulation particle= 7057 GeV/c^2
OPAL{0}> * Charge of simulation particle= 1.205e-15 [C]
Ippl{0}> Creating new file for data: AmrRingCyclotron.stat
Ippl{0}> Creating new file for data: AmrRingCyclotron.lbal
Ippl{0}> Creating new file for data: AmrRingCyclotron.grid
Ippl{0}> Creating new file for data: AmrRingCyclotron.mem
Ippl{0}> Creating new file for data: AmrRingCyclotron-bunch-0000.smb
Ippl{0}> Creating new file for data: AmrRingCyclotron-bunch-0001.smb
OPAL{0}> * **********************************************************************************
OPAL{0}> * Selected Tracking Method == CYCLOTRON-T, NEW TRACK
OPAL{0}> * **********************************************************************************
OPAL{0}> * Number of neighbour bunches= 2
OPAL{0}> * DT = 1e-12
OPAL{0}> * MAXSTEPS = 1080
OPAL{0}> * Phase space dump frequency = 4500
OPAL{0}> * Statistics dump frequency = 1 w.r.t. the time step.
OPAL{0}> * **********************************************************************************
OPAL{0}> Use 'GAMMA' for binning.
OPAL{0}> ***---------------------------- MULTI-BUNCHES MULTI-ENERGY-BINS MODE ----------------------------***
OPAL{0}> Read in: Bin=0 Particles Num=131072
OPAL{0}> FORCE mode: The multi bunches will be injected consecutively
OPAL{0}> after each revolution, until get "TURNS" bunches.
OPAL{0}>
OPAL{0}> * ************* D I S T R I B U T I O N ********************************************
OPAL{0}> *
OPAL{0}> * Number of particles: 32768
OPAL{0}> *
OPAL{0}> * Distribution type: GAUSS
OPAL{0}> *
OPAL{0}> * SIGMAX = 0.002 [m]
OPAL{0}> * SIGMAY = 0.002 [m]
OPAL{0}> * SIGMAZ = 0.002 [m]
OPAL{0}> * SIGMAPX = 1.46e-08 [Beta Gamma]
OPAL{0}> * SIGMAPY = 1.46e-08 [Beta Gamma]
OPAL{0}> * SIGMAPZ = 0 [Beta Gamma]
OPAL{0}> * AVRGPZ = 0 [Beta Gamma]
OPAL{0}> * CORRX = 0
OPAL{0}> * CORRY = 0
OPAL{0}> * CORRZ = 0
OPAL{0}> * R61 = 0
OPAL{0}> * R62 = 0
OPAL{0}> * R51 = 0
OPAL{0}> * R52 = 0
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.002 A
OPAL{0}> * FREQUENCY 50.65 MHz
OPAL{0}> * CHARGE +e * 1
OPAL{0}> * REST MASS 0.9383 GeV
OPAL{0}> * MOMENTUM 0.3746
OPAL{0}> * NPART 3.277e+04
OPAL{0}> * **********************************************************************************
OPAL{0}>
OPAL{0}> * ************* F I E L D S O L V E R **********************************************
OPAL{0}> * FIELDSOLVER FS1
OPAL{0}> * TYPE AMR_MG
OPAL{0}> * N-PROCESSORS 4
OPAL{0}> * MX 32
OPAL{0}> * MY 32
OPAL{0}> * MT 32
OPAL{0}> * BBOXINCR 20
OPAL{0}> * AMR_MAXLEVEL 1
OPAL{0}> * AMR_REFX 2
OPAL{0}> * AMR_REFY 2
OPAL{0}> * AMR_REFZ 2
OPAL{0}> * AMR_MAXGRIDX 32
OPAL{0}> * AMR_MAXGRIDY 32
OPAL{0}> * AMR_MAXGRIDZ 32
OPAL{0}> * AMR_BFX 16
OPAL{0}> * AMR_BFY 16
OPAL{0}> * AMR_BFZ 16
OPAL{0}> * AMR_TAGGING CHARGE_DENSITY
OPAL{0}> * AMR_DENSITY 1e-09
OPAL{0}> * AMR_MAX_NUM_PART 1
OPAL{0}> * AMR_MIN_NUM_PART 1
OPAL{0}> * AMR_DENSITY 1e-09
OPAL{0}> * AMR_SCALING 0.75
OPAL{0}> * AMR_DOMAIN_RATIO ( 1 0.75 0.75 )
OPAL{0}> * ITSOLVER (AMR_MG) SA
OPAL{0}> * AMR_MG_PREC NONE
OPAL{0}> * AMR_MG_REBALANCE 1
OPAL{0}> * AMR_MG_REUSE RAP
OPAL{0}> * AMR_MG_SMOOTHER GS
OPAL{0}> * AMR_MG_NSWEEPS 12
OPAL{0}> * AMR_MG_INTERP PC
OPAL{0}> * AMR_MG_NORM LINF
OPAL{0}> * AMR_MG_TOL 1e-10
OPAL{0}> * AMR_MG_VERBOSE 0
OPAL{0}> * BCFFTX DIRICHLET
OPAL{0}> * BCFFTY DIRICHLET
OPAL{0}> * BCFFTT (deprec.) DIRICHLET
OPAL{0}> * XDIM parallel
OPAL{0}> * YDIM parallel
OPAL{0}> * Z(T)DIM parallel
OPAL{0}>
OPAL{0}> * **********************************************************************************
OPAL{0}>
OPAL{0}> * -------------------------- Adding Cyclotron ---------------------------- *
OPAL{0}>
OPAL{0}> * Bunch global starting position:
OPAL{0}> * RINIT = 2130 [mm]
OPAL{0}> * PHIINIT = 0 [deg]
OPAL{0}> * ZINIT = 0 [mm]
OPAL{0}>
OPAL{0}> * Bunch global starting momenta:
OPAL{0}> * Initial gamma = 1.077
OPAL{0}> * Initial beta = 0.3708
OPAL{0}> * Reference total momentum (beta * gamma) = 399.2 [MCU]
OPAL{0}> * Reference azimuthal momentum (Pt) = 399.2 [MCU]
OPAL{0}> * Reference radial momentum (Pr) = -0.174 [MCU]
OPAL{0}> * Reference axial momentum (Pz) = 0 [MCU]
OPAL{0}>
OPAL{0}> * 8-fold field symmetry
OPAL{0}> * Field map file name = s03av.nar
OPAL{0}> * Type of cyclotron = RING
OPAL{0}> * Radial aperture = 0 ... 10 [m]
OPAL{0}> * Vertical aperture = -10 ... 10 [m]
OPAL{0}> * Number of trimcoils = 0
OPAL{0}> * Harmonic number h = 6
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.3333 [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}> * ------------------------------ RFCavity ------------------------------
OPAL{0}> * Minimal radius of cavity= 1900 [mm]
OPAL{0}> * Maximal radius of cavity= 4500 [mm]
OPAL{0}> * RF frequency (2*pi*f)= 3.182e+08 [rad/s]
OPAL{0}> * RF Field map file name= Cav1.dat
OPAL{0}> * Cavity azimuth position= 35 [deg]
OPAL{0}> * Cavity gap width= 220 [mm]
OPAL{0}> * Cavity Shift distance= 416 [mm]
OPAL{0}> * Initial RF phase (t=0)= 139.4 [deg]
OPAL{0}> * Read cavity voltage profile data
OPAL{0}> * Cavity voltage data read successfully!
OPAL{0}> * ------------------------------ RFCavity ------------------------------
OPAL{0}> * Minimal radius of cavity= 1900 [mm]
OPAL{0}> * Maximal radius of cavity= 4500 [mm]
OPAL{0}> * RF frequency (2*pi*f)= 3.182e+08 [rad/s]
OPAL{0}> * RF Field map file name= Cav1.dat
OPAL{0}> * Cavity azimuth position= 125 [deg]
OPAL{0}> * Cavity gap width= 220 [mm]
OPAL{0}> * Cavity Shift distance= 416 [mm]
OPAL{0}> * Initial RF phase (t=0)= 319.4 [deg]
OPAL{0}> * Read cavity voltage profile data
OPAL{0}> * Cavity voltage data read successfully!
OPAL{0}> * ------------------------------ RFCavity ------------------------------
OPAL{0}> * Minimal radius of cavity= 1900 [mm]
OPAL{0}> * Maximal radius of cavity= 4500 [mm]
OPAL{0}> * RF frequency (2*pi*f)= 9.547e+08 [rad/s]
OPAL{0}> * RF Field map file name= Cav3.dat
OPAL{0}> * Cavity azimuth position= 170 [deg]
OPAL{0}> * Cavity gap width= 250 [mm]
OPAL{0}> * Cavity Shift distance= 452 [mm]
OPAL{0}> * Initial RF phase (t=0)= 149.4 [deg]
OPAL{0}> * Read cavity voltage profile data
OPAL{0}> * Cavity voltage data read successfully!
OPAL{0}> * ------------------------------ RFCavity ------------------------------
OPAL{0}> * Minimal radius of cavity= 1900 [mm]
OPAL{0}> * Maximal radius of cavity= 4500 [mm]
OPAL{0}> * RF frequency (2*pi*f)= 3.182e+08 [rad/s]
OPAL{0}> * RF Field map file name= Cav1.dat
OPAL{0}> * Cavity azimuth position= 215 [deg]
OPAL{0}> * Cavity gap width= 220 [mm]
OPAL{0}> * Cavity Shift distance= 416 [mm]
OPAL{0}> * Initial RF phase (t=0)= 139.4 [deg]
OPAL{0}> * Read cavity voltage profile data
OPAL{0}> * Cavity voltage data read successfully!
OPAL{0}> * ------------------------------ RFCavity ------------------------------
OPAL{0}> * Minimal radius of cavity= 1900 [mm]
OPAL{0}> * Maximal radius of cavity= 4500 [mm]
OPAL{0}> * RF frequency (2*pi*f)= 3.182e+08 [rad/s]
OPAL{0}> * RF Field map file name= Cav1.dat
OPAL{0}> * Cavity azimuth position= 305 [deg]
OPAL{0}> * Cavity gap width= 220 [mm]
OPAL{0}> * Cavity Shift distance= 416 [mm]
OPAL{0}> * Initial RF phase (t=0)= 319.4 [deg]
OPAL{0}> * Read cavity voltage profile data
OPAL{0}> * Cavity voltage data read successfully!
OPAL{0}> * ---------------------------------------------------
OPAL{0}> * The selected Beam line elements are :
OPAL{0}> * -> Cyclotron
OPAL{0}> * -> RFCavity
OPAL{0}> * -> RFCavity
OPAL{0}> * -> RFCavity
OPAL{0}> * -> RFCavity
OPAL{0}> * -> RFCavity
OPAL{0}> * ---------------------------------------------------
OPAL{0}> * 4th order Runge-Kutta integrator
OPAL{0}>
OPAL{0}> * Store beam to H5 file for multibunch simulation ... Done.
OPAL{0}> * Before repartition: min local particle number 8120 max local particle number: 8255
OPAL{0}> * After repartition: min local particle number 8120 max local particle number: 8255
OPAL{0}>
OPAL{0}> * *********************** Bunch information in local frame: ************************
OPAL{0}> * ************** B U N C H *********************************************************
OPAL{0}> * NP = 32768
OPAL{0}> * Qtot = 39.487 [pC] Qi = 1.205 [fC]
OPAL{0}> * Ekin = 72.000 [MeV] dEkin = 1.833 [eV]
OPAL{0}> * rmax = ( 5.92953 , 5.98517 , 5.96754 ) [mm]
OPAL{0}> * rmin = ( -5.98308 , -6.00790 , -5.98716 ) [mm]
OPAL{0}> * rms beam size = ( 1.94725 , 1.95178 , 1.98421 ) [mm]
OPAL{0}> * rms momenta = ( 1.42075e-08 , 5.26836e-09 , 0.00000e+00 ) [beta gamma]
OPAL{0}> * mean position = ( 0.00000 , -0.00000 , -0.00000 ) [um]
OPAL{0}> * mean momenta = ( -1.66158e-17 , 3.99202e-01 , 0.00000e+00 ) [beta gamma]
OPAL{0}> * rms emittance = ( 6.93021e-11 , 2.57581e-11 , 0.00000e+00 ) (not normalized)
OPAL{0}> * rms correlation = ( 2.01014e-03 , 3.99906e-04 , 0.00000e+00 )
OPAL{0}> * hr = ( 600.78992 , 450.59244 , 450.59244 ) [um]
OPAL{0}> * dh = 2.00000e+01 [%]
OPAL{0}> * t = 0.000 [fs] dT = 54.843 [ps]
OPAL{0}> * spos = 0.000 [um]
OPAL{0}> * **********************************************************************************
OPAL{0}>
OPAL{0}>
OPAL{0}> * Phase space dump 0 (local frame) at integration step 0
OPAL{0}> * T = 0 ns, Live Particles: 32768
OPAL{0}> * E = 72 MeV, beta * gamma = 0.3992
OPAL{0}> * Bunch position: R = 2130 mm, Theta = 0.00027079 Deg, Z = 0.0017547 mm
OPAL{0}> * Local Azimuth = 0.024973 Deg, Local Elevation = 0 Deg
OPAL{0}>
OPAL{0}> * *********************** Bunch information in global frame: ***********************
OPAL{0}> * ************** B U N C H *********************************************************
OPAL{0}> * NP = 32768
OPAL{0}> * Qtot = 39.487 [pC] Qi = 1.205 [fC]
OPAL{0}> * Ekin = 72.000 [MeV] dEkin = 8.978 [eV]
OPAL{0}> * rmax = ( 2.13595 , 0.00600 , 0.00597 ) [m]
OPAL{0}> * rmin = ( 2.12404 , -0.00600 , -0.00599 ) [m]
OPAL{0}> * rms beam size = ( 1.94724 , 1.95179 , 1.98421 ) [mm]
OPAL{0}> * rms momenta = ( 1.42075e-08 , 2.58096e-08 , 0.00000e+00 ) [beta gamma]
OPAL{0}> * mean position = ( 2.13002 , 0.00001 , 0.00000 ) [m]
OPAL{0}> * mean momenta = ( -1.74000e-04 , 3.99202e-01 , 0.00000e+00 ) [beta gamma]
OPAL{0}> * rms emittance = ( 6.93019e-11 , 1.26189e-10 , 0.00000e+00 ) (not normalized)
OPAL{0}> * rms correlation = ( 2.01319e-03 , 7.99693e-05 , 0.00000e+00 )
OPAL{0}> * hr = ( 600.78992 , 450.59244 , 450.59244 ) [um]
OPAL{0}> * dh = 2.00000e+01 [%]
OPAL{0}> * t = 0.000 [fs] dT = 54.843 [ps]
OPAL{0}> * spos = 0.000 [um]
OPAL{0}> * **********************************************************************************
OPAL{0}>
OPAL{0}> * Beginning of this run is at t = 0 [ns]
OPAL{0}> * The time step is set to dt = 0.32906 [ns]
OPAL{0}> * MBM: Time interval between neighbour bunches is set to 118.46[ns]
OPAL{0}> * MBM: The particles energy bin reset frequency is set to 100
OPAL{0}> * Single particle trajectory dump frequency is set to 1000000
OPAL{0}> * The frequency to solve space charge fields is set to 2
OPAL{0}> * The repartition frequency is set to 5
OPAL{0}>
OPAL{0}> * --------------------------------- Start tracking ------------------------------------ *
******* WARNING *******
No repartitioning necessary: partitions were left unchanged by the repartitioner
******* WARNING *******
Only one partition: Skip call to the repartitioner.
--------------------------------------------------------------------------------
--- Multigrid Summary ---
--------------------------------------------------------------------------------
Number of levels = 3
Operator complexity = 1.55
Smoother complexity = 1.71
Cycle type = V
level rows nnz nnz/row c ratio procs
0 32768 223232 6.81 4
1 3952 117206 29.66 8.29 4
2 118 5066 42.93 33.49 1
Smoother (level 0) both : "Ifpack2::Chebyshev": {Initialized: true, Computed: true, "Ifpack2::Details::Chebyshev":{degree: 3, lambdaMax: 1.81374, alpha: 20, lambdaMin: 0.0906868, boost factor: 1.1}, Global matrix dimensions: [32768, 32768], Global nnz: 223232}
Smoother (level 1) both : "Ifpack2::Chebyshev": {Initialized: true, Computed: true, "Ifpack2::Details::Chebyshev":{degree: 3, lambdaMax: 1.22302, alpha: 20, lambdaMin: 0.0611509, boost factor: 1.1}, Global matrix dimensions: [3952, 3952], Global nnz: 117206}
Smoother (level 2) pre : "Ifpack2::Relaxation": {Initialized: true, Computed: true, Type: Gauss-Seidel, sweeps: 4, damping factor: 1, use l1: 1, l1 eta: 1.5, Global matrix dimensions: [118, 118], Global nnz: 5066}
Smoother (level 2) post : no smoother
OPAL{0}>
OPAL{0}> *** Finished turn 1, Total number of live particles: 32768
OPAL{0}>
OPAL{0}> * Phase space dump 1 (local frame) at integration step 360
OPAL{0}> * T = 118.46 ns, Live Particles: 32768
OPAL{0}> * E = 72 MeV, beta * gamma = 0.40585
OPAL{0}> * Bunch position: R = 2154.5 mm, Theta = -0.10697 Deg, Z = 0.0014297 mm
OPAL{0}> * Local Azimuth = -0.36698 Deg, Local Elevation = 2.3584e-05 Deg
OPAL{0}> * MBM: Injecting a new bunch ...
OPAL{0}>
OPAL{0}> * Read beam from H5 file for multibunch simulation ... Done.
Ippl{0}> Before reset Bin:
Ippl{0}> Bin 0 : particle number = 131072 gamma = 0.269805
Ippl{0}> minimal beta*gamma = 0.399202
Ippl{0}> After reset Bin:
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
OPAL{0}> * MBM: Bunch 2 injected, total particle number = 65536
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
OPAL{0}> * Initialization of all levels
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 0
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
OPAL{0}> * Initialization done.
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07798
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Before reset Bin:
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> minimal beta*gamma = 0.400842
Ippl{0}> After reset Bin:
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862
OPAL{0}> * Start regriding:
OPAL{0}> * Old finest level: 1
OPAL{0}[2]> * Start tagging of level 0
OPAL{0}[2]> * Finished tagging of level 0
OPAL{0}> * New finest level: 1
OPAL{0}> * Finished regriding
Ippl{0}> Bin 0 : particle number = 65536 gamma = 1.07862