# HG changeset patch
# User Antonio Pino Robles <data.script93@gmail.com>
# Date 1432830767 -7200
# Node ID af2b7695f1c41a49147a1f7d5436e93d6d9edb5b
# Parent  6f8c572f27fe9a2ec231dc452c6e7ff7ca2c8ef3
gallery.m: clean and vectorize part of qmult.

* scripts/special-matrix/gallery.m: make qmult() follow Octave guidelines, vectorize last for loop in qmult().

diff -r 6f8c572f27fe -r af2b7695f1c4 scripts/special-matrix/gallery.m
--- a/scripts/special-matrix/gallery.m	Thu May 28 00:08:47 2015 +0200
+++ b/scripts/special-matrix/gallery.m	Thu May 28 18:32:47 2015 +0200
@@ -2015,7 +2015,7 @@
     error ("gallery: 1 to 2 arguments are required for pei matrix.");
   elseif (! isnumeric (n) || ! isscalar (n) || fix (n) != n)
     error ("gallery: N must be an integer for pei matrix.");
-  elseif (! isnumeric (alpha) || ! isscalar (alpha)) # change w to alpha 
+  elseif (! isnumeric (alpha) || ! isscalar (alpha))
     error ("gallery: ALPHA must be a scalar for pei matrix.");
   endif
 
@@ -2876,7 +2876,7 @@
   #  Handle scalar A.
   if (max (n,m) == 1)
      n = A;
-     A = eye(n);
+     A = eye (n);
   endif
   
   d = zeros (n);
@@ -2896,10 +2896,7 @@
   endfor
   
   # Tidy up signs.
-  for i=1:n-1
-      A(i,:) = d(i) * A(i,:);
-  endfor
-  A(n,:) = A(n,:) * sign (randn);
+  A = [d'; sign(randn)] .* A;
   B = A;
 endfunction