[petsc-users] Add a matrix to MatNest twice

[Thomas Witkowski]

Am 21.02.2012 13:52, schrieb Jed Brown:
> On Tue, Feb 21, 2012 at 01:25, Thomas Witkowski 
> <thomas.witkowski at tu-dresden.de 
> <mailto:thomas.witkowski at tu-dresden.de>> wrote:
>
>     Is it possible to have the same matrix of type MatMPIAIJ twice in
>     a MatNest?
>
>
> That should be fine, let us know if something doesn't work.
Seems not to work correctly. I made some small changes in 
src/ksp/ksp/examples/tests/ex22.c:

21c21
<   np = 2;
---
 >   np = 3;
62c62
<   tmp[0][0] = A11;
---
 >   tmp[0][0] = A12;

When running with 2 threads, it ends with a segmentation violation. The 
stack is as follows:

[0]PETSC ERROR: [0] MatStashScatterEnd_Private line 109 
src/mat/utils/matstash.c
[0]PETSC ERROR: [0] MatAssemblyEnd_MPIAIJ line 641 
src/mat/impls/aij/mpi/mpiaij.c
[0]PETSC ERROR: [0] MatAssemblyEnd line 4934 src/mat/interface/matrix.c
[0]PETSC ERROR: [0] MatAssemblyEnd_Nest line 228 
src/mat/impls/nest/matnest.c
[0]PETSC ERROR: [0] MatAssemblyEnd line 4934 src/mat/interface/matrix.c
[0]PETSC ERROR: [0] test_solve line 17 src/ksp/ksp/examples/tests/ex22.c

Can you also reproduce this problem?

Thomas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120221/f4628770/attachment.htm>