Picard Solver

[Matthew Knepley]

On Mon, Mar 2, 2009 at 11:42 PM, Dave Lee <dave at vpac.org> wrote:

> Thats ok, thanks Matt. Just so i'm clear - if i call SNESCreate_Picard(),
> that will set all the other SNES functions to their _Picard equivalents,
> then i can just use the regular SNES function calls (without having to set
> the build F or build J functions) to solve my problem. This is right yeah?


Just do everything the same way you always use SNES, but pass PETSC_NULL for
the Jacobian function, and give
-snes_type picard.

  Matt


>
> Dave.
> - Show quoted text -
>
>
> ----- Original Message -----
> From: "Matthew Knepley" <knepley at gmail.com>
> To: "PETSc users list" <petsc-users at mcs.anl.gov>
> Sent: Tuesday, March 3, 2009 4:31:30 PM GMT +10:00 Canberra / Melbourne /
> Sydney
> Subject: Re: Picard Solver
>
> I am responsible for the Picard implementation, and also for the lack of
> documentation. It has only quadratic line search now. I can give an
> example of custom line search if yo want.
>
> Matt
>
>
> On Mon, Mar 2, 2009 at 8:43 PM, Dave Lee < dave at vpac.org > wrote:
>
>
> Does PETSc v. 3.0.0 include a fully implemented Picard non-linear solver?
> I'm just wondering because there seems to be some functionality in the
> src/snes/impls/picard directory, but there's very little documentation on
> any of this in the online index pages, and nothing in the version 3.0.0
> manual...
>
> Cheers, Dave.
>
>
>
> --
> 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
>



-- 
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/20090302/cf575d59/attachment.htm>