|
OR/MS Today - October 2007 GUPOR: Great Unsolved Problems in O.R. (Editor's note: This is another in a series of essays on the topic, "Great Unsolved Problems in O.R.," based on a special track of sessions organized by Saul I. Gass and Arjang A. Assad and presented at the 2006 INFORMS Fall Conference in Pittsburgh.) Computational Global Optimization: State-of-the-Art and Perspectives By János D. Pintér Nonlinear phenomena and processes are ubiquitous in both natural and man-made systems. Consequently, nonlinear descriptive and optimization models are of relevance across a vast range of studies in the natural sciences, engineering, economics and finances. Decision support models (strategic management, tactical optimization and real-time control) that incorporate an underlying nonlinear system description frequently possess multiple local and global optima. The objective of global optimization (GO) is to find the "absolutely best" solution of nonlinear optimization models under such circumstances. The GO paradigm subsumes and significantly extends the traditional (local) optimization framework in which an essentially convex model structure and/or the knowledge of a "sufficiently good" initial solution estimate is postulated. As is well known, local optimization strategies could (and often will) fail to find the global solution in multimodal problems. The continuous global optimization model can be stated as: min f(x) x∈D:={x: 1≤x≤u, g(x)≤0}: here x∈Rn the decision vector to optimize; l, u are finite component-wise vector bounds of x; f(x) is the continuous scalar objective function, and g(x) is the m-vector of continuous constraint functions. This model formally covers not only linear and convex nonlinear programming but as a highly special case also the entire class of combinatorial optimization models. This fact immediately hints at the tremendous potential difficulty of the GO model class. The theoretical and numerical complexity of GO models could increase at a massively exponential rate as the number of decision variables and/or of constraints increases. Due to its theoretical challenges and practical relevance, GO has become an important area of research. As of today, well over a hundred (English or other) textbooks and an increasing number of Web sites are devoted partly or completely to the subject. The most important GO model types and solution approaches (mostly exact, but also several prominent heuristic methods) are discussed in detail by the "Handbook of Global Optimization" volumes, edited by Horst and Pardalos (1995) and by Pardalos and Romeijn (2002). We also refer to the topical Webs site of Neumaier (2006), with numerous links to other useful information sources. The early development in GO (approximately until the late 1980s) was primarily devoted to model classifications and structural studies, and to the analysis of algorithmic frameworks to handle such problems with proven deterministic or stochastic global convergence properties. Several of the most important, broadly applicable GO model classes are: concave minimization (over a convex or nonconvex set), general complementarity problems, general quadratic programming, DC programming (all model functions can be represented as the difference of two convex functions), Lipschitz optimization (all functions are Lipschitz-continuous) or the most general case defined by merely continuous functions. These model classes are not disjoint; in fact, their hierarchy can be easily established. With respect to GO strategies, we mention here only a few prominent approaches. These include, in many variants, deterministically convergent branch-and-bound methods and stochastically convergent sequential search methods. The various heuristic approaches include adaptive neighborhood searches (simulated annealing and others), population-based (evolutionary and other) search procedures, as well as methods specifically targeted to handle GO models with computationally expensive functions. The key theoretical developments have been followed more recently by GO software implementations. As of today, professionally developed and supported software products are available for compiler platforms (C and FORTRAN), spreadsheets, optimization modeling languages (to our knowledge, currently for AIMMS, GAMS, LINGO and MPL), and for integrated scientific-technical computing systems (such as Maple, Mathematica and MATLAB). These model development and corresponding solver environments are aimed at meeting the needs of arguably different types of users. Key categories include educational users (instructors and students); research scientists, engineers, consultants and other practitioners (possibly, but not necessarily equipped with an in-depth optimization related background); optimization experts; software application developers; and other "power users." (The user categories listed are not necessarily disjoint.) The pros and cons of the individual software products (in terms of their hardware and software demands, ease of usage, model prototyping options, detailed code development and maintenance features, optimization model verification and processing tools, availability of solver options and other auxiliary tools, program execution speed, overall level of system integration, quality of related documentation and support, customization options and communication with end users) make the corresponding modeling and solver approaches more or less attractive for various user groups. In addition to professional software developments, non-commercial software projects are also devoted to solving various GO problem types. Both Mittelmann (2006) and Neumaier (2006) provide more extensive information on non-commercial systems. In recent decades, global optimization gradually has become an established discipline that is now taught worldwide at leading academic institutions. GO software is also increasingly applied in various contexts, including education, research, industrial and consulting practice. The currently available professional software implementations are routinely used to handle models with tens, hundreds and sometimes even thousands of variables and constraints. We recall the potential numerical difficulty of GO model instances; if one is interested in a guaranteed high-quality solution, then the necessary runtime could become minutes, hours, days or more, even on today's high-performance computers. Although one can expect further significant speed-up due to both further algorithm and software development and to progress in hardware/software technology, the theoretically exponential "curse of dimensionality" associated with the subject of GO will always remain. In the most general terms, global optimization technology is well-suited to analyze and solve complex nonlinear models in applied mathematics and the natural sciences; advanced engineering (acoustic, aerospace, chemical, control, electrical, environmental, process, telecommunications); biotechnology, medical and pharmaceutical studies; econometrics and financial modeling; and other areas. For GO application examples, case studies and further perspectives, see e.g. Corliss and Kearfott (1999); Edgar, Himmelblau and Lasdon (2001); Floudas et al. (1999); Gao, Ogden and Stavroulakis (2001); Grossmann (1996); Neumaier (2006); Nowak (2005); Papalambros and Wilde (2000); Pardalos, Shalloway and Xue (1996); Pintér (1996, 2002, 2006a and b); Schittkowski (2002); Tawarmalani and Sahinidis (2002); Zabinsky (2003).
OR/MS Today copyright © 2007 by the Institute for Operations Research and the Management Sciences. All rights reserved. Lionheart Publishing, Inc. 506 Roswell Rd., Suite 220, Marietta, GA 30060 USA Phone: 770-431-0867 | Fax: 770-432-6969 E-mail: lpi@lionhrtpub.com URL: http://www.lionhrtpub.com Web Site © Copyright 2007 by Lionheart Publishing, Inc. All rights reserved. |