Mercurial > hg > octave-lyh

diff liboctave/Sparse.cc @ 9177:39be2c4531c8

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix sparse indexing bug
author Jaroslav Hajek <highegg@gmail.com>
date Tue, 05 May 2009 15:18:55 +0200
parents 542015fada9e
children c6edba80dfae
line wrap: on
line diff
--- a/liboctave/Sparse.cc
+++ b/liboctave/Sparse.cc
@@ -2715,13 +2715,17 @@
 		      i++;
 		      while (ii < nc && c_lhs.cidx(ii+1) <= i)
 			ii++;
+                      tmp.xridx (kk++) = 0;
 		    }
 		  else
 		    {
 		      while (ic <= jj)
 			tmp.xcidx (ic++) = kk;
 		      if (scalar_non_zero)
-			tmp.xdata (kk) = scalar;
+                        {
+                          tmp.xdata (kk) = scalar;
+                          tmp.xridx (kk++) = 0;
+                        }
 		      if (ii == jj)
 			{
 			  i++;
@@ -2732,7 +2736,6 @@
 		      if (j < n)
 			jj = lhs_idx.elem(j);
 		    }
-		  tmp.xridx (kk++) = 0;
 		}
 
 	      for (octave_idx_type iidx = ic; iidx < max_idx+1; iidx++)