# HG changeset patch
# User Jaroslav Hajek <highegg@gmail.com>
# Date 1262619536 -3600
# Node ID 38600f8cba837ae46a4f778af6b3c83811995339
# Parent  8f7f325fa6785759ecc0929a934ba796bd3a37d7
optimize strmatch

diff -r 8f7f325fa678 -r 38600f8cba83 scripts/ChangeLog
--- a/scripts/ChangeLog	Mon Jan 04 12:32:07 2010 +0100
+++ b/scripts/ChangeLog	Mon Jan 04 16:38:56 2010 +0100
@@ -1,3 +1,7 @@
+2010-01-04  Jaroslav Hajek  <highegg@gmail.com>
+
+	* strings/strmatch.m: Optimize.
+
 2010-01-04  Jaroslav Hajek  <highegg@gmail.com>
 
 	* strings/strrep.m: Remove.
diff -r 8f7f325fa678 -r 38600f8cba83 scripts/strings/strmatch.m
--- a/scripts/strings/strmatch.m	Mon Jan 04 12:32:07 2010 +0100
+++ b/scripts/strings/strmatch.m	Mon Jan 04 16:38:56 2010 +0100
@@ -1,5 +1,6 @@
 ## Copyright (C) 2000, 2005, 2006, 2007, 2009 Paul Kienzle
 ## Copyright (C) 2003 Alois Schloegl
+## Copyright (C) 2010 VZLU Prague
 ##
 ## This file is part of Octave.
 ##
@@ -22,8 +23,9 @@
 ## Return indices of entries of @var{a} that match the string @var{s}.
 ## The second argument @var{a} may be a string matrix or a cell array of
 ## strings.  If the third argument @code{"exact"} is not given, then
-## @var{s} only needs to match @var{a} up to the length of @var{s}.  Nul
-## characters match blanks.  Results are returned as a column vector. 
+## @var{s} only needs to match @var{a} up to the length of @var{s}.
+## Trailing whitespace is ignored.
+## Results are returned as a column vector. 
 ## For example:
 ##
 ## @example
@@ -50,38 +52,50 @@
     print_usage ();
   endif
 
-  [nr, nc] = size (A);
-  nel = numel (A);
+  if (! ischar (s))
+    error ("strmatch: first argument must be a string");
+  endif
+  
+  ## Truncate trailing whitespace.
+  s = strtrimr (s);
+
+  len = length (s);
+
+  exact = nargin == 3 && ischar (exact) && strcmp (exact, "exact");
+
   if (iscell (A))
-    match = zeros (nel, 1);
-    if (nargin > 2)
-      for k = 1:nel
-	match(k) = strcmp (s, A{k}); 
-      endfor
-    else
-      for k = 1:nel
-	match(k) = strncmp (s, A{k}, length (s));
-      endfor
+    idx = find (strncmp (s, A, len));
+    if (exact)
+      ## We can't just use strcmp, because we need to ignore whitespace.
+      B = cellfun (@strtrimr, A(idx), "UniformOutput", false);
+      idx = idx (strcmp (s, B));
     endif
-    idx = find (match);
-  elseif (length (s) > nc)
-    idx = [];
-  else
-    if (nargin == 3 && length(s) < nc)
-      s(1,nc) = " ";
+  elseif (ischar (A))
+    [nr, nc] = size (A);
+    if (len > nc)
+      idx = [];
+    else
+      match = all (bsxfun (@eq, A(:,1:len), s), 2);
+      if (exact)
+        AA = A(:,len+1:nc);
+        match &= all (AA == "\0" | AA == " ", 2);
+      endif
+      idx = find (match);
     endif
-    s(s == 0) = " ";
-    A(A == 0) = " ";
-    match = s(ones(size(A,1),1),:) == A(:,1:length(s));
-    if (length (s) == 1)
-      idx = find (match);
-    else
-      idx = find (all (match')');
-    endif
+  else
+    error ("strmatch: second argument must be a string or cell array of strings");
   endif
     
 endfunction 
 
+## Removes nuls and blanks from the end of the array
+function s = strtrimr (s)
+  i = find (s == "\0" | s == " ", 1, "last");
+  if (i)
+    s = s(1:i);
+  endif
+endfunction
+
 %!error <Invalid call to strmatch> strmatch();
 %!error <Invalid call to strmatch> strmatch("a", "aaa", "exact", 1);
 %!assert (strmatch("a", {"aaa", "bab", "bbb"}), 1);