[Octave-bug-tracker] [bug #51241] Infinite loop in normest1

[Rik]

Follow-up Comment #9, bug #51241 (project octave):

I don't understand the algorithm, but here is a dump of the variables.


debug> i
i =  2
debug> imax
imax =  3
debug> p
p = 0
debug> S
S =

   1   1  -1
  -1   1   1
  -1  -1   1

debug> Sold
Sold =

   1   1  -1
   1  -1  -1
   1   1   1



The issue is this expression


any (abs (S(:,i)' * Sold) == n)


which is always true.  I wonder if the code should read


if (expression)
  ## i-th column of S parallel to previous one
  do stuff...
endif
i++;


rather than


if (expression)
  ## i-th column of S parallel to previous one
  do stuff...
else
  i++;
endif




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?51241>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/