Backport "Fix structure test for strided matrices"

stdlib/LinearAlgebra/src/dense.jl

@@ -1341,25 +1341,26 @@ This returns a `5×5 Bidiagonal{Float64}`, which can now be passed to other line
 (e.g. eigensolvers) which will use specialized methods for `Bidiagonal` types.
 """
 function factorize(A::AbstractMatrix{T}) where T
+    require_one_based_indexing(A)
     m, n = size(A)
     if m == n
         if m == 1 return A[1] end
         utri    = true
         utri1   = true
         herm    = true
         sym     = true
-        for j = 1:n-1, i = j+1:m
-            if utri1
+        for j = 1:n, i = j:m
+            if (j < n) && (i > j) && utri1 # indices are off-diagonal
                 if A[i,j] != 0
                     utri1 = i == j + 1
                     utri = false
                 end
             end
             if sym
-                sym &= A[i,j] == A[j,i]
+                sym &= A[i,j] == transpose(A[j,i])
             end
             if herm
-                herm &= A[i,j] == conj(A[j,i])
+                herm &= A[i,j] == adjoint(A[j,i])
             end
             if !(utri1|herm|sym) break end
         end
@@ -1372,7 +1373,7 @@ function factorize(A::AbstractMatrix{T}) where T
         if ltri1
             for i = 1:n-1
                 if A[i,i+1] != 0
-                    ltri &= false
+                    ltri = false
                     break
                 end
             end

stdlib/LinearAlgebra/test/dense.jl

@@ -1234,4 +1234,10 @@ Base.:+(x::TypeWithZero, ::TypeWithoutZero) = x
     @test diagm(0 => [TypeWithoutZero()]) isa Matrix{TypeWithZero}
 end
 
+@testset "structure of dense matrices" begin
+    # A is neither triangular nor symmetric/Hermitian
+    A = [1 im 2; -im 0  3; 2 3 im]
+    @test factorize(A) isa LU{ComplexF64, Matrix{ComplexF64}, Vector{Int}}
+end
+
 end # module TestDense