[petsc-users] Inertia computation fails

Jose E. Roman jroman at dsic.upv.es
Thu Mar 5 07:58:51 CST 2020


Are you using the EPSKrylovSchurSetDetectZeros() option? As in this example https://slepc.upv.es/documentation/current/src/eps/examples/tutorials/ex25.c.html
If so, then the explanation is probably that one of the endpoints of your interval coincides with an eigenvalue. Try with a slightly different interval.

Jose


> El 5 mar 2020, a las 14:51, Satish Balay via petsc-users <petsc-users at mcs.anl.gov> escribió:
> 
> The address here should be petsc-users - not petsc-users-bounces
> 
> balay at sb /home/balay/git-repo/slepc (master=)
> $ git grep 'Inertia computation fails'
> src/eps/impls/krylov/krylovschur/ks-slice.c:          if (zeros) SETERRQ1(((PetscObject)eps)->comm,PETSC_ERR_CONV_FAILED,"Inertia computation fails in %g",sr->int1);
> src/eps/impls/krylov/krylovschur/ks-slice.c:      if (zeros) SETERRQ1(((PetscObject)eps)->comm,PETSC_ERR_CONV_FAILED,"Inertia computation fails in %g",newShift);
> src/pep/impls/krylov/stoar/qslice.c:        if (!nconv) SETERRQ1(((PetscObject)pep)->comm,PETSC_ERR_CONV_FAILED,"Inertia computation fails in %g",nzshift);
> src/pep/impls/krylov/stoar/qslice.c:      if (zeros) SETERRQ1(((PetscObject)pep)->comm,PETSC_ERR_CONV_FAILED,"Inertia computation fails in %g",newShift);
> 
> So the message is likely coming from slepc
> 
> Satish
> 
> ---------- Forwarded message ----------
> Date: Thu, 05 Mar 2020 14:43:13 +0100
> From: Perceval Desforges <perceval.desforges at polytechnique.edu>
> To: petsc-users <petsc-users-bounces at mcs.anl.gov>
> Subject: Inertia computation fails
> 
> Dear PETSc developpers, I am using SLEPC and MUMPS to calculate the eigenvalues of a real
> symmetric matrix in an interval. I have come upon a crash and I was
> unable to find any documentation on the error I got. The error is:  Inertia computation fails in 2.19893 Is this a slepc, a Petsc or a mumps problem? Thanks, Perceval



More information about the petsc-users mailing list