<div dir="ltr">I noticed that the current next is broken in MatZeroRows<br><br>The following patch solves the problem<br><br>[szampini@node362 petsc]$ git diff src/mat/impls/aij/mpi/mpiaij.c<br>diff --git a/src/mat/impls/aij/mpi/mpiaij.c b/src/mat/impls/aij/mpi/mpiaij.c<br>
index 820d52e..41f89e9 100644<br>--- a/src/mat/impls/aij/mpi/mpiaij.c<br>+++ b/src/mat/impls/aij/mpi/mpiaij.c<br>@@ -778,6 +778,7 @@ PetscErrorCode MatZeroRows_MPIAIJ(Mat A,PetscInt N,const PetscInt rows[],PetscSc<br>   ierr = PetscSFSetGraph(sf, n, N, NULL, PETSC_OWN_POINTER, rrows, PETSC_OWN_POINTER);CHKERRQ(ierr);<br>
   /* Collect flags for rows to be zeroed */<br>   ierr = PetscSFReduceBegin(sf, MPIU_INT, rows, lrows, MPI_LOR);CHKERRQ(ierr);<br>+  ierr = PetscSFReduceEnd(sf, MPIU_INT, rows, lrows, MPI_LOR);CHKERRQ(ierr);<br>   ierr = PetscSFDestroy(&sf);CHKERRQ(ierr);<br>
   /* Compress and put in row numbers */<br>   for (r = 0; r < n; ++r) if (lrows[r]) lrows[len++] = r;<br><br clear="all"><div><br>-- <br><div dir="ltr"><pre cols="72"><font><span style="font-family:garamond,serif"><span style="font-family:garamond,serif"></span>Ph. D. Stefano Zampini<br>
CINECA  
SuperComputing Applications and Innovations Department - SCAI
Via dei Tizii, 6 00185 Roma - ITALY
------------------------------------------------------------------------------------------------------------------------
Email: <a href="mailto:s.zampini@cineca.it" target="_blank"><span>s.zampini@cineca.it</span></a>
SkypeID: stefano.zampini
GoogleTalk: <a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>
Tel: <a value="+390644486720">+39 06.44486.707<br></a></span></font><a value="+390644486720"><font><span style="font-family:garamond,serif">------------------------------------------------------------------------------------------------------------------------
</span></font><br></a></pre></div>
</div></div>