File: SNESNRICHARDSON.html

package info (click to toggle)
petsc 3.4.2.dfsg1-8.1
  • links: PTS, VCS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 129,104 kB
  • ctags: 516,422
  • sloc: ansic: 395,939; cpp: 47,201; python: 34,788; makefile: 17,193; fortran: 16,251; f90: 1,592; objc: 954; sh: 822; xml: 621; java: 381; lisp: 293; csh: 241
file content (35 lines) | stat: -rw-r--r-- 2,757 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML3.2 EN">
<HTML>
<HEAD> <link rel="canonical" href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESNRICHARDSON.html" />
<META NAME="GENERATOR" CONTENT="DOCTEXT">
<TITLE>SNESNRICHARDSON</TITLE>
</HEAD>
<BODY BGCOLOR="FFFFFF">
   <div id="version" align=right><b>petsc-3.4.2 2013-07-02</b></div>
<A NAME="SNESNRICHARDSON"><H1>SNESNRICHARDSON</H1></A>
Richardson nonlinear solver that uses successive substitutions, also sometimes known as Picard iteration. 
<P>
<H3><FONT COLOR="#CC3333">Options Database</FONT></H3>
<TABLE border="0" cellpadding="0" cellspacing="0">
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-snes_linesearch_type &lt;l2,cp,basic&gt; Line search type.</B></TD><TD>- -   -snes_linesearch_damping&lt;1.0&gt; Damping for the line search.
</TD></TR>
<P>
Notes: If no inner nonlinear preconditioner is provided then solves F(x) - b = 0 using x^{n+1} = x^{n} - lambda
(F(x^n) - b) where lambda is obtained either <A HREF="../SNES/SNESLineSearchSetDamping.html#SNESLineSearchSetDamping">SNESLineSearchSetDamping</A>(), -snes_damping or a line search.  If
an inner nonlinear preconditioner is provided (either with -npc_snes_type or <A HREF="../SNES/SNESSetPC.html#SNESSetPC">SNESSetPC</A>()) then the inner
solver is called an initial solution x^n and the nonlinear Richardson uses x^{n+1} = x^{n} + lambda d^{n}
where d^{n} = \hat{x}^{n} - x^{n} where \hat{x}^{n} is the solution returned from the inner solver.
<P>
The update, especially without inner nonlinear preconditioner, may be ill-scaled.  If using the basic
linesearch, one may have to scale the update with -snes_linesearch_damping
<P>
This uses no derivative information thus will be much slower then Newton's method obtained with -snes_type ls
<P>
<H3><FONT COLOR="#CC3333">See Also</FONT></H3>
  <A HREF="../SNES/SNESCreate.html#SNESCreate">SNESCreate</A>(), <A HREF="../SNES/SNES.html#SNES">SNES</A>, <A HREF="../SNES/SNESSetType.html#SNESSetType">SNESSetType</A>(), <A HREF="../SNES/SNESNEWTONLS.html#SNESNEWTONLS">SNESNEWTONLS</A>, <A HREF="../SNES/SNESNEWTONTR.html#SNESNEWTONTR">SNESNEWTONTR</A>, <A HREF="../SNES/SNESNGMRES.html#SNESNGMRES">SNESNGMRES</A>, <A HREF="../SNES/SNESQN.html#SNESQN">SNESQN</A>, <A HREF="../SNES/SNESNCG.html#SNESNCG">SNESNCG</A>
<BR><P><B><P><B><FONT COLOR="#CC3333">Level:</FONT></B>beginner
<BR><FONT COLOR="#CC3333">Location:</FONT></B><A HREF="../../../src/snes/impls/richardson/snesrichardson.c.html#SNESNRICHARDSON">src/snes/impls/richardson/snesrichardson.c</A>
<BR><A HREF="./index.html">Index of all SNES routines</A>
<BR><A HREF="../../index.html">Table of Contents for all manual pages</A>
<BR><A HREF="../singleindex.html">Index of all manual pages</A>
</BODY></HTML>