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 147<br />

By (10.9) and taking µ ≥ √ 2, D ≥ 2, β(F ti+1 , X i ) ≤ (1 − ɛ 1 )α/2.<br />

Therefore,<br />

d proj (N(F ti+1 , X i ), z ti+1 ) ≤ 2 1 − ɛ 1 1 − α<br />

1<br />

D 3/2 µ ψ(α) α2 r 1 (α)<br />

1 − (1 − ɛ 1 )α/2<br />

using (10.13).<br />

Putting all together,<br />

L(f t , z t ; t i , t i+1 ) ≥<br />

2<br />

D 3/2√ n ×<br />

×<br />

(ɛ 2−a 0) √ 1−ɛ 2 1<br />

(1+ɛ 1)<br />

− a 0 r 0 (a 0 ) − (1 − ɛ 1 )<br />

1−α<br />

ψ(α)(1−(1−ɛ 1)α/2) α2 r 1 (α)<br />

1 + ɛ 1 + π(1 − ɛ 1 )αr 0 (α)/ √ D<br />

The final bound was obtained numerically, assuming D ≥ 2. We<br />

check computationally that<br />

2<br />

(ɛ 2−a 0) √ 1−ɛ 2 1<br />

(1+ɛ 1)<br />

− a 0 r 0 (a 0 ) − (1 − ɛ 1 )<br />

1−α<br />

ψ(α)(1−(1−ɛ 1)α/2) α2 r 1 (α)<br />

1 + ɛ 1 + π(1 − ɛ 1 )αr 0 (α)/ √ 2<br />

≥ C −1<br />

(10.18)<br />

Proof of Theorem 10.5. Suppose the algorithm terminates. We claim<br />

that for each t i , X i is a (β, µ, a 0 )-certified approximate zero of F ti ,<br />

and that its associated zero is z ti . This is true by hypothesis when<br />

i = 0. Therefore, assume this is true up to a certain i.<br />

Recall that β(F, X) scales as ‖X‖. In particular,<br />

β(F ti+1 , X i+1 ) = β(F t i+1<br />

, N(F ti+1 , X i ))<br />

‖N(F ti+1 , X i )‖<br />

By (10.9) again, β(F ti+1 , X i ) ≤ (1 − ɛ 1 )α/2.<br />

We apply (10.14) to obtain that<br />

D 3/2<br />

2 β(F s, X i+1 )µ(f s , [N(F s , X i )]) ≤ a 0 .<br />

≤ β(F t i+1<br />

, N(F ti+1 , X i ))<br />

.<br />

1 − β(F ti+1 , X i )<br />

From (10.11), X i is an approximate zero of the second kind for<br />

F s , s ∈ [t i , t i+1 ]. Since both α(F s , X i ) and β(F s , X i ) are bounded

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

Saved successfully!

Ooh no, something went wrong!