Compare commits

...

2 Commits

Author SHA1 Message Date
colin 57e12cbe32 app: stacked_cores: 50-xx: explore some more runs
got the amplification up to a bit over 0.3...
2022-10-26 08:08:15 -07:00
colin 6af2d1d9e3 app: stacked_cores: complete some more 50-xx runs 2022-10-25 15:29:33 -07:00
3 changed files with 209 additions and 26 deletions

View File

@ -191,16 +191,41 @@ _48xx_study = from_params(
_49xx_study = from_params(
[
SimParams47(5, 1, 400, "1e10"),
# TODO
SimParams50(5e2, 4e4, 4000, 200, 5, 1, 400, "2e10"),
# y(0)=0.09, y(1)=0.17, slope0=0.07. "best" comparable 47-xx sim
# SimParams47(5, 1, 400, "1e10"),
# y(0)=0.81, y(1)=0.99, slope0=0.50 to 0.30
SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "2e10"),
# y(0)=0.60, y(1)=0.89, slope0=0.40 to 0.29
SimParams50(2e3, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# y(0)=0.73, y(1)=0.97, slope0=0.44 to 0.32
SimParams50(2e3, 2e4, 2000, 100, 5, 1, 400, "2e10"),
# y(0)=0.62, y(1)=0.80, slope0=0.23 to 0.22
SimParams50(5e2, 1e4, 1000, 50, 5, 1, 400, "1e10"),
# y(0)=0.65, y(1)=0.80, slope0>0.18
SimParams50(5e2, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# too few measurements
# SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "5e9"),
# SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "2e10"),
# SimParams50(2e3, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# SimParams50(5e3, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# SimParams50(1e4, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# y(0)=0.63, y(1)=0.89, slope0>0.27
SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# y(0)=0.46, y(1)=0.66, slope0>0.23
SimParams50(2e3, 2e4, 2000, 100, 5, 1, 400, "5e9"),
# y(0)=0.50, y(1)=0.74, slope0>0.25
SimParams50(5e3, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# y(0)=0.61, y(1)=0.86, slope0>0.25
SimParams50(5e3, 2e4, 2000, 100, 5, 1, 400, "2e10"),
# y(0)=0.39, y(1)=0.48, slope0>0.11
# SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "5e9"),
# y(0)=0.40, y(1)=0.56, slope0>0.15
# SimParams50(1e4, 2e4, 2000, 100, 5, 1, 400, "1e10"),
# y(0)=0.51, y(1)=0.68, slope0>0.16
# SimParams50(1e4, 2e4, 2000, 100, 5, 1, 400, "2e10"),
# y(0)=0.30, y(1)=0.40, slope0>0.09
# SimParams50(2e4, 2e4, 2000, 100, 5, 1, 400, "1e10"),
]
)

View File

@ -431,9 +431,15 @@ sims = [
(SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "5e9"), None, 20000),
(SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "1e10"), None, 20000),
(SimParams50(1e3, 2e4, 2000, 100, 5, 1, 400, "2e10"), None, 20000),
(SimParams50(2e3, 2e4, 2000, 100, 5, 1, 400, "5e9"), None, 20000),
(SimParams50(2e3, 2e4, 2000, 100, 5, 1, 400, "1e10"), None, 20000),
(SimParams50(2e3, 2e4, 2000, 100, 5, 1, 400, "2e10"), None, 20000),
(SimParams50(5e3, 2e4, 2000, 100, 5, 1, 400, "1e10"), None, 20000),
(SimParams50(5e3, 2e4, 2000, 100, 5, 1, 400, "2e10"), None, 20000),
(SimParams50(1e4, 2e4, 2000, 100, 5, 1, 400, "1e10"), None, 20000),
(SimParams50(1e4, 2e4, 2000, 100, 5, 1, 400, "2e10"), None, 20000),
(SimParams50(2e4, 2e4, 2000, 100, 5, 1, 400, "1e10"), None, 20000),
(SimParams50(2e4, 2e4, 4000, 200, 5, 1, 400, "2e10"), None, 20000),
]
measurements = { real.machine_name: [human, norm, None, None] for (real, human, norm) in sims }
@ -4902,21 +4908,151 @@ set_meas("50-0.0004rad-1000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupli
Piecewise(
[
[ -17587, 10360 ], # -1.000, M1=-11417
[ -17437, 10523 ], # -0.350, M1=-11275
[ -17413, 10547 ], # -0.300, M1=-11253
[ -17379, 10577 ], # -0.250, M1=-11227
[ -17331, 10616 ], # -0.200, M1=-11192
[ -17227, 10685 ], # -0.150, M1=-11129
[ -17045, 10777 ], # -0.130, M1=-11046
[ -16825, 10886 ], # -0.120, M1=-10946
[ -16402, 11093 ], # -0.110, M1=-10756
[ -15273, 11624 ], # -0.100, M1=-10256
[ -12326, 12652 ], # -0.090, M1=-9237
[ -9209, 13601 ], # -0.080, M1=-8208
[ -6063, 14439 ], # -0.070, M1=-7208
[ -2898, 15175 ], # -0.060, M1=-6231
[ 274, 15787 ], # -0.050, M1=-5270
[ 3417, 16146 ], # -0.040, M1=-4324
[ 6421, 16260 ], # -0.030, M1=-3416
[ 8677, 16324 ], # -0.020, M1=-2721
[ 9570, 16350 ], # -0.010, M1=-2443
[ 10229, 16369 ], # 0.000, M1=-2240
[ 14759, 16483 ], # 0.100, M1= -875
[ 16335, 16527 ], # 0.200, M1= -354
[ 16771, 16545 ], # 1.000, M1= -120
]
)
""")
set_meas("50-0.0004rad-2000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-5e9-drive", """
Piecewise(
[
[ -17402, -1620 ], # -1.000, M1=-13338
[ 4620, 3431 ], # 0.000, M1=-6623
[ 16729, 5549 ], # 1.000, M1=-3374
]
)
""")
set_meas("50-0.0004rad-2000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-1e10-drive", """
Piecewise(
[
[ -17555, 2961 ], # -1.000, M1=-12179
[ -17243, 3193 ], # -0.200, M1=-11947
[ -14966, 4108 ], # -0.130, M1=-11033
[ -12700, 4794 ], # -0.120, M1=-10348
[ -10348, 5479 ], # -0.110, M1=-9664
[ -5593, 6821 ], # -0.090, M1=-8300
[ -3210, 7517 ], # -0.080, M1=-7588
[ -846, 8215 ], # -0.070, M1=-6872
[ 1466, 8878 ], # -0.060, M1=-6186
[ 7967, 10735 ], # 0.000, M1=-4212
[ 16755, 13090 ], # 1.000, M1=-1516
]
)
""")
set_meas("50-0.0004rad-2000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-2e10-drive", """
Piecewise(
[
[ -17730, 7307 ], # -1.000, M1=-11351
[ -17214, 7709 ], # -0.090, M1=-10960
[ -16931, 7834 ], # -0.080, M1=-10838
[ -16177, 8165 ], # -0.070, M1=-10510
[ -14947, 8667 ], # -0.065, M1=-10009
[ -12688, 9434 ], # -0.060, M1=-9233
[ -7938, 10941 ], # -0.050, M1=-7656
[ -3128, 12304 ], # -0.040, M1=-6135
[ 1707, 13513 ], # -0.030, M1=-4660
[ 6447, 14550 ], # -0.020, M1=-3240
[ 9775, 15181 ], # -0.010, M1=-2241
[ 10966, 15386 ], # 0.000, M1=-1880
[ 14177, 15832 ], # 0.050, M1= -937
[ 16107, 16041 ], # 0.100, M1= -310
[ 16789, 16128 ], # 1.000, M1= 65
]
)
""")
set_meas("50-0.0004rad-5000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-1e10-drive", """
Piecewise(
[
[ -17874, -336 ], # -1.000, M1=-12353
[ -17709, -216 ], # -0.350, M1=-12199
[ -17685, -199 ], # -0.300, M1=-12177
[ -17658, -180 ], # -0.250, M1=-12153
[ -17621, -156 ], # -0.200, M1=-12122
[ -17563, -121 ], # -0.150, M1=-12078
[ -17436, -58 ], # -0.100, M1=-11997
[ -8261, 2237 ], # -0.050, M1=-9027
[ 8890, 6091 ], # 0.000, M1=-4131
[ 16207, 7754 ], # 0.100, M1=-2005
[ 16491, 7861 ], # 0.200, M1=-1869
[ 16625, 7929 ], # 0.350, M1=-1781
[ 16811, 8050 ], # 1.000, M1=-1626
]
)
""")
set_meas("50-0.0004rad-5000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-2e10-drive", """
Piecewise(
[
[ -18082, 3583 ], # -1.000, M1=-11466
[ -17630, 3909 ], # -0.070, M1=-11064
[ 10894, 10687 ], # 0.000, M1=-2354
[ 16869, 12160 ], # 1.000, M1= -258
]
)
""")
set_meas("50-0.0004rad-10000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-1e10-drive", """
Piecewise(
[
[ -18342, -3636 ], # -1.000, M1=-12916
[ -18115, -3533 ], # -0.350, M1=-12746
[ -18084, -3519 ], # -0.300, M1=-12723
[ -18051, -3503 ], # -0.250, M1=-12695
[ -18014, -3484 ], # -0.200, M1=-12662
[ -17961, -3460 ], # -0.150, M1=-12623
[ -17887, -3424 ], # -0.100, M1=-12563
[ -17632, -3334 ], # -0.050, M1=-12415
[ 8753, 735 ], # 0.000, M1=-5331
[ 16488, 1959 ], # 0.100, M1=-3107
[ 16640, 2024 ], # 0.200, M1=-2992
[ 16754, 2075 ], # 0.350, M1=-2899
[ 17009, 2174 ], # 1.000, M1=-2717
]
)
""")
set_meas("50-0.0004rad-10000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-2e10-drive", """
Piecewise(
[
[ -18603, -61 ], # -1.000, M1=-11774
[ 10164, 4652 ], # 0.000, M1=-3814
[ 17128, 6016 ], # 1.000, M1=-1449
]
)
""")
set_meas("50-0.0004rad-20000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-1e10-drive", """
Piecewise(
[
[ -19216, -6910 ], # -1.000, M1=-13886
[ 7843, -4486 ], # 0.000, M1=-7707
[ 17714, -3467 ], # 1.000, M1=-5023
]
)
""")
if __name__ == '__main__': main()

View File

@ -5430,9 +5430,15 @@ fn main() {
for init_set in [
&[
// targeted
// -0.22,
// -0.18,
// -0.27,
// -0.05,
// -0.02,
// -0.03,
// -0.01,
// -0.04,
// -0.06,
// -0.065,
// -0.055,
// -0.045,
][..],
&[
// establish the domain/range
@ -5444,30 +5450,46 @@ fn main() {
][..],
&[
-0.20,
-0.35,
-0.25,
-0.30,
-0.15,
-0.10,
0.20,
0.10,
][..],
&[
-0.05,
0.35,
0.05,
-0.15,
-0.25,
-0.35,
0.20,
-0.30,
][..],
] {
for (ctl_cond, coupling_cond, clock_duration, clock_decay, coupling_loops, s0_loops, s_major, cur_flt) in [
(1e4, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
(5e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
(5e2, 4e4, ps(4000), ps(200), 5, 1, um(400), 2e10),
// worth pursuing in more detail:
// y(0)=0.61, y(1)=0.86, slope0>0.25
(5e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 2e10),
// y(0)=0.63, y(1)=0.89, slope0>0.27
(1e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
// y(0)=0.81, y(1)=0.99, slope0=0.50 to 0.30
(1e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 2e10),
// y(0)=0.73, y(1)=0.97, slope0=0.44 to 0.32
(2e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 2e10),
// y(0)=0.60, y(1)=0.89, slope0=0.40 to 0.29
(2e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
(1e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 5e9),
(1e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
(1e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 2e10),
// (1e4, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
// (5e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
// (2e4, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
// (2e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 5e9),
// (1e4, 2e4, ps(2000), ps(100), 5, 1, um(400), 2e10),
(5e2, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
// y(0)=0.62, slope(0)=0.23, y(1)=0.80
(5e2, 1e4, ps(1000), ps(50), 5, 1, um(400), 1e10),
// (1e3, 2e4, ps(2000), ps(100), 5, 1, um(400), 5e9),
// (5e2, 2e4, ps(2000), ps(100), 5, 1, um(400), 1e10),
// // y(0)=0.62, slope(0)=0.23, y(1)=0.80
// (5e2, 1e4, ps(1000), ps(50), 5, 1, um(400), 1e10),
// old 47-xx runs
// (5e2, 5e3, ps(1000), ps(50), 5, 1, um(400), 1e10),