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.

136 [CH. 10: HOMOTOPY<br />

form<br />

for<br />

x i+1 = N proj (F ti , x i ),<br />

F t = (1 − t)F 0 + tF 1 , 0 = t 0 ≤ t i ≤ t τ = 1.<br />

The major difficulty was finding an adequate starting pair (F 0 , x 0 ).<br />

Only the existence of such a pair was known, without any clue on how<br />

to find one in polynomial time.<br />

A minor difficulty was the choice of the t i . This can be done<br />

by trial and error. By doing so, there is no guarantee that one is<br />

approximating an actual continuous solution path F t (x t ) ≡ 0. This is<br />

trouble when attempting to find all the roots of a polynomial system,<br />

or when investigating the corresponding Galois group.<br />

In 2006, Carlos Beltrán and Luis Miguel Pardo demonstrated in<br />

his doctoral thesis [6, 11] the existence of a good ‘questor set’ from<br />

which an adequate random pair (F 0 , x 0 ) could be drawn with a good<br />

probability.<br />

A randomized algorithm is said to be of Vegas type if it returns an<br />

answer with probability 1 − ɛ for some ɛ, and the answer it returns is<br />

always correct. This is by opposition to Monte-Carlo type algorithms,<br />

that would return a correct answer with probability 1 − ɛ.<br />

Theorem 10.2 (Beltrán and Pardo). Let ɛ > 0. Then there is a<br />

Vegas type algorithm such that, given n, d = d 1 , . . . d n and a random<br />

F 1 ∈ (H d , dH d ), finds with probability 1 − ɛ an approximate zero X<br />

for F 1 , within expected time O(N 5 ɛ −2 ), where N = dim H d is the<br />

input size.<br />

This result and its proof was greatly improved in subsequent papers<br />

by Beltrán and Pardo such as [13]. The running time was reduced<br />

to<br />

E(τ) = C(max d i ) 3/2 nN<br />

homotopy steps.<br />

In another development, Peter Bürgisser and Felipe Cucker gave a<br />

deterministic algorithm for solving random systems within expected<br />

O(log log N)<br />

E(τ) = N

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

Saved successfully!

Ooh no, something went wrong!