21.06.2015 Views

Nonlinear Equations - UFRJ

Nonlinear Equations - UFRJ

Nonlinear Equations - UFRJ

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

[SEC. 10.1: HOMOTOPY ALGORITHM 137<br />

homotopy steps. They pointed out that this solves Smale’s 17 th problem<br />

for the ‘case’ max d i ≤ n 1<br />

1+ɛ while the ‘case’ max d i ≥ n 1+ɛ<br />

follows from resultant based algorithms such as [67]. When<br />

Smale’s 17 th problem is still open.<br />

n 1<br />

1+ɛ ≤ max di ≤ n 1+ɛ ,<br />

Another recent advance are ‘condition-length’ based algorithms.<br />

While previous algorithm have a complexity bound in terms of the<br />

line integral of µ(F t , z t ) 2 in P(H d ), condition-length algorithms (suggested<br />

in [14,69] and developed in [7,31] have a complexity bound in<br />

terms of a geometric invariant, the condition length. This allows to<br />

reduce Smale’s 17 th problem (Open Problem 1.11) to a ‘variational’<br />

problem.<br />

In the rest of this chapter, I will give a simplified version of the<br />

algorithm in [31], together with its complexity analysis. Then, I will<br />

discuss how to use this algorithm to obtain results analogous to those<br />

of [13] and [25]. In the last section, I will review some recent results<br />

on the geometry of the condition metric.<br />

10.1 Homotopy algorithm<br />

Let d = (d 1 , . . . , d n ) be fixed, and set D = max d i . Recall that<br />

H d is the space of homogeneous polynomial systems in n variables of<br />

degree d 1 , . . . , d n . We want to find solutions z ∈ P n , and those will be<br />

represented by elements of C n+1 \{0}. We keep the convention of the<br />

previous chapter, where we set Z for a representative of z. However,<br />

we will prefer representatives with norm one whenever possible.<br />

We will consider an affine path in H d given by<br />

F t = (1 − t)F 0 + tF 1<br />

where F 0 and F 1 are scaled such that<br />

‖F 0 ‖ = 1 F 0 ⊥ F 1 − F 0 (10.1)<br />

with an extra bound,<br />

‖F 1 − F 0 ‖ ≤ 1. (10.2)

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

Saved successfully!

Ooh no, something went wrong!