Predictor-corrector algorithm for solvingP *(κ)-matrix LCP from arbitrary positive starting points

A new predictor-corrector algorithm is proposed for solvingP _*(κ)-matrix linear complementarity problems. If the problem is solvable, then the algorithm converges from an arbitrary positive starting point (x ⁰,s ⁰). The computational complexity of the algorithm depends on the quality of the starting point. If the starting point is feasible or close to being feasible, it has -iteration complexity, whereρ ₀ is the ratio of the smallest and average coordinate ofX ⁰ s ⁰. With appropriate initialization, a modified version of the algorithm terminates in O((1+κ)²(n/ρ ₀)L) steps either by finding a solution or by determining that the problem has no solution in a predetermined, arbitrarily large, region. The algorithm is quadratically convergent for problems having a strictly complementary solution. We also propose an extension of a recent algorithm of Mizuno toP _*(κ)-matrix linear complementarity problems such that it can start from arbitrary positive points and has superlinear convergence without a strictly complementary condition. The work of this author was supported in part by NSF, Grant DMS 9305760 and by an Oberman fellowship from the University of Iowa Center for Advanced Studies.     

Existence and Construction of H 1-Splines of Class C k on a Three Directional Mesh

Let be the uniform triangulation generated by the usual three directional mesh of the plane and let H ₁ be the regular hexagon formed by the six triangles of surrounding the origin. We study the space of piecewise polynomial functions in C ^k (R ²) with support H ₁ having a sufficiently high degree n, which are invariant with respect to the group of symmetries of H ₁ and whose sum of integer translates is constant. Such splines are called H ₁-splines. We first compute the dimension of this space in function of n and k. Then we prove the existence of a unique H ₁-spline of minimal degree for any fixed k0. Finally, we describe an algorithm computing the Bernstein–Bézier coefficients of this spline.     

Linear best approximation using a class ofk-majorl p norms

Consideration is given to problems of linear best approximation using a variant of the usuall _p norms referred to ask-majorl _p norms, for the cases when 1<p<. The underlying problem is the minimization of a piecewise smooth function. Best approximations are characterized, and a descent algorithm is developed.     

An infeasible-interior-point predictor-corrector algorithm for theP *-geometric LCP

AP _*-geometric linear complementarity problem (P _*GP) as a generalization of the monotone geometric linear complementarity problem is introduced. In particular, it contains the monotone standard linear complementarity problem and the horizontal linear complementarity problem. Linear and quadratic programming problems can be expressed in a “natural” way (i.e., without any change of variables) asP _*GP. It is shown that the algorithm of Mizunoet al. [6] can be extended to solve theP _*GP. The extended algorithm is globally convergent and its computational complexity depends on the quality of the starting points. The algorithm is quadratically convergent for problems having a strictly complementary solution. The work of F. A. Potra was supported in part by NSF Grant DMS 9305760     

<Emphasis Type="Italic">L</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript>-mixed intersection bodies

In this paper the author first introduce a new concept of L _p -dual mixed volumes of star bodies which extends the classical dual mixed volumes. Moreover, we extend the notions of L _p intersection body to L _p -mixed intersection body. Inequalities for L _p -dual mixed volumes of L _p -mixed intersection bodies are established and the results established here provide new estimates for these type of inequalities. This work was supported by the Natural Science Foundation of Zhejiang Province of China (Grant No. Y605065) and the Foundation of the Education Department of Zhejiang Province of China (Grant No. 20050392)     

Solving ℓ1 Regularization Problems With Piecewise Linear Losses

This article is concerned with the computational aspect of ?₁ regularization problems with a certain class of piecewise linear loss functions. The problem of computing the ?₁ regularization path for a piecewise linear loss can be formalized as a parametric linear programming problem. We propose an efficient implementation method of the parametric simplex algorithm for such a problem. We also conduct a simulation study to investigate the behavior of the number of “breakpoints” of the regularization path when both the number of observations and the number of explanatory variables vary. Our method is also applicable to the computation of the regularization path for a piecewise linear loss and the blockwise ?_∞ penalty. This article has supplementary material online.     

Exact L 2-norm plane separation

We consider the problem of separating two sets of points in an n-dimensional real space with a (hyper)plane that minimizes the sum of L _p -norm distances to the plane of points lying on the wrong side of it. Despite recent progress, practical techniques for the exact solution of cases other than the L ₁ and L _∞-norm were unavailable. We propose and implement a new approach, based on non-convex quadratic programming, for the exact solution of the L ₂-norm case. We solve in reasonable computing times artificial problems of up to 20000 points (in 6 dimensions) and 13 dimensions (with 2000 points). We also observe that, for difficult real-life instances from the UCI Repository, computation times are substantially reduced by incorporating heuristic results in the exact solution process. Finally, we compare the classification performance of the planes obtained for the L ₁, L ₂ and L _∞ formulations. It appears that, despite the fact that L ₂ formulation is computationally more expensive, it does not give significantly better results than the L ₁ and L _∞ formulations.     

Embedding into l ∞2 Is Easy, Embedding into l ∞3 Is?NP-Complete

We give a new algorithm for enumerating all possible embeddings of a metric space (i.e., the distances between every pair within a set of n points) into ℝ² Cartesian space preserving their l _∞ (or l ₁) metric distances. Its expected time is (i.e., within a poly-log of the size of the input) beating the previous algorithm. In contrast, we prove that detecting l _∞³ embeddings is NP-complete. The problem is also NP-complete within l ₁² or l _∞² with the added constraint that the locations of two of the points are given or alternatively that the two dimensions are curved into a three-dimensional sphere. We also refute a compaction theorem by giving a metric space that cannot be embedded in l _∞³; however, it can be embedded if any single point is removed. This research is partially supported by NSERC grants. I would like to thank Steven Watson for his extensive help on this paper.     

The L p -curvature images of convex bodies and L p -projection bodies

Associated with the L _p -curvature image defined by Lutwak, some inequalities for extended mixed p-affine surface areas of convex bodies and the support functions of L _p -projection bodies are established. As a natural extension of a result due to Lutwak, an L _p -type affine isoperimetric inequality, whose special cases are L _p -Busemann-Petty centroid inequality and L _p -affine projection inequality, respectively, is established. Some L _p -mixed volume inequalities involving L _p -projection bodies are also established.     

The decomposition of Kv into K2×K5''''s

The decomposition of the complete graph Kv into Kr×Kc's, the products of Kr and Kc,is originated from the use of DNA library screening. In this paper, we consider the case where r=2 and c = 5, and show that such a decomposition exists if and only if v ≡ 1 (mod 25).     

A new multivariable <Subscript><Emphasis Type="BoldItalic">6</Emphasis></Subscript><Emphasis Type="BoldItalic">ψ</Emphasis><Subscript><Emphasis Type="BoldItalic">6</Emphasis></Subscript> summation formula

By multidimensional matrix inversion, combined with an A _r extension of Jackson’s ₈ φ ₇ summation formula by Milne, a new multivariable ₈ φ ₇ summation is derived. By a polynomial argument this ₈ φ ₇ summation is transformed to another multivariable ₈ φ ₇ summation which, by taking a suitable limit, is reduced to a new multivariable extension of the nonterminating ₆ φ ₅ summation. The latter is then extended, by analytic continuation, to a new multivariable extension of Bailey’s very-well-poised ₆ ψ ₆ summation formula. Partly supported by FWF Austrian Science Fund grants P17563-N13, and S9607 (the second is part of the Austrian National Research Network “Analytic Combinatorics and Probabilistic Number Theory”).     

Global Linear and Quadratic One-step Smoothing Newton Method for P 0-LCP

We propose a new smoothing Newton method for solving the P ₀-matrix linear complementarity problem (P ₀-LCP) based on CHKS smoothing function. Our algorithm solves only one linear system of equations and performs only one line search per iteration. It is shown to converge to a P ₀-LCP solution globally linearly and locally quadratically without the strict complementarity assumption at the solution. To the best of author's knowledge, this is the first one-step smoothing Newton method to possess both global linear and local quadratic convergence. Preliminary numerical results indicate that the proposed algorithm is promising.     

Some complete metrics on spaces of fuzzy subsets

The classes of the L_p,∞- and L_p-metrics play an important role to develop a probability theory in fuzzy sample spaces. All of these metrics are known to be separable, but not complete. The classes are closely related as for each L_p,∞-metric there exists some L_p-metric which induces the same topology. This paper deals with the completion of the L_p,∞- and L_p-metrics. We can also show that the relationship between the classes of L_p,∞- and L_p-metrics still holds for the obtained respective classes of their completions.     

The Smallest Degree Sum That Yields Potentially Kr＋1 - K3-Graphic Sequences

Let a（Kr,＋1 - K3,n） be the smallest even integer such that each n-term graphic sequence п= （d1,d2,…dn） with term sum σ（п） = d1 ＋ d2 ＋…＋ dn 〉 σ（Kr＋1 -K3,n） has a realization containing Kr＋1 - K3 as a subgraph, where Kr＋1 -K3 is a graph obtained from a complete graph Kr＋1 by deleting three edges which form a triangle. In this paper, we determine the value σ（Kr＋1 - K3,n） for r ≥ 3 and n ≥ 3r＋ 5.     

<Emphasis Type="Italic">∂</Emphasis>-operations and <Emphasis Type="Italic">H</Emphasis><Subscript><Emphasis Type="Italic">V</Emphasis></Subscript>-fields

The hyperoperations, called theta-operations （δ）, are motivated from the usual property, which the derivative has on the derivation of a product of functions. Using any map on a set, one can define δ-operations. In this paper, we continue our study on the δ-operations on groupoids, rings, fields and vector spaces or on the corresponding hyperstructures. Using δ-operations one obtains, mainly, Hwstructures, which form the largest class of the hyperstructures. For representation theory of hyperstructures, by hypermatrices, one needs special Hv-rings or Hy-fields, so these hyperstructures can be used. Moreover, we study the relation of these δ-structures with other classes of hyperstructures, especially with the Hv-structures.     

The complexity of anl p method for discrete Tchebycheff approximation in exact arithmetic

A discrete Tchebycheff problem is approximated by a sequence ofl _p norm problems. This is an algorithm to be used if a large number of variables is involved as is the case in the approximation of a function of several variables by a polynomial. The complexity of this procedure is investigated and a lower bound for the number of steps to reach ε-optimality is established. Supported by NIH Grant RR01243 at the University of Washington     

Localization of <Emphasis Type="Italic">LM</Emphasis><Subscript><Emphasis Type="Italic">n</Emphasis></Subscript>-algebras

The aim of this paper is to define the localization LM _n -algebra of an LM _n —algebra L with respect to a topology F on L; in Section 5 we prove that the maximal LM _n -algebra of fractions (defined in [3]) and the LM _n -algebra of fractions relative to an Λ—closed system (defined in Section 2) are LM _n -algebras of localization.     

F_2+vF_2上的循环码

本文研究了环F_2+vF_2上的循环码.利用标准形生成元集刻画了环F_2+vF_2上的循环码的代数结构,证明了环F_2+vF_2上的每一个非零的循环码均有唯一的标准形生成元集,进而得到了每一个循环码均是由一个多项式生成的.     

Weak Equality and Exact Sequences in H v -modules

The largest class of multivalued systems satisfying the module-like axioms is the H_v-module. H_v-modules first were introduced by Vougiouklis. In this paper we define weak equality between two subsets of an H_v-module and introduced the notion of exact sequences of H_v-modules. Also some results on the weak equality and exact sequences are given.     

Absolute summability of infinite series

It is shown in [4] that if a normal matrix,A satisfies some conditions then |C,1| _k summability implies |A| _k summability wherek≥1. In the present paper, we consider the converse implication.