Working version of Mgje

Author: Omnikar

lib.ua

@@ -37,14 +37,30 @@ GaussEliminate‼ ↚ ∧(
 MDet ← ◌⍣GaussEliminate‼(⍥¯𝄐𝄐/≠)𝄐𝄐× ⊙0⊙1
 
 # Perform Gauss-Jordan elimination
-Mgje ← ⍣GaussEliminate‼(|)(⍤"Elimination cannot be completed"0)
+Mgje ← +0⍥(
+  ↻₁⍣(⍜⊙⊢÷⟜(-⊞×⊙⊸⊢⊂0÷)°⊂⊸≡⊏⊢)◌⊚≠0⊸⊢
+)⊸⧻ ⊏⍏⊸=0
+# Mgje ← ⍣GaussEliminate‼(|)(⍤"Elimination cannot be completed"0)
+
 # Matrix Inverse
 MInv ← |1 ⌅(⍜⍉(↘÷2⊸⧻) ⍣Mgje(⍤"Matrix has no inverse"0) ˜≡⊂˙⊞=°⊏|MInv)
 
 ┌─╴test
   ⍤⤙≍ 5 ⁅₃ MDet [1_¯3_¯1 1_¯1_¯2 1_0_0]
   ⍤⤙≍ [0_0_1 ¯0.4_0.2_0.2 0.2_¯0.6_0.4] ⁅₃ MInv [1_¯3_¯1 1_¯1_¯2 1_0_0]
   ⍤⤙≍ ˙⊞=⇡3 ⁅₃ Mgje [1_¯3_¯1 1_¯1_¯2 1_0_0]
+
+  Mgje [1_5_5
+        4_5_¯10]
+  ⍤⤙≍ [1_0_¯5
+       0_1_2]
+
+  ⁅₉ Mgje [0_0_2_¯4_¯5
+           0_1_¯1_1_3
+           0_6_0_¯6_5]
+  ⍤⤙≍ [0_1_0_¯1_0
+       0_0_1_¯2_0
+       0_0_0_0_1]
 └─╴
 
 MSquarify ↚ ⍜△[.√⊢]