Use generic methods for IsDiagonalMat, Is{Upper,Lower}TriangularMat

gap-system · Jan 6, 2019 · e1a1fb6 · e1a1fb6
1 parent 0ec0771
commit e1a1fb6
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 65 deletions.
diff --git a/lib/matobjplist.gi b/lib/matobjplist.gi
@@ -1160,68 +1160,6 @@ InstallMethod( Randomize, "for a mutable plist matrix",
     od;
   end );
 
-InstallMethod( IsDiagonalMat, "for a plist matrix",
-  [ IsPlistMatrixRep ],
-  function( m )
-    local i,j,l,n;
-    if m![RLPOS] <> Length(m![ROWSPOS]) then
-        ErrorNoReturn("IsDiagonalMat: Matrix not square");
-    fi;
-    n := m![RLPOS];
-    l := m![ROWSPOS];
-    for i in [1..n] do
-        for j in [1..i-1] do
-            if not IsZero(l[i,j]) then
-                return false;
-            fi;
-        od;
-        for j in [i+1..n] do
-            if not IsZero(l[i,j]) then
-                return false;
-            fi;
-        od;
-    od;
-    return true;
-  end );
-
-InstallMethod( IsLowerTriangularMat, "for a plist matrix",
-  [ IsPlistMatrixRep ],
-  function( m )
-    local i,j,l,n;
-    if m![RLPOS] <> Length(m![ROWSPOS]) then
-        ErrorNoReturn("IsLowerTriangularMat: Matrix not square");
-    fi;
-    n := m![RLPOS];
-    l := m![ROWSPOS];
-    for i in [1..n] do
-        for j in [i+1..n] do
-            if not IsZero(l[i,j]) then
-                return false;
-            fi;
-        od;
-    od;
-    return true;
-  end );
-
-InstallMethod( IsUpperTriangularMat, "for a plist matrix",
-  [ IsPlistMatrixRep ],
-  function( m )
-    local i,j,l,n;
-    if m![RLPOS] <> Length(m![ROWSPOS]) then
-        ErrorNoReturn("IsUpperTriangularMat: Matrix not square");
-    fi;
-    n := m![RLPOS];
-    l := m![ROWSPOS];
-    for i in [1..n] do
-        for j in [1..i-1] do
-            if not IsZero(l[i,j]) then
-                return false;
-            fi;
-        od;
-    od;
-    return true;
-  end );
-
 InstallMethod( TransposedMatMutable, "for a plist matrix",
   [ IsPlistMatrixRep ],
   function( m )

diff --git a/lib/matrix.gi b/lib/matrix.gi
@@ -165,7 +165,9 @@ InstallOtherMethod( IsDiagonalMat,
     [ IsMatrix ],
     function( mat )
    local  i, j,z;
-    if IsEmpty(mat) then return true;fi;
+    if NrRows(mat) = 0 then
+        return true;
+    fi;
     z:=ZeroOfBaseDomain(mat);
     for i  in [ 1 .. NrRows( mat ) ]  do
         for j  in [ 1 .. NrCols( mat ) ]  do
@@ -189,7 +191,9 @@ InstallOtherMethod( IsUpperTriangularMat,
     [ IsMatrix ],
     function( mat )
     local  i, j,z;
-    if IsEmpty(mat) then return true;fi;
+    if NrRows(mat) = 0 then
+        return true;
+    fi;
     z:=ZeroOfBaseDomain(mat);
     for i  in [ 1 .. NrRows( mat ) ]  do
         for j  in [ 1 .. i-1]  do
@@ -210,7 +214,9 @@ InstallOtherMethod( IsLowerTriangularMat,
     [ IsMatrix ],
     function( mat )
     local  i, j,z;
-    if IsEmpty(mat) then return true;fi;
+    if NrRows(mat) = 0 then
+        return true;
+    fi;
     z:=ZeroOfBaseDomain(mat);
     for i  in [ 1 .. NrRows( mat ) ]  do
         for j  in [ i+1 .. NrCols( mat ) ]  do