for
ProvedGaussianMatrix.specNorm_pinvR_sqc-identitydeterministiclinear-algebrapseudoinversespectral-norm
Let be a real matrix and let be its right pseudoinverse, where the inverse of a singular matrix is taken to be . Then, in the spectral norm,
Write , which is symmetric. By the -identity we get . If is singular, both sides are . Together with this gives , which turns tail bounds for into tail bounds for .
Formalization Note. pinvR A = Aᵀ * (A * Aᵀ)⁻¹ uses Mathlib's total inverse, and specNorm is the operator norm. The identity holds unconditionally.
Preamble
import Definitions.Def_GaussianMatrix_basic open MeasureTheory ProbabilityTheory open scoped Matrix
Formal statement
namespace GaussianMatrix
theorem specNorm_pinvR_sq {r k : ℕ} (A : Matrix (Fin r) (Fin k) ℝ) :
specNorm (pinvR A) ^ 2 = specNorm (A * Aᵀ)⁻¹ := by sorry
end GaussianMatrixSource
standard fact: C*-identity ‖B‖² = ‖BᵀB‖ for the spectral norm, applied to B = Aᵀ(AAᵀ)⁻¹, for which BᵀB = (AAᵀ)⁻¹ (e.g. N. Halko, P.-G. Martinsson, J. A. Tropp, SIAM Review 53 (2011), §A.2, ‖G†‖ = σ_min(G)⁻¹).