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.     

Interior penalty method for the indefinite time-harmonic Maxwell equations

Summary In this paper, we introduce and analyze the interior penalty discontinuous Galerkin method for the numerical discretization of the indefinite time-harmonic Maxwell equations in the high-frequency regime. Based on suitable duality arguments, we derive a-priori error bounds in the energy norm and the L²-norm. In particular, the error in the energy norm is shown to converge with the optimal order (h^min{s,}) with respect to the mesh size h, the polynomial degree , and the regularity exponent s of the analytical solution. Under additional regularity assumptions, the L²-error is shown to converge with the optimal order (h⁺¹). The theoretical results are confirmed in a series of numerical experiments.Supported by the EPSRC (Grant GR/R76615).Supported by the Swiss National Science Foundation under project 21-068126.02.Supported in part by the Natural Sciences and Engineering Council of Canada.     

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)     

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.     

Invertible modules for commutative -algebras with residue fields

The aim of this note is to understand under which conditions invertible modules over a commutative -algebra in the sense of Elmendorf, Kriz, Mandell & May give rise to elements in the algebraic Picard group of invertible graded modules over the coefficient ring by taking homotopy groups. If a connective commutative -algebra R has coherent localizations for every maximal ideal , then for every invertible R-module U, U_*=π_*U is an invertible graded R_*-module. In some non-connective cases we can carry the result over under the additional assumption that the commutative -algebra has ‘residue fields’ for all maximal ideals if the global dimension of R_* is small or if R is 2-periodic with underlying Noetherian complete local regular ring R₀. We apply these results to finite abelian Galois extensions of Lubin-Tate spectra.     

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.     

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.     

Rate of convergence of regularization procedures and finite element approximations for the total variation flow

Summary We derive rates of convergence for regularization procedures (characterized by a parameter ) and finite element approximations of the total variation flow, which arises from image processing, geometric analysis and materials sciences. Practically useful error estimates, which depend on only in low polynomial orders, are established for the proposed fully discrete finite element approximations. As a result, scaling laws which relate mesh parameters to the regularization parameter are also obtained. Numerical experiments are provided to validate the theoretical results and show efficiency of the proposed numerical methods.     

A numerical algorithm for the nonlinear Kirchhoff string equation

The initial boundary value problem is considered for the dynamic string equation . Its solution is found by means of an algorithm, the constituent parts of which are the Galerkin method, the modified Crank-Nicolson difference scheme used to perform approximation with respect to spatial and time variables, and also a Picard type iteration process for solving the system of nonlinear equations obtained by discretization. Errors of the three parts of the algorithm are estimated and, as a result, its total error estimate is obtained.     

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     

Fast Runge-Kutta approximation of inhomogeneous parabolic equations

The result after N steps of an implicit Runge-Kutta time discretization of an inhomogeneous linear parabolic differential equation is computed, up to accuracy ɛ, by solving only linear systems of equations. We derive, analyse, and numerically illustrate this fast algorithm.     

Asymptotically exact functional error estimators based on superconvergent gradient recovery

The use of dual/adjoint problems for approximating functionals of solutions of PDEs with great accuracy or to merely drive a goal-oriented adaptive refinement scheme has become well-accepted, and it continues to be an active area of research. The traditional approach involves dual residual weighting (DRW). In this work we present two new functional error estimators and give conditions under which we can expect them to be asymptotically exact. The first is of DRW type and is derived for meshes in which most triangles satisfy an -approximate parallelogram property. The second functional estimator involves dual error estimate weighting (DEW) using any superconvergent gradient recovery technique for the primal and dual solutions. Several experiments are done which demonstrate the asymptotic exactness of a DEW estimator which uses a gradient recovery scheme proposed by Bank and Xu, and the effectiveness of refinement done with respect to the corresponding local error indicators. Resubmitted to Numerische Mathematik, June 30, 2005, with changes suggested by referees.     

Hybrid cross approximation of integral operators

The efficient treatment of dense matrices arising, e.g., from the finite element discretisation of integral operators requires special compression techniques. In this article we use the -matrix representation that approximates the dense stiffness matrix in admissible blocks (corresponding to subdomains where the underlying kernel function is smooth) by low-rank matrices. The low-rank matrices are assembled by a new hybrid algorithm (HCA) that has the same proven convergence as standard interpolation but also the same efficiency as the (heuristic) adaptive cross approximation (ACA).     

On a class of weak Tchebycheff systems

In this paper we study the approximation power, the existence of a normalized B-basis and the structure of a degree-raising process for spaces of the formrequiring suitable assumptions on the functions u and v. The results about degree raising are detailed for special spaces of this form which have been recently introduced in the area of CAGD.     

A schlichtness theorem for envelopes of holomorphy

Let Ω be a domain in . We prove the following theorem. If the envelope of holomorphy of Ω is schlicht over Ω, then the envelope is in fact schlicht. We provide examples showing that the conclusion of the theorem does not hold in , n>2. Additionally, we show that the theorem cannot be generalized to provide information about domains in whose envelopes are multiply sheeted.     

Besov regularity and new error estimates for finite volume approximations of the p-laplacian

Summary. In [1], we have constructed a family of finite volume schemes on rectangular meshes for the p-laplacian and we proved error estimates in case the exact solution lies in W^2,p. Actually, W^2,p is not a natural space for solutions of the p-laplacian in the case p>2. Indeed, for general L^p data it can be shown that the solution only belongs to the Besov space In this paper, we prove Besov kind a priori estimates on the approximate solution for any data in L^p. We then obtain new error estimates for such solutions in the case of uniform meshes     

Holomorphic dynamics near germs of singular curves

Let M be a two dimensional complex manifold, p ∈ M and a germ of holomorphic foliation of M at p. Let be a germ of an irreducible, possibly singular, curve at p in M which is a separatrix for . We prove that if the Camacho-Sad-Suwa index Ind then there exists another separatrix for at p. A similar result is proved for the existence of parabolic curves for germs of holomorphic diffeomorphisms near a curve of fixed points.     

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 an eigenvalue problem related to the critical exponent

In this paper we study the eigenvalue problemwhere is a smooth bounded domain, and u is a positive solution of the problemsuch thatwhere S is the best Sobolev constant for the embedding of H¹₀() into L²*(), We prove several estimates for the eigenvalues _i, of (I), i=2,..,N+2 and some qualitative properties of the corresponding eigenfunctions.Supported by M.I.U.R., project Variational methods and nonlinear differential equations.     

Interior regularity criteria in weak spaces for the Navier-Stokes equations

We study the interior regularity of weak solutions of the incompressible Navier-Stokes equations in ×(0,T), where and 0<T<. The local boundedness of a weak solution u is proved under the assumption that is sufficiently small for some (r,s) with and 3r<. Our result extends the well-known criteria of Serrin (1962), Struwe (1988) and Takahashi (1990) to the weak space-time spaces.Mathematics Subject Classification (2000): 35Q30, 76N10