Nonlinear Equations - UFRJ
Nonlinear Equations - UFRJ
Nonlinear Equations - UFRJ
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
[SEC. 10.2: PROOF OF THEOREM 10.5 141<br />
1. If the algorithm terminates, then X is a (β, µ, a 0 ) certified approximate<br />
zero for F 1 .<br />
2. If the algorithm terminates, and z 0 denotes the zero of F 0 associated<br />
to X 0 , then z 1 is the zero of F 1 associated to X where<br />
f t (z t ) ≡ 0 is a continuous path.<br />
3. There is a constant C < 16.26 such that if the condition length<br />
L(f t , z t ; 0, 1) is finite, then the algorithm always terminates after<br />
at most<br />
1 + Cn 1/2 D 3/2 L(f t , z t ; 0, 1) (10.4)<br />
steps.<br />
The actual theorem in [31] is stronger, because the algorithm<br />
thereby allows for approximations instead of exact calculations. It<br />
is more general, as the path does not need to be linear. Also, it is<br />
worded in terms of the projective Newton operator N proj . This is<br />
why the constants are different. But the important feature of the<br />
theorem is an explicit step bound in terms of the condition length,<br />
and this is reproduced here.<br />
Remark 10.6. We can easily bound<br />
L(f t , z t ; 0, 1) ≤<br />
∫ 1<br />
0<br />
‖ḟt‖ ft µ(f t , z t ) 2 dt<br />
and recover the complexity analysis of previously known algorithms.<br />
Remark 10.7. The factor on √ n in the complexity bound comes from<br />
the approximation of µ by µ F . It can be removed at some cost. The<br />
price to pay is a more complicated subroutine for norm estimation,<br />
and a harder complexity analysis.<br />
10.2 Proof of Theorem 10.5<br />
Towards the proof of Theorem 10.5, we need five technical Lemmas.<br />
For the geometric insight, see figure 10.1.