Documentation

MIPRE.Background.QLD.Complete

Completing the pair measurement and its evaluated marginals (lem:qld-global-complete, #

lem:qld-global-sandwich)

Stage 4 ends with the interface SimulPair (lem:qld-simultaneous): a projective measurement with outcomes pairs (g_X, g_Z) of polynomials in m variables, whose evaluated marginals are consistent with the opposite party's expanded point measurements. This file builds it from a global measurement G with outcomes in LowIndDegPoly F (4m) d (the polynomials on the padded point) and the estimates of stage 4b, abstractly: for a projective G of the first player and two projective families X x, Z z of the second.

The completion #

A good outcome g = α g_X(x) + β g_Z(z) (IsGood) is relabelled by the pair (g_X, g_Z), read off its two coefficient polynomials through LowIndDegPoly.blockPoly (a polynomial reading only the x block, as a polynomial in m variables); every other outcome is sent to the fixed pair (0, 0). This is a coarse-graining of G, so it is a complete projective measurement (pairMeas, isPVM_pairMeas): the paper's completion by the orthogonal complement R, whose weight is the weight of the non-good outcomes, bounded by lem:qld-global-separate.

The marginals #

From the products estimate ∑_g E_u ‖(G_g ⊗ (1 - B_u(g(u)))) Φ‖² ≤ Δ, the triangle inequality and Cauchy--Schwarz give ∑_g E_u ‖(G_g ⊗ B_u(g(u))) Φ‖² ≥ 1 - 2√Δ (one_sub_two_sqrt_le_sum_snorm_sq). For a good g and the order X_a Z_b, the fibre analysis of lem:qld-global-separate bounds E_u ‖(G_g ⊗ B_u(g(u))) Φ‖² by 2/q · ⟨G_g⟩ plus the diagonal weight E_z ⟨G_g ⊗ Z_{g_Z(z)}(z)⟩, and the non-good outcomes carry weight at most δ_G; hence the Z marginal E_z ∑_g ⟨G_g ⊗ Z_{g_Z(z)}(z)⟩ ≥ 1 - 2√Δ - 2/q - δ_G (marg_Z_ge), and the order Z_b X_a gives the X marginal. The consistency form (inconsistency_evalMarg_Z_le, _X_le) is the same statement read through inconsistency. The instance of SimulPair is in PaddedLIDT.lean.

Stated in a bipartite model #

Stated in a bipartite model (Phase 5 of planning/mipco-track.md). The state is that of a bipartite model K (the model in which the global measurements live), the global measurement is a POVM G in the first player's algebra with IsPVMIn G.op, the families X x, Z z are projective families of the second player's algebra, a Born probability is K.bornProb, and a squared norm K.snorm (K.πA (G.op g) * K.πB B) ^ 2. The completed pair measurement is the coarse-graining G.map (pairOf hd) in the same algebra, so it needs no transport to the shape of SimulPair: the matrix associativity reindexing of the doubly extended registers (isPVM_reindex, isPVM_reindex_povm, reindex_aOp_aOp, POVM.aOp_aOp_reindex) and the matrix projective-measurement packaging (ProjectiveMeasurement.isPVM_M) are gone.

The pair of a good outcome #

theorem MIPRE.QLD.isLinAB_of_isGood {F : Type u_1} [Field F] {m d : ℕ} [NeZero m] {g : LIDT.LowIndDegPoly} (hg : IsGood g) :
theorem MIPRE.QLD.notMem_xSet_iff {m : ℕ} [NeZero m] {k : Fin (4 * m)} :
k ∉ xSet ↔ ∀ (j : Fin m), xIdx m j ≠ k
theorem MIPRE.QLD.notMem_zSet_iff {m : ℕ} [NeZero m] {k : Fin (4 * m)} :
k ∉ zSet ↔ ∀ (j : Fin m), zIdx m j ≠ k
theorem MIPRE.QLD.not_depOutside_gA_of_isGood {F : Type u_1} [Field F] {m d : ℕ} [NeZero m] (hd : 1 ≤ d) {g : LIDT.LowIndDegPoly} (hg : IsGood g) :

The α coefficient of a good outcome reads only the x block.

theorem MIPRE.QLD.not_depOutside_gB_of_isGood {F : Type u_1} [Field F] {m d : ℕ} [NeZero m] (hd : 1 ≤ d) {g : LIDT.LowIndDegPoly} (hg : IsGood g) :

The β coefficient of a good outcome reads only the z block.

noncomputable def MIPRE.QLD.gX {F : Type u_1} [Field F] {m d : ℕ} [NeZero m] (hd : 1 ≤ d) (g : LIDT.LowIndDegPoly) :

The X polynomial of an outcome: its α coefficient read on the x block.

Equations
Instances For
    noncomputable def MIPRE.QLD.gZ {F : Type u_1} [Field F] {m d : ℕ} [NeZero m] (hd : 1 ≤ d) (g : LIDT.LowIndDegPoly) :

    The Z polynomial of an outcome: its β coefficient read on the z block.

    Equations
    Instances For
      theorem MIPRE.QLD.eval_gX {F : Type u_1} [Field F] {m d : ℕ} [NeZero m] (hd : 1 ≤ d) {g : LIDT.LowIndDegPoly} (hg : IsGood g) (u : LIDT.Point F (4 * m)) :
      (gX hd g).eval (xBlk u) = (gA hd g).eval u
      theorem MIPRE.QLD.eval_gZ {F : Type u_1} [Field F] {m d : ℕ} [NeZero m] (hd : 1 ≤ d) {g : LIDT.LowIndDegPoly} (hg : IsGood g) (u : LIDT.Point F (4 * m)) :
      (gZ hd g).eval (zBlk u) = (gB hd g).eval u
      noncomputable def MIPRE.QLD.pairOf {F : Type u_1} [Field F] [DecidableEq F] {m d : ℕ} [NeZero m] (hd : 1 ≤ d) (g : LIDT.LowIndDegPoly) :
      PolyPair F m d

      The relabelling: a good outcome becomes its pair (g_X, g_Z), every other outcome the fixed pair (0, 0).

      Equations
      Instances For
        theorem MIPRE.QLD.pairOf_of_isGood {F : Type u_1} [Field F] [DecidableEq F] {m d : ℕ} [NeZero m] (hd : 1 ≤ d) {g : LIDT.LowIndDegPoly} (hg : IsGood g) :
        pairOf hd g = (gX hd g, gZ hd g)
        noncomputable def MIPRE.QLD.pairMeas {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m d : ℕ} [NeZero m] {R : Type u_2} [Ring R] [StarRing R] [PartialOrder R] [StarOrderedRing R] (hd : 1 ≤ d) (G : POVMIn LIDT.LowIndDegPoly R) :
        POVMIn (PolyPair F m d) R

        The completed pair measurement: the global measurement coarse-grained by the relabelling. The non-good outcomes are absorbed into the pair (0, 0): the paper's complement R, added to a fixed legal outcome.

        Equations
        Instances For
          theorem MIPRE.QLD.pairMeas_mats {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m d : ℕ} [NeZero m] {R : Type u_2} [Ring R] [StarRing R] [PartialOrder R] [StarOrderedRing R] (hd : 1 ≤ d) (G : POVMIn LIDT.LowIndDegPoly R) (p : PolyPair F m d) :
          (pairMeas hd G).op p = ∑ g : LIDT.LowIndDegPoly with pairOf hd g = p, G.op g
          theorem MIPRE.QLD.isPVM_pairMeas {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m d : ℕ} [NeZero m] {R : Type u_2} [Ring R] [StarRing R] [PartialOrder R] [StarOrderedRing R] (hd : 1 ≤ d) {G : POVMIn LIDT.LowIndDegPoly R} (hG : IsPVMIn G.op) :

          The completed pair measurement is projective.

          theorem MIPRE.QLD.evalMarg_pairMeas {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m d : ℕ} [NeZero m] {R : Type u_2} [Ring R] [StarRing R] [PartialOrder R] [StarOrderedRing R] (hd : 1 ≤ d) (G : POVMIn LIDT.LowIndDegPoly R) (W : Bas) (u : LIDT.Point F m) :
          evalMarg (pairMeas hd G) W u = POVMIn.map (fun (g : LIDT.LowIndDegPoly) => (PolyPair.proj W (pairOf hd g)).eval u) G

          The evaluated marginal of the completed measurement is the global measurement coarse-grained by the evaluated component of the pair.

          The inconsistency of a coarse-grained global measurement #

          theorem MIPRE.QLD.inconsistency_map_eq {F : Type u_1} [Fintype F] [DecidableEq F] {n d : ℕ} {𝒞 : Type u_2} {𝒜 : Type u_3} {ℬ : Type u_4} [Ring 𝒞] [StarRing 𝒞] [Algebra ℂ 𝒞] [Ring 𝒜] [StarRing 𝒜] [Algebra ℂ 𝒜] [Ring ℬ] [StarRing ℬ] [Algebra ℂ ℬ] [PartialOrder 𝒜] [StarOrderedRing 𝒜] [PartialOrder ℬ] {Q : Type u_5} [Fintype Q] (μ : Q → ℝ) (hμ : ∑ q : Q, μ q = 1) {K : BipartiteModel 𝒞 𝒜 ℬ} (hK : ‖K.ψ‖ = 1) (G : POVMIn LIDT.LowIndDegPoly 𝒜) (f : Q → LIDT.LowIndDegPoly → F) (P : Q → POVMIn F ℬ) :
          K.inconsistency μ (fun (q : Q) => POVMIn.map (f q) G) P = 1 - ∑ q : Q, μ q * ∑ g : LIDT.LowIndDegPoly, K.bornProb (G.op g) ((P q).op (f q g))

          The inconsistency of a coarse-graining of a global measurement, against a family of POVMs indexed by the same questions: one minus the average agreement, outcome by outcome of the original measurement. inconsistency_evalPOVM_eq is the case of evaluation at the point.

          The block marginals of the uniform padded point #

          theorem MIPRE.QLD.sum_uniform_zBlk {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m : ℕ} [NeZero m] (f : LIDT.Point F m → ℝ) :
          ∑ u : LIDT.Point F (4 * m), uniform (LIDT.Point F (4 * m)) u * f (zBlk u) = ∑ z : LIDT.Point F m, uniform (LIDT.Point F m) z * f z

          The z block of a uniform padded point is uniform.

          theorem MIPRE.QLD.sum_uniform_xBlk {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m : ℕ} [NeZero m] (f : LIDT.Point F m → ℝ) :
          ∑ u : LIDT.Point F (4 * m), uniform (LIDT.Point F (4 * m)) u * f (xBlk u) = ∑ x : LIDT.Point F m, uniform (LIDT.Point F m) x * f x

          The x block of a uniform padded point is uniform.

          The evaluated marginals #

          theorem MIPRE.QLD.sq_sub_two_mul_le_sq {a b c : ℝ} (ha : 0 ≤ a) (hb : 0 ≤ b) (hc : 0 ≤ c) (h : a ≤ c + b) :
          a ^ 2 - 2 * (a * b) ≤ c ^ 2

          a ≤ c + b for nonnegative reals gives a² - 2ab ≤ c².

          theorem MIPRE.QLD.one_sub_two_sqrt_le_sum_snorm_sq {F : Type u_1} [Field F] [Fintype F] {m d : ℕ} {𝒞 : Type u_2} {𝒜 : Type u_3} {ℬ : Type u_4} [Ring 𝒞] [StarRing 𝒞] [Algebra ℂ 𝒞] [Ring 𝒜] [StarRing 𝒜] [Algebra ℂ 𝒜] [Ring ℬ] [StarRing ℬ] [Algebra ℂ ℬ] [PartialOrder 𝒜] {K : BipartiteModel 𝒞 𝒜 ℬ} (hK : ‖K.ψ‖ = 1) (G : POVMIn LIDT.LowIndDegPoly 𝒜) (hG : IsPVMIn G.op) (B : LIDT.Point F (4 * m) → F → ℬ) {Δ : ℝ} (hprod : ∑ u : LIDT.Point F (4 * m), uniform (LIDT.Point F (4 * m)) u * ∑ g : LIDT.LowIndDegPoly, K.snorm (K.πA (G.op g) * K.πB (1 - B u (g.eval u))) ^ 2 ≤ Δ) :
          1 - 2 * √Δ ≤ ∑ u : LIDT.Point F (4 * m), uniform (LIDT.Point F (4 * m)) u * ∑ g : LIDT.LowIndDegPoly, K.snorm (K.πA (G.op g) * K.πB (B u (g.eval u))) ^ 2

          The products estimate bounds the weight through B from below: from ∑_g E_u ‖(G_g ⊗ (1 - B_u(g(u)))) Φ‖² ≤ Δ, the triangle inequality and Cauchy--Schwarz give ∑_g E_u ‖(G_g ⊗ B_u(g(u))) Φ‖² ≥ 1 - 2√Δ.

          theorem MIPRE.QLD.marg_Z_ge {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m d : ℕ} [NeZero m] {𝒞 : Type u_2} {𝒜 : Type u_3} {ℬ : Type u_4} [Ring 𝒞] [StarRing 𝒞] [Algebra ℂ 𝒞] [Ring 𝒜] [StarRing 𝒜] [Algebra ℂ 𝒜] [Ring ℬ] [StarRing ℬ] [Algebra ℂ ℬ] [PartialOrder 𝒜] [StarOrderedRing 𝒜] [PartialOrder ℬ] [StarOrderedRing ℬ] {X Z : LIDT.Point F m → F → ℬ} (hd : 1 ≤ d) {K : BipartiteModel 𝒞 𝒜 ℬ} (hK : ‖K.ψ‖ = 1) (G : POVMIn LIDT.LowIndDegPoly 𝒜) (hG : IsPVMIn G.op) (hX : ∀ (x : LIDT.Point F m), IsPVMIn (X x)) (hZ : ∀ (z : LIDT.Point F m), IsPVMIn (Z z)) {Δ δG : ℝ} (hprod : ∑ u : LIDT.Point F (4 * m), uniform (LIDT.Point F (4 * m)) u * ∑ g : LIDT.LowIndDegPoly, K.snorm (K.πA (G.op g) * K.πB (1 - ordComb ordXZ X Z u (g.eval u))) ^ 2 ≤ Δ) (hbad : ∑ g : LIDT.LowIndDegPoly with ¬IsGood g, K.bornProb (G.op g) 1 ≤ δG) :
          1 - 2 * √Δ - 2 / ↑(Fintype.card F) - δG ≤ ∑ z : LIDT.Point F m, uniform (LIDT.Point F m) z * ∑ g : LIDT.LowIndDegPoly, K.bornProb (G.op g) (Z z ((PolyPair.proj Bas.Z (pairOf hd g)).eval z))

          The Z marginal of the completed measurement tracks the second player's Z point measurement: from the products estimate for the order X_a Z_b and the weight δ_G of the non-good outcomes, E_z ∑_g ⟨G_g ⊗ Z_{g_Z(z)}(z)⟩ ≥ 1 - 2√Δ - 2/q - δ_G, where g_Z is the second component of the relabelled outcome.

          theorem MIPRE.QLD.marg_X_ge {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m d : ℕ} [NeZero m] {𝒞 : Type u_2} {𝒜 : Type u_3} {ℬ : Type u_4} [Ring 𝒞] [StarRing 𝒞] [Algebra ℂ 𝒞] [Ring 𝒜] [StarRing 𝒜] [Algebra ℂ 𝒜] [Ring ℬ] [StarRing ℬ] [Algebra ℂ ℬ] [PartialOrder 𝒜] [StarOrderedRing 𝒜] [PartialOrder ℬ] [StarOrderedRing ℬ] {X Z : LIDT.Point F m → F → ℬ} (hd : 1 ≤ d) {K : BipartiteModel 𝒞 𝒜 ℬ} (hK : ‖K.ψ‖ = 1) (G : POVMIn LIDT.LowIndDegPoly 𝒜) (hG : IsPVMIn G.op) (hX : ∀ (x : LIDT.Point F m), IsPVMIn (X x)) (hZ : ∀ (z : LIDT.Point F m), IsPVMIn (Z z)) {Δ δG : ℝ} (hprod : ∑ u : LIDT.Point F (4 * m), uniform (LIDT.Point F (4 * m)) u * ∑ g : LIDT.LowIndDegPoly, K.snorm (K.πA (G.op g) * K.πB (1 - ordComb ordZX X Z u (g.eval u))) ^ 2 ≤ Δ) (hbad : ∑ g : LIDT.LowIndDegPoly with ¬IsGood g, K.bornProb (G.op g) 1 ≤ δG) :
          1 - 2 * √Δ - 2 / ↑(Fintype.card F) - δG ≤ ∑ x : LIDT.Point F m, uniform (LIDT.Point F m) x * ∑ g : LIDT.LowIndDegPoly, K.bornProb (G.op g) (X x ((PolyPair.proj Bas.X (pairOf hd g)).eval x))

          The X marginal, from the order Z_b X_a: the mirror image.

          theorem MIPRE.QLD.inconsistency_evalMarg_Z_le {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m d : ℕ} [NeZero m] {𝒞 : Type u_2} {𝒜 : Type u_3} {ℬ : Type u_4} [Ring 𝒞] [StarRing 𝒞] [Algebra ℂ 𝒞] [Ring 𝒜] [StarRing 𝒜] [Algebra ℂ 𝒜] [Ring ℬ] [StarRing ℬ] [Algebra ℂ ℬ] [PartialOrder 𝒜] [StarOrderedRing 𝒜] [PartialOrder ℬ] [StarOrderedRing ℬ] {X Z : LIDT.Point F m → F → ℬ} (hd : 1 ≤ d) {K : BipartiteModel 𝒞 𝒜 ℬ} (hK : ‖K.ψ‖ = 1) (G : POVMIn LIDT.LowIndDegPoly 𝒜) (hG : IsPVMIn G.op) (hX : ∀ (x : LIDT.Point F m), IsPVMIn (X x)) (hZ : ∀ (z : LIDT.Point F m), IsPVMIn (Z z)) (N : LIDT.Point F m → POVMIn F ℬ) (hN : ∀ (z : LIDT.Point F m) (a : F), (N z).op a = Z z a) {Δ δG : ℝ} (hprod : ∑ u : LIDT.Point F (4 * m), uniform (LIDT.Point F (4 * m)) u * ∑ g : LIDT.LowIndDegPoly, K.snorm (K.πA (G.op g) * K.πB (1 - ordComb ordXZ X Z u (g.eval u))) ^ 2 ≤ Δ) (hbad : ∑ g : LIDT.LowIndDegPoly with ¬IsGood g, K.bornProb (G.op g) 1 ≤ δG) :
          K.inconsistency (uniform (LIDT.Point F m)) (fun (z : LIDT.Point F m) => evalMarg (pairMeas hd G) Bas.Z z) N ≤ 2 * √Δ + 2 / ↑(Fintype.card F) + δG

          lem:qld-global-sandwich, the Z marginal in consistency form: the evaluated Z marginal of the completed measurement against a family N whose elements are the Z z a.

          theorem MIPRE.QLD.inconsistency_evalMarg_X_le {F : Type u_1} [Field F] [Fintype F] [DecidableEq F] {m d : ℕ} [NeZero m] {𝒞 : Type u_2} {𝒜 : Type u_3} {ℬ : Type u_4} [Ring 𝒞] [StarRing 𝒞] [Algebra ℂ 𝒞] [Ring 𝒜] [StarRing 𝒜] [Algebra ℂ 𝒜] [Ring ℬ] [StarRing ℬ] [Algebra ℂ ℬ] [PartialOrder 𝒜] [StarOrderedRing 𝒜] [PartialOrder ℬ] [StarOrderedRing ℬ] {X Z : LIDT.Point F m → F → ℬ} (hd : 1 ≤ d) {K : BipartiteModel 𝒞 𝒜 ℬ} (hK : ‖K.ψ‖ = 1) (G : POVMIn LIDT.LowIndDegPoly 𝒜) (hG : IsPVMIn G.op) (hX : ∀ (x : LIDT.Point F m), IsPVMIn (X x)) (hZ : ∀ (z : LIDT.Point F m), IsPVMIn (Z z)) (N : LIDT.Point F m → POVMIn F ℬ) (hN : ∀ (x : LIDT.Point F m) (a : F), (N x).op a = X x a) {Δ δG : ℝ} (hprod : ∑ u : LIDT.Point F (4 * m), uniform (LIDT.Point F (4 * m)) u * ∑ g : LIDT.LowIndDegPoly, K.snorm (K.πA (G.op g) * K.πB (1 - ordComb ordZX X Z u (g.eval u))) ^ 2 ≤ Δ) (hbad : ∑ g : LIDT.LowIndDegPoly with ¬IsGood g, K.bornProb (G.op g) 1 ≤ δG) :
          K.inconsistency (uniform (LIDT.Point F m)) (fun (x : LIDT.Point F m) => evalMarg (pairMeas hd G) Bas.X x) N ≤ 2 * √Δ + 2 / ↑(Fintype.card F) + δG

          lem:qld-global-sandwich, the X marginal in consistency form.