--- liboctave/dMatrix.cc.orig	2008-04-23 19:05:43.000000000 +0200
+++ liboctave/dMatrix.cc	2008-04-23 19:05:49.000000000 +0200
@@ -2592,17 +2592,17 @@
   for (octave_idx_type i = 0; i < nc; i++)
     iperm(i) = i;  // 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<octave_idx_type> (dpermute(i)) - 1;
       octave_idx_type tmp = iperm(i);
-      iperm(i) = iperm (swapidx);
+      iperm(i) = iperm(swapidx);
       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<octave_idx_type> (dpermute(i)) - 1;
       octave_idx_type tmp = iperm(i);