# Multi-objective optimization

Plot of objectives when maximizing return and minimizing risk in financial portfolios (Pareto optimal points in red)

Multi-objective optimization (also known as multi-objective programming, vector optimization, multicriteria optimization, multiattribute optimization or Pareto optimization) is an area of multiple criteria decision making, that is concerned with mathematical optimization problems involving more than one objective function to be optimized simultaneously. Multi-objective optimization has been applied in many fields of science, including engineering, economics and logistics (see the section on applications for detailed examples) where optimal decisions need to be taken in the presence of trade-offs between two or more conflicting objectives. Minimizing weight while maximizing the strength of a particular component, and maximizing performance whilst minimizing fuel consumption and emission of pollutants of a vehicle are examples of multi-objective optimization problems involving two and three objectives, respectively. In practical problems, there can be more than three objectives.

For a nontrivial multi-objective optimization problem, there does not exist a single solution that simultaneously optimizes each objective. In that case, the objective functions are said to be conflicting, and there exists a (possibly infinite number of) Pareto optimal solutions. A solution is called nondominated, Pareto optimal, Pareto efficient or noninferior, if none of the objective functions can be improved in value without degrading some of the other objective values. Without additional subjective preference information, all Pareto optimal solutions are considered equally good (as vectors cannot be ordered completely). Researchers study multi-objective optimization problems from different viewpoints and, thus, there exist different solution philosophies and goals when setting and solving them. The goal may be to find a representative set of Pareto optimal solutions, and/or quantify the trade-offs in satisfying the different objectives, and/or finding a single solution that satisfies the subjective preferences of a human decision maker (DM).

## Introduction

A multi-objective optimization problem is an optimization problem that involves multiple objective functions.[1][2] In mathematical terms, a multi-objective optimization problem can be formulated as

\begin{align} \min &\left(f_1(x), f_2(x),\ldots, f_k(x) \right) \\ \text{s.t. } &x\in X, \end{align}

where the integer $k\geq 2$ is the number of objectives and the set $X$ is the feasible set of decision vectors. The feasible set is typically defined by some constraint functions. In addition, the vector-valued objective function is often defined as

$f:X\to\mathbb R^k, \ f(x)= (f_1(x),\ldots,f_k(x))^T$. If some objective function is to be maximized, it is equivalent to minimize its negative. The image of $X$ is denoted by $Y \in \mathbb R^k$

An element $x^*\in X$ is called a feasible solution or a feasible decision. A vector $z^* := f(x^*)\in \mathbb R^k$ for a feasible solution $x^*$ is called an objective vector or an outcome. In multi-objective optimization, there does not typically exist a feasible solution that minimizes all objective functions simultaneously. Therefore, attention is paid to Pareto optimal solutions; that is, solutions that cannot be improved in any of the objectives without degrading at least one of the other objectives. In mathematical terms, a feasible solution $x^1\in X$ is said to (Pareto) dominate another solution $x^2\in X$, if

1. $f_i(x^1)\leq f_i(x^2)$ for all indices $i \in \left\{ {1,2,\dots,k } \right\}$ and
2. $f_j(x^1) < f_j(x^2)$ for at least one index $j \in \left\{ {1,2,\dots,k } \right\}$.

A solution $x^1\in X$ (and the corresponding outcome $f(x^*)$) is called Pareto optimal, if there does not exist another solution that dominates it. The set of Pareto optimal outcomes is often called the Pareto front or Pareto boundary.

The Pareto front of a multi-objective optimization problem is bounded by a so-called nadir objective vector $z^{nad}$ and an ideal objective vector $z^{ideal}$, if these are finite. The nadir objective vector is defined as

$z^{nad}_i= \sup_{x\in X\text{ is Pareto optimal}} f_i(x) \text{ for all } i=1,\ldots,k$

and the ideal objective vector as

$z^{ideal}_i=\inf{f_i(x)}\text{ for all } i=1,\ldots,k.$

In other words, the components of a nadir and an ideal objective vector define upper and lower bounds for the objective function values of Pareto optimal solutions, respectively. In practice, the nadir objective vector can only be approximated as, typically, the whole Pareto optimal set is unknown.

## Examples of multi-objective optimization applications

### Economics

In economics, many problems involve multiple objectives along with constraints on what combinations of those objectives are attainable. For example, consumer's demand for various goods is determined by the process of maximization of the utilities derived from those goods, subject to a constraint based on how much income is available to spend on those goods and on the prices of those goods. This constraint allows more of one good to be purchased only at the sacrifice of consuming less of another good; therefore, the various objectives (more consumption of each good is preferred) are in conflict with each other. A common method for analyzing such a problem is to use a graph of indifference curves, representing preferences, and a budget constraint, representing the trade-offs that the consumer is faced with.

Another example involves the production possibilities frontier, which specifies what combinations of various types of goods can be produced by a society with certain amounts of various resources. The frontier specifies the trade-offs that the society is faced with — if the society is fully utilizing its resources, more of one good can be produced only at the expense of producing less of another good. A society must then use some process to choose among the possibilities on the frontier.

Macroeconomic policy-making is a context requiring multi-objective optimization. Typically a central bank must choose a stance for monetary policy that balances competing objectives — low inflation, low unemployment, low balance of trade deficit, etc. To do this, the central bank uses a model of the economy that quantitatively describes the various causal linkages in the economy; it simulates the model repeatedly under various possible stances of monetary policy, in order to obtain a menu of possible predicted outcomes for the various variables of interest. Then in principle it can use an aggregate objective function to rate the alternative sets of predicted outcomes, although in practice central banks use a non-quantitative, judgement-based, process for ranking the alternatives and making the policy choice.

### Finance

In finance, a common problem is to choose a portfolio when there are two conflicting objectives — the desire to have the expected value of portfolio returns be as high as possible, and the desire to have risk, measured by the standard deviation of portfolio returns, be as low as possible. This problem is often represented by a graph in which the efficient frontier shows the best combinations of risk and expected return that are available, and in which indifference curves show the investor's preferences for various risk-expected return combinations. The problem of optimizing a function of the expected value (first moment) and the standard deviation (square root of the second moment) of portfolio return is called a two-moment decision model.

### Optimal control

In engineering and economics, many problems involve multiple objectives which are not describable as the-more-the-better or the-less-the-better; instead, there is an ideal target value for each objective, and the desire is to get as close as possible to the desired value of each objective. For example, one might want to adjust a rocket's fuel usage and orientation so that it arrives both at a specified place and at a specified time; or one might want to conduct open market operations so that both the inflation rate and the unemployment rate are as close as possible to their desired values.

Often such problems are subject to linear equality constraints that prevent all objectives from being simultaneously perfectly met, especially when the number of controllable variables is less than the number of objectives and when the presence of random shocks generates uncertainty. Commonly a multi-objective quadratic objective function is used, with the cost associated with an objective rising quadratically with the distance of the objective from its ideal value. Since these problems typically involve adjusting the controlled variables at various points in time and/or evaluating the objectives at various points in time, intertemporal optimization techniques are employed.

### Optimal design

Product and process design can be largely improved using modern modeling, simulation and optimization techniques. The key question in optimal design is the measure of what is good or desirable about a design. Before looking for optimal designs it is important to identify characteristics which contribute the most to the overall value of the design. A good design typically involves multiple criteria/objectives such as capital cost/investment, operating cost, profit, quality and/or recovery of the product, efficiency, process safety, operation time etc. Therefore, in practical applications, the performance of process and product design is often measured with respect to multiple objectives. These objectives typically are conflicting, i.e., achieving the optimal value for one objective requires some compromise on one or more of other objectives.

For example, in paper industry when designing a paper mill, one can seek to decrease the amount of capital invested in a paper mill and enhance the quality of paper simultaneously. If the design of a paper mill is defined by large storage volumes and paper quality is defined by quality parameters, then the problem of optimal design of a paper mill can include objectives such as: i) minimization of expected variation of those quality parameter from their nominal values, ii) minimization of expected time of breaks and iii) minimization of investment cost of storage volumes. Here, maximum volume of towers are design variables. This example of optimal design of a paper mill is a simplification of the model used in.[3]

### Radio resource management

The purpose of radio resource management is to satisfy the data rates that are requested by the users of a cellular network.[4] The main resources are time intervals, frequency blocks, and transmit powers. Each user has its own objective function that, for example, can represent some combination of the data rate, latency, and energy efficiency. These objectives are conflicting since the frequency resources are very scarce, thus there is a need for tight spatial frequency reuse which causes immense inter-user interference if not properly controlled. Multi-user MIMO techniques are nowadays used to reduce the interference by adaptive precoding. The network operator would like to both bring great coverage and high data rates, thus the operator would like to find a Pareto optimal solution that balance the total network data throughput and the user fairness in an appropriate subjective manner.

Radio resource management is often solved by scalarization; that is, selection of a network utility function that tries to balance throughput and user fairness. The choice of utility function has a large impact on the computational complexity of the resulting single-objective optimization problem.[4] For example, the common utility of weighted sum rate gives an NP-hard problem with a complexity that scales exponentially with the number of users, while the weighted max-min fairness utility results in a quasi-convex optimization problem with only a polynomial scaling with the number of users.[5]

## Solving a multi-objective optimization problem

As there usually exist multiple Pareto optimal solutions for multi-objective optimization problems, what it means to solve such a problem is not as straightforward as it is for a conventional single-objective optimization problem. Therefore, different researchers have defined the term "solving a multi-objective optimization problem" in various ways. This section summarizes some of them and the contexts in which they are used. Many methods convert the original problem with multiple objectives into a single-objective optimization problem. This is called a scalarized problem. If scalarization is done carefully, Pareto optimality of the solutions obtained can be guaranteed.

Solving a multi-objective optimization problem is sometimes understood as approximating or computing all or a representative set of Pareto optimal solutions.[6][7]

When decision making is emphasized, the objective of solving a multi-objective optimization problem is referred to supporting a decision maker in finding the most preferred Pareto optimal solution according to his/her subjective preferences.[1][8] The underlying assumption is that one solution to the problem must be identified to be implemented in practice. Here, a human decision maker (DM) plays an important role. The DM is expected to be an expert in the problem domain.

The most preferred solution can be found using different philosophies. Multi-objective optimization methods can be divided into four classes. [2] In so-called no preference methods, no DM is expected to be available, but a neutral compromise solution is identified without preference information.[1] The other classes are so-called a priori, a posteriori and interactive methods and they all involve preference information from the DM in different ways.

In a priori methods, preference information is first asked from the DM and then a solution best satisfying these preferences is found. In a posteriori methods, a representative set of Pareto optimal solutions is first found and then the DM must choose one of them. In interactive methods, the decision maker is allowed to iteratively search for the most preferred solution. In each iteration of the interactive method, the DM is shown Pareto optimal solution(s) and describes how the solution(s) could be improved. The information given by the decision maker is then taken into account while generating new Pareto optimal solution(s) for the DM to study in the next iteration. In this way, the DM learns about the feasibility of his/her wishes and can concentrate on solutions that are interesting to him/her. The DM may stop the search whenever he/she wants to. More information and examples of different methods in the four classes are given in the following sections.

## Scalarizing multi-objective optimization problems

Scalarizing a multi-objective optimization problem means formulating a single-objective optimization problem such that optimal solutions to the single-objective optimization problem are Pareto optimal solutions to the multi-objective optimization problem.[2] In addition, it is often required that every Pareto optimal solution can be reached with some parameters of the scalarization.[2] With different parameters for the scalarization, different Pareto optimal solutions are produced. Very well-known examples are the so-called linear scalarization

$\min_{x\in X} \sum_{i=1}^k w_if_i(x),$

where the weights of the objectives $w_i>0$ are the parameters of the scalarization, and the $\epsilon$-constraint method (see, e.g.,[1])

$\begin{array}{ll} \min & f_j(x)\\ \text{s.t. }&x \in S\\ &f_i(x)\leq \epsilon_j \text{ for }i\in\{1,\ldots,k\}\setminus\{j\}, \end{array}$

where upper bounds $\epsilon_j$ are parameters as above and $f_j$ is the objective to be minimized.

Portfolio optimization is often conducted in terms of mean-variance analysis. In this context, the efficient set is a subset of the portfolios parametrized by the portfolio mean return $\mu_P$ in the problem of choosing portfolio shares so as to minimize the portfolio's variance of return $\sigma_P$ subject to a given value of $\mu_P$; see Mutual fund separation theorem for details. Alternatively, the efficient set can be specified by choosing the portfolio shares so as to maximize the function $\mu_P - b \sigma_P$; the set of efficient portfolios consists of the solutions as b ranges from zero to infinity.

## No-preference methods

Multi-objective optimization methods that do not require any preference information to be explicitly articulated by a decision maker can be classified as no-preference methods.[2] A well-known example is the method of global criterion,[9] in which a scalarized problem of the form

\begin{align} \min&\|f(x)-z^{ideal}\|\\ \text{s.t. }&x\in X \end{align}

is solved. In the above problem, $\|\cdot\|$ can be any $L_p$ norm, with common choices including $L_1$, $L_2$ and $L_\infty$.[1] The method of global criterion is sensitive to the scaling of the objective functions, and thus, it is recommended that the objectives are normalized into a uniform, dimensionless scale.[1][8]

## A priori methods

A priori methods require that sufficient preference information is expressed before the solution process.[2] Well-known examples of a priori methods include the utility function method, lexicographic method, and goal programming.

In the utility function method, it is assumed that the decision maker's utility function is available. A mapping $u\colon Y\rightarrow\mathbb{R}$ is a utility function if for all $\mathbf{y}^1,\mathbf{y}^2\in Y$ it holds that $u(\mathbf{y}^1)>u(\mathbf{y}^2)$ if the decision maker prefers $\mathbf{y}^1$ to $\mathbf{y}^2$, and $u(\mathbf{y}^1)=u(\mathbf{y}^2)$ if the decision maker is indifferent between $\mathbf{y}^1$ and $\mathbf{y}^2$. The utility function specifies an ordering of the decision vectors (recall that vectors can be ordered in many different ways). Once $u$ is obtained, it suffices to solve

$\max\;u(\mathbf{f}(\mathbf{x}))\text{ subject to }\mathbf{x}\in X,$

but in practice it is very difficult to construct a utility function that would accurately represent the decision maker's preferences[1] - particularly since the Pareto front is unknown before the optimization begins.

Lexicographic method assumes that the objectives can be ranked in the order of importance. We can assume, without loss of generality, that the objective functions are in the order of importance so that $f_1$ is the most important and $f_k$ the least important to the decision maker. The lexicographic method consists of solving a sequence of single-objective optimization problems of the form

\begin{align} \min&f_l(\mathbf{x})\\ \text{s.t. }&f_j(\mathbf{x})\leq\mathbf{y}^*_j,\;j=1,\dotsc,l-1,\\ &\mathbf{x}\in X, \end{align}

where $\mathbf{y}^*_j$ is the optimal value of the above problem with $l=j$. Thus, $\mathbf{y}^*_1:=\min\{f_1(\mathbf{x})\mid\mathbf{x}\in X\}$ and each new problem of the form in the above problem in the sequence adds one new constraint as $l$ goes from $1$ to $k$.

## A posteriori methods

A posteriori methods aim at producing all the Pareto optimal solutions or a representative subset of the Pareto optimal solutions. Well-known examples are the Normal Boundary Intersection (NBI),[10] Modified Normal Boundary Intersection (NBIm),[11] Normal Constraint (NC),[12][13] Successive Pareto Optimization (SPO)[14] and Directed Search Domain (DSD)[15] methods that solve the multi-objective optimization problem by constructing several scalarizations. The solution to each scalarization yields a Pareto optimal solution, whether locally or globally. The scalarizations of the NBI, NBIm, NC and DSD methods are constructed with the target of obtaining evenly distributed Pareto points that give a good evenly distributed approximation of the real set of Pareto points.

Evolutionary algorithms are popular approaches to generating Pareto optimal solutions to a multi-objective optimization problem. Currently, most evolutionary multi-objective optimization (EMO) algorithms apply Pareto-based ranking schemes. Evolutionary algorithms such as the Non-dominated Sorting Genetic Algorithm-II (NSGA-II) and Strength Pareto Evolutionary Algorithm 2 (SPEA-2) have become standard approaches, although some schemes based on particle swarm optimization and simulated annealing[16] are significant. The main advantage of evolutionary algorithms, when applied to solve multi-objective optimization problems, is the fact that they typically generate sets of solutions, allowing computation of an approximation of the entire Pareto front. The main disadvantage of evolutionary algorithms is their lower speed and the Pareto optimality of the solutions cannot be guaranteed. It is only known that none of the generated solutions dominates the others.

Other a posteriori methods are:

## Interactive methods

In interactive methods, the solution process is iterative and the decision maker continuously interacts with the method when searching for the most preferred solution (see e.g.,,[1][21]). In other words, the decision maker is expected to express preferences at each iteration in order to get Pareto optimal solutions that are of interest to him/her and learn what kind of solutions are attainable. The following steps are commonly present in interactive methods:[21]

1. initialize (e.g., calculate ideal and approximated nadir objective vectors and show them to the decision maker)
2. generate a Pareto optimal starting point (by using e.g. some no-preference method or solution given by the decision maker)
3. ask for preference information from the decision maker (e.g., aspiration levels or number of new solutions to be generated)
4. generate new Pareto optimal solution(s) according to the preferences and show it/them and possibly some other information about the problem to the decision maker
5. if several solutions were generated, ask the decision maker to select the best solution so far
6. stop, if the decision maker wants to; otherwise, go to step 3).

Above, aspiration levels refer to desirable objective function values forming a reference point. Instead of mathematical convergence that is often used as a stopping criterion in mathematical optimization methods, a psychological convergence is emphasized in interactive methods. Generally speaking, a method is terminated when the decision maker is confident that (s)he has found the most preferred solution available.

Different interactive methods involve different types of preference information. For example, three types can be identified: methods based on 1) trade-off information, 2) reference points and 3) classification of objective functions.[21] On the other hand, a fourth type of generating a small sample of solutions is included in[22] and.[23] An example of interactive method utilizing trade-off information is the Zionts-Wallenius method,[24] where the decision maker is shown several objective trade-offs at each iteration, and (s)he is expected to say whether (s)he likes, dislikes or is indifferent with respect to each trade-off. In reference point based methods (see e.g.,,[25][26]), the decision maker is expected at each iteration to specify a reference point consisting of desired values for each objective and a corresponding Pareto optimal solution(s) is then computed and shown to him/her for analysis. In classification based interactive methods, the decision maker is assumed to give preferences in the form of classifying objectives at the current Pareto optimal solution into different classes indicating how the values of the objectives should be changed to get a more preferred solution. Then, the classification information given is taken into account when new (more preferred) Pareto optimal solution(s) are computed. In the satisficing trade-off method (STOM)[27] three classes are used: objectives whose values 1) should be improved, 2) can be relaxed, and 3) are acceptable as such. In the NIMBUS method,[28][29] two additional classes are also used: objectives whose values 4) should be improved until a given bound and 5) can be relaxed until a given bound.

## Hybrid methods

Different hybrid methods exist, but here we consider hybridizing MCDM (multi-criteria decision making) and EMO (evolutionary multi-objective optimization). A hybrid algorithm in the context of multi-objective optimization is a combination of algorithms/approaches from these two fields (see e.g.,[21]). Hybrid algorithms of EMO and MCDM are mainly used to overcome shortcomings by utilizing strengths. Several types of hybrid algorithms have been proposed in the literature, e.g., incorporating MCDM approaches into EMO algorithms as a local search operator and to lead a DM to the most preferred solution(s) etc. A local search operator is mainly used to enhance the rate of convergence of EMO algorithms.

The roots for hybrid multi-objective optimization can be traced to the first Dagstuhl seminar organized in November 2004 (see, here). Here some of the best minds in EMO (Professor Kalyanmoy Deb, Professor Jürgen Branke etc.) and MCDM (Professor Kaisa Miettinen, Professor Ralph E. Steuer etc.) realized the potential in combining ideas and approaches of MCDM and EMO fields to prepare hybrids of them. Subsequently many more Dagstuhl seminars have been arranged to foster collaboration. Recently, hybrid multi-objective optimization has become an important theme in several international conferences in the area of EMO and MCDM (see e.g.,[30] and.[31])

## Visualization of the Pareto front

Visualization of the Pareto front is one of the a posteriori preference techniques of multi-objective optimization. The a posteriori preference techniques (see, for example,[1]) provide an important class of multi-objective optimization techniques. Usually the a posteriori preference techniques include four steps: (1) computer approximates the Pareto front, i.e. the Pareto optimal set in the objective space; (2) the decision maker studies the Pareto front approximation; (3) the decision maker identifies the preferred point at the Pareto front; (4) computer provides the Pareto optimal decision, which output coincides with the objective point identified by the decision maker. From the point of view of the decision maker, the second step of the a posteriori preference techniques is the most complicated one. There are two main approaches to informing the decision maker. First, a number of points of the Pareto front can be provided in the form of a list (interesting discussion and references are given in[32]) or using Heatmaps.[33] Alternative idea consists in visualizing the Pareto front.

### Visualization in bi-objective problems: tradeoff curve

In the case of bi-objective problems, informing the decision maker concerning the Pareto front is usually carried out by its visualization: the Pareto front, often named the tradeoff curve in this case, can be drawn at the objective plane. The tradeoff curve gives full information on objective values and on objective tradeoffs, which inform how improving one objective is related to deteriorating the second one while moving along the tradeoff curve. The decision maker takes this information into account while specifying the preferred Pareto optimal objective point. The idea to approximate and visualize the Pareto front was introduced for linear bi-objective decision problems by S.Gass and T.Saaty.[34] This idea was developed and applied in environmental problems by J.L. Cohon.[35] A review of methods for approximating the Pareto front for various decision problems with a small number of objectives (mainly, two) is provided in.[36]

### Visualization in high-order multi-objective optimization problems

There are two generic ideas how to visualize the Pareto front in high-order multi-objective decision problems (problems with more than two objectives). One of them, which is applicable in the case of a relatively small number of objective points that represent the Pareto front, is based on using the visualization techniques developed in statistics (various diagrams, etc – see the corresponding subsection below). The second idea proposes the display of bi-objective cross-sections (slices) of the Pareto front. It was introduced by W.S. Meisel in 1973[37] who argued that such slices inform the decision maker on objective tradeoffs. The figures that display a series of bi-objective slices of the Pareto front for three-objective problems are known as the decision maps. They give a clear picture of tradeoffs between three criteria. Disadvantages of such an approach are related to two following facts. First, the computational procedures for constructing the bi-objective slices of the Pareto front are not stable since the Pareto front is usually not stable. Secondly, it is applicable in the case of only three objectives. In the 1980s, the idea W.S. Meisel of implemented in a different form – in the form of the Interactive Decision Maps (IDM) technique.[38]

### Interactive decision maps technique

The IDM technique is based on approximating the Edgeworth-Pareto Hull (EPH) of the feasible objective set, that is, the feasible objective set broadened by the objective points dominated by it. Alternatively, this set is known as Free Disposal Hull. It is important that the EPH has the same Pareto front as the feasible objective set, but the bi-objective slices of the EPH look much simpler. The frontiers of bi-objective slices of the EPH contain the slices of the Pareto front. It is important that, in contrast to the Pareto front itself, the EPH is usually stable in respect to disturbances of data. The IDM technique applies fast on-line display of bi-objective slices of the EPH approximated in advance.

Since the bi-objective slices of the EPH for two selected objectives are extending (or shrinking) monotonically, while the value of one of the other objectives (the “third” objective) changes monotonically, the frontiers of the slices of the EPH, for which the values only of the “third” objective changes, do not intersect. This is why a figure with superimposed bi-objective slices of the EPH looks like an ordinary topographical map and is named the decision map, too. To study influence of the other (fourth, fifth, etc.) objectives, one can use animation of the decision maps. Such animation is possible due to the preliminary approximating the EPH. Alternatively, one can study various collections of snap-shots of the animation. Computers can visualize the Pareto front in the form of decision maps for linear and nonlinear decision problems for three to about eight objectives. Computer networks are able to bring, for example, Java applets that display graphs of the Pareto fronts on request. Real-life applications of the IDM technique are described in.[38]

### Illustration of the IDM technique

The above figure represents a gray scale copy of a color computer display for a real-life water quality problem[38] involving five objectives. The decision map consists of four superimposed bi-objective differently colored slices. A palette shows the relation between the values of the “third” objective and colors. Two scroll-bars are related to the values of the fourth and the fifth objectives.

A movement of a scroll-bar results in a change of the decision map. One can move the slider manually. However, the most effective form of displaying information to the DM is based on an automatic movement of the slider, that is, on a gradual increment (or decrement) in the constraint imposed on the value of an objective. A fast replacement of the decision maps offers the effect of animation. Because any reasonable number of scroll-bars can be located on the display, one can explore the influence of the fourth, the fifth (and maybe even the sixth and the seventh etc.) objectives on the decision map.

### Approximating the EPH

The EPH must be approximated in the IDM technique before the decision maps are displayed. Methods for approximating the EPH depend on the convexity properties of the EPH. Approximation methods are typically based either on approximation of the EPH by a convex polyhedral set or on approximation of the EPH by a large but finite number of domination cones in objective space with vertices that are close to the Pareto front. The first form can be applied only in the convex problems, while the second form is universal and can be used in general nonlinear problems[38]

### Approximation and visualization in the case of convex EPH

The EPH approximated by a polyhedral set is described by a system of a finite number of linear inequalities, which must be constructed by the approximation technique. Mathematical theory of optimal polyhedral approximation of convex bodies was developed during recently, and its results can be applied for developing the effective methods for approximating the EPH (see details in[38]). A large number of bi-objective slices of such approximations can be computed and displayed in the form of a decision map in several seconds.

### Point-wise approximation of the Pareto front and its visualization

An EPH approximation by a large but finite number of domination cones can be constructed on the basis of any point-wise approximation of the Pareto front, which can be found by using a broad range of techniques from classic single-objective optimization methods[1] up to modern evolutionary methods[39] Hybrid methods for approximating the EPH based on combination of classic and evolutionary methods can be used, too.[40] The bi-objective slices of such approximation can be computed very fast as well. Application of these methods results in decision maps that look fairly understandable if the number of approximating points is sufficiently large.

### Search for the preferred decision

In the IDM technique, search for the preferred decision is based on identification of a preferred Pareto optimal objective point (feasible goal). Decision maps help the user to identify the goal directly at a tradeoff curve drawn at the computer display. Then, a Pareto optimal decision associated with the goal is found automatically.

Detailed discussion of the Pareto front visualization problems is provided in.[41]

## Multi-objective optimization software

Name (alphabetically) License Brief info
1000Minds free for academic software for Multi-Criteria Decision-Making, prioritisation and resource allocation. Internet-based.
Accord by Robust Decisions implementing the Bayesian Team Support technique;no longer commercially availble but still used as a support tool
BENSOLVE GPL Free MatLab implementation of Benson's algorithm to solve linear vector optimization problems
Decisionarium free for academic global space for decision support (for academic use)
D-Sight proprietary visual and interactive tool for multicriteria decision aid problems based on the PROMETHEE methods and Multi-Attribute Utility Theory
GUIMOO LGPL Graphical User Interface for Multi Objective Optimization from INRIA
IDS proprietary Intelligent Decision System for Multiple Criteria Decision Analysis
IDSS Software free for non-profit activities MCDM software of the Laboratory of Intelligent Decision Support Systems (University of Poznan, Poland)
IND-NIMBUS proprietary implementation of the interactive NIMBUS method that can be connected with different simulation and modelling tools
interalg BSD solver with specifiable accuracy from OpenOpt - free universal cross-platform numerical optimization framework written in Python language using NumPy arrays, see its MOP page and other problems involved.
IRIS and VIP proprietary IRIS - Interactive Robustness analysis and parameters' Inference softward for multicriteria Sorting problems and VIP - Variable Interdependent Parameters Analysis software
jMetal LGPL Java-based framework for multi-objecive optimization with metaheuristics.
MACBETH for MCDA proprietary Measuring Attractiveness by a Categorical Based Evaluation TecHnique in MultiCriteria Decison Aid
MakeItRational proprietary AHP based decision software
Collection of Multiple Criteria Decision Support Software different by Dr. Roland Weistroffer
NIMBUS free for academic for solving nonlinear (and even nondifferentiable) multiobjective optimization problems in an interactive way. Operates via the Internet.
ParadisEO-MOEO CeCill module specifically devoted to multiobjective optimization in ParadisEO, software framework for the design and implementation of metaheuristics, hybrid methods as well as parallel and distributed models from INRIA
PROMETHEE-GAIA free for academic
PaGMO / PyGMO free Parallel Global Multiobjective Optimizer (and its Python alter ego PyGMO) offers a user-friendly access to a wide array of global and local optimization algorithms and problems. The main purpose of the software is to provide a parallelization engine common to all algorithms through the 'generalized island model'. Initially developed within the European Space Agency, the code was intended to help the automated design of interplanetary trajectories and spacecraft transfers in general. The user can implement his own problem and algorithm both in C++ and in Python.
MCDA software by Quartzstar Ltd.: OnBalance for evaluation decisions and HiPriority for resource allocation
RGDB Graphic tool that helps to select preferable rows from relational databases
TransparentChoice proprietary Strategic decision-making software, MCDM software
VISA proprietary Web based Multi-Criteria Decision Making Software.

Weistroffer et al. have written a book chapter[42] on multi-objective optimization software.

## References

1. Kaisa Miettinen (1999). Nonlinear Multiobjective Optimization. Springer. ISBN 978-0-7923-8278-2. Retrieved 29 May 2012.
2. Ching-Lai Hwang; Abu Syed Md Masud (1979). Multiple objective decision making, methods and applications: a state-of-the-art survey. Springer-Verlag. ISBN 978-0-387-09111-2. Retrieved 29 May 2012.
3. ^ Ropponen, A.; Ritala, R.; Pistikopoulos, E. N. (2011). "Optimization issues of the broke management system in papermaking". Computers & Chemical Engineering 35 (11): 2510. doi:10.1016/j.compchemeng.2010.12.012. edit
4. ^ a b E. Björnson and E. Jorswieck, Optimal Resource Allocation in Coordinated Multi-Cell Systems, Foundations and Trends in Communications and Information Theory, vol. 9, no. 2-3, pp. 113-381, 2013.
5. ^ Z.-Q. Luo and S. Zhang, Dynamic spectrum management: Complexity and duality, IEEE Journal of Selected Topics in Signal Processing, vol. 2, no. 1, pp. 57–73, 2008.
6. ^ Matthias Ehrgott (1 June 2005). Multicriteria Optimization. Birkhäuser. ISBN 978-3-540-21398-7. Retrieved 29 May 2012.
7. ^ Carlos A. Coello Coello; Gary B. Lamont; David A. Van Veldhuisen (2007). Evolutionary Algorithms for Solving Multi-Objective Problems. Springer. ISBN 978-0-387-36797-2. Retrieved 1 November 2012.
8. ^ a b Jürgen Branke; Kalyanmoy Deb; Kaisa Miettinen; Roman Slowinski (21 November 2008). Multiobjective Optimization: Interactive and Evolutionary Approaches. Springer. ISBN 978-3-540-88907-6. Retrieved 1 November 2012.
9. ^ Zeleny, M. (1973), "Compromise Programming", in Cochrane, J.L.; Zeleny, M., Multiple Criteria Decision Making, University of South Carolina Press, Columbia, pp. 262–301
10. ^ Das, I.; Dennis, J. E. (1998). "Normal-Boundary Intersection: A New Method for Generating the Pareto Surface in Nonlinear Multicriteria Optimization Problems". SIAM Journal on Optimization 8 (3): 631. doi:10.1137/S1052623496307510. edit
11. ^ S. Motta, Renato; Afonso, Silvana M. B., Lyra, Paulo R. M. (8 January 2012). "A modified NBI and NC method for the solution of N-multiobjective optimization problems". Structural and Multidisciplinary Optimization. doi:10.1007/s00158-011-0729-5.
12. ^ Messac, A.; Ismail-Yahaya, A.; Mattson, C.A. (2003). "The normalized normal constraint method for generating the Pareto frontier". Structural and multidisciplinary optimization 25 (2): 86–98.
13. ^ Messac, A.; Mattson, C. A. (2004). "Normal constraint method with guarantee of even representation of complete Pareto frontier". AIAA journal 42 (10): 2101–2111.
14. ^ Mueller-Gritschneder, Daniel; Graeb, Helmut; Schlichtmann, Ulf (2009). "A Successive Approach to Compute the Bounded Pareto Front of Practical Multiobjective Optimization Problems". SIAM Journal on Optimization 20 (2): 915–934.
15. ^ Erfani, Tohid; Utyuzhnikov, Sergei V. (2011). "Directed Search Domain: A Method for Even Generation of Pareto Frontier in Multiobjective Optimization" (pdf). Journal of Engineering Optimization 43 (5): 1–18. Retrieved October 17, 2011.
16. ^ Suman, B.; Kumar, P. (2006). "A survey of simulated annealing as a tool for single and multiobjective optimization". Journal of the Operational Research Society 57 (10): 1143–1160. doi:10.1057/palgrave.jors.2602068.
17. ^ Craft, D.; Halabi, T.; Shih, H.; Bortfeld, T. (2006). "Approximating convex Pareto surfaces in multiobjective radiotherapy planning". Medical Physics 33 (9): 3399–3407.
18. ^ Beume, N.; Naujoks, B.; Emmerich, M. (2007). "SMS-EMOA: Multiobjective selection based on dominated hypervolume". European Journal of Operational Research 181 (3): 1653. doi:10.1016/j.ejor.2006.08.008. edit
19. ^ Battiti, Roberto; Mauro Brunato; Franco Mascia (2008). Reactive Search and Intelligent Optimization. Springer Verlag. ISBN 978-0-387-09623-0.
20. ^ Battiti, Roberto; Mauro Brunato (2011). Reactive Business Intelligence. From Data to Models to Insight.. Trento, Italy: Reactive Search Srl. ISBN 978-88-905795-0-9.
21. ^ a b c d Miettinen, K.; Ruiz, F.; Wierzbicki, A. P. (2008). "Introduction to Multiobjective Optimization: Interactive Approaches". Multiobjective Optimization. Lecture Notes in Computer Science 5252. p. 27. doi:10.1007/978-3-540-88908-3_2. ISBN 978-3-540-88907-6. edit
22. ^ Luque, M.; Ruiz, F.; Miettinen, K. (2008). "Global formulation for interactive multiobjective optimization". OR Spectrum 33: 27. doi:10.1007/s00291-008-0154-3. edit
23. ^ Ruiz, F.; Luque, M.; Miettinen, K. (2011). "Improving the computational efficiency in a global formulation (GLIDE) for interactive multiobjective optimization". Annals of Operations Research 197: 47. doi:10.1007/s10479-010-0831-x. edit
24. ^ Zionts, S.; Wallenius, J. (1976). "An Interactive Programming Method for Solving the Multiple Criteria Problem". Management Science 22 (6): 652. doi:10.1287/mnsc.22.6.652. edit
25. ^ Wierzbicki, A. P. (1986). "On the completeness and constructiveness of parametric characterizations to vector optimization problems". OR Spektrum 8 (2): 73–78. doi:10.1007/BF01719738. edit
26. ^ Andrzej P. Wierzbicki; Marek Makowski; Jaap Wessels (31 May 2000). Model-Based Decision Support Methodology with Environmental Applications. Springer. ISBN 978-0-7923-6327-9. Retrieved 17 September 2012.
27. ^ Nakayama, H.; Sawaragi, Y. (1984), "Satisficing Trade-Off Method for Multiobjective Programming", in Grauer, M.; Wierzbicki, A. P., Interactive Decision Analysis, Springer-Verlag Berlin, Heidelberg, pp. 113–122
28. ^ Miettinen, K.; Mäkelä, M. M. (1995). "Interactive bundle-based method for nondifferentiable multiobjeective optimization: Nimbus§". Optimization 34 (3): 231. doi:10.1080/02331939508844109. edit
29. ^ Miettinen, K.; Mäkelä, M. M. (2006). "Synchronous approach in interactive multiobjective optimization". European Journal of Operational Research 170 (3): 909. doi:10.1016/j.ejor.2004.07.052. edit
30. ^ Sindhya, K.; Ruiz, A. B.; Miettinen, K. (2011). "A Preference Based Interactive Evolutionary Algorithm for Multi-objective Optimization: PIE". Evolutionary Multi-Criterion Optimization. Lecture Notes in Computer Science 6576. p. 212. doi:10.1007/978-3-642-19893-9_15. ISBN 978-3-642-19892-2. edit
31. ^ Sindhya, K.; Deb, K.; Miettinen, K. (2008). "A Local Search Based Evolutionary Multi-objective Optimization Approach for Fast and Accurate Convergence". Parallel Problem Solving from Nature – PPSN X. Lecture Notes in Computer Science 5199. p. 815. doi:10.1007/978-3-540-87700-4_81. ISBN 978-3-540-87699-1. edit
32. ^ Benson, Harold P.; Sayin, Serpil (1997). "Towards finding global representations of the efficient set in multiple objective mathematical programming". Naval Research Logistics 44 (1): 47–67. doi:10.1002/(SICI)1520-6750(199702)44:1<47::AID-NAV3>3.0.CO;2-M. ISSN 0894-069X.
33. ^ Pryke, Andy; Sanaz Mostaghim, Alireza Nazemi (2007). "Heatmap Visualisation of Population Based Multi Objective Algorithms". Evolutionary Multi-Criterion Optimization: 361–375.
34. ^ Gass, Saul; Saaty, Thomas (1955). "The computational algorithm for the parametric objective function". Naval Research Logistics Quarterly 2 (1-2): 39–45. doi:10.1002/nav.3800020106. ISSN 0028-1441.
35. ^ Jared L. Cohon (13 January 2004). Multiobjective Programming and Planning. Courier Dover Publications. ISBN 978-0-486-43263-2. Retrieved 29 May 2012.
36. ^ Ruzika, S.; Wiecek, M. M. (2005). "Approximation Methods in Multiobjective Programming". Journal of Optimization Theory and Applications 126 (3): 473–501. doi:10.1007/s10957-005-5494-4. ISSN 0022-3239.
37. ^ Meisel, W. L. (1973), "Tradeoff decision in multiple criteria decision making", in J. L. Cochrane; M. Zeleny, Multiple Criteria Decision Making (S.C. University of Columbia): 461–476
38. A. V. Lotov; V. A. Bushenkov; G. K. Kamenev (29 February 2004). Interactive Decision Maps: Approximation and Visualization of Pareto Frontier. Springer. ISBN 978-1-4020-7631-2. Retrieved 29 May 2012.
39. ^ Kalyanmoy Deb (23 March 2009). Multi-Objective Optimization Using Evolutionary Algorithms. John Wiley & Sons. ISBN 978-0-470-74361-4. Retrieved 1 November 2012.
40. ^ Berezkin, V. E.; Kamenev, G. K.; Lotov, A. V. (2006). "Hybrid adaptive methods for approximating a nonconvex multidimensional Pareto frontier". Computational Mathematics and Mathematical Physics 46 (11): 1918. doi:10.1134/S096554250611008X. edit
41. ^ Lotov, A. V.; Miettinen, K. (2008). "Visualizing the Pareto Frontier". Multiobjective Optimization. Lecture Notes in Computer Science 5252. p. 213. doi:10.1007/978-3-540-88908-3_9. ISBN 978-3-540-88907-6. edit
42. ^ Weistroffer, H. R.; Smith, C. H.; Narula, S. C. (2005). "Multiple Criteria Decision Support Software". Multiple Criteria Decision Analysis: State of the Art Surveys. International Series in Operations Research & Management Science 78. p. 989. doi:10.1007/0-387-23081-5_24. ISBN 0-387-23067-X. edit