On the augmented system approach to sparse least-squares problems

Summary We study the augmented system approach for the solution of sparse linear least-squares problems. It is well known that this method has better numerical properties than the method based on the normal equations. We use recent work by Arioli et al. (1988) to introduce error bounds and estimates for the components of the solution of the augmented system. In particular, we find that, using iterative refinement, we obtain a very robust algorithm and our estimates of the error are accurate and cheap to compute. The final error and all our error estimates are much better than the classical or Skeel's error analysis (1979) indicates. Moreover, we prove that our error estimates are independent of the row scaling of the augmented system and we analyze the influence of the Björck scaling (1967) on these estimates. We illustrate this with runs both on large-scale practical problems and contrived examples, comparing the numerical behaviour of the augmented systems approach with a code using the normal equations. These experiments show that while the augmented system approach with iterative refinement can sometimes be less efficient than the normal equations approach, it is comparable or better when the least-squares matrix has a full row, and is, in any case, much more stable and robust.This author was visiting Harwell and was funded by a grant from the Italian National Council of Research (CNR), Istituto di Elaborazione dell'Informazione-CNR, via S. Maria 46, I-56100 Pisa, ItalyThis author was visiting Harwell from Faculty of Mathematics and Computer Science of the University of Amsterdam     

增算子的不动点和耦合不动点

Under some general continuous and compact conditions, the existence problems of fiked points andd coupled fixed points for increasing operators are studied. an application, we utilize the results obtained to study the existence of solutions for differential inclusions in Banach spaces.     

Multiderivations of Coxeter arrangements

Let V be an ℓ-dimensional Euclidean space. Let G⊂O(V) be a finite irreducible orthogonal reflection group. Let ? be the corresponding Coxeter arrangement. Let S be the algebra of polynomial functions on V. For H∈? choose α _H ∈V ^* such that H=ker(α _H ). For each nonnegative integer m, define the derivation module D ^(m) (?)={θ∈Der _S |θ(α _H )∈Sα ^m _H }. The module is known to be a free S-module of rank ℓ by K. Saito (1975) for m=1 and L. Solomon-H. Terao (1998) for m=2. The main result of this paper is that this is the case for all m. Moreover we explicitly construct a basis for D ^(m) (?). Their degrees are all equal to mh/2 (when m is even) or are equal to ((m−1)h/2)+m _i (1≤i≤ℓ) (when m is odd). Here m ₁≤···≤m _ℓ are the exponents of G and h=m _ℓ+1 is the Coxeter number. The construction heavily uses the primitive derivation D which plays a central role in the theory of flat generators by K. Saito (or equivalently the Frobenius manifold structure for the orbit space of G). Some new results concerning the primitive derivation D are obtained in the course of proof of the main result. Oblatum 27-XI-2001 & 4-XII-2001?Published online: 18 February 2002     

The two-dimensional quantum Euclidean algebra

The algebra dual to Woronowicz's deformation of the two-dimensional Euclidean group is constructed. The same algebra is obtained from SU _q (2) via contraction on both the group and algebra levels.     

Properties of holomorphic Wiener functions —skeleton,contraction, and local Taylor expansion

Summary We show that each holomorphic Wiener function has a skeleton which is intrinsic from several viewpoints. In particular, we study the topological aspects of the skeletons by using the local Taylor expansion for holomorphic Wiener functions.Supported in part by the Grant-in-Aid for Science Research 03740120 Min. Education     

Caratheodory interpolation kernels

Kernels over the unit disk for which a version of Carathéodory interpolation is true are characterized in a simple computationally verifiable manner.     

Two Hadamard matrices of order 956 of Goethals-Seidel type

No Hadamard matrices of order 956 appear in the literature. In this note we construct two such matrices. They are both of Goethals-Seidel type. As a consequence of this and of a result of Seberry and Yamada, we conclude that there exist 4-Williamson type matrices of order 1913.This work was supported by NSERC of Canada Grants A-5285 and EQP0107539.     

Graphical operations on projective spaces

Motivated by the work of Crapo and Rota [6] on the lifting of a projective complex, we introduce a class of invariant operations associated to integral-weighted graphs, which we call graphical operations. Such operations generalize the sixth harmonic of a quadranguler set on a projective line. We determine the expansion of the graphical operations in terms of multi-linear bracket polynomials in a Grassmann-Cayley algebra. Reducibility and compositions of such invariant operations are also investigated with a number of examples.Supported by Courant Instructorship, New York University.     

The dual Brunn-Minkowski theory for bounded Borel sets: Dual affine quermassintegrals and inequalities

This paper develops a significant extension of E. Lutwak's dual Brunn-Minkowski theory, originally applicable only to star-shaped sets, to the class of bounded Borel sets. The focus is on expressions and inequalities involving chord-power integrals, random simplex integrals, and dual affine quermassintegrals. New inequalities obtained include those of isoperimetric and Brunn-Minkowski type. A new generalization of the well-known Busemann intersection inequality is also proved. Particular attention is given to precise equality conditions, which require results stating that a bounded Borel set, almost all of whose sections of a fixed dimension are essentially convex, is itself essentially convex.