Denis: Do you use shift-and-invert method for solving eigenvalue problem? If so, the linear problems would be extremely ill-conditioned, for which the direct solver, such LU or Cholesky are usually the only working option.
You may run your petsc/slepc code with option '-ksp_monitor' to observe convergence behavior. Hong After running in debug mode it seems that the GAMG solver indeed did not > converge, however throwing the error leads to SIGABRT (backtrace and frames > are below). > It is still very suspicious why would solving for (unchanged) mass matrix > wouldn't converge inside SLEPc's spectral transformation. > > p.s. valgrind takes enormous amount of time on this problem, > will try to leave it over the weekend. > > Denis. > > =============== > Program received signal SIGABRT, Aborted. > 0x00007fffea87fcc9 in __GI_raise (sig=sig@entry=6) > at ../nptl/sysdeps/unix/sysv/linux/raise.c:56 > 56 ../nptl/sysdeps/unix/sysv/linux/raise.c: No such file or directory. > (gdb) bt > #0 0x00007fffea87fcc9 in __GI_raise (sig=sig@entry=6) > at ../nptl/sysdeps/unix/sysv/linux/raise.c:56 > #1 0x00007fffea8830d8 in __GI_abort () at abort.c:89 > #2 0x00007fffeb790c91 in PetscTraceBackErrorHandler (comm=0x2a09bd0, > line=798, fun=0x7fffed0e24b9 <__func__.20043> "KSPSolve", > file=0x7fffed0e1620 > "/home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/ksp/ksp/interface/itfunc.c", > n=91, p=PETSC_ERROR_INITIAL, > mess=0x7fffffffac30 "KSPSolve has not converged", ctx=0x0) > at > /home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/sys/error/errtrace.c:243 > #3 0x00007fffeb78b8b9 in PetscError (comm=0x2a09bd0, line=798, > func=0x7fffed0e24b9 <__func__.20043> "KSPSolve", > file=0x7fffed0e1620 > "/home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/ksp/ksp/interface/itfunc.c", > n=91, p=PETSC_ERROR_INITIAL, > mess=0x7fffed0e1e7a "KSPSolve has not converged") > at > /home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/sys/error/err.c:377 > #4 0x00007fffec75e1e7 in KSPSolve (ksp=0x367227d0, b=0x35b285c0, > x=0x35d89250) > at > /home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/ksp/ksp/interface/itfunc.c:798 > #5 0x00007fffe32a8657 in STMatSolve (st=0x3672d820, b=0x35b285c0, > x=0x35d89250) > at > /home/davydden/.hashdist/tmp/slepc-22nb32nbgvhx/src/sys/classes/st/interface/stsles.c:166 > ---Type <return> to continue, or q <return> to quit---q > Quit > (gdb) f 5 > #5 0x00007fffe32a8657 in STMatSolve (st=0x3672d820, b=0x35b285c0, > x=0x35d89250) > at > /home/davydden/.hashdist/tmp/slepc-22nb32nbgvhx/src/sys/classes/st/interface/stsles.c:166 > 166 ierr = KSPSolve(st->ksp,b,x);CHKERRQ(ierr); > (gdb) f 4 > #4 0x00007fffec75e1e7 in KSPSolve (ksp=0x367227d0, b=0x35b285c0, > x=0x35d89250) > at > /home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/ksp/ksp/interface/itfunc.c:798 > 798 if (ksp->errorifnotconverged && ksp->reason < 0) > SETERRQ(comm,PETSC_ERR_NOT_CONVERGED,"KSPSolve has not converged"); > (gdb) f 3 > #3 0x00007fffeb78b8b9 in PetscError (comm=0x2a09bd0, line=798, > func=0x7fffed0e24b9 <__func__.20043> "KSPSolve", > file=0x7fffed0e1620 > "/home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/ksp/ksp/interface/itfunc.c", > n=91, p=PETSC_ERROR_INITIAL, > mess=0x7fffed0e1e7a "KSPSolve has not converged") > at > /home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/sys/error/err.c:377 > 377 if (!eh) ierr = > PetscTraceBackErrorHandler(comm,line,func,file,n,p,lbuf,0); > (gdb) f 2 > #2 0x00007fffeb790c91 in PetscTraceBackErrorHandler (comm=0x2a09bd0, > line=798, fun=0x7fffed0e24b9 <__func__.20043> "KSPSolve", > file=0x7fffed0e1620 > "/home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/ksp/ksp/interface/itfunc.c", > n=91, p=PETSC_ERROR_INITIAL, > mess=0x7fffffffac30 "KSPSolve has not converged", ctx=0x0) > at > /home/davydden/.hashdist/tmp/petsc-hujktg3j6hq7/src/sys/error/errtrace.c:243 > 243 abort(); > (gdb) f 1 > #1 0x00007fffea8830d8 in __GI_abort () at abort.c:89 > 89 abort.c: No such file or directory. > >
