Definitions for Tight bounds for private communication over bosonic Gaussian channels based on teleportation simulation with optimal finite resources
DefinitionQAlgorithms_Papers_Independent_Intakethird100_P17The objects the formalized statements of Tight bounds for private communication over bosonic Gaussian channels based on teleportation simulation with optimal finite resources (arXiv:1808.00608) are stated in. Each was checked against the paper by a blind read-back and a moderator.
GaussianInConvention: Appendix A, PDF p.7, (A1): rescale the pinned [q,p]=2i coordinates. Covariance here is the centered half anticommutator, and the vacuum is I/(2κ). Finite photon energy is explicit; no moment identity is a field.quadratureKernel: Appendix A, PDF p.7, before (A1): actual quadrature coefficients on the finite-particle core. The momentum has +i on creation and -i on annihilation.polynomialKernelProduct: Appendix A, PDF p.7, (A9): composition of finite-band polynomial kernels. The theorem asserts convergence of the internal sums, not a default tsum value.PolynomialProductSummable: Domain guard for the unbounded polynomial product on its common core.quadraticKernel: Appendix A, PDF p.7, (A7): coefficients of the genuine q^T A q.densityCoefficients: Actual Fock matrix coefficients of a bounded oscillator density.PolynomialPairingSummable: Absolute convergence of the density-weighted polynomial trace pairing.polynomialExpectation: Appendix A, PDF p.7, (A8)--(A12): actual polynomial expectation. Only the summable domain is used, and the source theorem concludes that domain.quadraticVariance: Appendix A, PDF p.7, (A12): variance formed from actual first/second moments.complexCovariance: Entrywise real-to-complex conversion of a finite covariance matrix.StrictGibbsCovariance: Appendix A, PDF p.7, finite Gibbs domain of (A2)--(A4). Strict uncertainty is equivalent to all symplectic eigenvalues >1/(2κ). No pure-mode log or inverse-coth endpoint is assigned a finite value.gibbsSeriesTerm: Appendix A, PDF p.7, (A4): analytic arcoth series in inverse powers. The strict Gibbs domain has spectral radius of the inverse <1; convergence is concluded in each theorem using this definition. No normality of sigma*iOmega is assumed, and no covariance commutativity is required.kappaGibbsMatrix: Appendix A, PDF p.7, (A4): G=2iOmega arcoth(2κ sigma iOmega). The result is real on the stated physical domain.gaussianSigma: Appendix A, PDF p.7, (A19): the covariance expression, independently of operator entropy. Strict uncertainty makes the determinant positive real.oscillatorLogCutoff: Appendix A, PDF pp.7--8: bounded support-log regularization in BITS. At faithful but not bounded-below densities the unbounded logarithm is obtained only through the explicitly convergent density-weighted limits below.crossEntropyCutoff: Appendix A, PDF p.7, (A18): independently computed regularized cross entropy.oscillatorCrossEntropy: The actual cross entropy on its convergent finite-log-moment domain. Each identity below also states convergence, excluding off-domain limUnder defaults.oscillatorRelativeEntropy: Independently defined Umegaki relative entropy, in bits, retaining +infinity. The accepted finite-projection definition is used, not the Gaussian formula.relativeVarianceCutoff: Appendix A, PDF p.7, (A14)--(A15): square of the ACTUAL centered operator log difference at a bounded cutoff, without assuming the density operators commute.oscillatorRelativeVariance: Relative entropy variance on the domain where the density-weighted squared support-log difference converges. Convergence is part of the theorem.gaussianVarianceExpression: Appendix A, PDF p.8, final unnumbered formula. The Omega term is positive as printed: replacing the last Omega by its transpose would reverse its sign.
import Mathlib
import Definitions.Def_QAlgorithms_Papers_Independent_Claimed50_P30
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P13
import Definitions.Def_QAlgorithms_Papers_Independent_Next46_P17
/-!
arXiv:1808.00608v4, Appendix A, PDF pp.7--8.
Actual oscillator densities use the pinned complete Fock Weyl system. All
paper-local definitions have bodies; no entropy or moment identity is assumed.
Operational nodes and source ambiguities are recorded individually in HARD.md.
Compilation and independent moderation are pending the separate handoff.
-/
noncomputable section
open Filter
open scoped BigOperators ENNReal ComplexOrder Topology Matrix
attribute [local instance] Classical.propDecidable Classical.decEq
namespace QAlgorithms.Papers.Independent.Intakethird100.P17
open QAlgorithms.Papers.Independent.Claimed50.P30
(ModeOperator Quadrature PhaseVector numberKet trace TraceSummable
FockWeylSystem GaussianMoments canonicalOmega)
open QAlgorithms.Papers.Independent.Dependency50.P13 (IsTraceClass relativeEntropy)
open QAlgorithms.Papers.Independent.Next46.P17 (raisingCoefficient)
/-- Appendix A, PDF p.7, (A1): rescale the pinned [q,p]=2i coordinates.
Covariance here is the centered half anticommutator, and the vacuum is I/(2κ).
Finite photon energy is explicit; no moment identity is a field. -/
def GaussianInConvention {n : ℕ} (W : FockWeylSystem n) (kappa : ℝ)
(rho : ModeOperator n) (u : PhaseVector n)
(sigma : Matrix (Quadrature n) (Quadrature n) ℝ) : Prop :=
GaussianMoments W rho (Real.sqrt (2 * kappa) • u) ((2 * kappa) • sigma) ∧
Summable (fun j : Fin n → ℕ =>
(∑ m, (j m : ℝ)) * (inner ℂ (numberKet j) (rho (numberKet j))).re)
/-- Appendix A, PDF p.7, before (A1): actual quadrature coefficients on the
finite-particle core. The momentum has +i on creation and -i on annihilation. -/
def quadratureKernel {n : ℕ} (kappa : ℝ) (a : Quadrature n)
(i j : Fin n → ℕ) : ℂ :=
if a.2 = 0 then
(raisingCoefficient a.1 i j + raisingCoefficient a.1 j i) /
(Real.sqrt (2 * kappa) : ℂ)
else Complex.I *
(raisingCoefficient a.1 i j - raisingCoefficient a.1 j i) /
(Real.sqrt (2 * kappa) : ℂ)
/-- Appendix A, PDF p.7, (A9): composition of finite-band polynomial kernels.
The theorem asserts convergence of the internal sums, not a default tsum value. -/
def polynomialKernelProduct {n : ℕ}
(A B : Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) :
Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ :=
fun i j => ∑' k, A i k * B k j
/-- Domain guard for the unbounded polynomial product on its common core. -/
def PolynomialProductSummable {n : ℕ}
(A B : Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) : Prop :=
∀ i j, Summable (fun k => A i k * B k j)
/-- Appendix A, PDF p.7, (A7): coefficients of the genuine q^T A q. -/
def quadraticKernel {n : ℕ} (kappa : ℝ)
(A : Matrix (Quadrature n) (Quadrature n) ℝ) :
Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ :=
fun i j => ∑ a, ∑ b, (A a b : ℂ) *
polynomialKernelProduct (quadratureKernel kappa a) (quadratureKernel kappa b) i j
/-- Actual Fock matrix coefficients of a bounded oscillator density. -/
def densityCoefficients {n : ℕ} (rho : ModeOperator n) :
Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ :=
fun i j => rho (numberKet j) i
/-- Absolute convergence of the density-weighted polynomial trace pairing. -/
def PolynomialPairingSummable {n : ℕ} (rho : ModeOperator n)
(A : Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) : Prop :=
Summable (fun z : (Fin n → ℕ) × (Fin n → ℕ) =>
densityCoefficients rho z.1 z.2 * A z.2 z.1)
/-- Appendix A, PDF p.7, (A8)--(A12): actual polynomial expectation.
Only the summable domain is used, and the source theorem concludes that domain. -/
def polynomialExpectation {n : ℕ} (rho : ModeOperator n)
(A : Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) : ℝ :=
(∑' z : (Fin n → ℕ) × (Fin n → ℕ),
densityCoefficients rho z.1 z.2 * A z.2 z.1).re
/-- Appendix A, PDF p.7, (A12): variance formed from actual first/second moments. -/
def quadraticVariance {n : ℕ} (rho : ModeOperator n) (kappa : ℝ)
(A : Matrix (Quadrature n) (Quadrature n) ℝ) : ℝ :=
let Q := quadraticKernel kappa A
polynomialExpectation rho (polynomialKernelProduct Q Q) - polynomialExpectation rho Q ^ 2
/-- Entrywise real-to-complex conversion of a finite covariance matrix. -/
def complexCovariance {n : ℕ} (sigma : Matrix (Quadrature n) (Quadrature n) ℝ) :
Matrix (Quadrature n) (Quadrature n) ℂ := fun i j => (sigma i j : ℂ)
/-- Appendix A, PDF p.7, finite Gibbs domain of (A2)--(A4).
Strict uncertainty is equivalent to all symplectic eigenvalues >1/(2κ).
No pure-mode log or inverse-coth endpoint is assigned a finite value. -/
def StrictGibbsCovariance {n : ℕ} (kappa : ℝ)
(sigma : Matrix (Quadrature n) (Quadrature n) ℝ) : Prop :=
sigma.transpose = sigma ∧ sigma.PosDef ∧
((kappa : ℂ) • complexCovariance sigma +
(Complex.I / 2) • complexCovariance (canonicalOmega n)).PosDef
/-- Appendix A, PDF p.7, (A4): analytic arcoth series in inverse powers.
The strict Gibbs domain has spectral radius of the inverse <1; convergence is
concluded in each theorem using this definition. No normality of sigma*iOmega
is assumed, and no covariance commutativity is required. -/
def gibbsSeriesTerm {n : ℕ} (kappa : ℝ)
(sigma : Matrix (Quadrature n) (Quadrature n) ℝ) (k : ℕ) :
Matrix (Quadrature n) (Quadrature n) ℂ :=
let B := ((2 * kappa : ℝ) : ℂ) •
(complexCovariance sigma * (Complex.I • complexCovariance (canonicalOmega n)))
(((2 * k + 1 : ℕ) : ℂ)⁻¹) • (B⁻¹) ^ (2 * k + 1)
/-- Appendix A, PDF p.7, (A4): G=2iOmega arcoth(2κ sigma iOmega).
The result is real on the stated physical domain. -/
def kappaGibbsMatrix {n : ℕ} (kappa : ℝ)
(sigma : Matrix (Quadrature n) (Quadrature n) ℝ) :
Matrix (Quadrature n) (Quadrature n) ℝ :=
let G := (2 * Complex.I) •
(complexCovariance (canonicalOmega n) * ∑' k, gibbsSeriesTerm kappa sigma k)
fun i j => (G i j).re
/-- Appendix A, PDF p.7, (A19): the covariance expression, independently of
operator entropy. Strict uncertainty makes the determinant positive real. -/
def gaussianSigma {n : ℕ} (kappa : ℝ)
(sigma1 sigma2 : Matrix (Quadrature n) (Quadrature n) ℝ)
(delta : PhaseVector n) : ℝ :=
let G := kappaGibbsMatrix kappa sigma2
let Zsquare := ((kappa : ℂ) • complexCovariance sigma2 +
(Complex.I / 2) • complexCovariance (canonicalOmega n)).det.re
(Real.log Zsquare + kappa * (sigma1 * G).trace +
kappa * ∑ i, ∑ j, delta i * G i j * delta j) / (2 * Real.log 2)
/-- Appendix A, PDF pp.7--8: bounded support-log regularization in BITS.
At faithful but not bounded-below densities the unbounded logarithm is obtained
only through the explicitly convergent density-weighted limits below. -/
def oscillatorLogCutoff {n : ℕ} (rho : ModeOperator n) (k : ℕ) : ModeOperator n := by
letI : ContinuousFunctionalCalculus ℝ (ModeOperator n) IsSelfAdjoint :=
IsSelfAdjoint.instContinuousFunctionalCalculus (A := ModeOperator n)
exact cfc (fun x : ℝ => Real.log (max x (1 / ((k : ℝ) + 1))) / Real.log 2) rho
/-- Appendix A, PDF p.7, (A18): independently computed regularized cross entropy. -/
def crossEntropyCutoff {n : ℕ} (rho sigma : ModeOperator n) (k : ℕ) : ℝ :=
-(trace (rho.comp (oscillatorLogCutoff sigma k))).re
/-- The actual cross entropy on its convergent finite-log-moment domain.
Each identity below also states convergence, excluding off-domain limUnder defaults. -/
def oscillatorCrossEntropy {n : ℕ} (rho sigma : ModeOperator n) : ℝ :=
Filter.limUnder atTop (crossEntropyCutoff rho sigma)
/-- Independently defined Umegaki relative entropy, in bits, retaining +infinity.
The accepted finite-projection definition is used, not the Gaussian formula. -/
def oscillatorRelativeEntropy {n : ℕ} (rho sigma : ModeOperator n)
(hrho : IsTraceClass rho) (hsigma : IsTraceClass sigma) : ℝ≥0∞ :=
relativeEntropy ⟨rho, hrho⟩ ⟨sigma, hsigma⟩ / ENNReal.ofReal (Real.log 2)
/-- Appendix A, PDF p.7, (A14)--(A15): square of the ACTUAL centered operator
log difference at a bounded cutoff, without assuming the density operators commute. -/
def relativeVarianceCutoff {n : ℕ} (rho sigma : ModeOperator n)
(mean : ℝ) (k : ℕ) : ℝ :=
let L := oscillatorLogCutoff rho k - oscillatorLogCutoff sigma k -
(mean : ℂ) • (1 : ModeOperator n)
(trace (rho.comp (L.comp L))).re
/-- Relative entropy variance on the domain where the density-weighted squared
support-log difference converges. Convergence is part of the theorem. -/
def oscillatorRelativeVariance {n : ℕ} (rho sigma : ModeOperator n) (mean : ℝ) : ℝ :=
Filter.limUnder atTop (relativeVarianceCutoff rho sigma mean)
/-- Appendix A, PDF p.8, final unnumbered formula. The Omega term is positive
as printed: replacing the last Omega by its transpose would reverse its sign. -/
def gaussianVarianceExpression {n : ℕ} (kappa : ℝ)
(sigma1 sigma2 : Matrix (Quadrature n) (Quadrature n) ℝ)
(delta : PhaseVector n) : ℝ :=
let G1 := kappaGibbsMatrix kappa sigma1
let G2 := kappaGibbsMatrix kappa sigma2
let D := G1 - G2
let Omega := canonicalOmega n
let B := (8 * kappa ^ 2) • (G2 * sigma1 * G2)
(4 * kappa ^ 2 * (sigma1 * D * sigma1 * D).trace +
(D * Omega * D * Omega).trace + ∑ i, ∑ j, delta i * B i j * delta j) /
(2 * (2 * Real.log 2) ^ 2)
example {n : ℕ} (m : Fin n) (j : Fin n → ℕ) :
raisingCoefficient m (Function.update j m (j m + 1)) j =
(Real.sqrt ((j m : ℝ) + 1) : ℂ) := by
simp [raisingCoefficient]
example {n : ℕ} (kappa : ℝ) (m : Fin n) (i j : Fin n → ℕ) :
quadratureKernel kappa (m, 1) i j = Complex.I *
(raisingCoefficient m i j - raisingCoefficient m j i) /
(Real.sqrt (2 * kappa) : ℂ) := by
simp [quadratureKernel]
example {n : ℕ} (rho : ModeOperator n) (k : ℕ) :
relativeVarianceCutoff rho rho 0 k = 0 := by
simp [relativeVarianceCutoff, trace,
QAlgorithms.Papers.Independent.Dependency50.P13.operatorTrace]
example {n : ℕ} (kappa : ℝ) (sigma : Matrix (Quadrature n) (Quadrature n) ℝ) :
gaussianVarianceExpression kappa sigma sigma 0 = 0 := by
simp [gaussianVarianceExpression]
end QAlgorithms.Papers.Independent.Intakethird100.P17
Read-back
What the Lean code literally says, in plain math · gpt-6.1-sol
QAlgorithms.Papers.Independent.Intakethird100.P17.GaussianInConvention
For any natural number , including , put , , and , and let have coordinate at and elsewhere. The argument is a family of bounded unitary operators , indexed by , satisfying
where, for and ,
For any real , bounded complex-linear operator , vector , and real matrix indexed by , this predicate means all of the following: there exist sequences with such that for every ; and for every ; the family is summable and has sum ; the real matrix is symmetric and positive definite, and is positive semidefinite, where
for every , the family is summable and
and the real family
is summable. There is no explicit sign restriction on : the real square root is when ; when and , the required positive definiteness of cannot hold. When , is a singleton, is empty, and the last summability condition concerns a single zero term.
QAlgorithms.Papers.Independent.Intakethird100.P17.PolynomialPairingSummable
For any natural number , let , , and be the coordinate unit vector at . For any bounded complex-linear operator and arbitrary complex matrix , this predicate means that the single family indexed by pairs
is summable in , equivalently
No positivity or normalization of , or symmetry of , is required. When , the indexing set is a singleton and this predicate holds automatically.
QAlgorithms.Papers.Independent.Intakethird100.P17.PolynomialProductSummable
For any natural number , with , and arbitrary complex matrices , this predicate means that for every , the family
is summable in , equivalently . This is an entrywise requirement and imposes no uniform bound over . When , is a singleton and the predicate holds for every .
QAlgorithms.Papers.Independent.Intakethird100.P17.StrictGibbsCovariance
For any natural number , real , and real matrix indexed by , this predicate means that , that is positive definite over , and that the complex matrix
is positive definite over . Here the subscript means entrywise embedding of a real matrix into , and
Positive definiteness requires the respective Hermitian quadratic forms to be strictly positive on every nonzero vector. There is no separate sign hypothesis on . When , both index sets are empty and these conditions hold vacuously for the unique empty matrix, for every .
QAlgorithms.Papers.Independent.Intakethird100.P17.complexCovariance
For any natural number , including , and arbitrary real matrix indexed by , this definition gives the complex matrix with entries
It imposes no symmetry, positivity, or covariance condition. When , it is the unique matrix with empty row and column index sets.
QAlgorithms.Papers.Independent.Intakethird100.P17.crossEntropyCutoff
For every natural number , bounded complex linear operators on , and natural number , this defines
where is the standard unit vector at , and is the real continuous functional calculus for self-adjoint operators, with value zero if is not self-adjoint. No density, trace-class, or summability hypotheses are required. The displayed infinite sum is total and has value zero if its complex-valued summands are not summable. The cases and are included.
QAlgorithms.Papers.Independent.Intakethird100.P17.densityCoefficients
For any natural number , let , let , and let denote the vector with coordinate at and elsewhere. For any bounded complex-linear operator , this definition gives the complex matrix
No positivity, trace-class, or normalization assumption on is imposed. When , is a singleton, so this matrix has one entry.
QAlgorithms.Papers.Independent.Intakethird100.P17.gaussianSigma
For every natural number , real number , real matrices indexed by , and real vector indexed by , this defines
where , , all its other entries are zero, and real matrices inside complex expressions are embedded entrywise. Here
No covariance, positivity, invertibility, or summability conditions are imposed. Singular matrix inverses and nonsummable matrix series have value zero. The real logarithm is total, with and for negative , so need not be positive. The case is included and gives zero.
QAlgorithms.Papers.Independent.Intakethird100.P17.gaussianVarianceExpression
For every natural number , real number , real matrices indexed by , and real vector indexed by , let , , with all other entries zero. Regarding real matrices as complex matrices entrywise, set
The definition returns the real number
No positivity, symmetry, invertibility, convergence, or sign conditions on the result are required. Singular matrix inverses and nonsummable matrix series have value zero. Both and are allowed; for , the result is zero.
QAlgorithms.Papers.Independent.Intakethird100.P17.gibbsSeriesTerm
For any natural number , real , arbitrary real matrix indexed by , and natural number , including , embed entrywise into the complex matrix , and define by
With , this definition gives the complex matrix
The exponent is a positive natural number, and its reciprocal is taken in . No symmetry, positivity, strict Gibbs condition, or invertibility hypothesis is imposed. The matrix inverse is total: for a singular , , hence . When , the result is the unique empty matrix.
QAlgorithms.Papers.Independent.Intakethird100.P17.kappaGibbsMatrix
For every natural number , every real number , and every real matrix indexed by , this defines a real matrix . Let , , and let all other entries of be zero; regard real matrices as complex matrices entrywise. Then
There are no positivity, symmetry, invertibility, or convergence assumptions. Matrix inversion is total and gives the zero matrix for a singular matrix; the infinite matrix sum is the total mathematical sum, which is zero when the series is not summable. In particular, is allowed. When , the result is the unique matrix on the empty index set.
QAlgorithms.Papers.Independent.Intakethird100.P17.oscillatorCrossEntropy
For every natural number and bounded complex linear operators on , this defines the real value selected by the total filter-limit operation, as through natural numbers, for the sequence
Here is the standard number basis, indexed by all functions , and uses the real continuous functional calculus for self-adjoint operators, giving zero outside that condition. Nonsummable complex sums have value zero. No convergence, self-adjointness, density, or trace-class hypothesis is imposed: if the real sequence converges, the selected value is its limit, but the definition itself does not assert convergence. The case is included.
QAlgorithms.Papers.Independent.Intakethird100.P17.oscillatorLogCutoff
For every natural number , every bounded complex linear operator on , and every natural number , this defines the operator obtained by applying the real continuous functional calculus for self-adjoint operators to with the scalar function
No self-adjointness, positivity, normalization, or trace-class assumption on is supplied; the total functional-calculus operation gives zero when its self-adjointness condition fails. The scalar logarithm always has a positive argument here. Both , where the lower cutoff is , and , where is one-dimensional, are included.
QAlgorithms.Papers.Independent.Intakethird100.P17.oscillatorRelativeEntropy
For every natural number and bounded complex linear operators on , supplied with proofs that each admits a representation
for some sequences , this defines an element of . For every natural number and every orthonormal family in , form
Assign this pair the value
where and use the real continuous functional calculus for self-adjoint matrices, including its total default outside the required conditions. The defined value is
with the denominator interpreted in . No positivity or trace-one assumption is imposed on either operator. Empty orthonormal families , and the case , are included.
QAlgorithms.Papers.Independent.Intakethird100.P17.oscillatorRelativeVariance
For every natural number , bounded complex linear operators on , and real number , this defines the real value selected by the total filter-limit operation, as , for
Here is the standard number basis indexed by ; the real functional calculus gives zero outside self-adjoint operators, and a nonsummable complex trace series has value zero. No convergence, density, trace-class, or positivity hypotheses are imposed, and is arbitrary. If the displayed real sequence converges, the definition gives its limit; no such convergence is asserted by the definition. The case is included.
QAlgorithms.Papers.Independent.Intakethird100.P17.polynomialExpectation
For any natural number , let , , and be the coordinate unit vector at . For any bounded complex-linear operator and arbitrary complex matrix , this definition assigns the real number
The sum is over pairs, rather than an iterated sum. There is no summability, positivity, normalization, or symmetry hypothesis; a nonsummable complex family has total sum , so the defined value is then . When , the expression is the real part of the product of the sole operator coefficient and the sole entry of .
QAlgorithms.Papers.Independent.Intakethird100.P17.polynomialKernelProduct
For any natural number , let . For arbitrary complex matrices , their defined kernel product is the matrix
No summability hypothesis is required: each sum is the total infinite-sum operation, which returns if its complex-valued family is not summable. When , is a singleton and this is ordinary multiplication of the two sole entries.
QAlgorithms.Papers.Independent.Intakethird100.P17.quadraticKernel
For any natural number , real , and arbitrary real matrix indexed by , this definition gives a complex matrix on by
Here, writing , let be if is obtained from by increasing coordinate by , and otherwise; then
All real factors are embedded in . No symmetry or positivity of , or positivity of , is assumed. Infinite sums use the total sum operation, returning for a nonsummable family. For , division by the zero real square root makes every , and hence , zero. For , the outer sums are empty and the resulting matrix is zero on the singleton .
QAlgorithms.Papers.Independent.Intakethird100.P17.quadraticVariance
For any natural number , real , bounded complex-linear operator on , where , and arbitrary real matrix indexed by , let be the coordinate unit vector at . Define when is obtained by increasing coordinate of by , and otherwise, and set
The defined real number is
No summability, normalization, positivity, or symmetry assumptions are imposed, and no nonnegativity of this value is asserted. Each infinite sum returns if its family is not summable. Real square roots of nonpositive numbers and division by zero have value ; consequently the expression is when . It is also when , because then is empty and the matrix is zero.
QAlgorithms.Papers.Independent.Intakethird100.P17.quadratureKernel
For any natural number , real , quadrature index , and occupation indices , define
The complex-valued kernel is
with real quantities embedded in . No positivity assumption on is imposed. If , its real square-root denominator is , and complex division by makes the kernel . For , there is no quadrature index , so there are no such evaluations.
QAlgorithms.Papers.Independent.Intakethird100.P17.relativeVarianceCutoff
For every natural number , bounded complex linear operators on , real number , and natural number , this defines
The vectors are the standard unit vectors indexed by all . Each functional calculus is the real continuous functional calculus for self-adjoint operators and gives zero when its self-adjointness condition fails; a nonsummable complex trace series has value zero. There are no density, trace-class, positivity, summability, or restrictions on , and is not required to be an expectation. The product is . The cases and are included.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.