Mercurial > octave-antonio

changeset 20170:af2b7695f1c4 draft default tip @

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
gallery.m: clean and vectorize part of qmult. * scripts/special-matrix/gallery.m: make qmult() follow Octave guidelines, vectorize last for loop in qmult().
author Antonio Pino Robles <data.script93@gmail.com>
date Thu, 28 May 2015 18:32:47 +0200
parents 6f8c572f27fe
children
files scripts/special-matrix/gallery.m
diffstat 1 files changed, 3 insertions(+), 6 deletions(-) [+]
line wrap: on
line diff
--- 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