handle non-square case like LA

JuliaArrays · Jul 2, 2024 · 3847275 · 3847275
1 parent ea7594e
commit 3847275
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/src/inv.jl b/src/inv.jl
@@ -83,4 +83,10 @@ end
     end
 end
 
-inv(LUp::LU) = LUp.U \ (LUp.L \ typeof(LUp.L.data)(I)[LUp.p,:])
+function inv(LUp::LU)
+    if !(LUp.L isa LowerTriangular)
+        checksquare(LUp.L)
+        checksquare(LUp.U)
+    end
+    LUp.U \ (LUp.L \ typeof(parent(LUp.L))(I)[LUp.p,:])
+end