Documentation

MIPRE.Background.LIDT.MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerIdentities.ProjectorApprox

Section 5 — P-Q approximation for QXP layers #

The final comparison estimates between the projective family P produced from XHat and the original projective Q layer.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.qxpMixedCrossSum_eq_sqrt {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (data : QXPLayerData Outcome ι) :
a : Outcome, (Xa data a).conjTranspose * data.xHat * Pa data a = CFC.sqrt (QTotal data.qLayer)

The mixed X--Xhat summation appearing in the proof of lem:P-Q-approx.

After rewriting Q_a and P_a through the matrices X, Xhat, and the projective measurement T, the sum of the mixed terms collapses to Xᴴ * Xhat, hence to sqrt Q by lem:X-times-X-hat.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.pQApprox {Outcome : Type u_1} {ι : Type u_2} [Fintype ι] [DecidableEq ι] [Fintype Outcome] (ψ : QuantumState ι) (A : Measurement Outcome ι) (ζ : Error) (data : QXPLayerData Outcome ι) ( : ψ.IsNormalized) ( : 0 ζ) (hζ_small : ζ 1 / 4) :

P is close to Q (lem:P-Q-approx).

The final internal comparison in the paper's repair step is derived from the primitive X/XHat/P identities in QXPLayerData, the rank-reduction witness, and the standard small-error hypotheses. No closeness bound is stored inside QXPLayerData; the proof below follows the paper's expansion through squaredDifference, qAlmostProjective, and sqrtQCompleteness.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.pQApprox_ofRankReductionSigmaRangeAndSvdIdentities {Outcome : Type u_1} [Fintype Outcome] [DecidableEq Outcome] {ι : Type u_2} [Fintype ι] [DecidableEq ι] (ψ : QuantumState ι) (A : Measurement Outcome ι) (ζ : Error) {qLayer : QLayerData Outcome ι} (hRank : RankReductionWitness ψ A ζ qLayer) [Nonempty (FiniteHilbertSpace.sigmaFinCarrier fun (a : Outcome) => Matrix.rank (qLayer.q.outcome a))] (xHat : Matrix (sigmaRangeCarrier qLayer.q) ι ) (xHat_coisometry : xHat * xHat.conjTranspose = 1) (xHat_mixed : (sigmaFinRangeEmbedding qLayer.q.outcome ).conjTranspose * xHat = CFC.sqrt (QTotal qLayer)) ( : ψ.IsNormalized) ( : 0 ζ) (hζ_small : ζ 1 / 4) :
∃ (data : QXPLayerData Outcome ι) (hq : data.qLayer = sigmaRangeQLayer qLayer.q), (hq data.x = have this := sigmaFinRangeEmbedding qLayer.q.outcome ; this) (hq data.xHat = have this := xHat; this) SDDOpRel ψ (uniformDistribution Unit) (constOpFamily data.qLayer.q) (constOpFamily (PFamily data)) (30 * zetaQuarterRoot ζ)

Apply lem:P-Q-approx to the canonical sigma-space QXP layer obtained from a rank-reduction witness.

The theorem keeps the SVD/polar data for Xhat explicit, but removes the remaining bookkeeping needed to use pQApprox: the rank-reduction witness is transported to the sigma-space layer, and the resulting QXPLayerData is the canonical one built from sigmaFinRangeEmbedding.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.pQApprox_ofRankReductionSigmaRangePositiveGram {Outcome : Type uOutcome} [Fintype Outcome] [DecidableEq Outcome] {ι : Type uι} [Fintype ι] [DecidableEq ι] (ψ : QuantumState ι) (A : Measurement Outcome ι) (ζ : Error) {qLayer : QLayerData Outcome ι} (hRank : RankReductionWitness ψ A ζ qLayer) [Nonempty (FiniteHilbertSpace.sigmaFinCarrier fun (a : Outcome) => Matrix.rank (qLayer.q.outcome a))] ( : ψ.IsNormalized) ( : 0 ζ) (hζ_small : ζ 1 / 4) :
∃ (xHat : Matrix (sigmaRangeCarrier qLayer.q) ι ), xHat * xHat.conjTranspose = 1 (sigmaFinRangeEmbedding qLayer.q.outcome ).conjTranspose * xHat = CFC.sqrt (QTotal qLayer) ∃ (data : QXPLayerData Outcome ι) (hq : data.qLayer = sigmaRangeQLayer qLayer.q), (hq data.x = have this := sigmaFinRangeEmbedding qLayer.q.outcome ; this) (hq data.xHat = have this := xHat; this) SDDOpRel ψ (uniformDistribution Unit) (constOpFamily data.qLayer.q) (constOpFamily (PFamily data)) (30 * zetaQuarterRoot ζ)

Apply lem:P-Q-approx to the positive-Gram sigma-space QXP layer.

This is the constructor-facing form of the local Q -> X -> Xhat -> P stage. The rank-reduction witness supplies the sigma-space X; the positive-Gram polar construction supplies Xhat; and the theorem concludes the paper's P-versus-Q approximation without any external SVD data.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.pQApprox_ofRankReductionSigmaRangePositiveGram_with_x_coisometry {Outcome : Type uOutcome} [Fintype Outcome] [DecidableEq Outcome] {ι : Type uι} [Fintype ι] [DecidableEq ι] (ψ : QuantumState ι) (A : Measurement Outcome ι) (ζ : Error) {qLayer : QLayerData Outcome ι} (hRank : RankReductionWitness ψ A ζ qLayer) (hsum_le_one : a : Outcome, qLayer.q.outcome a 1) [Nonempty (FiniteHilbertSpace.sigmaFinCarrier fun (a : Outcome) => Matrix.rank (qLayer.q.outcome a))] ( : ψ.IsNormalized) ( : 0 ζ) (hζ_small : ζ 1 / 4) :
∃ (xHat : Matrix (sigmaRangeCarrier qLayer.q) ι ), xHat * xHat.conjTranspose = 1 (sigmaFinRangeEmbedding qLayer.q.outcome ).conjTranspose * xHat = CFC.sqrt (QTotal qLayer) ∃ (data : QXPLayerData Outcome ι) (hq : data.qLayer = sigmaRangeQLayer qLayer.q), (hq data.x = have this := sigmaFinRangeEmbedding qLayer.q.outcome ; this) (hq data.xHat = have this := xHat; this) data.x * data.x.conjTranspose = 1 SDDOpRel ψ (uniformDistribution Unit) (constOpFamily data.qLayer.q) (constOpFamily (PFamily data)) (30 * zetaQuarterRoot ζ)

Apply lem:P-Q-approx to the positive-Gram sigma-space QXP layer, and also record coisometry of the original sigma embedding X.

The additional hypothesis is the subnormalization of the projective Q family. Under this hypothesis the range basis vectors chosen for distinct outcomes are orthogonal, so the finite sigma-range embedding has orthonormal rows. This is the construction-level coisometry condition later used to preserve the fresh option-completion row and hence to obtain the QXP-internal comparison Q_none ≤ P_none. An additional source-to-Q comparison is still required to recover (optionCompletion A).outcome none ≤ P.outcome none; see docs/reports/issue-1642-restrictsome-residual-domination-obstruction.md.

theorem MIPStarRE.LDT.MakingMeasurementsProjective.pQApprox_ofRankReductionSigmaRangeAndRectangularSvdSquareRootUnitaryGroup {Outcome : Type uOutcome} [Fintype Outcome] [DecidableEq Outcome] {ι : Type uι} [Fintype ι] [DecidableEq ι] (ψ : QuantumState ι) (A : Measurement Outcome ι) (ζ : Error) {qLayer : QLayerData Outcome ι} (hRank : RankReductionWitness ψ A ζ qLayer) [Nonempty (FiniteHilbertSpace.sigmaFinCarrier fun (a : Outcome) => Matrix.rank (qLayer.q.outcome a))] (U : (Matrix.unitaryGroup (sigmaRangeCarrier qLayer.q) )) (V : (Matrix.unitaryGroup ι )) (S Iro : Matrix (sigmaRangeCarrier qLayer.q) ι ) (hIro : Iro * Iro.conjTranspose = 1) (hx : sigmaFinRangeEmbedding qLayer.q.outcome = U * S * (↑V).conjTranspose) (hMiddle_nonneg : 0 V * (S.conjTranspose * Iro) * (↑V).conjTranspose) (hMiddle_sq : V * (S.conjTranspose * Iro) * (↑V).conjTranspose * (V * (S.conjTranspose * Iro) * (↑V).conjTranspose) = QTotal qLayer) ( : ψ.IsNormalized) ( : 0 ζ) (hζ_small : ζ 1 / 4) :
∃ (data : QXPLayerData Outcome ι) (hq : data.qLayer = sigmaRangeQLayer qLayer.q), (hq data.x = have this := sigmaFinRangeEmbedding qLayer.q.outcome ; this) (hq data.xHat = have this := U * Iro * (↑V).conjTranspose; this) SDDOpRel ψ (uniformDistribution Unit) (constOpFamily data.qLayer.q) (constOpFamily (PFamily data)) (30 * zetaQuarterRoot ζ)

Apply lem:P-Q-approx to unitary-group rectangular SVD data and the positive-square characterization of the middle factor.

The square SVD factors are represented as elements of Matrix.unitaryGroup; hence the unitarity laws are not separate hypotheses.