Stefano Zampini <[email protected]> writes: > I noticed that the current next is broken in MatZeroRows > > The following patch solves the problem
Thanks, applied in 'knepley/fix-mat-zero-rows'. > [szampini@node362 petsc]$ git diff src/mat/impls/aij/mpi/mpiaij.c > diff --git a/src/mat/impls/aij/mpi/mpiaij.c b/src/mat/impls/aij/mpi/mpiaij.c > index 820d52e..41f89e9 100644 > --- a/src/mat/impls/aij/mpi/mpiaij.c > +++ b/src/mat/impls/aij/mpi/mpiaij.c > @@ -778,6 +778,7 @@ PetscErrorCode MatZeroRows_MPIAIJ(Mat A,PetscInt > N,const PetscInt rows[],PetscSc > ierr = PetscSFSetGraph(sf, n, N, NULL, PETSC_OWN_POINTER, rrows, > PETSC_OWN_POINTER);CHKERRQ(ierr); > /* Collect flags for rows to be zeroed */ > ierr = PetscSFReduceBegin(sf, MPIU_INT, rows, lrows, > MPI_LOR);CHKERRQ(ierr); > + ierr = PetscSFReduceEnd(sf, MPIU_INT, rows, lrows, > MPI_LOR);CHKERRQ(ierr); FWIW, your mail client is messing up line wrapping so this change has to be applied manually rather than using "git apply" or "git am". It's easy enough to do by hand for a one-liner, but it's best not to send multi-line patches this way. "git format-patch" or "git send-email" are good ways to send patches by mail.
pgpj83iCORAUh.pgp
Description: PGP signature
