<div dir="ltr"><div><div><div>Ok. I used -st_type sinvert only, <br><br></div>I am comparing it eigs after normalization. Even now the eigenvectors differ. <br><br></div>Whether I should use a different eps solver ? <br><br>
</div>Pls let me know. <br><br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Mar 23, 2014 at 9:46 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 23/03/2014, a las 12:56, venkatesh g escribió:<br>
<div class=""><br>
> Hi,<br>
><br>
> I am using SLEPC example EX7.C, I am using complex binary matrices A and B written from matlab using Petscbinarywrite.m.<br>
><br>
> I ran the example using "./ex7 -f1 A -f2 B -eps_type krylovschur -st_type sinvert -evecs out.mat -eps_smallest_magnitude"<br>
<br>
</div>You should not use -st_type sinvert together with -eps_smallest_magnitude. Read the documentation.<br>
<div class=""><br>
><br>
> My plotted my eigenvectors using Petscbinaryread and comparing them to MATLAB's eigs command.<br>
><br>
> It is wrong. It is right only if I use real matrices A and B.<br>
><br>
> I compiled Petsc for complex scalars. How to get the correct eigenvectors for complex matrices, please note my eigenvalue is right but not my eigenvector.<br>
<br>
</div>How do you know the eigenvectors are wrong? Did you check the residuals? If you want to compare the eigenvectors to the output of eigs, you may need to normalize them in the same way, for instance with v=v/max(v).<br>
<span class="HOEnZb"><font color="#888888"><br>
Jose<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
><br>
> Any help is greatly appreciated! Kindly help me out.<br>
><br>
> cheers,<br>
> Venkatesh<br>
> Sr. Research Fellow<br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div>