From 408943c4e0e0bb54d110073ba89c003c31e82142 Mon Sep 17 00:00:00 2001 From: "Viral B. Shah" Date: Sat, 13 Dec 2025 15:33:05 -0500 Subject: [PATCH] Fix #1496 --- test/blas.jl | 13 ++++++++----- 1 file changed, 8 insertions(+), 5 deletions(-) diff --git a/test/blas.jl b/test/blas.jl index b7f4a03a..6cda7fbe 100644 --- a/test/blas.jl +++ b/test/blas.jl @@ -721,15 +721,18 @@ end @test BLAS.her!('L', real(elty(2)), x, A) isa WrappedArray{elty,2} @test A == WrappedArray(elty[5 2+2im; 11+3im 20]) # Level 3 - A = WrappedArray(elty[1+im 2+2im; 3+3im 4+4im]) + # Hermitian matrices require real diagonal elements + A = WrappedArray(elty[1 2+2im; 2-2im 4]) B = WrappedArray(elty[1+im 2+2im; 3+3im 4+4im]) - C = WrappedArray(elty[1+im 2+2im; 3+3im 4+4im]) + C = WrappedArray(elty[1 2+2im; 2-2im 4]) @test BLAS.hemm!('L', 'U', elty(2), A, B, elty(1), C) isa WrappedArray{elty,2} - @test C == WrappedArray([3+27im 6+38im; 35+27im 52+36im]) + @test C == WrappedArray([3+26im 6+38im; 34+22im 52+32im]) + C = WrappedArray(elty[1 2+2im; 2-2im 4]) # reset C to Hermitian @test BLAS.herk!('U', 'N', real(elty(2)), A, real(elty(1)), C) isa WrappedArray{elty,2} - @test C == WrappedArray([23 50+38im; 35+27im 152]) + @test C == WrappedArray([19 22+22im; 2-2im 52]) + C = WrappedArray(elty[1 2+2im; 2-2im 4]) # reset C to Hermitian @test BLAS.her2k!('U', 'N', elty(2), A, B, real(elty(1)), C) isa WrappedArray{elty,2} - @test C == WrappedArray([63 138+38im; 35+27im 352]) + @test C == WrappedArray([37 56+20im; 2-2im 68]) end end