20.05.2014 Views

link to my thesis

link to my thesis

link to my thesis

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.

Summary<br />

The problem of computing the solutions of a system of multivariate polynomial equations<br />

can be approached by the Stetter-Möller matrix method which casts the problem<br />

in<strong>to</strong> a large eigenvalue problem. This Stetter-Möller matrix method forms the starting<br />

point for the development of computational procedures for the two main applications<br />

addressed in this <strong>thesis</strong>:<br />

• The global optimization of a multivariate polynomial, described in Part II of<br />

this <strong>thesis</strong>, and<br />

• the H 2 model-order reduction problem, described in Part III of this <strong>thesis</strong>.<br />

Part I of this <strong>thesis</strong> provides an introduction in the background of algebraic geometry<br />

and an overview of various methods <strong>to</strong> solve systems of polynomial equations.<br />

In Chapter 4 a global optimization method is worked out which computes the<br />

global minimum of a Minkowski dominated polynomial. The Stetter-Möller matrix<br />

method transforms the problem of finding solutions <strong>to</strong> the system of first-order conditions<br />

of this polynomial in<strong>to</strong> an eigenvalue problem. This method is described in<br />

[48], [50] and [81]. A drawback of this approach is that the matrices involved in this<br />

eigenvalue problem are usually very large.<br />

The research question which plays a central role in Part II of this <strong>thesis</strong> is formulated<br />

as follows: How <strong>to</strong> improve the efficiency of the Stetter-Möller matrix method,<br />

applied <strong>to</strong> the global optimization of a multivariate polynomial, by means of an nDsystems<br />

approach? The efficiency of this method is improved in this <strong>thesis</strong> by using a<br />

matrix-free implementation of the matrix-vec<strong>to</strong>r products and by using an iterative<br />

eigenvalue solver instead of a direct eigenvalue solver.<br />

The matrix-free implementation is achieved by the development and implementation<br />

of an nD-system of difference equations as described in Chapter 5. This yields a<br />

routine that computes the action of the involved large and sparse matrix on a given<br />

257

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

Saved successfully!

Ooh no, something went wrong!