Nur, do you have PETSc configured without hypre by any chance? step-55 runs correctly for me.
On Fri, Aug 11, 2017 at 7:24 AM, Nur Fadel <nurai...@gmail.com> wrote: > Hi all, > I have installed Deal II 8.5 and trunk versions on Daint (the CSCS machine). > When I run the example step-5 with PETSc I have the following error: > Running using PETSc. > Cycle 0: > Number of degrees of freedom: 659 (578+81) > +---------------------------------------------+------------+------------+ > | Total wallclock time elapsed since start | 0.354s | | > | | | | > | Section | no. calls | wall time | % of total | > +---------------------------------+-----------+------------+------------+ > | assembly | 1 | 0.14s | 40% | > | setup | 1 | 0.0873s | 25% | > | solve | 1 | 0.365s | 1e+02% | > +---------------------------------+-----------+------------+------------+ > ERROR: Uncaught exception in MPI_InitFinalize on proc 0. Skipping > MPI_Finalize() to avoid a deadlock. > ---------------------------------------------------- > Exception on processing: > -------------------------------------------------------- > An error occurred in line <66> of file > </project/csstaff/nfadel/dealii/source/lac/petsc_precondition.cc> in > function > void > dealii::PETScWrappers::PreconditionerBase::vmult(dealii::PETScWrappers::VectorBase&, > const dealii::PETScWrappers::VectorBase&) const > The violated condition was: > pc != nullptr > Additional information: (none) > -------------------------------------------------------- > Aborting! > ---------------------------------------------------- > [0]PETSC ERROR: #1 PetscCommDestroy() line 214 in > /project/csstaff/nfadel/petsc-vanilla/src/sys/objects/tagm.c > [0]PETSC ERROR: #2 PetscHeaderDestroy_Private() line 116 in > /project/csstaff/nfadel/petsc-vanilla/src/sys/objects/inherit.c > > on github I found the issue "step-55 PETSc fails #4150" already solved but > it seems not the same problem. > > I reproduce the crash both in serial and parallel, while trilinos works > fine. > > The error is caught when Deal II checks if the pointer to preconditioner > object's pc is a NULL. > AssertThrow (pc != nullptr, StandardExceptions::ExcInvalidState ()); > > It happens when the minres solver does: > precondition.vmult (v,*u[1]); > > Do you have any hint how to fix it? > > > Thank you in advance! > Nur > > > -- > The deal.II project is located at > https://urldefense.proofpoint.com/v2/url?u=http-3A__www.dealii.org_&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=naKdiC5MpUUon2ztgdpg6EXMaY_Y4siuEIMCo2REU-U&s=bkBgAEUtvkgVku1tMVte7BoHucoYxCvhYeXIy_OHMzU&e= > > For mailing list/forum options, see > https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_forum_dealii-3Fhl-3Den&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=naKdiC5MpUUon2ztgdpg6EXMaY_Y4siuEIMCo2REU-U&s=GdBkwqIdCUvJTldnuwEtrzs_3XawG9py6IWg74MwFdI&e= > > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to dealii+unsubscr...@googlegroups.com. > For more options, visit > https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_optout&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=naKdiC5MpUUon2ztgdpg6EXMaY_Y4siuEIMCo2REU-U&s=iXMtDNGI9SX47AwdcwsJ6UoDB5XU5Ea9cnzidyZIZaE&e= > . -- Timo Heister http://www.math.clemson.edu/~heister/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
