app: stacked_cores: 50-xx: complete a run with high *control* conductivity, and schedule a few more

in fact, M2 is initialized improperly: this actually acts as an
(overpowered) single-clock-cycle inverter.

crates/applications/stacked_cores/scripts/plot_inverters.py

@@ -17,6 +17,10 @@ class SimParamsCascaded(SimParams):
     def human_name(self) -> str:
         return f"Cascade: {self.p1.human_name} -> {self.p2.human_name}"
 
+from_params = lambda l: [
+    (p, get_meas(p)) for p in l if get_meas(p)
+]
+
 
 # plot pre-40xx sims
 for (name, curve) in [
@@ -93,7 +97,7 @@ of_interest = []
 #     ]
 # ]
 
-of_interest += [(p, get_meas(p)) for p in
+of_interest += from_params(
     [
         # SimParams41(4,  3, 400, "2e10"),
         # SimParams41(4,  3, 400, "4e10"),
@@ -126,12 +130,12 @@ of_interest += [(p, get_meas(p)) for p in
         # SimParams41(9,  1, 400, "5e9"),
         # SimParams41(9,  1, 400, "1e10"),
         # SimParams41(9,  1, 400, "2e10"),
-    ] if get_meas(p)
-]
+    ]
+)
 
 all_viable_inverters = filter_meas(viable_inverter=True);
 # inverters with steepest starting slope
-inverters_with_steepest_slope0 = [(p, get_meas(p)) for p in
+inverters_with_steepest_slope0 = from_params(
     [
         SimParams48(5e2, 4e4, 4000, 200, 9, 1, 400, "1e10"),
         # SimParams48(5e2, 2e4, 2000, 100, 9, 1, 400, "1e10"),
@@ -142,9 +146,10 @@ inverters_with_steepest_slope0 = [(p, get_meas(p)) for p in
         SimParams41(16, 2, 800,  "1e10"),
         SimParams41(12, 2, 600,  "1e10"),
     ]
-]
+)
+_47xx_all = filter_meas(run="47")
 _48xx_all = filter_meas(run="48")
-_48xx_study = [(p, get_meas(p)) for p in
+_48xx_study = from_params(
     [
         # T0: y0=0.62, slope0=1.4 until x=0.20
         SimParams48(5e2, 1e4, 1000, 50,  9, 1, 400, "1e10"),
@@ -182,7 +187,22 @@ _48xx_study = [(p, get_meas(p)) for p in
         # y0=0.62, slope0=1.3 until x=0.20
         SimParams48(5e2, 1e5, 10000, 500, 9, 1, 400, "5e9"),
     ]
-]
+)
+
+_49xx_study = from_params(
+    [
+        SimParams47(5,  1, 400,  "1e10"),
+        SimParams50(5e2, 1e4, 1000, 50,   5, 1, 400, "1e10"),
+        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"),
+    ]
+)
 
 
 
@@ -201,7 +221,9 @@ _48xx_study = [(p, get_meas(p)) for p in
 # of_interest += filter_meas(run="48", coupling_cond=1e4, drive=5e9)
 
 # of_interest += inverters_with_steepest_slope0
-of_interest += _48xx_study
+# of_interest += _47xx_all
+# of_interest += _48xx_study
+of_interest += _49xx_study
 
 # plot cascaded inverter -> buffer
 # for (inv_p, inv_curve) in filter_meas(is_inverter=True):

crates/applications/stacked_cores/scripts/stacked_cores_40xx_db.py

@@ -171,7 +171,10 @@ class SimParams47(SimParams):
     def t_mid(self) -> float:
         return 4e-9
 
-class SimParams48(SimParams):
+class SimParamsV3(SimParams):
+    """
+    SimParams where the conductivities and clock length are configurable
+    """
     def __init__(
         self,
         control_cond: int,
@@ -206,6 +209,7 @@ class SimParams48(SimParams):
         name = super().human_name
         return f"{name} ({self.clock_length_ps}ps {self.control_cond}/{self.coupling_cond} cond)"
 
+class SimParams48(SimParamsV3):
     @property
     def run(self) -> str:
         return "48"
@@ -214,6 +218,23 @@ class SimParams48(SimParams):
     def is_inverter(self) -> bool:
         return False
 
+class SimParams50(SimParamsV3):
+    @property
+    def run(self) -> str:
+        return "50"
+
+    @property
+    def is_inverter(self) -> bool:
+        return False
+
+    @property
+    def t_last(self) -> float:
+        return 5e-12 * self.clock_length_ps
+
+    @property
+    def t_mid(self) -> float:
+        return 4e-12 * self.clock_length_ps
+
 sims = [
     # params, human friendly db name, normalization
     (SimParams40(6,  1,  400, "5e10"), "fwd_40_6_3_1_5e10",        17000),
@@ -404,6 +425,15 @@ sims = [
     (SimParams48(5e2, 1e5,10000, 500, 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),
+
+    (SimParams50(5e2, 1e4, 1000, 50,  5, 1, 400, "1e10"), None, 20000),
+    (SimParams50(5e2, 2e4, 2000, 100, 5, 1, 400, "1e10"), None, 20000),
+    (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, "1e10"), None, 20000),
+    (SimParams50(5e3, 2e4, 2000, 100, 5, 1, 400, "1e10"), None, 20000),
+    (SimParams50(1e4, 2e4, 2000, 100, 5, 1, 400, "1e10"), None, 20000),
 ]
 
 measurements = { real.machine_name: [human, norm, None, None] for (real, human, norm) in sims }
@@ -4819,4 +4849,74 @@ Piecewise(
 )
 """)
 
+set_meas("50-0.0004rad-500ctl_cond-10000coupling_cond-1000ps-50ps-5ctl-5coupling-3_1_winding-1e10-drive", """
+Piecewise(
+    [
+        [ -17425,   3923 ],  # -1.000, M1=-12660
+        [ -17001,   4090 ],  # -0.350, M1=-12447
+        [ -16429,   4261 ],  # -0.300, M1=-12228
+        [ -13120,   5011 ],  # -0.250, M1=-11235
+        [  -9690,   5770 ],  # -0.220, M1=-10232
+        [  -7385,   6255 ],  # -0.200, M1=-9575
+        [  -1872,   7376 ],  # -0.150, M1=-7948
+        [    802,   7871 ],  # -0.100, M1=-7164
+        [   5084,   8592 ],  #  0.000, M1=-5944
+        [   9278,   9209 ],  #  0.100, M1=-4806
+        [  13045,   9699 ],  #  0.200, M1=-3809
+        [  16741,  10142 ],  #  1.000, M1=-2780
+    ]
+)
+""")
+
+set_meas("50-0.0004rad-500ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-1e10-drive", """
+Piecewise(
+    [
+        [ -17349,   4910 ],  # -1.000, M1=-12580
+        [   4848,   8967 ],  #  0.000, M1=-5823
+        [  16730,  10287 ],  #  1.000, M1=-2565
+    ]
+)
+""")
+
+set_meas("50-0.0004rad-1000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-5e9-drive", """
+Piecewise(
+    [
+        [ -17260,  -3784 ],  # -1.000, M1=-13550
+        [   1294,  -1671 ],  #  0.000, M1=-8291
+        [  16709,   -676 ],  #  1.000, M1=-4803
+    ]
+)
+""")
+
+set_meas("50-0.0004rad-1000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-1e10-drive", """
+Piecewise(
+    [
+        [ -17427,   4252 ],  # -1.000, M1=-12305
+        [   6604,  10894 ],  #  0.000, M1=-4913
+        [  16740,  13291 ],  #  1.000, M1=-1934
+    ]
+)
+""")
+
+set_meas("50-0.0004rad-1000ctl_cond-20000coupling_cond-2000ps-100ps-5ctl-5coupling-3_1_winding-2e10-drive", """
+Piecewise(
+    [
+        [ -17587,  10360 ],  # -1.000, M1=-11417
+        [  10229,  16369 ],  #  0.000, M1=-2240
+        [  16771,  16545 ],  #  1.000, M1= -120
+    ]
+)
+""")
+
+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
+        [   7967,  10735 ],  #  0.000, M1=-4212
+        [  16755,  13090 ],  #  1.000, M1=-1516
+    ]
+)
+""")
+
 if __name__ == '__main__': main()

crates/applications/stacked_cores/src/main.rs

@@ -1170,6 +1170,25 @@ fn drive_map_3stack_with_init_buff_cascade(amp0: f32) -> [[ClockState; 3]; 5] {
     ]
 }
 
+#[allow(unused)]
+fn drive_map_3stack_multi_clock_buffer(amp0: f32) -> [[ClockState; 3]; 5] {
+    use ClockState as C;
+    // amplitudes are inverted from what you would expect.
+    // hold(-1) puts the core into a positive M
+    [
+        // init S0; charge S1 negative, S2 negative
+        [C::hold(-amp0),    C::hold_high(),    C::hold_high()],
+        // let the cores settle; open S1 for receive
+        [C::release(-amp0), C::release_high(), C::hold_high()],
+        // write S0 -> S1 (pushes S1 from neg -> pos)
+        [C::hold_high(),    C::float(),        C::hold_high()],
+        // open S2 for receive
+        [C::hold_high(),    C::float(),        C::release_high()],
+        // write S1 -> S2 by (pushes S2 neg -> low)
+        [C::hold_high(),    C::hold_high(),     C::float()],
+    ]
+}
+
 #[allow(unused)]
 fn drive_map_4stack_with_init4(amp0: f32, amp1: f32, amp2: f32, amp3: f32) -> [[ClockState; 4]; 3] {
     use ClockState as C;
@@ -5103,7 +5122,7 @@ fn main() {
         }
     }
 
-    if true {
+    if false {
         for init_set in [
             &[
                 // targeted for (5e2, 5e3, ps(1000), ps(50),  9,  1, um(400),  1e10)
@@ -5407,6 +5426,95 @@ fn main() {
             }
         }
     }
+    if true {
+        for init_set in [
+            &[
+                // targeted
+                // -0.22,
+                // -0.18,
+                // -0.27,
+            ][..],
+            &[
+                // establish the domain/range
+                 1.00,
+                -1.00,
+            ][..],
+            &[
+                 0.00,
+            ][..],
+            &[
+                -0.20,
+                -0.35,
+                -0.25,
+                -0.30,
+                -0.15,
+
+                -0.10,
+                 0.20,
+                 0.10,
+                -0.05,
+                 0.35,
+            ][..],
+        ] {
+            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),
+                (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),
+
+                (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),
+                // (5e2, 5e3, ps(1000), ps(50), 7,  1, um(600),  8e9),
+                // (5e2, 5e3, ps(1000), ps(50), 7,  1, um(600),  1e10),
+                // (5e2, 5e3, ps(1000), ps(50), 7,  1, um(600),  2e10),
+                // (5e2, 5e3, ps(1000), ps(50), 7,  1, um(600),  5e9),
+                // (5e2, 5e3, ps(1000), ps(50), 7,  1, um(600),  3e10),
+                // (5e2, 5e3, ps(1000), ps(50), 5,  1, um(400),  5e9),
+                // (5e2, 5e3, ps(1000), ps(50), 4,  2, um(600),  2e10),
+                // (5e2, 5e3, ps(1000), ps(50), 4,  2, um(600),  5e10),
+                // (5e2, 5e3, ps(1000), ps(50), 5,  1, um(400),  8e9),
+                // (5e2, 5e3, ps(1000), ps(50), 5,  1, um(400),  2e10),
+                // (5e2, 5e3, ps(1000), ps(50), 5,  1, um(400),  4e10),
+            ] {
+                for &init_flt in init_set {
+                    // coupling loops (M0 -> M1) + (M1 -> M2) + control slots
+                    let slots_per_asym = 2*s0_loops;
+                    let net_slots = 2*slots_per_asym + 1;
+                    let mut params = params_v2
+                        .with_clock_phase_duration(clock_duration)
+                        .with_clock_decay(clock_decay)
+                        .with_ctl_conductivity(ctl_cond)
+                        .with_coupling_conductivity(coupling_cond)
+                        .with_s_major(s_major)
+                        .with_coupling_loops(coupling_loops)
+                        .with_input_magnitude(cur_flt)
+                        // control loops
+                        .with_coupling(0, 0, 0, net_slots, CouplingMethod::Control)
+                        .with_coupling(1, 1, 0, net_slots, CouplingMethod::Control)
+                        .with_coupling(2, 2, 0, net_slots, CouplingMethod::Control)
+                    ;
+                    params = couple_asymmetric_buffer(&params, 0 /* sender core */, s0_loops, 1 /* slot offset */, net_slots);
+                    params = couple_asymmetric_buffer(&params, 1 /* sender core */, s0_loops, 1 + slots_per_asym /* slot offset */, net_slots);
+
+                    let name = asymmetric_inverter_name_v3(
+                        &params, "50", coupling_loops /* ctl loops */, coupling_loops, 2*s0_loops + 1, init_flt
+                    );
+                    run_sim(
+                        &name,
+                        drive_map_3stack_multi_clock_buffer(init_flt),
+                        params,
+                    );
+                }
+            }
+        }
+    }
 }