diff --git a/stdlib/LinearAlgebra/test/hessenberg.jl b/stdlib/LinearAlgebra/test/hessenberg.jl index 84fb08237362f..94ed21c091fcd 100644 --- a/stdlib/LinearAlgebra/test/hessenberg.jl +++ b/stdlib/LinearAlgebra/test/hessenberg.jl @@ -46,7 +46,6 @@ let n = 10 complex.(Areal, Aimg) : Areal) A = herm ? Hermitian(A_ + A_') : A_ - Am = Matrix(A) # workaround for #32001 H = hessenberg(A) eltyh = eltype(H) @@ -67,7 +66,7 @@ let n = 10 @test convert(Array, 2 * H) ≈ 2 * A ≈ convert(Array, H * 2) @test convert(Array, H + 2I) ≈ A + 2I ≈ convert(Array, 2I + H) - @test convert(Array, H + (2+4im)I) ≈ Am + (2+4im)I ≈ convert(Array, (2+4im)I + H) + @test convert(Array, H + (2+4im)I) ≈ A + (2+4im)I ≈ convert(Array, (2+4im)I + H) @test convert(Array, H - 2I) ≈ A - 2I ≈ -convert(Array, 2I - H) @test convert(Array, -H) == -convert(Array, H) @test convert(Array, 2*(H + (2+4im)I)) ≈ 2Am + (4+8im)I @@ -77,17 +76,17 @@ let n = 10 @test H \ b ≈ A \ b ≈ H \ complex(b) @test H \ B ≈ A \ B ≈ H \ complex(B) @test (H - I) \ B ≈ (A - I) \ B - @test (H - (3+4im)I) \ B ≈ (Am - (3+4im)I) \ B + @test (H - (3+4im)I) \ B ≈ (A - (3+4im)I) \ B @test b' / H ≈ b' / A ≈ complex.(b') / H @test B' / H ≈ B' / A ≈ complex(B') / H @test B' / (H - I) ≈ B' / (A - I) - @test B' / (H - (3+4im)I) ≈ B' / (Am - (3+4im)I) - @test (H - (3+4im)I)' \ B ≈ (Am - (3+4im)I)' \ B - @test B' / (H - (3+4im)I)' ≈ B' / (Am - (3+4im)I)' + @test B' / (H - (3+4im)I) ≈ B' / (A - (3+4im)I) + @test (H - (3+4im)I)' \ B ≈ (A - (3+4im)I)' \ B + @test B' / (H - (3+4im)I)' ≈ B' / (A - (3+4im)I)' for shift in (0,1,3+4im) - @test det(H + shift*I) ≈ det(Am + shift*I) - @test logabsdet(H + shift*I) ≅ logabsdet(Am + shift*I) + @test det(H + shift*I) ≈ det(A + shift*I) + @test logabsdet(H + shift*I) ≅ logabsdet(A + shift*I) end end end