Re: [deal.II] crash with step-55 with PETSc both in serial and parallel

2017-08-31 Thread Nur Fadel
Sorry, for the late answer, 
I haven't setup Hypre, this is the problem!
Thanks!

Nur

Il giorno domenica 13 agosto 2017 18:25:14 UTC+2, Timo Heister ha scritto:
>
> 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 <nura...@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 
> >  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_=DwIBaQ=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw=naKdiC5MpUUon2ztgdpg6EXMaY_Y4siuEIMCo2REU-U=bkBgAEUtvkgVku1tMVte7BoHucoYxCvhYeXIy_OHMzU=
>  
> > For mailing list/forum options, see 
> > 
> https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_forum_dealii-3Fhl-3Den=DwIBaQ=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw=naKdiC5MpUUon2ztgdpg6EXMaY_Y4siuEIMCo2REU-U=GdBkwqIdCUvJTldnuwEtrzs_3XawG9py6IWg74MwFdI=
>  
> > --- 
> > 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+un...@googlegroups.com . 
> > For more options, visit 
> https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_optout=DwIBaQ=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw=naKdiC5MpUUon2ztgdpg6EXMaY_Y4siuEIMCo2REU-U=iXMtDNGI9SX47AwdcwsJ6UoDB5XU5Ea9cnzidyZIZaE=
>  
> . 
>
>
>
> -- 
> 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.


[deal.II] crash with step-55 with PETSc both in serial and parallel

2017-08-11 Thread Nur Fadel
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 
 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 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.