Quadruple and octuple layer potentials in two dimensions I: Analytical apparatus

A detailed analysis is presented of all pseudo-differential operators of orders up to 2 encountered in classical potential theory in two dimensions. Each of the operators under investigation turns out to be a sum of one or more of standard operators (second derivative, derivative of the Hilbert transform, etc.), and an integral operator with smooth kernel. This classification leads to an extremely simple analysis of spectra of such operators, and simplifies the design of procedures for their numerical evaluation. In a sequel to this paper, the obtained apparatus will be used to construct stable discretizations of arbitrarily high order for a variety of boundary value problems for elliptic partial differential equations.     

On probabilistic constraints induced by rectangular sets and multivariate normal distributions

In this paper, we consider optimization problems under probabilistic constraints which are defined by two-sided inequalities for the underlying normally distributed random vector. As a main step for an algorithmic solution of such problems, we prove a derivative formula for (normal) probabilities of rectangles as functions of their lower or upper bounds. This formula allows to reduce the calculus of such derivatives to the calculus of (normal) probabilities of rectangles themselves thus generalizing a similar well-known statement for multivariate normal distribution functions. As an application, we consider a problem from water reservoir management. One of the outcomes of the problem solution is that the (still frequently encountered) use of simple individual probabilistic constraints can completely fail. By contrast, the (more difficult) use of joint probabilistic constraints, which heavily depends on the derivative formula mentioned before, yields very reasonable and robust solutions over the whole time horizon considered.     

Genetic algorithm-based subproblem solution procedures for a modified shifting bottleneck heuristic for complex job shops

In this paper, we consider a modified shifting bottleneck heuristic for complex job shops. The considered job shop environment contains parallel batching machines, machines with sequence-dependent setup times and reentrant process flows. Semiconductor wafer fabrication facilities (Wafer Fabs) are typical examples for manufacturing systems with these characteristics. Our primary performance measure is total weighted tardiness (TWT). The shifting bottleneck heuristic uses a disjunctive graph to decompose the overall scheduling into scheduling problems for single tool groups. The scheduling algorithms for these scheduling problems are called subproblem solution procedures (SSPs). In previous research, only subproblem solution procedures based on dispatching rules have been considered. In this paper, we are interested in how much we can gain in terms of TWT if we apply more sophisticated subproblem solution procedures like genetic algorithms for parallel machine scheduling. We conduct simulation experiments in a dynamic job shop environment in order to assess the performance of the suggested subproblem solution procedures. It turns out that using near to optimal subproblem solution procedures leads in many situations to improved results compared to dispatching-based subproblem solution procedures.     

Necessary and Sufficient Conditions for the Interchange Between Infimum and the Symbol of Integration

We make a study of various notions of decomposability for subsets of measurable functions in relation with the interchange results between infimum and integration. For this we introduce the notions of serial decomposability and of decomposability relatively to an integrand. A characterization of closed serially decomposable subsets of the Lebesgue spaces L _p is given. The second notion of decomposability introduced is characteristic for the interchange property studied. Many examples are presented. The links are made with R. T. Rockafellar’s decomposability, F. Hiai, H. Umegaki’s decomposability, G. Bouchitté and M. Valadier’s stability and normal decomposability introduced by O. Anza Hafsa and J.-P. Mandallena. As applications we obtain exact lower bounds for minimization problems of integral functionals on normally decomposable spaces (spaces of continuous functions for example), and for the minimization of a class of functionals of the Calculus of Variations.     

An approximate numerical model for radiation problems

A variety of engineering problems can be successfully solved by coupling finite element and boundary element procedures. Approximate boundary elements, which can be used when dealing with radiation problems in unlimited domains are presented. They are simple to implement and can be easily inserted in existing frontal solution packages. Numerical examples are also reported.     

New formulation and relaxation to solve a concave-cost network flow problem

In this paper we study a minimum cost, multicommodity network flow problem in which the total cost is piecewise linear, concave of the total flow along the arcs. Specifically, the problem can be defined as follows. Given a directed network, a set of pairs of communicating nodes and a set of available capacity ranges and their corresponding variable and fixed cost components for each arc, the problem is to select for each arc a range and identify a path for each commodity between its source and destination nodes so as to minimize the total costs. We also extend the problem to the case of piecewise nonlinear, concave cost function. New mathematical programming formulations of the problems are presented. Efficient solution procedures based on Lagrangean relaxations of the problems are developed. Extensive computational results across a variety of networks are reported. These results indicate that the solution procedures are effective for a wide range of traffic loads and different cost structures. They also show that this work represents an improvement over previous work made by other authors. This improvement is the result of the introduction of the new formulations of the problems and their relaxations.     

Context-specific independence in graphical log-linear models

Log-linear models are the popular workhorses of analyzing contingency tables. A log-linear parameterization of an interaction model can be more expressive than a direct parameterization based on probabilities, leading to a powerful way of defining restrictions derived from marginal, conditional and context-specific independence. However, parameter estimation is often simpler under a direct parameterization, provided that the model enjoys certain decomposability properties. Here we introduce a cyclical projection algorithm for obtaining maximum likelihood estimates of log-linear parameters under an arbitrary context-specific graphical log-linear model, which needs not satisfy criteria of decomposability. We illustrate that lifting the restriction of decomposability makes the models more expressive, such that additional context-specific independencies embedded in real data can be identified. It is also shown how a context-specific graphical model can correspond to a non-hierarchical log-linear parameterization with a concise interpretation. This observation can pave way to further development of non-hierarchical log-linear models, which have been largely neglected due to their believed lack of interpretability.     

Effective and simple EOQ-like solutions for stochastic demand periodic review systems

In this paper we examine a periodic review system under stochastic demand with variable stockout costs. The optimal values for cycle length and amount of safety stock are difficult to obtain because one of the First Order Conditions does not have a closed form solution. However, by using a Taylor series expansion to approximate part of the cost function, we produce a simple cost function structure which is similar to that of deterministic models.We argue that this simple structure is also beneficial to promote the solution in other problems where coordination of cycles is required. To illustrate, we use the joint replenishment problem for multiple items under stochastic demand and suggest simple and efficient solution procedures.     

Layering strategies for creating exploitable structure in linear and integer programs

The strategy of subdividing optimization problems into layers by splitting variables into multiple copies has proved useful as a method for inducing exploitable structure in a variety of applications, particularly those involving embedded pure and generalized networks. A framework is proposed in this paper which leads to new relaxation and restriction methods for linear and integer programming based on our extension of this strategy. This framework underscores the use of constructions that lead to stronger relaxations and more flexible strategies than previous applications. Our results establish the equivalence of all layered Lagrangeans formed by parameterizing the equal value requirement of copied variables for different choices of the principal layers. It is further shown that these Lagrangeans dominate traditional Lagrangeans based on incorporating non-principal layers into the objective function. In addition a means for exploiting the layered Lagrangeans is provided by generating subgradients based on a simple averaging calculation. Finally, we show how this new layering strategy can be augmented by an integrated relaxation/restriction procedure, and indicate variations that can be employed to particular advantage in a parallel processing environment. Preliminary computational results on fifteen real world zero-one personnel assignment problems, comparing two layering approaches with five procedures previously found best for those problems, are encouraging. One of the layering strategies tested dominated all non-layering procedures in terms of both quality and solution time.This research was supported in part by the Office of Naval Research Contract N00014-78-C-0222 with the Center for Business Decision Analysis and by the US Department of Agriculture Contract 51-3142-4020 with Management Science Software Systems.     

Well-posedness for a cauchy problem associated to time dependent free boundaries with nonlocal leading terms

We study a class of free boundary problems, where the normal velocity of the interface is proportional to the derivative of the solution of an elliptic PDE; we give a simple, explicit criterion for the well-posedness of the linearized Cauchy problem. The method is then applied to two classical problems; the stefan problem and the Muskat problem.     

Flowshop scheduling research after five decades

Since Johnson’s seminal paper in 1954, flowshop scheduling problems have received considerable research attention over the last fifty years. As a result, several optimization and heuristic solution procedures are available to solve a variety of flowshop scheduling problems. This paper provides a brief glimpse into the evolution of flowshop scheduling problems and possible approaches for their solution over the last fifty years. It briefly introduces the current flowshop problems being solved and the approaches being taken to solve (optimally or approximately) them. The paper concludes with some fruitful directions for future research.     

A Neural Network Application in Personnel Scheduling

In this paper, we present and evaluate a neural network model for solving a typical personnel-scheduling problem, i.e. an airport ground staff rostering problem. Personnel scheduling problems are widely found in servicing and manufacturing industries. The inherent complexity of personnel scheduling problems has normally resulted in the development of integer programming-based models and various heuristic solution procedures. The neural network approach has been admitted as a promising alternative to solving a variety of combinatorial optimization problems. While few works relate neural network to applications of personnel scheduling problems, there is great theoretical and practical value in exploring the potential of this area. In this paper, we introduce a neural network model following a relatively new modeling approach to solve a real rostering case. We show how to convert a mixed integer programming formulation to a neural network model. We also provide the experiment results comparing the neural network method with three popular heuristics, i.e. simulated annealing, Tabu search and genetic algorithm. The computational study reveals some potential of neural networks in solving personnel scheduling problems.     

Implicit algorithms for solving the Cauchy problem for the equations of the dynamics of mechanical systems

It is shown that in the numerical solution of the Cauchy problem for systems of second-order ordinary differential equations, when solved for the highest-order derivative, it is possible to construct simple and economical implicit computational algorithms for step-by-step integration without using laborious iterative procedures based on processes of the Newton-Raphson iterative type. The initial problem must first be transformed to a new argument — the length of its integral curve. Such a transformation is carried out using an equation relating the initial parameter of the problem to the length of the integral curve. The linear acceleration method is used as an example to demonstrate the procedure of constructing an implicit algorithm using simple iterations for the numerical solution of the transformed Cauchy problem. Propositions concerning the computational properties of the iterative process are formulated and proved. Explicit estimates are given for an integration stepsize that guarantees the convergence of the simple iterations. The efficacy of the proposed procedure is demonstrated by the numerical solution of three problems. A comparative analysis is carried out of the numerical solutions obtained with and without parametrization of the initial problems in these three settings. As a qualitative test the problem of the celestial mechanics of the “Pleiades” is considered. The second example is devoted to modelling the non-linear dynamics of an elastic flexible rod fixed at one end as a cantilever and coiled in its initial (static) state into a ring by a bending moment. The third example demonstrates the numerical solution of the problem of the “unfolding” of a mechanical system consisting of three flexible rods with given control input.     

Subgradient Methods for Sharp Weakly Convex Functions

Subgradient methods converge linearly on a convex function that grows sharply away from its solution set. In this work, we show that the same is true for sharp functions that are only weakly convex, provided that the subgradient methods are initialized within a fixed tube around the solution set. A variety of statistical and signal processing tasks come equipped with good initialization and provably lead to formulations that are both weakly convex and sharp. Therefore, in such settings, subgradient methods can serve as inexpensive local search procedures. We illustrate the proposed techniques on phase retrieval and covariance estimation problems.     

Absalom: Balancing assembly lines with assignment restrictions

Assembly line balancing problems (ALBPs) arise whenever an assembly line is configured, redesigned or adjusted. An ALBP consists of distributing the total workload for manufacturing products among the work stations along the line. On the one hand, research has focussed on developing effective and fast solution methods for exactly solving the simple assembly line balancing problem (SALBP). On the other hand, a number of real-world extensions of SALBP have been introduced but solved with straight-forward and simple heuristics in many cases. Therefore, there is a lack of procedures for exactly solving such generalized ALBP.In this paper, we show how to extend the well-known solution procedure Salome [Scholl, A., Klein, R., 1997. Salome: A bidirectional branch-and-bound procedure for assembly line balancing. Informs J. Comput. 9 319–334], which is able to solve even large SALBP instances in a very effective manner, to a problem extension with different types of assignment restrictions (called ARALBP). The extended procedure, referred to as Absalom, employs a favorable branching scheme, an arsenal of bounding rules and a variety of logical tests using ideas from constraint programming.Computational experiments show that Absalom is a very promising exact solution approach although the additional assignment restrictions complicate the problem considerably and necessitate a relaxation of some components of Salome.     

Heuristic solutions and confidence intervals for the multicovering problem

The multicovering problem can be expressed as: Minimize CX subject to AX ⩾, b, X ϵ {0, 1}, where A is a zero-one matrix and b is a vector of positive integers. This mathematical model has many applications to scheduling and location problems. Large examples of such problems arise in industrial, commercial and military settings, and their size frequently exceeds the limits of computational tractability. For this important problem, we examine a variety of simple heuristic approaches which can be applied when optimal solutions are not available. The approximate solutions thus generated are used to construct confidence intervals for the unknown value of the optimal solution. A large-scale computational study for randomly generated problmes suggests that these intervals are both very narrow and very likely to contain the optimal solution value. A study of 10 very large real-world problems further supports the success of our methodology and the quality of the approximate solutions found.     

The stop-and-drop problem in nonprofit food distribution networks

In this paper, we introduce the stop-and-drop problem (SDRP), a new variant of location-routing problems, that is mostly applicable to nonprofit food distribution networks. In these distribution problems, there is a central warehouse that contains food items to be delivered to agencies serving the people in need. The food is delivered by trucks to multiple sites in the service area and partner agencies travel to these sites to pick up their food. The tactical decision problem in this setting involves how to jointly select a set of delivery sites, assign agencies to these sites, and schedule routes for the delivery vehicles. The problem is modeled as an integrated mixed-integer program for which we delineate a two-phase sequential solution approach. We also propose two Benders decomposition-based solution procedures, namely a linear programming relaxation based Benders implementation and a logic-based Benders decomposition heuristic. We show through a set of realistic problem instances that given a fixed time limit, these decomposition based methods perform better than both the standard branch-and-bound solution and the two-phase approach. The general problem and the realistic instances used in the computational study are motivated by interactions with food banks in southeastern United States.     

Solving non-linear portfolio optimization problems with the primal-dual interior point method

Stochastic programming is recognized as a powerful tool to help decision making under uncertainty in financial planning. The deterministic equivalent formulations of these stochastic programs have huge dimensions even for moderate numbers of assets, time stages and scenarios per time stage. So far models treated by mathematical programming approaches have been limited to simple linear or quadratic models due to the inability of currently available solvers to solve NLP problems of typical sizes. However stochastic programming problems are highly structured. The key to the efficient solution of such problems is therefore the ability to exploit their structure. Interior point methods are well-suited to the solution of very large non-linear optimization problems. In this paper we exploit this feature and show how portfolio optimization problems with sizes measured in millions of constraints and decision variables, featuring constraints on semi-variance, skewness or non-linear utility functions in the objective, can be solved with the state-of-the-art solver.     

Solving neutral delay differential equations with state-dependent delays

In this paper we consider a class of neutral delay differential equations with state dependent delays. For such equations the possible discontinuity in the derivative of the solution at the initial point may propagate along the integration interval giving rise to subsequent points, called “breaking points”, where the solution derivative is still discontinuous. As a consequence, in a right neighbourhood of each such point we have to face a Cauchy problem where the equation has a discontinuous right-hand side. In this case the existence and the uniqueness of the solution is no longer guaranteed to the right of such points and hence the solution of the neutral equation may either cease to exist or bifurcate. After illustrating why uniqueness and existence of the solution is no longer guaranteed for general state-dependent problems and showing a possible way to detect these occurrences automatically, we explain how to generalize/regularize the problem in order to suitably extend the solution beyond the breaking point. This is important, for example, when exploring numerically the presence of possible periodic orbits.     

Tutorial on Surrogate Constraint Approaches for Optimization in Graphs

Surrogate constraint methods have been embedded in a variety of mathematical programming applications over the past thirty years, yet their potential uses and underlying principles remain incompletely understood by a large segment of the optimization community. In a number of significant domains of combinatorial optimization, researchers have produced solution strategies without recognizing that they can be derived as special instances of surrogate constraint methods. Once the connection to surrogate constraint ideas is exposed, additional ways to exploit this framework become visible, frequently offering opportunities for improvement.We provide a tutorial on surrogate constraint approaches for optimization in graphs, illustrating the key ideas by reference to independent set and graph coloring problems, including constructions for weighted independent sets which have applications to associated covering and weighted maximum clique problems. In these settings, the surrogate constraints can be generated relative to well-known packing and covering formulations that are convenient for exposing key notions. The surrogate constraint approaches yield widely used heuristics for identifying independent sets as simple special cases, and also afford previously unidentified heuristics that have greater power in these settings. Our tutorial also shows how the use of surrogate constraints can be placed within the context of vocabulary building strategies for independent set and coloring problems, providing a framework for applying surrogate constraints that can be used in other applications.At a higher level, we show how to make use of surrogate constraint information, together with specialized algorithms for solving associated sub-problems, to obtain stronger objective function bounds and improved choice rules for heuristic or exact methods. The theorems that support these developments yield further strategies for exploiting surrogate constraint relaxations, both in graph optimization and integer programming generally.