Define factorize(Adjoint) to make e.g. inv(Adjoint) work (in most cases)

Author: andreasnoack

stdlib/LinearAlgebra/src/dense.jl

@@ -1222,6 +1222,7 @@ function factorize(A::StridedMatrix{T}) where T
     end
     qrfact(A, Val(true))
 end
+factorize(A::Adjoint) = adjoint(factorize(parent(A)))
 
 ## Moore-Penrose pseudoinverse
 

stdlib/LinearAlgebra/test/dense.jl

@@ -847,4 +847,13 @@ end
     end
 end
 
+@testset "inverse of Adjoint" begin
+    n = 5
+    A = randn(n, n)
+    @test inv(A')*A' ≈ I
+    B = complex.(A, randn(n, n))
+    B = B + transpose(B)
+    @test_broken inv(B')*B' ≈ I
+end
+
 end # module TestDense