+-- | Returns True if the given matrix is upper-triangular, and False
+-- otherwise.
+--
+-- Examples:
+--
+-- >>> let m = fromList [[1,0],[1,1]] :: Mat2 Int
+-- >>> is_upper_triangular m
+-- False
+--
+-- >>> let m = fromList [[1,2],[0,3]] :: Mat2 Int
+-- >>> is_upper_triangular m
+-- True
+--
+is_upper_triangular :: (Eq a, Ring.C a, Arity m, Arity n)
+ => Mat m n a -> Bool
+is_upper_triangular m =
+ and $ concat results
+ where
+ results = [[ test i j | i <- [0..(nrows m)-1]] | j <- [0..(ncols m)-1] ]
+
+ test :: Int -> Int -> Bool
+ test i j
+ | i <= j = True
+ | otherwise = m !!! (i,j) == 0
+
+
+-- | Returns True if the given matrix is lower-triangular, and False
+-- otherwise.
+--
+-- Examples:
+--
+-- >>> let m = fromList [[1,0],[1,1]] :: Mat2 Int
+-- >>> is_lower_triangular m
+-- True
+--
+-- >>> let m = fromList [[1,2],[0,3]] :: Mat2 Int
+-- >>> is_lower_triangular m
+-- False
+--
+is_lower_triangular :: (Eq a,
+ Ring.C a,
+ Arity m,
+ Arity n)
+ => Mat m n a
+ -> Bool
+is_lower_triangular = is_upper_triangular . transpose
+
+
+-- | Returns True if the given matrix is triangular, and False
+-- otherwise.
+--
+-- Examples:
+--
+-- >>> let m = fromList [[1,0],[1,1]] :: Mat2 Int
+-- >>> is_triangular m
+-- True
+--
+-- >>> let m = fromList [[1,2],[0,3]] :: Mat2 Int
+-- >>> is_triangular m
+-- True
+--
+-- >>> let m = fromList [[1,2],[3,4]] :: Mat2 Int
+-- >>> is_triangular m
+-- False
+--
+is_triangular :: (Eq a,
+ Ring.C a,
+ Arity m,
+ Arity n)
+ => Mat m n a
+ -> Bool
+is_triangular m = is_upper_triangular m || is_lower_triangular m
+
+
+-- | Return the (i,j)th minor of m.