[petsc-dev] broken TS example related to Jacobian and coloring etc.

Jed Brown jedbrown at mcs.anl.gov
Wed Aug 15 08:31:12 CDT 2012


On Wed, Aug 15, 2012 at 6:41 AM, Barry Smith <bsmith at mcs.anl.gov> wrote:

> I think it uses null in order to delay computing the matrix coloring etc
> from the time of SNESSetJacobian() to a later time.   I guess what could
> happen is it creates a n "empty" fdcoloring beast and passes that into
> SNESSetJacobian() and then fills it up later when it has the information it
> needs.


Right, I think I would prefer that. A side-effect is that the caller
setting coloring is now responsible for ownership of the context. but it is
for everything else anyway.

Lisandro asked for destructors for those contexts, but I'm not sure we need
a separate one for every context slot. Would one per SNES be enough?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120815/d98aae61/attachment.html>


More information about the petsc-dev mailing list