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>
