<div dir="ltr"><div><div><div><div><div><div>Ok. I am doing the following<br><br></div>1. Running "./ex7 -f1 A2 -f2 B2 -st_type sinvert -evecs VC"<br></div>2. In MATLAB I do "vc=PetscBinaryRead('VC');" <br>
<br></div><div>The output is one Eigenvector in VC which is 800x1.<br><br></div>And I did the normalization vc=vc/max(vc) and compared with the normalized original vector.. Also I must tell that if A2 and B2 are real, it works.<br>
</div><br></div>I also changed the line 179 of ex7.c like u said. <br><br></div>Pls let me know.<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Mar 24, 2014 at 2:24 PM, Jose E. Roman <span dir="ltr"><<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
El 24/03/2014, a las 05:59, venkatesh g escribió:<br>
<div class=""><br>
> Ok. I used -st_type sinvert only,<br>
><br>
> I am comparing it eigs after normalization. Even now the eigenvectors differ.<br>
><br>
> Whether I should use a different eps solver ?<br>
><br>
> Pls let me know.<br>
><br>
<br>
</div>It works for me, so you should provide more details about how you are loading the eigenvectors in Matlab and comparing them.<br>
<br>
Note that the complex version of the example also writes the xi vector (which should not) so you will have zero vectors interleaved with the eigenvectors. If you want to avoid this, make the following change (line 179 of ex7.c):<br>

<br>
#if !defined(PETSC_USE_COMPLEX)<br>
       if (!ishermitian) { ierr = VecView(xi,viewer);CHKERRQ(ierr); }<br>
#endif<br>
<span class="HOEnZb"><font color="#888888"><br>
Jose<br>
<br>
</font></span></blockquote></div><br></div>