21.06.2015 Views

Nonlinear Equations - UFRJ

Nonlinear Equations - UFRJ

Nonlinear Equations - UFRJ

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!