[petsc-users] Solving with SNES

[Matthew Knepley]

On Sat, Nov 12, 2011 at 4:20 PM, behzad baghapour <
behzad.baghapour at gmail.com> wrote:

> The state of my problem is set and saved in elements and faces. So I
> should pass and update them for residual and Jacobian calculations in each
> Newton step. It is better for me to use them since the element and face
> connectivities required in residual and jacobian matrix are saved there.
> How would be my best choice then???
>

I cannot understand what information you are talking about, so I can't make
any recommendation. If you can
clearly state, preferably with equations, what you are storing, we might be
able to help.

   Matt


> On Sat, Nov 12, 2011 at 6:04 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
>
>> On Sat, Nov 12, 2011 at 02:08, behzad baghapour <
>> behzad.baghapour at gmail.com> wrote:
>>
>>> I should pass a user-defined context, containing element and face local
>>> data obtained from PDE solution, to RHS and Jacobian evaluation of SNES
>>> process. Here I need to update my context each Newton Iteration. I may pass
>>> my context into a monitor function set by SNESMonitorSet. In the monitor
>>> function, I cast into my context and update my field (elements and faces)
>>> for next evaluations in RHS and Jacobian routines in Newton process.
>>
>>
>> Are you updating external forcing information? The state vector is passed
>> in to the residual and Jacobian routines, you should just use it.
>>
> --
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111112/a073d605/attachment.htm>