Compare commits

...

4 Commits

3 changed files with 233 additions and 12 deletions

View File

@ -133,6 +133,8 @@ all_viable_inverters = filter_meas(viable_inverter=True);
# inverters with steepest starting slope
inverters_with_steepest_slope0 = [(p, get_meas(p)) for p in
[
SimParams48(5e2, 4e4, 4000, 200, 9, 1, 400, "1e10"),
# SimParams48(5e2, 2e4, 2000, 100, 9, 1, 400, "1e10"),
SimParams48(5e2, 1e4, 1000, 50, 9, 1, 400, "1e10"),
SimParams48(1e3, 1e4, 2000, 100, 9, 1, 400, "1e10"),
SimParams41(36, 1, 1200, "5e9"),
@ -144,23 +146,30 @@ inverters_with_steepest_slope0 = [(p, get_meas(p)) for p in
_48xx_all = filter_meas(run="48")
_48xx_study = [(p, get_meas(p)) for p in
[
# y0=0.62, slope0=1.4
# T0: y0=0.62, slope0=1.4 until x=0.20
SimParams48(5e2, 1e4, 1000, 50, 9, 1, 400, "1e10"),
# y0=0.55, slope0=1.2
# y0=0.55, slope0=1.2 until x=0.30
SimParams48(5e2, 5e3, 1000, 50, 9, 1, 400, "1e10"),
# # y0=0.60, slope0=0.96
# SimParams48(5e2, 5e3, 1000, 50, 9, 1, 400, "2e10"),
# y0=0.57, slope0=1.25
# y0=0.57, slope0=1.25 until x=0.30
SimParams48(1e3, 1e4, 2000, 100, 9, 1, 400, "1e10"),
# y0=0.47, slope0=1.1
# y0=0.47, slope0=1.1 until x=0.30
SimParams48(1e3, 1e4, 2000, 100, 9, 1, 400, "5e9"),
# y0=0.52, slope0=1.3, but narrower than the 1e10 variant
# y0=0.52, slope0=1.3 until x=0.20
SimParams48(5e2, 1e4, 1000, 50, 9, 1, 400, "5e9"),
# TODO:
# y0=0.67, slope0=1.6 until x=0.20
SimParams48(5e2, 2e4, 2000, 100, 9, 1, 400, "1e10"),
# y0=0.57, slope0=1.3 until x=0.20
SimParams48(5e2, 2e4, 2000, 100, 9, 1, 400, "5e9"),
# y0=0.70, slope0=1.7 until x=0.15
SimParams48(5e2, 4e4, 4000, 200, 9, 1, 400, "1e10"),
# TODO
SimParams48(5e2, 4e4, 2000, 100, 9, 1, 400, "5e9"),
# y0=0.64, slope0=1.4 until x=0.20
SimParams48(2e2, 1e4, 1000, 50, 9, 1, 400, "1e10"),
]
]
@ -208,5 +217,5 @@ for (params, curve) in of_interest[:10]:
if not params.is_inverter:
curve = curve.logically_inverted()
if params.is_inverter or True:
curve.plot_equilibrium(title = f"{params.human_name} equilibrium")
# if params.is_inverter or True:
# curve.plot_equilibrium(title = f"{params.human_name} equilibrium")

View File

@ -384,6 +384,7 @@ sims = [
# control conductivity, coupling cond, clock decay (ps), clock duration, coupling loops, asymmetric s0 loops, radius, drive current
(SimParams48(5e2, 5e3, 1000, 50, 9, 1, 400, "1e10"), None, 20000),
(SimParams48(5e2, 5e3, 1000, 50, 9, 1, 400, "2e10"), None, 20000),
(SimParams48(2e2, 1e4, 1000, 50, 9, 1, 400, "1e10"), None, 20000),
(SimParams48(5e2, 1e4, 1000, 50, 9, 1, 400, "3e9"), None, 20000),
(SimParams48(5e2, 1e4, 1000, 50, 9, 1, 400, "5e9"), None, 20000),
(SimParams48(5e2, 1e4, 1000, 50, 9, 1, 400, "1e10"), None, 20000),
@ -395,6 +396,8 @@ sims = [
(SimParams48(5e2, 2e4, 4000, 200, 9, 1, 400, "1e10"), None, 20000),
(SimParams48(2e3, 2e4, 4000, 200, 9, 1, 400, "1e9"), None, 20000),
(SimParams48(2e3, 2e4, 4000, 200, 9, 1, 400, "2e9"), None, 20000),
(SimParams48(5e2, 4e4, 2000, 100, 9, 1, 400, "5e9"), None, 20000),
(SimParams48(5e2, 4e4, 4000, 200, 9, 1, 400, "1e10"), None, 20000),
(SimParams48(5e3, 5e4,10000, 500, 9, 1, 400, "5e8"), None, 20000),
(SimParams48(5e3, 5e4,10000, 500, 9, 1, 400, "1e9"), None, 20000),
]
@ -4144,6 +4147,41 @@ Piecewise(
)
""")
set_meas("48-0.0004rad-200ctl_cond-10000coupling_cond-1000ps-50ps-9ctl-9coupling-3_1_winding-1e10-drive", """
Piecewise(
[
[ -18665, 2331 ], # -1.00
[ -18367, 2808 ], # -0.35
[ -18262, 2964 ], # -0.30
[ -18225, 3017 ], # -0.29
[ -18177, 3086 ], # -0.28
[ -18103, 3188 ], # -0.27
[ -17962, 3383 ], # -0.26
[ -17758, 3673 ], # -0.25
[ -17482, 4070 ], # -0.24
[ -17147, 4571 ], # -0.23
[ -16723, 5206 ], # -0.22
[ -16212, 5981 ], # -0.21
[ -15523, 7019 ], # -0.20
[ -14381, 8713 ], # -0.19
[ -12852, 10876 ], # -0.18
[ -11271, 12966 ], # -0.17
[ -9673, 14779 ], # -0.16
[ -8061, 16111 ], # -0.15
[ -6445, 16527 ], # -0.14
[ -4827, 16712 ], # -0.13
[ -3209, 16747 ], # -0.12
[ -1592, 16771 ], # -0.11
[ 23, 16792 ], # -0.10
[ 7948, 16834 ], # -0.05
[ 11334, 16838 ], # 0.00
[ 12855, 16839 ], # 0.10
[ 17046, 16842 ], # 0.35
[ 17352, 16841 ], # 1.00
]
)
""")
set_meas("48-0.0004rad-500ctl_cond-10000coupling_cond-1000ps-50ps-9ctl-9coupling-3_1_winding-3e9-drive", """
Piecewise(
[
@ -4254,11 +4292,95 @@ Piecewise(
)
""")
set_meas("48-0.0004rad-500ctl_cond-20000coupling_cond-2000ps-100ps-9ctl-9coupling-3_1_winding-5e9-drive", """
Piecewise(
[
[ -18385, 377 ], # -1.00
[ -17714, 1264 ], # -0.40
[ -16717, 2796 ], # -0.35
[ -16429, 3229 ], # -0.34
[ -16098, 3715 ], # -0.33
[ -15695, 4277 ], # -0.32
[ -15192, 4953 ], # -0.31
[ -14391, 5983 ], # -0.30
[ -13426, 7231 ], # -0.29
[ -12434, 8475 ], # -0.28
[ -11434, 9620 ], # -0.27
[ -10426, 10668 ], # -0.26
[ -9409, 11790 ], # -0.25
[ -8390, 12770 ], # -0.24
[ -7369, 13639 ], # -0.23
[ -6347, 14376 ], # -0.22
[ -4301, 15224 ], # -0.20
[ -1233, 15693 ], # -0.17
[ 811, 15853 ], # -0.15
[ 3872, 16043 ], # -0.12
[ 5900, 16128 ], # -0.10
[ 8860, 16230 ], # -0.07
[ 10290, 16275 ], # -0.05
[ 11468, 16310 ], # -0.02
[ 11672, 16316 ], # 0.00
[ 11988, 16323 ], # 0.05
[ 12093, 16326 ], # 0.07
[ 12246, 16329 ], # 0.10
[ 12346, 16331 ], # 0.12
[ 12750, 16340 ], # 0.20
[ 15902, 16419 ], # 0.35
[ 17215, 16453 ], # 1.00
]
)
""")
set_meas("48-0.0004rad-500ctl_cond-20000coupling_cond-2000ps-100ps-9ctl-9coupling-3_1_winding-1e10-drive", """
Piecewise(
[
[ -18644, 2982 ], # -1.00
[ -18484, 3261 ], # -0.50
[ -18419, 3377 ], # -0.40
[ -18380, 3451 ], # -0.35
[ -18362, 3488 ], # -0.33
[ -18343, 3528 ], # -0.31
[ -18332, 3550 ], # -0.30
[ -18319, 3573 ], # -0.29
[ -18286, 3631 ], # -0.27
[ -18264, 3667 ], # -0.26
[ -18239, 3710 ], # -0.25
[ -18209, 3760 ], # -0.24
[ -18164, 3831 ], # -0.23
[ -18103, 3929 ], # -0.22
[ -17989, 4104 ], # -0.21
[ -17769, 4450 ], # -0.20
[ -17442, 4975 ], # -0.19
[ -17018, 5680 ], # -0.18
[ -16449, 6621 ], # -0.17
[ -15674, 7898 ], # -0.16
[ -14299, 10123 ], # -0.15
[ -12360, 13165 ], # -0.14
[ -10352, 15693 ], # -0.13
[ -8317, 16251 ], # -0.12
[ -6271, 16512 ], # -0.11
[ -4221, 16752 ], # -0.10
[ 1920, 16832 ], # -0.07
[ 5972, 16841 ], # -0.05
[ 10908, 16845 ], # -0.02
[ 11760, 16845 ], # 0.00
[ 12557, 16845 ], # 0.05
[ 12782, 16845 ], # 0.07
[ 13133, 16845 ], # 0.10
[ 13402, 16846 ], # 0.12
[ 16662, 16848 ], # 0.20
[ 17065, 16847 ], # 0.35
[ 17348, 16846 ], # 1.00
]
)
""")
set_meas("48-0.0004rad-1000ctl_cond-10000coupling_cond-2000ps-100ps-9ctl-9coupling-3_1_winding-3e9-drive", """
Piecewise(
[
[ -16999, -3777 ], # -1.00
[ -13768, -1435 ], # -0.35
[ -9559, 1353 ], # -0.30
[ -5278, 3430 ], # -0.25
[ -4420, 3777 ], # -0.24
[ -3562, 4109 ], # -0.23
@ -4387,6 +4509,7 @@ Piecewise(
[
[ -16695, -12962 ], # -1.00
[ 189, -7932 ], # -0.35
[ 893, -7705 ], # -0.30
[ 1302, -7574 ], # -0.25
[ 1728, -7439 ], # -0.20
[ 1996, -7355 ], # -0.17
@ -4447,11 +4570,36 @@ Piecewise(
)
""")
set_meas("48-0.0004rad-500ctl_cond-40000coupling_cond-4000ps-200ps-9ctl-9coupling-3_1_winding-1e10-drive", """
Piecewise(
[
[ -18503, 3956 ], # -1.00
[ -18244, 4450 ], # -0.35
[ -18105, 4724 ], # -0.25
[ -17635, 5520 ], # -0.20
[ -16347, 7816 ], # -0.17
[ -15588, 9160 ], # -0.16
[ -14219, 11518 ], # -0.15
[ -12285, 14698 ], # -0.14
[ -10278, 16161 ], # -0.13
[ -4150, 16834 ], # -0.10
[ 6054, 16875 ], # -0.05
[ 12063, 16878 ], # 0.00
[ 12697, 16879 ], # 0.05
[ 13136, 16879 ], # 0.10
[ 16644, 16881 ], # 0.20
[ 17028, 16882 ], # 0.35
[ 17304, 16881 ], # 1.00
]
)
""")
set_meas("48-0.0004rad-5000ctl_cond-50000coupling_cond-10000ps-500ps-9ctl-9coupling-3_1_winding-5e8-drive", """
Piecewise(
[
[ -16608, -12761 ], # -1.00
[ 107, -7801 ], # -0.35
[ 923, -7540 ], # -0.30
[ 1350, -7404 ], # -0.25
[ 1793, -7265 ], # -0.20
[ 2069, -7179 ], # -0.17
@ -4478,6 +4626,7 @@ Piecewise(
[
[ -16845, -3852 ], # -1.00
[ -14379, -2012 ], # -0.35
[ -10136, 1101 ], # -0.30
[ -5768, 3505 ], # -0.25
[ -1389, 5290 ], # -0.20
[ 1219, 6176 ], # -0.17

View File

@ -5130,6 +5130,65 @@ fn main() {
// -0.03,
// -0.02,
// -0.01,
// targeted for (5e2, 4e4, ps(4000), ps(200), 9, 1, um(400), 1e10)
// -0.17,
// -0.13,
// -0.16,
// -0.14,
// -0.18,
// -0.12,
// -0.19,
// -0.11,
// -0.21,
// -0.22,
// -0.23,
// targeted for (5e2, 2e4, ps(2000), ps(100), 9, 1, um(400), 1e10)
// -0.16,
// -0.18,
// -0.19,
// -0.21,
// -0.22,
// -0.23,
// -0.24,
// -0.30,
// -0.14,
// -0.13,
// -0.12,
// -0.11,
// fallthrough
// targeted for (5e2, 2e4, ps(2000), ps(100), 9, 1, um(400), 5e9)
// -0.23,
// -0.27,
// -0.29,
// -0.31,
// -0.33,
// -0.40,
// -0.22,
// -0.24,
// -0.26,
// -0.28,
// -0.30,
// -0.32,
// -0.34,
// -0.19,
// -0.50,
// targeted for (2e2, 1e4, ps(1000), ps(50), 9, 1, um(400), 1e10)
// -0.19,
// -0.21,
// -0.18,
// -0.22,
// -0.17,
// -0.23,
// -0.16,
// -0.24,
// -0.26,
// -0.27,
// -0.28,
// -0.29,
// -0.14,
// -0.13,
// -0.12,
// -0.11,
][..],
&[
// establish the domain/range
@ -5143,8 +5202,10 @@ fn main() {
-0.05,
-0.35,
-0.25,
-0.15,
-0.30,
0.10,
0.20,
// 0.20,
0.35,
][..],
&[
@ -5191,6 +5252,7 @@ fn main() {
][..],
] {
for (ctl_cond, coupling_cond, clock_duration, clock_decay, coupling_loops, s0_loops, s_major, cur_flt) in [
// (2e2, 1e4, ps(1000), ps(50), 9, 1, um(400), 1e10), // complete. y0=0.64, slope0=1.4
// (5e2, 1e4, ps(1000), ps(50), 9, 1, um(400), 1e10), // complete. y0=0.62, slope0=1.4
// (5e2, 5e3, ps(1000), ps(50), 9, 1, um(400), 1e10), // complete. y0=0.55, slope0=1.2
// (1e3, 1e4, ps(2000), ps(100), 9, 1, um(400), 1e10), // incomplete. y0=0.57, slope0=1.2
@ -5201,9 +5263,10 @@ fn main() {
// (5e2, 1e4, ps(1000), ps(50), 9, 1, um(400), 3e9),
(5e2, 2e4, ps(2000), ps(100), 9, 1, um(400), 1e10), // incomplete
(5e2, 2e4, ps(2000), ps(100), 9, 1, um(400), 5e9), // incomplete
(5e2, 4e4, ps(4000), ps(200), 9, 1, um(400), 1e10), // incomplete
(5e2, 4e4, ps(2000), ps(100), 9, 1, um(400), 5e9),
// (5e2, 4e4, ps(4000), ps(200), 9, 1, um(400), 1e10), // complete. y0=0.70, slope0=1.4
// (5e2, 2e4, ps(2000), ps(100), 9, 1, um(400), 1e10), // complete. slope0=1.6
// (5e2, 2e4, ps(2000), ps(100), 9, 1, um(400), 5e9), // complete. slope0=1.3
(1e3, 1e4, ps(2000), ps(100), 9, 1, um(400), 3e9),
// (2e3, 2e4, ps(4000), ps(200), 9, 1, um(400), 2e9), // too low slope