<br><div class="gmail_quote">Hello,<br><br>I input a hessian matrix of size (5148 X 5148 size) into SLEPc to find it's eigenvalues (using EPS) as follows:<br><br><b>call MatCreate(PETSC_COMM_WORLD,A,ierr)<br>      call MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,3*natoms,3*natoms,ierr)<br>

      call MatSetFromOptions(A,ierr)<br>      call MatSetUp(A,ierr)<br><br>      do n=1,number_of_elements<br>              call MatSetValues(A,3*natoms,indices(n,1),3*natoms,indices(n,2),elements(n),INSERT_VALUES,ierr)<br>

      end do<br><br>      call MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY,ierr)<br>      call MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY,ierr)</b><br><br>But somehow, the eigenvalues do not seem to converge as I find the following output:<br>

<b><br> Number of iterations of the method:   1<br> Solution method: krylovschur<br> Number of requested eigenvalues: 1<br> Number of iterations of the method:   1<br> Solution method: krylovschur<br> Number of requested eigenvalues: 1<br>

 Stopping condition: tol=1.0000E-08, maxit= 648<br> Number of converged eigenpairs: 0<br></b><br>Can someone please tell me what's the problem, because this matrix for sure has eigenvalues as I have obtained before without using SLEPc?<br>

<br>Thanks & Regards,<br>Paul<br><br><br><br>
</div><br>