<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi, all<br>

<br>

I am trying to solve a few eigenvalues in a interval of a sparse Hermitian matrix (F) with the dimension as 2362. I know there are only 20 eigenvalues in the interval [-80.0,-70.0] for the matrix F. I choose the EPSARNOLDI method. The code looks like as follows.

 But after several hours running, EPS can not get converged eigenvalues. Therefore, may I ask if I did some wrong in the code? Or should I add more solver parameters for my problem? Thank you in advance.<br>

<br>

#include <finclude/slepcepsdef.h><br>

      use slepceps<br>

#if defined(PETSC_USE_FORTRAN_DATATYPES)<br>

      type(Mat) F<br>

      type(EPS) solver<br>

#else<br>

      Mat F<br>

      EPS solver<br>

#endif<br>

      ... ...<br>

      ! creates eigensolver context<br>

      call EPSCreate(PETSC_COMM_WORLD, solver, ierr)<br>

      ! sets operators<br>

      call EPSSetOperators(solver, F, PETSC_NULL_OBJECT, ierr)<br>

      ! sets solver parameters<br>

      call EPSSetProblemType(solver, EPS_HEP, ierr)<br>

      call EPSSetWhichEigenpairs(solver, EPS_ALL, ierr)<br>

      call EPSSetInterval(solver, -80.0, 70.0, ierr)<br>

      call EPSSetType(solver, EPSARNOLDI, ierr)<br>

      call EPSSetTolerances(solver, 1.0D-9, 90000, ierr)<br>

      ! solve the eigensystem<br>

      call EPSSolve(solver, ierr)<br>

<br>

Cheers<br>

<br>

Gao<br>

</div>

</body>

</html>