[petsc-dev] Creation of snes->vec_func

Jed Brown jedbrown at mcs.anl.gov
Thu Dec 22 10:20:47 CST 2011


On Thu, Dec 22, 2011 at 10:18, Barry Smith <bsmith at mcs.anl.gov> wrote:

> In some cases after certain rows/columns have been removed due to active
> constraints certain OTHER rows end up being identically zero. This is some
> hacky code to remove those rows (since they screw up the linear solve).
> Just leave the code.


No, it is *never* executed due to the if (0) statement! In addition to that
code never being executed, the call to MatFindNonzeroRows() with unused
return argument is an *unconditional* memory leak.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111222/da6e3c8e/attachment.html>


More information about the petsc-dev mailing list