Trace of a matrix

This file defines the trace of a matrix, the map sending a matrix to the sum of its diagonal entries.

See also LinearAlgebra.Trace for the trace of an endomorphism.

Tags

matrix, trace, diagonal

def Matrix.trace {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] (A : Matrix n n R) : R

The trace of a square matrix. For more bundled versions, see:

• Matrix.traceAddMonoidHom
• Matrix.traceLinearMap

@[simp]

theorem Matrix.trace_diagonal {R : Type u_6} [AddCommMonoid R] {o : Type u_8} [Fintype o] [DecidableEq o] (d : o → R) : (diagonal d).trace = ∑ i : o, d i

@[simp]

theorem Matrix.trace_zero (n : Type u_3) (R : Type u_6) [Fintype n] [AddCommMonoid R] : trace 0 = 0

@[simp]

theorem Matrix.trace_eq_zero_of_isEmpty {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] [IsEmpty n] (A : Matrix n n R) : A.trace = 0

@[simp]

theorem Matrix.trace_add {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] (A B : Matrix n n R) : (A + B).trace = A.trace + B.trace

@[simp]

theorem Matrix.trace_smul {n : Type u_3} {α : Type u_5} {R : Type u_6} [Fintype n] [AddCommMonoid R] [DistribSMul α R] (r : α) (A : Matrix n n R) : (r • A).trace = r • A.trace

@[simp]

theorem Matrix.trace_transpose {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] (A : Matrix n n R) : A.transpose.trace = A.trace

@[simp]

theorem Matrix.trace_conjTranspose {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] [StarAddMonoid R] (A : Matrix n n R) : A.conjTranspose.trace = star A.trace

def Matrix.traceAddMonoidHom (n : Type u_3) (R : Type u_6) [Fintype n] [AddCommMonoid R] : Matrix n n R →+ R

Matrix.trace as an AddMonoidHom

@[simp]

theorem Matrix.traceAddMonoidHom_apply (n : Type u_3) (R : Type u_6) [Fintype n] [AddCommMonoid R] (A : Matrix n n R) : (traceAddMonoidHom n R) A = A.trace

def Matrix.traceLinearMap (n : Type u_3) (α : Type u_5) (R : Type u_6) [Fintype n] [AddCommMonoid R] [Semiring α] [Module α R] : Matrix n n R →ₗ[α] R

Matrix.trace as a LinearMap

@[simp]

theorem Matrix.traceLinearMap_apply (n : Type u_3) (α : Type u_5) (R : Type u_6) [Fintype n] [AddCommMonoid R] [Semiring α] [Module α R] (A : Matrix n n R) : (traceLinearMap n α R) A = A.trace

@[simp]

theorem Matrix.trace_list_sum {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] (l : List (Matrix n n R)) : l.sum.trace = (List.map trace l).sum

@[simp]

theorem Matrix.trace_multiset_sum {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] (s : Multiset (Matrix n n R)) : s.sum.trace = (Multiset.map trace s).sum

@[simp]

theorem Matrix.trace_sum {ι : Type u_1} {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] (s : Finset ι) (f : ι → Matrix n n R) : (∑ i ∈ s, f i).trace = ∑ i ∈ s, (f i).trace

theorem AddMonoidHom.map_trace {n : Type u_3} {R : Type u_6} {S : Type u_7} [Fintype n] [AddCommMonoid R] [AddCommMonoid S] {F : Type u_8} [FunLike F R S] [AddMonoidHomClass F R S] (f : F) (A : Matrix n n R) : f A.trace = (A.map ⇑f).trace

theorem Matrix.trace_blockDiagonal {n : Type u_3} {p : Type u_4} {R : Type u_6} [Fintype n] [Fintype p] [AddCommMonoid R] [DecidableEq p] (M : p → Matrix n n R) : (blockDiagonal M).trace = ∑ i : p, (M i).trace

theorem Matrix.trace_blockDiagonal' {p : Type u_4} {R : Type u_6} [Fintype p] [AddCommMonoid R] [DecidableEq p] {m : p → Type u_8} [(i : p) → Fintype (m i)] (M : (i : p) → Matrix (m i) (m i) R) : (blockDiagonal' M).trace = ∑ i : p, (M i).trace

@[simp]

theorem Matrix.trace_sub {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommGroup R] (A B : Matrix n n R) : (A - B).trace = A.trace - B.trace

@[simp]

theorem Matrix.trace_neg {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommGroup R] (A : Matrix n n R) : (-A).trace = -A.trace

@[simp]

theorem Matrix.trace_one {n : Type u_3} {R : Type u_6} [Fintype n] [DecidableEq n] [AddCommMonoidWithOne R] : trace 1 = ↑(Fintype.card n)

@[simp]

theorem Matrix.trace_transpose_mul {m : Type u_2} {n : Type u_3} {R : Type u_6} [Fintype m] [Fintype n] [AddCommMonoid R] [Mul R] (A : Matrix m n R) (B : Matrix n m R) : (A.transpose * B.transpose).trace = (A * B).trace

theorem Matrix.trace_mul_comm {m : Type u_2} {n : Type u_3} {R : Type u_6} [Fintype m] [Fintype n] [AddCommMonoid R] [CommMagma R] (A : Matrix m n R) (B : Matrix n m R) : (A * B).trace = (B * A).trace

theorem Matrix.trace_mul_cycle {m : Type u_2} {n : Type u_3} {p : Type u_4} {R : Type u_6} [Fintype m] [Fintype n] [Fintype p] [NonUnitalCommSemiring R] (A : Matrix m n R) (B : Matrix n p R) (C : Matrix p m R) : (A * B * C).trace = (C * A * B).trace

theorem Matrix.trace_mul_cycle' {m : Type u_2} {n : Type u_3} {p : Type u_4} {R : Type u_6} [Fintype m] [Fintype n] [Fintype p] [NonUnitalCommSemiring R] (A : Matrix m n R) (B : Matrix n p R) (C : Matrix p m R) : (A * (B * C)).trace = (C * (A * B)).trace

@[simp]

theorem Matrix.trace_replicateCol_mul_replicateRow {n : Type u_3} {R : Type u_6} [Fintype n] {ι : Type u_8} [Unique ι] [NonUnitalNonAssocSemiring R] (a b : n → R) : (replicateCol ι a * replicateRow ι b).trace = a ⬝ᵥ b

@[simp]

theorem Matrix.trace_vecMulVec {n : Type u_3} {R : Type u_6} [Fintype n] [NonUnitalNonAssocSemiring R] (a b : n → R) : (vecMulVec a b).trace = a ⬝ᵥ b

theorem Matrix.trace_submatrix_succ {R : Type u_6} {n : ℕ} [AddCommMonoid R] (M : Matrix (Fin n.succ) (Fin n.succ) R) : M 0 0 + (M.submatrix Fin.succ Fin.succ).trace = M.trace

theorem Matrix.trace_units_conj {m : Type u_2} {R : Type u_6} [Fintype m] [DecidableEq m] [CommSemiring R] (M : (Matrix m m R)ˣ) (N : Matrix m m R) : (↑M * N * ↑M⁻¹).trace = N.trace

theorem Matrix.trace_units_conj' {m : Type u_2} {R : Type u_6} [Fintype m] [DecidableEq m] [CommSemiring R] (M : (Matrix m m R)ˣ) (N : Matrix m m R) : (↑M⁻¹ * N * ↑M).trace = N.trace

Special cases for Fin n for low values of n

@[simp]

theorem Matrix.trace_fin_zero {R : Type u_6} [AddCommMonoid R] (A : Matrix (Fin 0) (Fin 0) R) : A.trace = 0

theorem Matrix.trace_fin_one {R : Type u_6} [AddCommMonoid R] (A : Matrix (Fin 1) (Fin 1) R) : A.trace = A 0 0

theorem Matrix.trace_fin_two {R : Type u_6} [AddCommMonoid R] (A : Matrix (Fin 2) (Fin 2) R) : A.trace = A 0 0 + A 1 1

theorem Matrix.trace_fin_three {R : Type u_6} [AddCommMonoid R] (A : Matrix (Fin 3) (Fin 3) R) : A.trace = A 0 0 + A 1 1 + A 2 2

@[simp]

theorem Matrix.trace_fin_one_of {R : Type u_6} [AddCommMonoid R] (a : R) : !![a].trace = a

@[simp]

theorem Matrix.trace_fin_two_of {R : Type u_6} [AddCommMonoid R] (a b c d : R) : !![a, b; c, d].trace = a + d

@[simp]

theorem Matrix.trace_fin_three_of {R : Type u_6} [AddCommMonoid R] (a b c d e f g h i : R) : !![a, b, c; d, e, f; g, h, i].trace = a + e + i

@[simp]

theorem Matrix.trace_single_eq_of_ne {n : Type u_10} {α : Type u_12} [DecidableEq n] [Fintype n] [AddCommMonoid α] (i j : n) (c : α) (h : i ≠ j) : (single i j c).trace = 0

@[simp]

theorem Matrix.trace_single_eq_same {n : Type u_10} {α : Type u_12} [DecidableEq n] [Fintype n] [AddCommMonoid α] (i : n) (c : α) : (single i i c).trace = c

theorem Matrix.trace_single_mul {m : Type u_9} {n : Type u_10} {R : Type u_11} [DecidableEq m] [DecidableEq n] [Fintype n] [NonUnitalNonAssocSemiring R] [Fintype m] (i : n) (j : m) (a : R) (x : Matrix m n R) : (single i j a * x).trace = a • x j i

theorem Matrix.trace_mul_single {m : Type u_9} {n : Type u_10} {R : Type u_11} [DecidableEq m] [DecidableEq n] [Fintype n] [NonUnitalNonAssocSemiring R] [Fintype m] (x : Matrix m n R) (i : n) (j : m) (a : R) : (x * single i j a).trace = MulOpposite.op a • x j i

theorem Matrix.trace_surjective {n : Type u_3} {R : Type u_6} [Fintype n] [AddCommMonoid R] [Nonempty n] : Function.Surjective trace

theorem Matrix.ext_iff_trace_mul_left {m : Type u_2} {n : Type u_3} {R : Type u_6} [Fintype m] [Fintype n] [NonAssocSemiring R] {A B : Matrix m n R} : A = B ↔ ∀ (x : Matrix n m R), (x * A).trace = (x * B).trace

Matrices A and B are equal iff (x * A).trace = (x * B).trace for all x.

theorem Matrix.ext_iff_trace_mul_right {m : Type u_2} {n : Type u_3} {R : Type u_6} [Fintype m] [Fintype n] [NonAssocSemiring R] {A B : Matrix m n R} : A = B ↔ ∀ (x : Matrix n m R), (A * x).trace = (B * x).trace

Matrices A and B are equal iff (A * x).trace = (B * x).trace for all x.