From 098f2bbdec263152ee07a490cee49fd85dadea76 Mon Sep 17 00:00:00 2001
From: vijay gopal chilkuri <vijay.gopal.c@gmail.com>
Date: Wed, 21 Dec 2016 11:53:13 +0100
Subject: [PATCH] trying to read array elements

---
 src/ex1.c | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/src/ex1.c b/src/ex1.c
index af1e530..36f0a77 100644
--- a/src/ex1.c
+++ b/src/ex1.c
@@ -211,17 +211,19 @@ int main(int argc,char **argv)
 
   if (nconv>0) {
 
-    for (i=0;i<nev;i++) {
       /*
         Get converged eigenpairs: i-th eigenvalue is stored in kr (real part) and
         ki (imaginary part)
       */
-      ierr = EPSGetEigenpair(eps,i,&kr,&ki,xr,xi);CHKERRQ(ierr);
+      ierr = EPSGetEigenpair(eps,1,&kr,&ki,xr,xi);CHKERRQ(ierr);
 
 	  ierr = VecNorm(xr, NORM_2, &norm);CHKERRQ(ierr);
 	  PetscPrintf(PETSC_COMM_WORLD," Norm = %18f \n", (double)norm);
+  	  for (i=Istart; i<Iend; i+=getdata.nnz) {
+		  ierr = VecGetValue(xr, 1, i, value);
+		  PetscPrintf(PETSC_COMM_WORLD," Element # = %d Value = %18f \n", i, value[0]); 
+	  }
 
-    }
     ierr = PetscPrintf(PETSC_COMM_WORLD,"\n");CHKERRQ(ierr);
   }