Rational approximation to trigonometric operators

We consider the approximation of trigonometric operator functions that arise in the numerical solution of wave equations by trigonometric integrators. It is well known that Krylov subspace methods for matrix functions without exponential decay show superlinear convergence behavior if the number of steps is larger than the norm of the operator. Thus, Krylov approximations may fail to converge for unbounded operators. In this paper, we propose and analyze a rational Krylov subspace method which converges not only for finite element or finite difference approximations to differential operators but even for abstract, unbounded operators. In contrast to standard Krylov methods, the convergence will be independent of the norm of the operator and thus of its spatial discretization. We will discuss efficient implementations for finite element discretizations and illustrate our analysis with numerical experiments. AMS subject classification (2000)  65F10, 65L60, 65M60, 65N22     

Extreme eigenvalues of real symmetric Toeplitz matrices

We exploit the even and odd spectrum of real symmetric Toeplitz matrices for the computation of their extreme eigenvalues, which are obtained as the solutions of spectral, or secular, equations. We also present a concise convergence analysis for a method to solve these spectral equations, along with an efficient stopping rule, an error analysis, and extensive numerical results.

     

Characterization of standard embeddings between complex Grassmannians by means of varieties of minimal rational tangents

In 1993,Tsal proved that a proper holomorphic mapping f:Ω→Ω' from an irreducible bounded symmetric domainΩof rank≥2 into a bounded symmetric domainΩ' is necessarily totally geodesic provided that r':=rank(Ω')≤rank(Ω):= r,proving a conjecture of the author's motivated by Hermitian metric rigidity.As a first step in the proof,Tsai showed that df preserves almost everywhere the set of tangent vectors of rank 1.Identifying bounded symmetric domains as open subsets of their compact duals by means of the Borel embedding,this means that the germ of f at a general point preserves the varieties of minimal rational tangents(VMRTs). In another completely different direction Hwang-Mok established with very few exceptions the Cartan- Fubini extension priniciple for germs of local biholomorphisms between Fano manifolds of Picard num- ber 1,showing that the germ of map extends to a global biholomorphism provided that it preserves VMRTs.We propose to isolate the problem of characterization of special holomorphic embeddings between Fano manifolds of Picard number 1,especially in the case of classical manifolds such as ratio- nal homogeneous spaces of Picard number 1,by a non-equidimensional analogue of the Cartan-Fubini extension principle.As an illustration we show along this line that standard embeddings between com- plex Grassmann manifolds of rank≤2 can be characterized by the VMRT-preserving property and a non-degeneracy condition,giving a new proof of a result of Neretin's which on the one hand paves the way for far-reaching generalizations to the context of rational homogeneous spaces and more generally Fano manifolds of Picard number 1,on the other hand should be applicable to the study of proper holomorphic mappings between bounded domains carrying some form of geometric structures.     

High accuracy time and space transform method for advection–diffusion equation in an unbounded domain

A new numerical method called high accuracy time and space transform method (TSTM) is introduced to solve the advection–diffusion equation in an unbounded domain. By a spatial transform, the advection–diffusion equation in the unbounded domain Rⁿ is converted to one on the bounded domain [?1, 1]ⁿ, and the Laplace transform is applied to eliminate time dependency. The consequent boundary value problem is solved by collocation on Chebyshev points. To face the well‐known computational challenge represented by the numerical inversion of the Laplace transform, Talbot's method is applied, consisting of numerically integrating the Bromwich integral on a special contour by means of trapezoidal or midpoint rules. Numerical experiments illustrate that TSTM has exponential rate in time and space. Copyright © 2008 John Wiley & Sons, Ltd.     

Spectral methods for orthogonal rational functions

We present an operator theoretic approach to orthogonal rational functions based on the identification of a suitable matrix representation of the multiplication operator associated with the corresponding orthogonality measure. Two alternatives are discussed, leading to representations which are linear fractional transformations with matrix coefficients acting on infinite Hessenberg or five-diagonal unitary matrices. This approach permits us to recover the orthogonality measure throughout the spectral analysis of an infinite matrix depending uniquely on the poles and the parameters of the recurrence relation for the orthogonal rational functions. Besides, the zeros of the orthogonal and para-orthogonal rational functions are identified as the eigenvalues of matrix linear fractional transformations of finite Hessenberg or five-diagonal matrices. As an application we use operator perturbation theory results to obtain new relations between the support of the orthogonality measure and the location of the poles and parameters of the recurrence relation for the orthogonal rational functions.     

Points of Bounded Height on Equivariant Compactifications of Vector Groups, I

We prove asymptotic formulas for the number of rational points of bounded height on certain equivariant compactifications of the affine plane.     

一类Radial基插值的保凸基

本文为 Radial 插值找到了一组保凸基.     

Toric Rigid Spaces

This paper gives a method to construct rigid spaces, which is similar to the method used to construct toric schemes.