Re: [petsc-dev] Fwd: [petsc-users] Local size 7 not compatible with block size 13

2014-01-15 Thread Barry Smith
On Jan 15, 2014, at 4:44 PM, Matthew Knepley wrote: > On Wed, Jan 15, 2014 at 4:36 PM, Barry Smith wrote: > >This problem comes up for MPIAIJ matrices when the “off-diagona” matrix is > compressed to only the columns that have non zeros in them. When a > MatDuplicate() is done on this i

Re: [petsc-dev] Fwd: [petsc-users] Local size 7 not compatible with block size 13

2014-01-15 Thread Matthew Knepley
On Wed, Jan 15, 2014 at 4:36 PM, Barry Smith wrote: > >This problem comes up for MPIAIJ matrices when the “off-diagona” matrix > is compressed to only the columns that have non zeros in them. When a > MatDuplicate() is done on this it detects the inconsistency that the number > of local colu

[petsc-dev] Fwd: [petsc-users] Local size 7 not compatible with block size 13

2014-01-15 Thread Barry Smith
This problem comes up for MPIAIJ matrices when the “off-diagona” matrix is compressed to only the columns that have non zeros in them. When a MatDuplicate() is done on this it detects the inconsistency that the number of local columns may no longer be divisible by the block size. We could h

Re: [petsc-dev] Fwd: [petsc-users] MatAssemblyEnd hangs during a parallel calculation with PETSc>3.3

2014-01-15 Thread Barry Smith
On Jan 15, 2014, at 11:20 AM, Matthew Knepley wrote: > This is a bug in the setup phase for Mat. I did not put in the logic, so I > might be missing the point. Here is > my outline: > > 1) When MatSetPreallocation_*() is called, we set the > NEW_NONZERO_ALLOCATION_ERR option, but ONLY >

Re: [petsc-dev] Include what you use

2014-01-15 Thread Karl Rupp
Hi Aron, > Stupid GMail. Link: http://sylvestre.ledru.info/blog/2013/09/08/include-what-you-use-yet Cool, thanks for the link! Best regards, Karli

[petsc-dev] Include what you use

2014-01-15 Thread Aron Ahmadia
This on your radar?

Re: [petsc-dev] Include what you use

2014-01-15 Thread Aron Ahmadia
Stupid GMail. Link: http://sylvestre.ledru.info/blog/2013/09/08/include-what-you-use-yet A On Wed, Jan 15, 2014 at 1:14 PM, Aron Ahmadia wrote: > This on your radar? >

[petsc-dev] Fwd: [petsc-users] MatAssemblyEnd hangs during a parallel calculation with PETSc>3.3

2014-01-15 Thread Matthew Knepley
This is a bug in the setup phase for Mat. I did not put in the logic, so I might be missing the point. Here is my outline: 1) When MatSetPreallocation_*() is called, we set the NEW_NONZERO_ALLOCATION_ERR option, but ONLY if the user passed a nonzero nz or nnz. 2) During assembly, we ch

Re: [petsc-dev] [GPU] Crash on ex19 with mpirun -np 2 (optimized build)

2014-01-15 Thread Karl Rupp
Hi Pierre, > I tried to rebuild the optimized PETSc library by changing several options and ran: mpirun -np 2 ./ex19 -cuda_show_devices -dm_mat_type aijcusp -dm_vec_type cusp -ksp_type fgmres -ksp_view -log_summary -pc_type none -snes_monitor_short -snes_rtol 1.e-5 Options used: --with-pthread

Re: [petsc-dev] Comp-comm overlap with DMPlex

2014-01-15 Thread Matthew Knepley
On Wed, Jan 15, 2014 at 7:16 AM, Michael Lange wrote: > Hi, > > I am trying to implement a computation-communication overlap with DMPlex, > where local values that do not need to be sent are updated while the DM/SF > is broadcasting the ghost values (DMGlobalToLocal() or PetscSFBcast(), I've > tr

[petsc-dev] Comp-comm overlap with DMPlex

2014-01-15 Thread Michael Lange
Hi, I am trying to implement a computation-communication overlap with DMPlex, where local values that do not need to be sent are updated while the DM/SF is broadcasting the ghost values (DMGlobalToLocal() or PetscSFBcast(), I've tried both). The problem is that this only works if I force the