On Mon, Mar 2, 2009 at 11:42 PM, Dave Lee <span dir="ltr"><<a href="mailto:dave@vpac.org">dave@vpac.org</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
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?</blockquote>
<div><br>Just do everything the same way you always use SNES, but pass PETSC_NULL for the Jacobian function, and give<br>-snes_type picard.<br><br> Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<font color="#888888"><br>
Dave.<br>
</font><div><div><span id="q_11fcaee4b14b5e49_2" class="h4">- Show quoted text -</span></div><div class="h5"><br>
<br>
----- Original Message -----<br>
From: "Matthew Knepley" <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>><br>
To: "PETSc users list" <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>><br>
Sent: Tuesday, March 3, 2009 4:31:30 PM GMT +10:00 Canberra / Melbourne / Sydney<br>
Subject: Re: Picard Solver<br>
<br>
I am responsible for the Picard implementation, and also for the lack of<br>
documentation. It has only quadratic line search now. I can give an<br>
example of custom line search if yo want.<br>
<br>
Matt<br>
<br>
<br>
On Mon, Mar 2, 2009 at 8:43 PM, Dave Lee < <a href="mailto:dave@vpac.org">dave@vpac.org</a> > wrote:<br>
<br>
<br>
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...<br>
<br>
Cheers, Dave.<br>
<br>
<br>
<br>
--<br>
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>