9,871 research outputs found
A path following algorithm for the graph matching problem
We propose a convex-concave programming approach for the labeled weighted
graph matching problem. The convex-concave programming formulation is obtained
by rewriting the weighted graph matching problem as a least-square problem on
the set of permutation matrices and relaxing it to two different optimization
problems: a quadratic convex and a quadratic concave optimization problem on
the set of doubly stochastic matrices. The concave relaxation has the same
global minimum as the initial graph matching problem, but the search for its
global minimum is also a hard combinatorial problem. We therefore construct an
approximation of the concave problem solution by following a solution path of a
convex-concave problem obtained by linear interpolation of the convex and
concave formulations, starting from the convex relaxation. This method allows
to easily integrate the information on graph label similarities into the
optimization problem, and therefore to perform labeled weighted graph matching.
The algorithm is compared with some of the best performing graph matching
methods on four datasets: simulated graphs, QAPLib, retina vessel images and
handwritten chinese characters. In all cases, the results are competitive with
the state-of-the-art.Comment: 23 pages, 13 figures,typo correction, new results in sections 4,5,
Numerical Schubert calculus
We develop numerical homotopy algorithms for solving systems of polynomial
equations arising from the classical Schubert calculus. These homotopies are
optimal in that generically no paths diverge. For problems defined by
hypersurface Schubert conditions we give two algorithms based on extrinsic
deformations of the Grassmannian: one is derived from a Gr\"obner basis for the
Pl\"ucker ideal of the Grassmannian and the other from a SAGBI basis for its
projective coordinate ring. The more general case of special Schubert
conditions is solved by delicate intrinsic deformations, called Pieri
homotopies, which first arose in the study of enumerative geometry over the
real numbers. Computational results are presented and applications to control
theory are discussed.Comment: 24 pages, LaTeX 2e with 2 figures, used epsf.st
Automatic Frechet differentiation for the numerical solution of boundary-value problems
A new solver for nonlinear boundary-value problems (BVPs) in Matlab is presented, based on the Chebfun software system for representing functions and operators automatically as numerical objects. The solver implements Newton's method in function space, where instead of the usual Jacobian matrices, the derivatives involved are Frechet derivatives. A major novelty of this approach is the application of automatic differentiation (AD) techniques to compute the operator-valued Frechet derivatives in the continuous context. Other novelties include the use of anonymous functions and numbering of each variable to enable a recursive, delayed evaluation of derivatives with forward mode AD. The AD techniques are applied within a new Chebfun class called chebop which allows users to set up and solve nonlinear BVPs in a few lines of code, using the "nonlinear backslash" operator (\). This framework enables one to study the behaviour of Newton's method in function space
- …