diff_out_elements

Index: eth.c
===================================================================
--- eth.c	(revision 1435)
+++ eth.c	(working copy)
@@ -309,6 +309,8 @@
   PetscScalar element, value, vecval;
   for (sitestride1 = Istart; sitestride1 < Iend; sitestride1++)
     {
+      ierr = PetscPrintf(PETSC_COMM_WORLD,"\n");
+      CHKERRQ (ierr);
       sitestride1_vec[0] = sitestride1;
       for (sitestride2 = 0; sitestride2 < matsize; sitestride2++)
 	{
@@ -358,6 +360,11 @@
 	      colvalues[sitestride2] = element;
 	    }
 	}
+      CHKERRQ (ierr);
+      for(count=0;count<matsize;count++){
+	ierr = PetscPrintf(PETSC_COMM_WORLD,"[%d %d]: %2.6lf ",sitestride1, count,colvalues[count]); 
+	CHKERRQ (ierr);
+      }
       //Insert the array of colvalues to the sitestride1^th row of H   
       ierr =
 	MatSetValues (AVG_BDIBJ, 1, &sitestride1, matsize, idx, colvalues,