An algorithm for computing Jordan chains and inverting analytic matrix functions

We present an efficient algorithm for obtaining a canonical system of Jordan chains for an n × n regular analytic matrix function A(λ) that is singular at the origin. For any analytic vector function b(λ), we show that each term in the Laurent expansion of A(λ)⁻¹b(λ) may be obtained from the previous terms by solving an (n + d) × (n+d) linear system, where d is the order of the zero of det A(λ) at λ = 0. The matrix representing this linear system contains A(0) as a principal submatrix, which can be useful if A(0) is sparse. The last several iterations can be eliminated if left Jordan chains are computed in addition to right Jordan chains. The performance of the algorithm in floating point and exact (rational) arithmetic is reported for several test cases. The method is shown to be forward stable in floating point arithmetic.     

Banach-stone theorems and separating maps

LetC(X,E) andC(Y,F) denote the spaces of continuous functions on the Tihonov spacesX andY, taking values in the Banach spacesE andF, respectively. A linear mapH:C(X,E)→C(Y,F) isseparating iff(x)g(x)=0 for allx inX impliesHf(y)Hg(y)=0 for ally inY. Some automatic continuity properties and Banach-Stone type theorems (i.e., asserting that isometries must be of a certain form) for separating mapsH between spaces of real- and complex-valued functions have already been developed. The extension of such results to spaces of vector-valued functions is the general subject of this paper. We prove in Theorem 4.1, for example, for compactX andY, that a linear isometryH betweenC(X,E) andC(Y,F) is a “Banach-Stone” map if and only ifH is “biseparating (i.e,H andH ⁻¹ are separating). The Banach-Stone theorems of Jerison and Lau for vector-valued functions are then deduced in Corollaries 4.3 and 4.4 for the cases whenE andF or their topological duals, respectively, are strictly convex. Research supported by the Fundació Caixa Castelló, MI/25.043/92     

REML estimation for binary data in GLMMs

The restricted maximum likelihood (REML) procedure is useful for inferences about variance components in mixed linear models. However, its extension to hierarchical generalized linear models (HGLMs) is often hampered by analytically intractable integrals. Numerical integration such as Gauss-Hermite quadrature (GHQ) is generally not recommended when the dimensionality of the integral is high. With binary data various extensions of the REML method have been suggested, but they have had unsatisfactory biases in estimation. In this paper we propose a statistically and computationally efficient REML procedure for the analysis of binary data, which is applicable over a wide class of models and design structures. We propose a bias-correction method for models such as binary matched pairs and discuss how the REML estimating equations for mixed linear models can be modified to implement more general models.     

On the fast reduction of a quasiseparable matrix to Hessenberg and tridiagonal forms

In this paper we design a fast new algorithm for reducing an N × N quasiseparable matrix to upper Hessenberg form via a sequence of N − 2 unitary transformations. The new reduction is especially useful when it is followed by the QR algorithm to obtain a complete set of eigenvalues of the original matrix. In particular, it is shown that in a number of cases some recently devised fast adaptations of the QR method for quasiseparable matrices can benefit from using the proposed reduction as a preprocessing step, yielding lower cost and a simplification of implementation.     

Differential operators on Hilbert modular forms

We investigate differential operators and their compatibility with subgroups of SL₂n(R). In particular, we construct Rankin-Cohen brackets for Hilbert modular forms, and more generally, multilinear differential operators on the space of Hilbert modular forms. As an application, we explicitly determine the Rankin-Cohen bracket of a Hilbert-Eisenstein series and an arbitrary Hilbert modular form. We use this result to compute the Petersson inner product of such a bracket and a Hilbert modular cusp form.     

Engel Values in Residually Finite Groups

The following theorem is proved. Let n be a positive integer and q a power of a prime p. There exists a number m = m(n, q) depending only on n and q such that if G is any residually finite group satisfying the identity ([x _1,n y ₁] ⋯ [x _m,n y _m ])^q ≡ 1, then the verbal subgroup of G corresponding to the nth Engel word is locally finite.     

Ordered OLPS and CF Nth numerators

Definitions, theorems and examples are established for a general model of Laurent polynomial spaces and ordered orthogonal Laurent polynomial sequences, ordered with respect to ordered bases and orthogonal with respect to inner products ·=L°⊙ decomposed into transition functional ⊙ and strong moment functional, or, more generally, sample functional L couplings. Under this formulation that is shown to subsume those in the existing literature, new fundamental results are produced, including necessary and sufficient conditions for ordered OLPS to be sequences of nth numerators of continued fractions, in contrast to the classical result concerning nth denominators which is shown to hold only in special cases.     

The three space problem in topological groups

We study compact, countably compact, pseudocompact, and functionally bounded sets in extensions of topological groups. A property P is said to be a three space property if, for every topological group G and a closed invariant subgroup N of G, the fact that both groups N and G/N have P implies that G also has P. It is shown that if all compact (countably compact) subsets of the groups N and G/N are metrizable, then G has the same property. However, the result cannot be extended to pseudocompact subsets, a counterexample exists under p=c. Another example shows that extensions of groups do not preserve the classes of realcompact, Dieudonné complete and μ-spaces: one can find a pseudocompact, non-compact Abelian topological group G and an infinite, closed, realcompact subgroup N of G such that G/N is compact and all functionally bounded subsets of N are finite. Several examples given in the article destroy a number of tempting conjectures about extensions of topological groups.     

Grothendieck groups of unbounded complexes of finitely generated modules

Let Λ be a left Artinian ring, D⁺(mod Λ) (resp., D⁻(mod Λ), D(mod Λ)) the derived category of bounded below complexes (resp., bounded above complexes, unbounded complexes) of finitely generated left Λ-modules. We show that the Grothendieck groups K₀(D⁺(mod Λ)), K₀(D⁻(mod Λ)) and K₀(D(mod Λ)) are trivial. Received: 7 April 2005     

On the torsion of Jacobians of principal modular curves of level 3<Superscript><Emphasis Type="Italic">n</Emphasis></Superscript>

We demonstrate that the 3-power torsion points of the Jacobians of the principal modular curves X(3ⁿ) are fixed by the kernel of the canonical outer Galois representation of the pro-3 fundamental group of the projective line minus three points. The proof proceeds by demonstrating the curves in question satisfy a two-part criterion given by Anderson and Ihara. Two proofs of the second part of the criterion are provided; the first relies on a theorem of Shimura, while the second uses the moduli interpretation. Received: 30 September 2005