<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial">Hello <br>every
 one. I am Wei. I want to use slepc to solver my quadratic eigenvalue 
problem. and now I start with the tutorial ex16 and compare with polyeig
 function matlab. but the two give different results.<br>here is the results from slepc<br><br>Number of requested eigenvalues:  6<br>Stopping condition: tol= 1e-08 , maxit= 100<br>Number of converged approximate eigenpairs:  8<br><br>          k           ||(k^2M+Ck+K)x||/||kx|| <br>-------------------- -------------------------<br>-7.85046229415e-16 + 2.81636155327 J, 3.31272023702e-09<br>-7.85046229415e-16 + -2.81636155327 J, 3.31272023702e-09<br>6.28036983474e-16 + 2.79837285903 J, 2.31471998545e-09<br>6.28036983474e-16 + -2.79837285903 J, 2.31471998545e-09<br>-5.4953236053e-16 + 2.78026777804 J, 6.1734628482e-10<br>-5.4953236053e-16 + -2.78026777804 J, 6.1734628482e-10<br>7.85046229433e-17 + 2.76882293959 J, 8.38075971231e-10<br>7.85046229433e-17 + -2.76882293959 J, 8.38075971231e-10<br><br>and here is the results from polyeig in Matlab.<br>tol =<br><br> &nb
 
 sp;1.0000e-008<br><br><dl class="buddy"><dd>ee =<br><br>        0 + 2.4425i<br>        0 - 2.4425i<br>        0 + 2.4218i<br>        0 - 2.4218i<br>        0 + 2.3876i<br>        0 - 2.3876i<br>        0 + 2.3405i<br>        0 - 2.3405i<br>        0 + 2.2815i<br>        0 - 2.2815i<br>        0 + 2.2117i<br>        0 - 2.2117i<br>        0&nb
 
 sp;+ 2.1324i<br>        0 - 2.1324i<br>        0 + 2.0456i<br>        0 - 2.0456i<br>        0 + 1.4262i<br>        0 - 1.4262i<br>        0 + 1.4612i<br>        0 - 1.4612i<br>        0 + 1.5164i<br>        0 - 1.5164i<br>        0 + 1.5881i<br>        0 - 1.5881i<br>        0 + 1.9533i<br>        0 - 1.9533i<br>   &
 
 nbsp;    0 + 1.6717i<br>        0 - 1.6717i<br>        0 + 1.7631i<br>        0 - 1.7631i<br>        0 + 1.8581i<br>        0 - 1.8581i  <br></dd></dl><br>it
 is different a lot especially for imagine part even if both converged 
to 1e-8. I tried both linear solver and qarnoldi method.   could any one
 tell me what is that? is it error?<br><div><br><p>--------------------------------------- </p>
<p>Wei zhang<br></p>PhD<br></div></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span> 




</div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>