<div class="gmail_quote">On Fri, Sep 16, 2011 at 23:00, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div>It is the name for this generic concept:<div><br></div><div>  Solve F(u) = 0, by successive approximations u^{n+1} = F(u^n)</div></div></blockquote></div><br><div>You have defined F(u^n) = x_n - f(x_n). Why is this the One True Way to solve f(x) = 0?</div>

<div><br></div><div>When you were a child and learned fixed-point iteration for scalar problems, you learned to reformulate f(x) = 0 as x = F(x) so that the fixed point iteration was convergent. This step was crucial because for any given problem, there are lots of was to formulate the fixed point iteration and most of them don't work. A worked example for f(x) = x^3 - sin x is in these notes:</div>
<div><br></div><div><a href="http://pages.cs.wisc.edu/~amos/412/lecture-notes/lecture03.pdf">http://pages.cs.wisc.edu/~amos/412/lecture-notes/lecture03.pdf</a></div><div><br></div><div><br></div><div>I don't see why systems are any different. Why do you assume that when given f(x), the only mathematically pure formulation is as x = x - f(x)?</div>