Julia LinearAlgebra.qr!用法及代码示例

qr!(A, pivot = NoPivot(); blocksize)

当 A 是 StridedMatrix 的子类型时，qr! 与 qr 相同，但通过覆盖输入 A 而不是创建副本来节省空间。如果因式分解产生的数字不能由 A 的元素类型表示，则会引发 InexactError 异常，例如对于整数类型。

例子

julia> a = [1. 2.; 3. 4.]
2×2 Matrix{Float64}:
 1.0  2.0
 3.0  4.0

julia> qr!(a)
LinearAlgebra.QRCompactWY{Float64, Matrix{Float64}}
Q factor:
2×2 LinearAlgebra.QRCompactWYQ{Float64, Matrix{Float64}}:
 -0.316228  -0.948683
 -0.948683   0.316228
R factor:
2×2 Matrix{Float64}:
 -3.16228  -4.42719
  0.0      -0.632456

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

julia> qr!(a)
ERROR: InexactError: Int64(3.1622776601683795)
Stacktrace:
[...]