<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><span style="white-space: pre-wrap; font-size: 12pt;">>If the eigenvector is complex then of course the eigenvalue is complex as well </span><div><span style="white-space: pre-wrap; font-size: 12pt;">>(I assume your matrix is real non-symmetric). You have to get both the real and >imaginary parts of the eigenvalue.</span></div><div><span style="white-space: pre-wrap; font-size: 12pt;">></span><a href="http://www.grycap.upv.es/slepc/documentation/current/docs/manualpages/EPS/EPSGetEigenpair.html" style="white-space: pre-wrap; font-size: 12pt;">http://www.grycap.upv.es/slepc/documentation/current/docs/manualpages/EPS/EPSGetEigenpair.html</a></div><div><span style="white-space: pre-wrap; font-size: 12pt;">>An alternative is to do all the computation in complex arithmetic (configure --with-scalar-type=complex).</span></div><div><span style="font-size: 12pt; white-space: pre-wrap;">>Jose</span></div><div><span style="font-size: 12pt; white-space: pre-wrap;"><br></span></div><div><span style="font-size: 12pt; white-space: pre-wrap;">Hi Jose,</span></div><div><span style="font-size: 12pt; white-space: pre-wrap;"><br></span></div><div><span style="font-size: 12pt; white-space: pre-wrap;">Probably you misunderstood my question my problem is: </span><span style="white-space: pre-wrap;">ierr = EPSGetEigenvalue(eps,i,&eigen_r,&eigen_i);CHKERRQ(ierr);</span></div><div><span style="white-space: pre-wrap;">returns all eigenvalues for each iteration where there is no imaginary parts (and set to zeros in eigen_i as expected).</span></div><div><span style="white-space: pre-wrap;"><br></span></div><div><span style="white-space: pre-wrap;">I cross checked that with another algorithm which proves that all the eigenvalues are correct (no imaginary parts),but </span></div><div><span style="white-space: pre-wrap;">glitch is with trying to fetch the eigenvectors. Eigenvectors obtained from the code below are mostly filled with zeros </span></div><div><span style="white-space: pre-wrap;">and only the 11th row(Slepc returns 11 eigenvalue) has real and </span><span style="white-space: pre-wrap; font-size: 12pt;">imaginary part which is not possible.</span></div><div><br></div> <span style="white-space: pre-wrap;">for (j=0;j<test_int;j++)
{
ierr = EPSGetEigenvector(eps,j,&vec[0],&vec1[0]);CHKERRQ(ierr);
PetscPrintf(MPI_COMM_WORLD," %d || %4.8f || %4.8f \n",j,vec[j],vec1[j]);
}</span> <div><span style="white-space: pre-wrap;"> </span></div><div><span style="white-space: pre-wrap;">I also provide you the main code and required files.</span></div><div><span style="white-space: pre-wrap;"><br></span></div><div><span style="white-space: pre-wrap;">Your help will be appreciated,</span></div><div><span style="white-space: pre-wrap;"><br></span></div><div><span style="white-space: pre-wrap;">Regards,</span></div> </div></body>
</html>