# HG changeset patch
# User jwe
# Date 791321085 0
# Node ID fb8ca633c439c5897e75a07debadfa437316ba14
# Parent  d1a1608f1028427687e0577a916b5988069f11ed
[project @ 1995-01-28 19:24:32 by jwe]

diff -r d1a1608f1028 -r fb8ca633c439 scripts/general/tril.m
--- a/scripts/general/tril.m	Fri Jan 27 22:26:33 1995 +0000
+++ b/scripts/general/tril.m	Sat Jan 28 19:24:45 1995 +0000
@@ -27,24 +27,22 @@
 
   if (nargin > 0)
     [nr, nc] = size (x);
-    retval = x;
+    retval = zeros (nr, nc);
   endif
 
   if (nargin == 1)
     k = 0;
   elseif (nargin == 2)
-    max_nr_nc = max (nr, nc);
-    if ((k > 0 && k > nr - 1) || (k < 0 && k < 1 - nc))
+    if ((k > 0 && k > nc) || (k < 0 && k < -nr))
       error ("tril: requested diagonal out of range")
     endif
   else
     usage ("tril (x [, k])");
   endif
 
-  for i = 1:nr
-    for j = i+1-k:nc
-      retval (i, j) = 0.0;
-    endfor
+  for j = 1 : min (nc, nc + abs (k))
+    nr_limit = max (1, j-k);
+    retval (nr_limit:nr, j) = x (nr_limit:nr, j);
   endfor
 
 endfunction
diff -r d1a1608f1028 -r fb8ca633c439 scripts/general/triu.m
--- a/scripts/general/triu.m	Fri Jan 27 22:26:33 1995 +0000
+++ b/scripts/general/triu.m	Sat Jan 28 19:24:45 1995 +0000
@@ -27,24 +27,22 @@
 
   if (nargin > 0)
     [nr, nc] = size (x);
-    retval = x;
+    retval = zeros (nr, nc);
   endif
 
   if (nargin == 1)
     k = 0;
   elseif (nargin == 2)
-    max_nr_nc = max (nr, nc);
-    if ((k > 0 && k > nc - 1) || (k < 0 && k < 1 - nr))
+    if ((k > 0 && k > nc) || (k < 0 && k < -nr))
       error ("triu: requested diagonal out of range")
     endif
   else
     usage ("triu (x [, k])");
   endif
 
-  for j = 1:nc
-    for i = j+1-k:nr
-      retval (i, j) = 0.0;
-    endfor
+  for j = max (1, k) : nc
+    nr_limit = min (nr, j-k);
+    retval (1:nr_limit, j) = x (1:nr_limit, j);
   endfor
 
 endfunction