20,151 research outputs found
SDPNAL+: A Matlab software for semidefinite programming with bound constraints (version 1.0)
SDPNAL+ is a {\sc Matlab} software package that implements an augmented
Lagrangian based method to solve large scale semidefinite programming problems
with bound constraints. The implementation was initially based on a majorized
semismooth Newton-CG augmented Lagrangian method, here we designed it within an
inexact symmetric Gauss-Seidel based semi-proximal ADMM/ALM (alternating
direction method of multipliers/augmented Lagrangian method) framework for the
purpose of deriving simpler stopping conditions and closing the gap between the
practical implementation of the algorithm and the theoretical algorithm. The
basic code is written in {\sc Matlab}, but some subroutines in C language are
incorporated via Mex files. We also design a convenient interface for users to
input their SDP models into the solver. Numerous problems arising from
combinatorial optimization and binary integer quadratic programming problems
have been tested to evaluate the performance of the solver. Extensive numerical
experiments conducted in [Yang, Sun, and Toh, Mathematical Programming
Computation, 7 (2015), pp. 331--366] show that the proposed method is quite
efficient and robust, in that it is able to solve 98.9\% of the 745 test
instances of SDP problems arising from various applications to the accuracy of
in the relative KKT residual
Modeling and solving the multi-period inventory routing problem with constant demand rates
The inventory routing problem (IRP) is one of the challenging optimization problems in supply chain logistics. It combines inventory control and vehicle routing optimization. The main purpose of the IRP is to determine optimal delivery times and quantities to be delivered to customers, as well as optimal vehicle routes to distribute these quantities. The IRP is an underlying logistical optimization problem for supply chains implementing vendor-managed inventory (VMI) policies, in which the supplier takes responsibility for the management of the customers' inventory. In this paper, we consider a multi-period inventory routing problem assuming constant demand rates (MP-CIRP). The proposed model is formulated as a linear mixed-integer program and solved with a Lagrangian relaxation method. The solution obtained by the Lagrangian relaxation method is then used to generate a close to optimal feasible solution of the MP-CIRP by solving a series of assignment problems. The numerical experiments carried out so far show that the proposed Lagrangian relaxation approach nds quite good solutions for the MP-CIRP and in reasonable computation times
Post-Newtonian SPH calculations of binary neutron star coalescence. I. Method and first results
We present the first results from our Post-Newtonian (PN) Smoothed Particle
Hydrodynamics (SPH) code, which has been used to study the coalescence of
binary neutron star (NS) systems. The Lagrangian particle-based code
incorporates consistently all lowest-order (1PN) relativistic effects, as well
as gravitational radiation reaction, the lowest-order dissipative term in
general relativity. We test our code on sequences of single NS models of
varying compactness, and we discuss ways to make PN simulations more relevant
to realistic NS models. We also present a PN SPH relaxation procedure for
constructing equilibrium models of synchronized binaries, and we use these
equilibrium models as initial conditions for our dynamical calculations of
binary coalescence. Though unphysical, since tidal synchronization is not
expected in NS binaries, these initial conditions allow us to compare our PN
work with previous Newtonian results.
We compare calculations with and without 1PN effects, for NS with stiff
equations of state, modeled as polytropes with . We find that 1PN
effects can play a major role in the coalescence, accelerating the final
inspiral and causing a significant misalignment in the binary just prior to
final merging. In addition, the character of the gravitational wave signal is
altered dramatically, showing strong modulation of the exponentially decaying
waveform near the end of the merger. We also discuss briefly the implications
of our results for models of gamma-ray bursts at cosmological distances.Comment: RevTeX, 37 pages, 17 figures, to appear in Phys. Rev. D, minor
corrections onl
- …