Metric regularity theory lies in the very heart of variational analysis, a relatively new discipline whose appearance was, to a large extent, determined by the needs of modern optimization theory in which such phenomena as nondifferentiability and set-valued mappings naturally appear. The roots of the theory go back to such fundamental results of the classical analysis as the implicit function theorem, Sard theorem and some others. The paper offers a survey of the state of the art of some principal parts of the theory along with a variety of its applications in analysis and optimization.

In this paper, a primal-dual interior point method is proposed for general constrained optimization, which incorporated a penalty function and a kind of new identification technique of the active set. At each iteration, the proposed algorithm only needs to solve two or three reduced systems of linear equations with the same coefficient matrix. The size of systems of linear equations can be decreased due to the introduction of the working set, which is an estimate of the active set. The penalty parameter is automatically updated and the uniformly positive definiteness condition on the Hessian approximation of the Lagrangian is relaxed. The proposed algorithm possesses global and superlinear convergence under some mild conditions. Finally, some preliminary numerical results are reported.

This paper investigates the eigenmode optimization problem governed by the scalar Helmholtz equation in continuum system in which the computed eigenmode approaches the prescribed eigenmode in the whole domain. The first variation for the eigenmode optimization problem is evaluated by the quadratic penalty method, the adjoint variable method, and the formula based on sensitivity analysis. A penalty optimization algorithm is proposed, in which the density evolution is accomplished by introducing an artificial time term and solving an additional ordinary differential equation. The validity of the presented algorithm is confirmed by numerical results of the first and second eigenmode optimizations in 1D and 2D problems.

Sequential decision problems can often be modeled as Markov decision processes. Classical solution approaches assume that the parameters of the model are known. However, model parameters are usually estimated and uncertain in practice. As a result, managers are often interested in how estimation errors affect the optimal solution. In this paper we illustrate how sensitivity analysis can be performed directly for a Markov decision process with uncertain reward parameters using the Bellman equations. In particular, we consider problems involving (i) a single stationary parameter, (ii) multiple stationary parameters, and (iii) multiple nonstationary parameters. We illustrate the applicability of this work through a capacitated stochastic lot-sizing problem.

We study Taylor series expansions of stationary characteristics of general-state-space Markov chains. The elements of the Taylor series are explicitly calculated and a lower bound for the radius of convergence of the Taylor series is established. The analysis provided in this paper applies to the case where the stationary characteristic is given through an unbounded sample performance function such as the second moment of the stationary waiting time in a queueing system.

We consider a multicomponent maintenance system controlled by an age replacement policy: when one of the components fails, it is immediately replaced; all components older than a threshold age θ are preventively replaced. Costs are associated with each maintenance action, such as replacement after failure or preventive replacement. We derive a weak derivative estimator for the derivative of the cost performance with respect to θ. The technique is quite general and can be applied to many other threshold optimization problems in maintenance. The estimator is easy to implement and considerably increases the efficiency of a Robbins-Monro type of stochastic approximation algorithm. The paper is self-contained in the sense that it includes a proof of the correctness of the weak derivative estimation algorithm.