aware of p_ext

src/controller.rs

@@ -32,8 +32,12 @@ fn get_weight(pd: f64, pm: f64) -> (f64, f64) {
         w2 = w2 / wsum;
     }
 
-    // w1 = 0.5;
+    // w1 = 0.1;
     // w2 = 1.0 - w1;
+    // p is w1 / w2
+    // let p = 10.;
+    // let w1 = p / (1. + p);
+    // let w2 = 1. / (1. + p);
     log::info!("w1: {}, w2: {}", w1, w2);
 
     (w1, w2)
@@ -77,8 +81,9 @@ pub fn runner2(xn: usize, zm: usize, measure_rounds: usize, measure_model: Measu
     // let p1 = 0.01;
     let probs = &measure_model.success.probs;
     let p_ext = probs[1] * 4.;
-    log::debug!("p_ext: {}", p_ext);
-    let (w1, w2) = get_weight(p1 + p_ext, measure_model.failure.prob_sum);
+    log::info!("p_ext: {}", p_ext);
+    let (w1, w2) = get_weight(p1 + p_ext, measure_model.failure.prob_sum + p_ext);
+    // let (w1, w2) = get_weight(p1 + p_ext, measure_model.failure.prob_sum);
 
     let channel = NoiseChannel::new_depolarizing_channel(p1);
     let error_model = ErrorModel::new(channel, Some(measure_model));