From e1867ee5412c599b3f69389c62b4113392249bca Mon Sep 17 00:00:00 2001 From: colin Date: Sat, 22 Oct 2022 05:37:34 -0700 Subject: [PATCH] app: stacked_cores: 48-xx: complete a very low control-conductivity run (2e2) --- .../stacked_cores/scripts/plot_inverters.py | 7 +++ .../scripts/stacked_cores_40xx_db.py | 41 ++++++++++++++ crates/applications/stacked_cores/src/main.rs | 55 +++++++++++++------ 3 files changed, 86 insertions(+), 17 deletions(-) diff --git a/crates/applications/stacked_cores/scripts/plot_inverters.py b/crates/applications/stacked_cores/scripts/plot_inverters.py index a01ef6e..8cd98c6 100755 --- a/crates/applications/stacked_cores/scripts/plot_inverters.py +++ b/crates/applications/stacked_cores/scripts/plot_inverters.py @@ -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"), @@ -163,6 +165,11 @@ _48xx_study = [(p, get_meas(p)) for p in 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"), ] ] diff --git a/crates/applications/stacked_cores/scripts/stacked_cores_40xx_db.py b/crates/applications/stacked_cores/scripts/stacked_cores_40xx_db.py index 278ad02..d56772a 100755 --- a/crates/applications/stacked_cores/scripts/stacked_cores_40xx_db.py +++ b/crates/applications/stacked_cores/scripts/stacked_cores_40xx_db.py @@ -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,7 @@ 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), @@ -4145,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( [ @@ -4343,6 +4380,7 @@ 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 @@ -4471,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 @@ -4560,6 +4599,7 @@ 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 @@ -4586,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 diff --git a/crates/applications/stacked_cores/src/main.rs b/crates/applications/stacked_cores/src/main.rs index 46da865..efaeca8 100644 --- a/crates/applications/stacked_cores/src/main.rs +++ b/crates/applications/stacked_cores/src/main.rs @@ -5157,21 +5157,38 @@ fn main() { // -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, + // -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 @@ -5185,8 +5202,10 @@ fn main() { -0.05, -0.35, -0.25, + -0.15, + -0.30, 0.10, - 0.20, + // 0.20, 0.35, ][..], &[ @@ -5233,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 @@ -5243,7 +5263,8 @@ fn main() { // (5e2, 1e4, ps(1000), ps(50), 9, 1, um(400), 3e9), - // (5e2, 4e4, ps(4000), ps(200), 9, 1, um(400), 1e10), // complete. slope0=1.4 + (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