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

[Rik]

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

Within normest1, the problem is this section of code


          while (i <= imax)
            ## The maximum number of parallel columns of length n with
entries
            ## {-1,1} is 2^(n-1). Therefore, if the number of columns of S is
            ## greater than 2^(n-1), for sure some of them are parallel to
some
            ## columns of Sold. Don't even try to change them (i <= 2^(n-1)).
            ## Now, check if S(:,i) is parallel to any previous column of S
            p = (any (abs (S(:,i)' * S(:,1:i-1)) == n));
            if (p || (any (abs (S(:,i)' * Sold) == n)))
              ## i-th column of S parallel to a previous or to a
              ## column of Sold: change it.
              S(:,i) = sign (2*rand (n, 1)-1);
            else
              i++;
            endif
          endwhile


If the first branch of the if statement is taken, then the variable i is never
incremented.


    _______________________________________________________

Reply to this item at:

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

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