<div dir="ltr"><div>Hello,</div><div><br></div><div>I am trying to solve Schroginger equation with open boundary conditions. The problem boils down to a complex polynomial eigenvalue equation of order two with A0,A1,A2 matrices of the size 300 x 300.</div><div>Here, A0 and A2 are tridiagonal matrices while A1 has only single element at with imaginary value.</div><div>The complex eigenvalues are such that real part is around 1 while imaginary part is in the range of 1e-12 to 1e-15.</div><div> </div><div> PEP solver gives out 600 complex eigenvalues having 300 complex conjugate pairs.</div><div>However problem is that the real part of two states is same while complex part of the eigenvalues differs by a order of one.( i mean ~e-14 and ~e-15)</div><div><br></div><div>The solver gives eigenvalues of the pair the same when the imaginary values are more than 1e-10.</div><div><br></div><div>I tried with -pep_tol 1e-18. -pep_max_it 10000. However, answer does not change for eigenvalues with imaginary part of the order of ~1e-15.</div><div><br></div><div>other options -pep_smallest imaginary and changing the solver types also does not seem to work.</div><div><br></div><div>Is there some way one can improve answers I mean more precise value of very small imaginary part of eigenvalue? </div><div> </div><div>We have tried  -st_type sinvert -st_transform  as well as -pep_target and -rg_type ellipse(region filtering) but there still the problem persists.</div><div> Kindly suggest any way to improve precision for imaginary values below (1e-10)   </div><div><br></div><div><div class="gmail_signature"><div dir="ltr">Regards,<div>Ajay Manwani<br><div><br></div></div></div></div></div>
</div>