[petsc-users] Changing the jacobian matrix during SNES iterations

Matthew Knepley knepley at gmail.com
Wed Jul 18 12:59:48 CDT 2018

On Wed, Jul 18, 2018 at 1:31 PM David Knezevic <david.knezevic at akselos.com>

> I'm using SNES for a finite element contact solve, in which the sparsity
> pattern of the jacobian can change from one Newton iteration to the next
> (since the nodes on the contact surface move).
> In order to handle this I figured the best way would be to destroy the
> jacobian matrix and re-allocate it with a new sparsity pattern inside each
> call to FormJacobian, does that seem like a reasonable approach in this
> context?


> Also, I recall from an earlier discussion that this matrix re-allocation
> inside FormJacobian is supported by SNES, but I just wanted to confirm that?


> Also, I was wondering if there is any example where the matrix is
> re-allocated inside SNES iterations so that I can make sure that I do it
> correctly?

No, unfortunately. Contributions always welcome :)



> Best regards,
> David

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180718/f3d3de4e/attachment.html>

More information about the petsc-users mailing list