[petsc-dev] Broken ZeroRows in next

Jed Brown jedbrown at mcs.anl.gov
Tue Nov 5 13:28:27 CST 2013


Stefano Zampini <s.zampini at cineca.it> 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 at 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.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20131105/f7906987/attachment.sig>


More information about the petsc-dev mailing list