Julia LinearAlgebra.BunchKaufman用法及代码示例

BunchKaufman <: Factorization

对称或 Hermitian 矩阵 A 的 Bunch-Kaufman 分解的矩阵分解类型为 P'UDU'P 或 P'LDL'P ，具体取决于上(默认)或下三角形是否存储在 A 中。如果 A 是复对称的，则 U' 和 L' 分别表示非共轭转置，即 transpose(U) 和 transpose(L) 。这是 bunchkaufman 的返回类型，对应的矩阵分解函数。

如果 S::BunchKaufman 是分解对象，则可以根据给定的 S.uplo 和 S.p 通过 S.D 、 S.U 或 S.L 获得分量。

在给定 S.uplo 和 S.p 的情况下，迭代分解产生组件 S.D 、 S.U 或 S.L 。

例子

julia> A = [1 2; 2 3]
2×2 Matrix{Int64}:
 1  2
 2  3

julia> S = bunchkaufman(A) # A gets wrapped internally by Symmetric(A)
BunchKaufman{Float64, Matrix{Float64}}
D factor:
2×2 Tridiagonal{Float64, Vector{Float64}}:
 -0.333333  0.0
  0.0       3.0
U factor:
2×2 UnitUpperTriangular{Float64, Matrix{Float64}}:
 1.0  0.666667
  ⋅   1.0
permutation:
2-element Vector{Int64}:
 1
 2

julia> d, u, p = S; # destructuring via iteration

julia> d == S.D && u == S.U && p == S.p
true

julia> S = bunchkaufman(Symmetric(A, :L))
BunchKaufman{Float64, Matrix{Float64}}
D factor:
2×2 Tridiagonal{Float64, Vector{Float64}}:
 3.0   0.0
 0.0  -0.333333
L factor:
2×2 UnitLowerTriangular{Float64, Matrix{Float64}}:
 1.0        ⋅
 0.666667  1.0
permutation:
2-element Vector{Int64}:
 2
 1