# HG changeset patch
# User dbateman
# Date 1151790726 0
# Node ID 50d43cdbec80a005f2ab01294e7e83c810ddfc61
# Parent  565d0cd4d9d04f69e288eb1d9006036792f9fbca
[project @ 2006-07-01 21:52:05 by dbateman]

diff --git a/liboctave/CSparse.cc b/liboctave/CSparse.cc
--- a/liboctave/CSparse.cc
+++ b/liboctave/CSparse.cc
@@ -782,9 +782,9 @@
 		      if (typ == MatrixType::Upper)
 			colUp = cidx(j+1) - 1;
 		      else
-			colUp = cidx(j) - 1;
+			colUp = cidx(j);
 		      Complex pivot = data(colUp);
-		      if (pivot == 0. || colUp != j) 
+		      if (pivot == 0. || ridx(colUp) != j) 
 			{
 			  (*current_liboctave_error_handler) 
 			    ("division by zero");
@@ -810,9 +810,9 @@
 		  if (typ == MatrixType::Upper)
 		    colUp = cidx(i+1) - 1;
 		  else
-		    colUp = cidx(i) - 1;
+		    colUp = cidx(i);
 		  Complex pivot = data(colUp);
-		  if (pivot == 0. || colUp != i) 
+		  if (pivot == 0. || ridx(colUp) != i) 
 		    {
 		      (*current_liboctave_error_handler) ("division by zero");
 		      goto inverse_singular;
@@ -878,7 +878,7 @@
 		      if (typ == MatrixType::Permuted_Upper)
 			pivot = data(cidx(jidx+1) - 1);
 		      else
-			pivot = data(cidx(jidx) - 1);
+			pivot = data(cidx(jidx));
 		      if (pivot == 0.) 
 			{
 			  (*current_liboctave_error_handler) 
@@ -894,7 +894,7 @@
 		  if (typ == MatrixType::Permuted_Upper)
 		    colUp = cidx(perm[iidx]+1) - 1;
 		  else
-		    colUp = cidx(perm[iidx]) - 1;		  
+		    colUp = cidx(perm[iidx]);		  
 
   		  Complex pivot = data(colUp);
 		  if (pivot == 0.)
diff --git a/liboctave/dSparse.cc b/liboctave/dSparse.cc
--- a/liboctave/dSparse.cc
+++ b/liboctave/dSparse.cc
@@ -858,9 +858,9 @@
 		      if (typ == MatrixType::Upper)
 			colUp = cidx(j+1) - 1;
 		      else
-			colUp = cidx(j) - 1;
+			colUp = cidx(j);
 		      double pivot = data(colUp);
-		      if (pivot == 0. || colUp != j) 
+		      if (pivot == 0. || ridx(colUp) != j) 
 			{
 			  (*current_liboctave_error_handler) 
 			    ("division by zero");
@@ -886,9 +886,9 @@
 		  if (typ == MatrixType::Upper)
 		    colUp = cidx(i+1) - 1;
 		  else
-		    colUp = cidx(i) - 1;
+		    colUp = cidx(i);
 		  double pivot = data(colUp);
-		  if (pivot == 0. || colUp != i) 
+		  if (pivot == 0. || ridx(colUp) != i) 
 		    {
 		      (*current_liboctave_error_handler) ("division by zero");
 		      goto inverse_singular;
@@ -954,7 +954,7 @@
 		      if (typ == MatrixType::Permuted_Upper)
 			pivot = data(cidx(jidx+1) - 1);
 		      else
-			pivot = data(cidx(jidx) - 1);
+			pivot = data(cidx(jidx));
 		      if (pivot == 0.) 
 			{
 			  (*current_liboctave_error_handler) 
@@ -970,7 +970,7 @@
 		  if (typ == MatrixType::Permuted_Upper)
 		    colUp = cidx(perm[iidx]+1) - 1;
 		  else
-		    colUp = cidx(perm[iidx]) - 1;		  
+		    colUp = cidx(perm[iidx]);
 
 		  double pivot = data(colUp);
 		  if (pivot == 0.)