<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Am 21.02.2012 13:52, schrieb Jed Brown:
<blockquote
cite="mid:CAM9tzSn=G4R-+h+fETTuEtggNoLh7U11VxVSdD5q3sf7A0if+Q@mail.gmail.com"
type="cite">
<div class="gmail_quote">On Tue, Feb 21, 2012 at 01:25, Thomas
Witkowski <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:thomas.witkowski@tu-dresden.de">thomas.witkowski@tu-dresden.de</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
Is it possible to have the same matrix of type MatMPIAIJ twice
in a MatNest?<br>
</blockquote>
<div><br>
</div>
<div>That should be fine, let us know if something doesn't work.</div>
</div>
</blockquote>
Seems not to work correctly. I made some small changes in
src/ksp/ksp/examples/tests/ex22.c:<br>
<br>
21c21<br>
< np = 2;<br>
---<br>
> np = 3;<br>
62c62<br>
< tmp[0][0] = A11;<br>
---<br>
> tmp[0][0] = A12;<br>
<br>
When running with 2 threads, it ends with a segmentation violation.
The stack is as follows:<br>
<br>
[0]PETSC ERROR: [0] MatStashScatterEnd_Private line 109
src/mat/utils/matstash.c<br>
[0]PETSC ERROR: [0] MatAssemblyEnd_MPIAIJ line 641
src/mat/impls/aij/mpi/mpiaij.c<br>
[0]PETSC ERROR: [0] MatAssemblyEnd line 4934
src/mat/interface/matrix.c<br>
[0]PETSC ERROR: [0] MatAssemblyEnd_Nest line 228
src/mat/impls/nest/matnest.c<br>
[0]PETSC ERROR: [0] MatAssemblyEnd line 4934
src/mat/interface/matrix.c<br>
[0]PETSC ERROR: [0] test_solve line 17
src/ksp/ksp/examples/tests/ex22.c<br>
<br>
Can you also reproduce this problem?<br>
<br>
Thomas<br>
</body>
</html>