# HG changeset patch
# User jwe
# Date 1049749417 0
# Node ID 51cd16a73307a6e1bbf335c983a40af293411e6b
# Parent  fa80721998d68cbd517cc81f32c01ef0cdb4ab50
[project @ 2003-04-07 21:03:37 by jwe]

diff --git a/scripts/ChangeLog b/scripts/ChangeLog
--- a/scripts/ChangeLog
+++ b/scripts/ChangeLog
@@ -1,3 +1,8 @@
+2003-04-07  John W. Eaton  <jwe@bevo.che.wisc.edu>
+
+	* control/base/__bodquist__.m: Don't convert pdbig and fdbig to
+	column vectors.
+
 2003-03-24  Quentin Spencer <qspencer@ieee.org>
 
 	* linear-algebra/null.m: Handle empty matrix arg.
diff --git a/scripts/control/base/__bodquist__.m b/scripts/control/base/__bodquist__.m
--- a/scripts/control/base/__bodquist__.m
+++ b/scripts/control/base/__bodquist__.m
@@ -102,12 +102,12 @@
     dmag = 0.2;         # desired max change in magnitude
     while(pcnt)
       pd = abs(diff(phase));                    # phase variation
-      pdbig = vec(find(pd > dphase));
+      pdbig = find(pd > dphase);
 
       lp = length(f);  lp1 = lp-1;              # relative variation
       fd = abs(diff(f));
       fm = max(abs([f(1:lp1); f(2:lp)]));
-      fdbig = vec(find(fd > fm/10));
+      fdbig = find(fd > fm/10);
 
       bigpts = union(fdbig, pdbig);