<div dir="ltr"><div><div>Thanks a lot. I did what you told. <br></div><br>It is working now like a charm. <br><br></div>However, I have another query, which I will post separately. <br><br></div><div class="gmail_extra"><br>
<br><div class="gmail_quote">On Mon, Mar 24, 2014 at 3:21 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 10:28, venkatesh g escribió:<br>
<div class=""><br>
> Ok. I am doing the following<br>
><br>
> 1. Running "./ex7 -f1 A2 -f2 B2 -st_type sinvert -evecs VC"<br>
> 2. In MATLAB I do "vc=PetscBinaryRead('VC');"<br>
><br>
> The output is one Eigenvector in VC which is 800x1.<br>
><br>
> 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>
><br>
> I also changed the line 179 of ex7.c like u said.<br>
><br>
> Pls let me know.<br>
><br>
<br>
</div>You should do<br>
<br>
  vc=PetscBinaryRead('VC','complex',true);<br>
<br>
If you want to load several eigenvectors, then<br>
<br>
  vc=PetscBinaryRead('VC','complex',true,'cell',nconv);<br>
<br>
where nconv is the number of converged eigenpairs reported by SLEPc. In that case, the result is a cell array so vc{1} is the first eigenvector, vc{2} the second one, and so on.<br>
<span class="HOEnZb"><font color="#888888"><br>
Jose<br>
<br>
</font></span></blockquote></div><br></div>