Commit 8f09a6b0 authored by adelmann's avatar adelmann

not yet ready, do not use it in the regressiontests yet

parent f08ced4c
......@@ -384,6 +384,16 @@ RegressionTests/RingCyclotronMTS/reference/RingCyclotronMTS.out.md5 -text
RegressionTests/RingCyclotronMTS/reference/RingCyclotronMTS.stat -text
RegressionTests/RingCyclotronMTS/reference/RingCyclotronMTS.stat.md5 -text
RegressionTests/RingCyclotronMTS/s03av.nar -text
RegressionTests/RingCyclotronMatched/Cav1.dat -text
RegressionTests/RingCyclotronMatched/Cav3.dat -text
RegressionTests/RingCyclotronMatched/RingCyclotronMatched.in -text
RegressionTests/RingCyclotronMatched/RingCyclotronMatched.local -text
RegressionTests/RingCyclotronMatched/RingCyclotronMatched.rt -text
RegressionTests/RingCyclotronMatched/RingCyclotronMatched.sge -text
RegressionTests/RingCyclotronMatched/reference/RingCyclotronMatched.out -text
RegressionTests/RingCyclotronMatched/ring590_bfld.dat -text
RegressionTests/RingCyclotronMatched/ring_tc15_sector.dat -text
RegressionTests/RingCyclotronMatched/s03av.nar -text
RegressionTests/SAAMG-Test-1/Pipe_1m_10cm.h5 -text
RegressionTests/SAAMG-Test-1/SAAMG-Test-1.in -text
RegressionTests/SAAMG-Test-1/SAAMG-Test-1.local -text
......
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
OPTION, ECHO = FALSE;
OPTION, PSDUMPFREQ = 1;
OPTION, SPTDUMPFREQ = 1;
OPTION, PSDUMPEACHTURN = FALSE;
OPTION, PSDUMPLOCALFRAME = FALSE;
OPTION, REPARTFREQ = 20;
OPTION, CZERO = FALSE;
OPTION, TELL = TRUE;
Title,string="OPAL-cycl: test matched distribution PSI Ring";
Edes=.580;
gamma=(Edes+PMASS)/PMASS;
beta=sqrt(1-(1/gamma^2));
gambet=gamma*beta;
P0 = gamma*beta*PMASS;
brho = (PMASS*1.0e9*gambet) / CLIGHT;
//value,{gamma,brho,Edes,beta,gambet};
phi01=139.4281;
phi02=phi01+180.0;
phi04=phi01;
phi05=phi01+180.0;
phi03=phi01+10.0;
volt1st=0.9;
volt3rd=0.9*4.0*0.112;
turns = 0;
nstep = 0;
frequency=50.650;
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", FMLOWE=72.0, FMHIGHE=590.0;
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, DISTRIBUTION=GAUSSMATCHED, LINE=l1, FMAPFN="ring590_bfld.dat",
EX = 1.0e-06,
EY = 1.0e-06,
ET = 1.0e-06;
/*
Dist1:DISTRIBUTION, DISTRIBUTION=GAUSSMATCHED, LINE=l1, FMAPFN="ring_tc15_sector.dat",
EX = 1.0e-06,
EY = 1.0e-06,
ET = 1.0e-06;
*/
Fs1:FIELDSOLVER, FSTYPE=FFT, MX=32, MY=32, MT=16,
PARFFTX=true, PARFFTY=true, PARFFTT=true,
BCFFTX=open, BCFFTY=open, BCFFTT=open;
beam1: BEAM, PARTICLE=PROTON, pc=P0, NPART=1e5, BCURRENT=2.2E-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;
endtrack;
Stop;
#!/bin/bash
mpirun -np 4 $OPAL_EXE_PATH/opal RingCyclotronMatched.in --info 0 --warn 0 2>&1
"One turn RK-4 tracking in PSI Ring with space charge"
out "Ekin" last 1E-15
out "rms beam size" last 1E-04
out "rms emittance" last 1E-07
#!/bin/bash
#$ -cwd
#$ -j y
#$ -pe orte 4
#$ -N RingCyclotron-RT
#$ -v LD_LIBRARY_PATH,OPAL_EXE_PATH,OPENMPI,REG_TEST_DIR
MACHINE_FILE=$TMPDIR/machinefile
awk '/^felsim/ {print $1" slots="$2}' $PE_HOSTFILE > $MACHINE_FILE
cp $MACHINE_FILE machinefile.last
cd $REG_TEST_DIR
OPAL="$OPAL_EXE_PATH/opal RingCyclotronMatched.in --info 0 --warn 0 2>&1"
CMD="$OPENMPI/bin/mpirun -x LD_LIBRARY_PATH -machinefile $MACHINE_FILE -np $NSLOTS $OPAL "
$CMD
This diff is collapsed.
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment