Orthonormal rational function vectors

Summary In this paper, we develop a matrix framework to solve the problem of finding orthonormal rational function vectors with prescribed poles with respect to a certain discrete inner product that is defined by a set of data points and corresponding weight vectors w_i,j. Our algorithm for solving the problem is recursive, and it is of complexity If all data points are real or lie on the unit circle, then the complexity is reduced by an order of magnitude.     

Fast and stable QR eigenvalue algorithms for generalized companion matrices and secular equations

Summary We introduce a class of n×n structured matrices which includes three well-known classes of generalized companion matrices: tridiagonal plus rank-one matrices (comrade matrices), diagonal plus rank-one matrices and arrowhead matrices. Relying on the structure properties of , we show that if A then A=RQ , where A=QR is the QR decomposition of A. This allows one to implement the QR iteration for computing the eigenvalues and the eigenvectors of any A with O(n) arithmetic operations per iteration and with O(n) memory storage. This iteration, applied to generalized companion matrices, provides new O(n²) flops algorithms for computing polynomial zeros and for solving the associated (rational) secular equations. Numerical experiments confirm the effectiveness and the robustness of our approach.The results of this paper were presented at the Workshop on Nonlinear Approximations in Numerical Analysis, June 22 – 25, 2003, Moscow, Russia, at the Workshop on Operator Theory and Applications (IWOTA), June 24 – 27, 2003, Cagliari, Italy, at the Workshop on Numerical Linear Algebra at Universidad Carlos III in Leganes, June 16 – 17, 2003, Leganes, Spain, at the SIAM Conference on Applied Linear Algebra, July 15 – 19, 2003, Williamsburg, VA and in the Technical Report [8]. This work was partially supported by MIUR, grant number 2002014121, and by GNCS-INDAM. This work was supported by NSF Grant CCR 9732206 and PSC CUNY Awards 66406-0033 and 65393-0034.     

Multi-level spectral galerkin method for the navier-stokes problem I : spatial discretization

A multi-level spectral Galerkin method for the two-dimensional non-stationary Navier-Stokes equations is presented. The method proposed here is a multiscale method in which the fully nonlinear Navier-Stokes equations are solved only on a low-dimensional space subsequent approximations are generated on a succession of higher-dimensional spaces j=2, . . . ,J, by solving a linearized Navier-Stokes problem around the solution on the previous level. Error estimates depending on the kinematic viscosity 0<ν<1 are also presented for the J-level spectral Galerkin method. The optimal accuracy is achieved when We demonstrate theoretically that the J-level spectral Galerkin method is much more efficient than the standard one-level spectral Galerkin method on the highest-dimensional space . The work of this author was supported in part by the NSF of China 10371095, City University of Hong Kong Research Project 7001093 Hong Kong and the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. CityU 1084/02P)     

A branch-and-cut algorithm for the stochastic uncapacitated lot-sizing problem

This paper addresses a multi-stage stochastic integer programming formulation of the uncapacitated lot-sizing problem under uncertainty. We show that the classical (ℓ,S) inequalities for the deterministic lot-sizing polytope are also valid for the stochastic lot-sizing polytope. We then extend the (ℓ,S) inequalities to a general class of valid inequalities, called the inequalities, and we establish necessary and sufficient conditions which guarantee that the inequalities are facet-defining. A separation heuristic for inequalities is developed and incorporated into a branch-and-cut algorithm. A computational study verifies the usefulness of the inequalities as cuts. This research has been supported in part by the National Science Foundation under Award number DMII-0121495.     

Preconditioning Landweber iteration in Hilbert scales

In this paper we investigate convergence of Landweber iteration in Hilbert scales for linear and nonlinear inverse problems. As opposed to the usual application of Hilbert scales in the framework of regularization methods, we focus here on the case s≤0, which (for Tikhonov regularization) corresponds to regularization in a weaker norm. In this case, the Hilbert scale operator L^−2s appearing in the iteration acts as a preconditioner, which significantly reduces the number of iterations needed to match an appropriate stopping criterion. Additionally, we carry out our analysis under significantly relaxed conditions, i.e., we only require instead of which is the usual condition for regularization in Hilbert scales. The assumptions needed for our analysis are verified for several examples and numerical results are presented illustrating the theoretical ones. supported by the Austrian Science Foundation (FWF) under grant SFB/F013     

The Rees algebra of a positive normal affine semigroup ring

Let R be a positive normal affine semigroup ring of dimension d and let be the maximal homogeneous ideal of R. We show that the integral closure of is equal to for all n ∈ℕ with n ≥ d − 2. From this we derive that the Rees algebra R[t] is normal in case that d ≤ 3. If emb dim(R) = d + 1, we can give a necessary and sufficient condition for R[t] to be normal.     

On Cross-intersecting Families of Sets

A family of -element subsets and a family of k-element subsets of an n-element set are cross-intersecting if every set from has a nonempty intersection with every set from . We compare two previously established inequalities each related to the maximization of the product , and give a new and short proof for one of them. We also determine the maximum of for arbitrary positive weights ,_k.     

On topological Kadec norms

We present a topological analogue of the classic Kadec Renorming Theorem, as follows. Let be two separable metric topologies on the same set X. We prove that every point in X has an -neighbourhood basis consisting of sets that are -closed if and only if there exists a function φ: X→ℝ that is -lower semi-continuous and such that is the weakest topology on X that contains and that makes φ continuous. An immediate corollary is that the class of almost n-dimensional spaces consists precisely of the graphs of lower semi-continuous functions with at most n-dimensional domains.     

Summing norms of identities between unitary ideals

Based on abstract interpolation, we prove asymptotic formulae for the (F,2)-summing norm of inclusions id: , where E and F are two Banach sequence spaces. Here, stands for the unitary ideal of operators on the n-dimensional Hilbert space whose singular values belong to E, and for the Hilbert-Schmidt operators. Our results are noncommutative analogues of results due to Bennett and Carl, as well as their recent generalizations to Banach sequence spaces. As an application, we give lower and upper estimates for certain s-numbers of the embeddings id: and id: . In the concluding section, we finally consider mixing norms. The second named author was supported by KBN Grant 2 P03A 042 18.     

On the exact Hausdorff dimension of the set of Liouville numbers. II

Let denote the set of Liouville numbers. For a dimension function h, we write for the h-dimensional Hausdorff measure of . In previous work, the exact ``cut-point' at which the Hausdorff measure of drops from infinity to zero has been located for various classes of dimension functions h satisfying certain rather restrictive growth conditions. In the paper, we locate the exact ``cut-point' at which the Hausdorff measure of drops from infinity to zero for all dimension functions h. Namely, if h is a dimension function for which the function increases faster than any power function near 0, then , and if h is a dimension function for which the function increases slower than some power function near 0, then . This provides a complete characterization of all Hausdorff measures of without assuming anything about the dimension function h, and answers a question asked by R. D. Mauldin. We also show that if then does not have σ-finite measure. This answers another question asked by R. D. Mauldin. This work was done while Dave L. Renfro was at the Department of Mathematics at Central Michigan University.     

On the approximation order of extremal point methods for hyperbolic minimal energy problems

Summary. Let be an analytic Jordan curve in the unit disk We regard the hyperbolic minimal energy problem where () denotes the set of all probability measures on . There exist several extremal point discretizations of *, among others introduced by M. Tsuji (Tsuji points) or by K. Menke (hyperbolic Menke points). In the present article, it is proven that hyperbolic Menke points approach the images of roots of unity under a conformal map from onto geometrically fast if the number of points tends to infinity. This establishes a conjecture of K. Menke. In particular, explicit bounds for the approximation error are given. Finally, an effective method for the numerical determination of * providing a geometrically shrinking error bound is presented.Mathematics Subject Classification (1991): 30C85, 30E10, 31C20The notation Menke points has been introduced by D. Gaier.     

On the zero-distribution of Epstein zeta-functions

We prove that the mean value of the real parts of the nontrivial zeros of the Epstein zeta-function associated with a positive definite quadratic form in n variables is equal to . Furthermore, we show that Epstein zeta-functions in general have an asymmetric zero-distribution with respect to the critical line Re .     

Hausdorff measures, capacities and compact composition operators

It is shown that there exist analytic self-maps ϕ of the unit disc inducing compact composition operators on the Hardy space , 1 ≤ p < ∞ such that the Hausdorff dimension of the set is one; sharpening a classical result due to Schwartz. Moreover, the same holds in the weighted Dirichlet spaces with 0 < α < 1. As a consequence, we deduce that there exist symbols ϕ inducing compact composition operators on such that the α-capacity of E_ϕ is positive, which is no longer true for those just inducing Hilbert-Schmidt composition operators on . First author is partially supported by Plan Nacional I+D grant no. BFM2003-00034, and Gobierno de Aragón research group Análisis Matemático y Aplicaciones, ref. DGA E-64 . Second author is partially supported by Plan Nacional I+D grant no. BFM2002-00571 and Junta de Andalucía RNM-314.     

New shape functions for triangular p-FEM using integrated Jacobi polynomials

In this paper, the second order boundary value problem −∇·((x,y)∇u)=f is discretized by the Finite Element Method using piecewise polynomial functions of degree p on a triangular mesh. On the reference element, we define integrated Jacobi polynomials as interior ansatz functions. If is a constant function on each triangle and each triangle has straight edges, we prove that the element stiffness matrix has not more than nonzero matrix entries. An application for preconditioning is given. Numerical examples show the advantages of the proposed basis.     

On the rank of abelian varieties over function fields

Let be a smooth projective curve defined over a number field k, A/k() an abelian variety and (τ, B) the k()/k-trace of A. We estimate how the rank of A(k())/τB(k) varies when we take a finite geometrically abelian cover defined over k. This work was partially supported by CNPq research grant 304424/2003-0, Pronex 41.96.0830.00 and CNPq Edital Universal 470099/2003-8. I would like to thank Douglas Ulmer for comments on how to treat the case of arbitrary ramification, but the conductor prime to the ramification locus, in the case of elliptic fibrations. I would also like to thank Marc Hindry for comments on the inequality comparing the conductors of A and A'. Finally, I also thank the referee for his comments and criticisms.     

Estimates for stable hypersurfaces of prescribed F-mean curvature

We consider immersed hypersurfaces :Mⁿ→ℝⁿ⁺¹ with prescribed anisotropic mean curvature . Such hypersurfaces can be characterized as critical points of parametric functionals of the type with an elliptic Lagrangian F depending on normal directions and a smooth vectorfield Q satisfying . We establish curvature estimates for stable hypersurfaces of dimension n≤5, provided F is C³-close to the area integrand.     

2-universal -lattices over real quadratic fields

Let be a real quadratic field with m a square-free positive rational integer, and be the ring of integers in F. An -lattice L on a totally positive definite quadratic space V over F is called r-universal if L represents all totally positive definite -lattices l with rank r over . We prove that there exists no 2-universal -lattice over F with rank less than 6, and there exists a 2-universal -lattice over F with rank 6 if and only if m=2, 5. Moreover there exists only one 2-universal -lattice with rank 6, up to isometry, over .     

Mean-value theorems for multiplicative functions on Beurling’s generalized integers (I)

Halász’s general mean-value theorem for multiplicative functions on ℕ is classical in probabilistic number theory. We extend this theorem to functions f, defined on a set of generalized integers associated with a set of generalized primes in Beurling’s sense, which satisfies Halász’s conditions, in particular,Assume that the distribution function N(x) of satisfieswith γ>γ₀, where ρ₁<ρ₂<···<ρ_m are constants with ρ_m≥1 and A₁,···,A_m are real constants with A_m>0. Also, assume that the Chebyshev function ψ(x) of satisfieswith M>M₀. Then the asymptoticimplieswhere τ is a positive constant with τ≥1 and L(u) is a slowly oscillating function with |L(u)|=1.     

Fast convolution with radial kernels at nonequispaced knots

Summary. We develop a new algorithm for the fast evaluation of linear combinations of radial functions based on the recently developed fast Fourier transform at nonequispaced knots. For smooth kernels, e.g. the Gaussian, our algorithm requires arithmetic operations. In case of singular kernels an additional regularization procedure must be incorporated and the algorithm has the arithmetic complexity if either the points y_j or the points x_k are reasonably uniformly distributed. We prove error estimates to obtain clues about the choice of the involved parameters and present numerical examples for various singular and smooth kernels in two dimensions.Mathematics Subject Classification (2000): 65T40, 65T50, 65F30Revised version received December 3, 2003     

Partial regularity for stationary harmonic maps at a free boundary

Let and be smooth Riemannian manifolds, of the dimension n≥2 with nonempty boundary, and compact without boundary. We consider stationary harmonic maps u ∈ H¹(, ) with a free boundary condition of the type u(∂) ⊂ Γ, given a submanifold Γ⊂. We prove partial boundary regularity, namely (sing(u))=0, a result that was until now only known in the interior of the domain (see [B]). The key of the proof is a new lemma that allows an extension of u by a reflection construction. Once the partial regularity theorem is known, it is possible to reduce the dimension of the singular set further under additional assumptions on the target manifold and the submanifold Γ.