Documentation

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

Section 5 — Rectangular SVD constructors for Q/X/XHat/P data #

Rectangular-SVD matrix identities and constructors for the paper's Q/X/XHat/P intermediate layer.

def MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvdX {μ : Type u_1} {ι : Type u_2} [Fintype μ] [Fintype ι] (U : Matrix μ μ ) (V : Matrix ι ι ) (S : Matrix μ ι ) :
Matrix μ ι

The rectangular singular-value reconstruction of X.

If U and V are the square unitary factors and S is the rectangular diagonal singular-value matrix, this is the matrix U * S * V† appearing in the paper.

Equations
Instances For
    def MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvdXHat {μ : Type u_1} {ι : Type u_2} [Fintype μ] [Fintype ι] (U : Matrix μ μ ) (V : Matrix ι ι ) (Iro : Matrix μ ι ) :
    Matrix μ ι

    The rectangular-SVD candidate for Xhat.

    This is the paper's matrix U * I_{m × d} * V†, where Iro is the rectangular row-coisometry replacing the displayed rectangular identity matrix.

    Equations
    Instances For
      def MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvdDomainMiddle {μ : Type u_1} {ι : Type u_2} [Fintype μ] [Fintype ι] (V : Matrix ι ι ) (S Iro : Matrix μ ι ) :
      Matrix ι ι

      The domain-space middle operator in the product X† * Xhat.

      For rectangular-SVD data this is V * (S† * Iro) * V†, the operator on the original Hilbert space which is later identified with sqrt Q.

      Equations
      Instances For
        def MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvdRangeMiddle {μ : Type u_1} {ι : Type u_2} [Fintype μ] [Fintype ι] (U : Matrix μ μ ) (S Iro : Matrix μ ι ) :
        Matrix μ μ

        The row-space middle operator in the product X * Xhat†.

        For rectangular-SVD data this is U * (S * Iro†) * U†, the formal counterpart of the square matrix Σ_{m × m} in the paper.

        Equations
        Instances For
          theorem MIPStarRE.LDT.MakingMeasurementsProjective.unitaryGroup_mul_rectangular_coisometry {μ : Type u_1} {ι : Type u_2} [Fintype μ] [DecidableEq μ] [Fintype ι] (U : (Matrix.unitaryGroup μ )) (W : Matrix μ ι ) (hW : W * W.conjTranspose = 1) :
          U * W * (U * W).conjTranspose = 1

          Left multiplication by a unitary group element preserves rectangular row coisometries.

          If W W† = I, then (U W)(U W)† = I. This is the square-unitary part of the paper's calculation of lem:X-hat-squared.

          Right multiplication by the adjoint of a unitary group element preserves rectangular row coisometries.

          If W W† = I, then (W V†)(W V†)† = I. This is the right-square-unitary part of the paper's calculation of lem:X-hat-squared.

          theorem MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvd_xHat_coisometry_unitaryGroup {μ : Type u_1} {ι : Type u_2} [Fintype μ] [DecidableEq μ] [Fintype ι] [DecidableEq ι] (U : (Matrix.unitaryGroup μ )) (V : (Matrix.unitaryGroup ι )) (Iro : Matrix μ ι ) (hIro : Iro * Iro.conjTranspose = 1) :
          U * Iro * (↑V).conjTranspose * (U * Iro * (↑V).conjTranspose).conjTranspose = 1

          The row-coisometry identity for the rectangular SVD choice of Xhat, with the square factors represented as unitary group elements.

          theorem MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvd_xHat_mixed_raw {μ : Type u_1} {ι : Type u_2} [Fintype μ] [DecidableEq μ] [Fintype ι] (x : Matrix μ ι ) (U : Matrix μ μ ) (V : Matrix ι ι ) (S Iro : Matrix μ ι ) (hU_right : U.conjTranspose * U = 1) (hx : x = U * S * V.conjTranspose) :

          The mixed product obtained by multiplying the rectangular SVD formulae.

          This lemma contains only the matrix algebra. The spectral identification of the right hand side with a square root is supplied separately, since downstream constructors usually know the square root in the form CFC.sqrt Q.

          theorem MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvd_xHat_mixed_raw_unitaryGroup {μ : Type u_1} {ι : Type u_2} [Fintype μ] [DecidableEq μ] [Fintype ι] (x : Matrix μ ι ) (U : (Matrix.unitaryGroup μ )) (V : Matrix ι ι ) (S Iro : Matrix μ ι ) (hx : x = U * S * V.conjTranspose) :

          The mixed product obtained by multiplying the rectangular SVD formulae, with the left square factor represented as a unitary group element.

          theorem MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvd_x_mul_xHat_conjTranspose_raw {μ : Type u_1} {ι : Type u_2} [Fintype μ] [Fintype ι] [DecidableEq ι] (x : Matrix μ ι ) (U : Matrix μ μ ) (V : Matrix ι ι ) (S Iro : Matrix μ ι ) (hV_right : V.conjTranspose * V = 1) (hx : x = U * S * V.conjTranspose) :

          The first mixed product obtained from the rectangular SVD formulae.

          This is the first identity in the paper's lem:X-times-X-hat, written in the rectangular notation used by the formalization. If X = U * S * Vᴴ and Xhat = U * Iro * Vᴴ, then X * Xhatᴴ = U * (S * Iroᴴ) * Uᴴ. The middle factor S * Iroᴴ is the formal counterpart of the square matrix Σ_{m × m} appearing in the paper.

          theorem MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvd_x_mul_xHat_conjTranspose_raw_unitaryGroup {μ : Type u_1} {ι : Type u_2} [Fintype μ] [Fintype ι] [DecidableEq ι] (x : Matrix μ ι ) (U : Matrix μ μ ) (V : (Matrix.unitaryGroup ι )) (S Iro : Matrix μ ι ) (hx : x = U * S * (↑V).conjTranspose) :

          The first mixed product obtained from the rectangular SVD formulae, with the right square factor represented as a unitary group element.

          A positive operator whose square is Q is the CFC square root of Q.

          This is the uniqueness of the positive square root, stated in the matrix language used in the projectivization layer.

          theorem MIPStarRE.LDT.MakingMeasurementsProjective.rectangularSvd_middle_eq_sqrt_of_square {μ : Type u_1} {ι : Type u_2} [Fintype μ] [Fintype ι] [NonUnitalContinuousFunctionalCalculus (Matrix ι ι ) IsSelfAdjoint] (V : Matrix ι ι ) (S Iro : Matrix μ ι ) (Q : Matrix ι ι ) (hMiddle_nonneg : 0 V * (S.conjTranspose * Iro) * V.conjTranspose) (hMiddle_sq : V * (S.conjTranspose * Iro) * V.conjTranspose * (V * (S.conjTranspose * Iro) * V.conjTranspose) = Q) :

          The square-root identification for the middle factor in the rectangular SVD calculation.

          If the middle operator V * (Sᴴ * Iro) * Vᴴ is positive and its square is the target operator Q, then it is the positive square root of Q. This is the spectral input which turns the raw SVD calculation into the paper's identity X† Xhat = sqrt Q.

          The mixed rectangular SVD identity with the left square factor represented as a unitary group element and the target square root supplied as an external operator Q.

          theorem MIPStarRE.LDT.MakingMeasurementsProjective.exists_xHat_of_rectangularSvd_unitaryGroup {μ : Type u_1} {ι : Type u_2} [Fintype μ] [DecidableEq μ] [Fintype ι] [DecidableEq ι] (x : Matrix μ ι ) (U : (Matrix.unitaryGroup μ )) (V : (Matrix.unitaryGroup ι )) (S Iro : Matrix μ ι ) (Q : Matrix ι ι ) (hIro : Iro * Iro.conjTranspose = 1) (hx : x = U * S * (↑V).conjTranspose) (hSqrt : V * (S.conjTranspose * Iro) * (↑V).conjTranspose = CFC.sqrt Q) :
          ∃ (xHat : Matrix μ ι ), xHat = U * Iro * (↑V).conjTranspose xHat * xHat.conjTranspose = 1 x.conjTranspose * xHat = CFC.sqrt Q

          The rectangular SVD data determine a candidate Xhat and its two primitive identities, with the square factors represented as unitary group elements.

          noncomputable def MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerData.ofRankReductionAndSvdIdentities {Outcome : Type u_1} [Fintype Outcome] {ι : Type u_2} [Fintype ι] [DecidableEq ι] {ψ : QuantumState ι} {A : Measurement Outcome ι} {ζ : Error} {qLayer : QLayerData Outcome ι} (hRank : RankReductionWitness ψ A ζ qLayer) (x xHat : Matrix qLayer.auxSpace.carrier ι ) (qa_eq : ∀ (a : Outcome), qLayer.q.outcome a = x.conjTranspose * Ta qLayer a * x) (xHat_coisometry : xHat * xHat.conjTranspose = 1) (xHat_mixed : x.conjTranspose * xHat = CFC.sqrt (QTotal qLayer)) :
          QXPLayerData Outcome ι

          Assemble QXPLayerData from a rank-reduction witness and the SVD identities for Xhat.

          The rank-reduction witness supplies the projectivity of each Q_a and the identity ∑_a Q_a = Q. The remaining hypotheses are exactly the local matrix-decomposition data for Q_a = X† T_a X and the two SVD-derived identities for the chosen Xhat. Thus this constructor removes the rank-reduction fields from the caller's obligations. The rectangular SVD/polar decomposition that provides xHat is supplied by the sigma-range / rectangular polar-decomposition route, which provides the unitary and coisometry factors from the positive spectral subspace of Q. The end-to-end chain through the rounding-to-projectors, rank-reduction, and orthogonalization lemmas remains to be closed upstream.

          Equations
          • One or more equations did not get rendered due to their size.
          Instances For
            noncomputable def MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerData.ofRankReductionAndRectangularSvdUnitaryGroup {Outcome : Type u_1} [Fintype Outcome] {ι : Type u_2} [Fintype ι] [DecidableEq ι] {ψ : QuantumState ι} {A : Measurement Outcome ι} {ζ : Error} {qLayer : QLayerData Outcome ι} (hRank : RankReductionWitness ψ A ζ qLayer) (x : Matrix qLayer.auxSpace.carrier ι ) (U : (Matrix.unitaryGroup qLayer.auxSpace.carrier )) (V : (Matrix.unitaryGroup ι )) (S Iro : Matrix qLayer.auxSpace.carrier ι ) (qa_eq : ∀ (a : Outcome), qLayer.q.outcome a = x.conjTranspose * Ta qLayer a * x) (hIro : Iro * Iro.conjTranspose = 1) (hx : x = U * S * (↑V).conjTranspose) (hSqrt : V * (S.conjTranspose * Iro) * (↑V).conjTranspose = CFC.sqrt (QTotal qLayer)) :
            QXPLayerData Outcome ι

            Assemble QXPLayerData from rank-reduction data and rectangular SVD data whose square factors are represented as Mathlib unitary-group elements.

            The left and right unitarity laws are carried by the type of U and V; the only remaining rectangular law is the coisometry of Iro.

            Equations
            • One or more equations did not get rendered due to their size.
            Instances For
              noncomputable def MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerData.ofRankReductionAndRectangularSvdSquareRootUnitaryGroup {Outcome : Type u_1} [Fintype Outcome] {ι : Type u_2} [Fintype ι] [DecidableEq ι] {ψ : QuantumState ι} {A : Measurement Outcome ι} {ζ : Error} {qLayer : QLayerData Outcome ι} (hRank : RankReductionWitness ψ A ζ qLayer) (x : Matrix qLayer.auxSpace.carrier ι ) (U : (Matrix.unitaryGroup qLayer.auxSpace.carrier )) (V : (Matrix.unitaryGroup ι )) (S Iro : Matrix qLayer.auxSpace.carrier ι ) (qa_eq : ∀ (a : Outcome), qLayer.q.outcome a = x.conjTranspose * Ta qLayer a * x) (hIro : Iro * Iro.conjTranspose = 1) (hx : x = 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) :
              QXPLayerData Outcome ι

              Assemble QXPLayerData from unitary-group rectangular SVD data whose middle factor is characterized as the positive square root of Q.

              The unitarity of the square factors is represented by Matrix.unitaryGroup, while the square-root identification is supplied by positivity and the square equation for the middle factor.

              Equations
              • One or more equations did not get rendered due to their size.
              Instances For
                noncomputable def MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerData.ofSigmaRangeAndRectangularSvdUnitaryGroup {Outcome : Type uOutcome} [Fintype Outcome] [DecidableEq Outcome] {ι : Type uι} [Fintype ι] [DecidableEq ι] (q : OpFamily Outcome ι) (qa_projective : ∀ (a : Outcome), Quantum.IsProj (q.outcome a)) (q_sum_eq_total : a : Outcome, q.outcome a = q.total) [Nonempty (FiniteHilbertSpace.sigmaFinCarrier fun (a : Outcome) => Matrix.rank (q.outcome a))] (U : (Matrix.unitaryGroup (sigmaRangeCarrier q) )) (V : (Matrix.unitaryGroup ι )) (S Iro : Matrix (sigmaRangeCarrier q) ι ) (hIro : Iro * Iro.conjTranspose = 1) (hx : sigmaFinRangeEmbedding q.outcome qa_projective = U * S * (↑V).conjTranspose) (hSqrt : V * (S.conjTranspose * Iro) * (↑V).conjTranspose = CFC.sqrt q.total) :
                QXPLayerData Outcome ι

                Assemble the canonical sigma-space Q/X/Xhat/P layer from rectangular SVD data whose square factors are Matrix.unitaryGroup elements.

                The unitarity hypotheses for U and V are carried by their types.

                Equations
                • One or more equations did not get rendered due to their size.
                Instances For
                  noncomputable def MIPStarRE.LDT.MakingMeasurementsProjective.QXPLayerData.ofSigmaRangeAndRectangularSvdSquareRootUnitaryGroup {Outcome : Type uOutcome} [Fintype Outcome] [DecidableEq Outcome] {ι : Type uι} [Fintype ι] [DecidableEq ι] (q : OpFamily Outcome ι) (qa_projective : ∀ (a : Outcome), Quantum.IsProj (q.outcome a)) (q_sum_eq_total : a : Outcome, q.outcome a = q.total) [Nonempty (FiniteHilbertSpace.sigmaFinCarrier fun (a : Outcome) => Matrix.rank (q.outcome a))] (U : (Matrix.unitaryGroup (sigmaRangeCarrier q) )) (V : (Matrix.unitaryGroup ι )) (S Iro : Matrix (sigmaRangeCarrier q) ι ) (hIro : Iro * Iro.conjTranspose = 1) (hx : sigmaFinRangeEmbedding q.outcome qa_projective = 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) = q.total) :
                  QXPLayerData Outcome ι

                  Assemble the canonical sigma-space Q/X/Xhat/P layer from unitary-group rectangular SVD data whose middle factor is characterized as a positive square root.

                  Equations
                  • One or more equations did not get rendered due to their size.
                  Instances For
                    theorem MIPStarRE.LDT.MakingMeasurementsProjective.exists_qxpLayerData_ofRankReductionSigmaRangeAndRectangularSvdUnitaryGroup {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) (hSqrt : V * (S.conjTranspose * Iro) * (↑V).conjTranspose = 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 := U * Iro * (↑V).conjTranspose; this

                    Rank-reduction existence form for the canonical sigma-space QXP layer from unitary-group rectangular SVD data.

                    theorem MIPStarRE.LDT.MakingMeasurementsProjective.exists_qxpLayerData_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) :
                    ∃ (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

                    Rank-reduction existence form for the canonical sigma-space QXP layer from unitary-group rectangular SVD data whose middle factor is characterized as a positive square root.