Refine
Language
- English (16)
Keywords
- KKT-System (1)
- LU decomposition (1)
- TR1 Update (1)
- adjoint based updating (1)
- chromatic staff notation (1)
- convergence order (1)
- equally tempered scale (1)
- local convergence (1)
- major thirds (1)
- nonlinear equations (1)
Automatic, or algorithmic, differentiation addresses the need for the accurate
and efficient calculation of derivative values in scientific computing. To this
end procedural programs for the evaluation of problem-specific functions are
transformed into programs that also compute the required derivative values
at the same numerical arguments in floating point arithmetic. Disregarding
many important implementation issues, we examine in this article complexity
bounds and other more mathematical aspects of the program transformation
task sketched above.
In this paper we introduce a new method for the computation of KKT matrices that arise from solving constrained, nonlinear optimization problems. This method requires updating of null-space factorizations after a low rank modification. The update procedure has the advantage that it is significantly cheaper than a re-factorization of the system at each new iterate. This paper focuses on the cheap update of a rectangular LU decomposition after a rank-1 modification.
Two different procedures for updating the LU factorization are presented in detail and compared regarding their costs of computation and their stability. Moreover we will introduce an extension of these algorithms which further improves the computation time. This turns out to be an excellent alternative to algorithms based on orthogonal transformations.
Quasi-Newton methods based on least change secant updating
formulas that solve linear equations $Ax=b$ in $n=\dim(x)=\dim(b)$ steps
can be expected to solve corresponding smooth nonlinear
systems $n$-step quadratically, i.e. with an $r$-order
of $\rho = 2^{1/n} = 1 + 1/n +O(1/n^2)$. The best rate one can
possibly expect on general problems is given by the positive root
$\rho_n$ of $\rho^n(\rho -1)=1$, for which
$\rho_n-1 = \ln(n)/n + O(1/n^2)$. To show that this upper bound is
actually achieved one usually has to impose a priori some kind of
linear independence condition on the sequence of steps taken by the
quasi-Newton iteration in question. Without any such assumptions we
establish in this paper the convergence order $\rho_n$ for the
two-sided rank one formula proposed by Schlenkrich et al in \cite{SGW06}.
It requires the evaluation of adjoint vectors, is invariant with respect
to linear transformations on the variable domain and combines the
properties of bounded deterioration and heredity.
For the solution of nonlinear equation systems
quasi-Newton methods based on low-rank updates are of particular interest. We analyze a class
of TR1 update formulas to approximate the system Jacobian. The local q-superlinear convergence for nonlinear problems is proved for a particular subclass of updates. Moreover, we give an estimate of the r-order of convergence. Numerical results comparing the TR1 method to Newton's and other quasi-Newton methods atr presented.
Cross–derivatives are mixed partial derivatives that are obtained by differentiating at most
once in every coordinate direction. They are a computational tool in combinatorics and high–
dimensional integration. Here we present two methods of computing exact values of all cross–
derivatives at a given point both following the general philosophy of automatic differentiation.
Implementation details are discussed and numerical results given.
This paper proposes a new mathematical model for the open pit mine planning problem,
based on continuous functional analysis. The traditional models for this problem have been
constructed by using discrete 0-1 decision variables, giving rise to large-scale combinatorial
and Mixed Integer Programming (MIP) problems. Instead, we use a continuous approach
which allows for a refined imposition of slope constraints associated with geotechnical stability.
The model introduced here is posed in a suitable functional space, essentially the
real-valued functions that are Lipschitz continuous on a given two dimensional bounded region.
We derive existence results and investigate some qualitative properties of the solutions
Duality Results for Stationary Problems of Open Pit Mine Planning in a Continuous Function Framework
(2011)
Open Pit Mine Planning problems are usually considered in a Mixed Integer Programming
context. Characterizing each attainable profile by a continuous function yields a continuous
framework. It allows for a more detailed modeling of slope constraints and other material
properties of slanted layers. Although the resulting nonlinear programming problems are
in general non-convex and non-differentiable, they provide certain advantages as one can
directly compute sensitivities of optimal solutions w.r.t. small data perturbations. In this
work duality results are derived for the stationary problems of the continuous framework
employing an additional condition called convex-likeness.