Convex Hull of Planar H-Polyhedra

Suppose are planar (convex) H-polyhedra, that is, $A_i \in \mathbb{R}^{n_i \times 2}$ and $\vec{c}_i \in \mathbb{R}^{n_i}$. Let $P_i = \{\vec{x} \in \mathbb{R}^2 \mid A_i\vec{x} \leq \vec{c}_i \}$ and $n = n_1 + n_2$. We present an $O(n \log n)$ algorithm for calculating an H-polyhedron with the smallest $P = \{\vec{x} \in \mathbb{R}^2 \mid A\vec{x} \leq \vec{c} \}$ such that $P_1 \cup P_2 \subseteq P$

A Complete Characterization of the Gap between Convexity and SOS-Convexity

Our first contribution in this paper is to prove that three natural sum of squares (sos) based sufficient conditions for convexity of polynomials, via the definition of convexity, its first order characterization, and its second order characterization, are equivalent. These three equivalent algebraic conditions, henceforth referred to as sos-convexity, can be checked by semidefinite programming whereas deciding convexity is NP-hard. If we denote the set of convex and sos-convex polynomials in $n$ variables of degree $d$ with $\tilde{C}_{n,d}$ and $\tilde{\Sigma C}_{n,d}$ respectively, then our main contribution is to prove that $\tilde{C}_{n,d}=\tilde{\Sigma C}_{n,d}$ if and only if $n=1$ or $d=2$ or $(n,d)=(2,4)$. We also present a complete characterization for forms (homogeneous polynomials) except for the case $(n,d)=(3,4)$ which is joint work with G. Blekherman and is to be published elsewhere. Our result states that the set $C_{n,d}$ of convex forms in $n$ variables of degree $d$ equals the set $\Sigma C_{n,d}$ of sos-convex forms if and only if $n=2$ or $d=2$ or $(n,d)=(3,4)$. To prove these results, we present in particular explicit examples of polynomials in $\tilde{C}_{2,6}\setminus\tilde{\Sigma C}_{2,6}$ and $\tilde{C}_{3,4}\setminus\tilde{\Sigma C}_{3,4}$ and forms in $C_{3,6}\setminus\Sigma C_{3,6}$ and $C_{4,4}\setminus\Sigma C_{4,4}$, and a general procedure for constructing forms in $C_{n,d+2}\setminus\Sigma C_{n,d+2}$ from nonnegative but not sos forms in $n$ variables and degree $d$. Although for disparate reasons, the remarkable outcome is that convex polynomials (resp. forms) are sos-convex exactly in cases where nonnegative polynomials (resp. forms) are sums of squares, as characterized by Hilbert.Comment: 25 pages; minor editorial revisions made; formal certificates for computer assisted proofs of the paper added to arXi

Optimality conditions for linear copositive programming problems with isolated immobile indices

In the present paper, we apply our recent results on optimality for convex semi-infinite programming to a problem of linear copositive programming (LCP). We prove explicit optimality conditions that use concepts of immobile indices and their immobility orders and do not require the Slater constraint qualification to be satisfied. The only assumption that we impose here is that the set of immobile indices consists of isolated points and hence is finite. This assumption is weaker than the Slater condition; therefore, the optimality conditions obtained in the paper are more general when compared with those usually used in LCP. We present an example of a problem in which the new optimality conditions allow one to test the optimality of a given feasible solution while the known optimality conditions fail to do this. Further, we use the immobile indices to construct a pair of regularized dual copositive problems and show that regardless of whether the Slater condition is satisfied or not, the duality gap between the optimal values of these problems is zero. An example of a problem is presented for which the standard strict duality fails, but the duality gap obtained by using the regularized dual problem vanishes.publishe

The extension problem for partial Boolean structures in Quantum Mechanics

Alternative partial Boolean structures, implicit in the discussion of classical representability of sets of quantum mechanical predictions, are characterized, with definite general conclusions on the equivalence of the approaches going back to Bell and Kochen-Specker. An algebraic approach is presented, allowing for a discussion of partial classical extension, amounting to reduction of the number of contexts, classical representability arising as a special case. As a result, known techniques are generalized and some of the associated computational difficulties overcome. The implications on the discussion of Boole-Bell inequalities are indicated.Comment: A number of misprints have been corrected and some terminology changed in order to avoid possible ambiguitie

Quadrilateral-octagon coordinates for almost normal surfaces

Normal and almost normal surfaces are essential tools for algorithmic 3-manifold topology, but to use them requires exponentially slow enumeration algorithms in a high-dimensional vector space. The quadrilateral coordinates of Tollefson alleviate this problem considerably for normal surfaces, by reducing the dimension of this vector space from 7n to 3n (where n is the complexity of the underlying triangulation). Here we develop an analogous theory for octagonal almost normal surfaces, using quadrilateral and octagon coordinates to reduce this dimension from 10n to 6n. As an application, we show that quadrilateral-octagon coordinates can be used exclusively in the streamlined 3-sphere recognition algorithm of Jaco, Rubinstein and Thompson, reducing experimental running times by factors of thousands. We also introduce joint coordinates, a system with only 3n dimensions for octagonal almost normal surfaces that has appealing geometric properties.Comment: 34 pages, 20 figures; v2: Simplified the proof of Theorem 4.5 using cohomology, plus other minor changes; v3: Minor housekeepin

A Bichromatic Incidence Bound and an Application

We prove a new, tight upper bound on the number of incidences between points and hyperplanes in Euclidean d-space. Given n points, of which k are colored red, there are O_d(m^{2/3}k^{2/3}n^{(d-2)/3} + kn^{d-2} + m) incidences between the k red points and m hyperplanes spanned by all n points provided that m = \Omega(n^{d-2}). For the monochromatic case k = n, this was proved by Agarwal and Aronov. We use this incidence bound to prove that a set of n points, no more than n-k of which lie on any plane or two lines, spans \Omega(nk^2) planes. We also provide an infinite family of counterexamples to a conjecture of Purdy's on the number of hyperplanes spanned by a set of points in dimensions higher than 3, and present new conjectures not subject to the counterexample.Comment: 12 page

Lines, Circles, Planes and Spheres

Let $S$ be a set of $n$ points in $\mathbb{R}^3$, no three collinear and not all coplanar. If at most $n-k$ are coplanar and $n$ is sufficiently large, the total number of planes determined is at least $1 + k \binom{n-k}{2}-\binom{k}{2}(\frac{n-k}{2})$. For similar conditions and sufficiently large $n$, (inspired by the work of P. D. T. A. Elliott in \cite{Ell67}) we also show that the number of spheres determined by $n$ points is at least $1+\binom{n-1}{3}-t_3^{orchard}(n-1)$, and this bound is best possible under its hypothesis. (By $t_3^{orchard}(n)$, we are denoting the maximum number of three-point lines attainable by a configuration of $n$ points, no four collinear, in the plane, i.e., the classic Orchard Problem.) New lower bounds are also given for both lines and circles.Comment: 37 page

A Transfer Matrix for the Backbone Exponent of Two-Dimensional Percolation

Rephrasing the backbone of two-dimensional percolation as a monochromatic path crossing problem, we investigate the latter by a transfer matrix approach. Conformal invariance links the backbone dimension D_b to the highest eigenvalue of the transfer matrix T, and we obtain the result D_b=1.6431 \pm 0.0006. For a strip of width L, T is roughly of size 2^{3^L}, but we manage to reduce it to \sim L!. We find that the value of D_b is stable with respect to inclusion of additional ``blobs'' tangent to the backbone in a finite number of points.Comment: 19 page

Rescaled coordinate descent methods for linear programming

We propose two simple polynomial-time algorithms to find a positive solution to Ax=0Ax=0 . Both algorithms iterate between coordinate descent steps similar to von Neumann’s algorithm, and rescaling steps. In both cases, either the updating step leads to a substantial decrease in the norm, or we can infer that the condition measure is small and rescale in order to improve the geometry. We also show how the algorithms can be extended to find a solution of maximum support for the system Ax=0Ax=0 , x≥0x≥0 . This is an extended abstract. The missing proofs will be provided in the full version

Localization via fractional moments for models on Z\mathbb{Z} with single-site potentials of finite support

One of the fundamental results in the theory of localization for discrete Schr\"odinger operators with random potentials is the exponential decay of Green's function and the absence of continuous spectrum. In this paper we provide a new variant of these results for one-dimensional alloy-type potentials with finitely supported sign-changing single-site potentials using the fractional moment method.Comment: LaTeX-file, 26 pages with 2 LaTeX figure