--- liboctave/CMatrix.cc.orig 2008-04-23 19:05:38.000000000 +0200 +++ liboctave/CMatrix.cc 2008-04-23 19:05:55.000000000 +0200 @@ -2967,8 +2967,8 @@ for (octave_idx_type i = 0; i < nc; i++) iperm(i) = i; // initialize to identity permutation - // leading permutations in forward order - for (octave_idx_type i = 0; i < (ilo-1); i++) + // trailing permutations must be done in reverse order + for (octave_idx_type i = nc - 1; i >= ihi; i--) { octave_idx_type swapidx = static_cast (dpermute(i)) - 1; octave_idx_type tmp = iperm(i); @@ -2976,8 +2976,8 @@ iperm(swapidx) = tmp; } - // trailing permutations must be done in reverse order - for (octave_idx_type i = nc - 1; i >= ihi; i--) + // leading permutations in forward order + for (octave_idx_type i = 0; i < (ilo-1); i++) { octave_idx_type swapidx = static_cast (dpermute(i)) - 1; octave_idx_type tmp = iperm(i);