[petsc-dev] MatZeroRows(), MatZeroRowsIS() have new arguments

Jed Brown jed at 59A2.org
Wed Oct 27 08:21:01 CDT 2010

On Wed, Oct 27, 2010 at 14:47, Barry Smith <bsmith at mcs.anl.gov> wrote:

> I've considered it but have not seen a specific need. Do you have examples
> where it really matters?

The natural one would be Dirichlet boundary conditions for
(non-Galerkin) geometric multigrid with either variable coefficients or
non-uniform mesh size.  I still prefer modifying the element matrices before
assembly (which also maintains symmetry), but if you want to use MatZeroRows
to impose the boundary condition, you would also need a non-constant

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20101027/27bd5698/attachment.html>

More information about the petsc-dev mailing list