On application of the Lanczos method to solution of some partial differential equations

Let A be a square symmetric n × n matrix, φ be a vector from ⁿ, and f be a function defined on the spectral interval of A. The problem of computation of the vector u = f(A)φ arises very often in mathematical physics.

We propose the following method to compute u. First, perform m steps of the Lanczos method with A and φ. Define the spectral Lanczos decomposition method (SLDM) solution as u_m = φ Qf(H)e₁, where Q is the n × m matrix of the m Lanczos vectors and H is the m × m tridiagonal symmetric matrix of the Lanczos method. We obtain estimates for u − u_m that are stable in the presence of computer round-off errors when using the simple Lanczos method.

We concentrate on computation of exp(− tA)φ, when A is nonnegative definite. Error estimates for this special case show superconvergence of the SLDM solution. Sample computational results are given for the two-dimensional equation of heat conduction. These results show that computational costs are reduced by a factor between 3 and 90 compared to the most efficient explicit time-stepping schemes. Finally, we consider application of SLDM to hyperbolic and elliptic equations.     

Full quadrature sums for pth powers of polynomials with Freud weights

In this paper, we provide a solution of the quadrature sum problem of R. Askey for a class of Freud weights. Let r> 0, b (− ∞, 2]. We establish a full quadrature sum estimate

1 p < ∞, for every polynomial P of degree at most n + rn^1/3, where W² is a Freud weight such as exp(−¦x¦), > 1, λ_jn are the Christoffel numbers, x_jn are the zeros of the orthonormal polynomials for the weight W², and C is independent of n and P. We also prove a generalisation, and that such an estimate is not possible for polynomials P of degree M = m(n) if m(n) = n + ξ_nn^1/3, where ξ_n → ∞ as n → ∞. Previous estimates could sum only over those x_jn with ¦x_jn¦ σx_1n, some fixed 0 < σ < 1.     

A characterization on n-critical economical generalized tic-tac-toe games

There is a so called generalized tic-tac-toe game playing on a finite set X with winning sets A₁, A₂,…, A_m. Two players, F and S, take in turn a previous untaken vertex of X, with F going first. The one who takes all the vertices of some winning set first wins the game. Erd s and Selfridge proved that if |A₁|=|A₂|==|A_m|=n and m<2ⁿ⁻¹, then the game is a draw. This result is best possible in the sense that once m=2ⁿ⁻¹, then there is a family A₁, A₂,…, A_m so that F can win. In this paper we characterize all those sets A₁,…, A_2n−1 so that F can win in exactly n moves. We also get similar result in the biased games.     

A note on linearization of actions of finitely semisimple Hopf algebras on local algebras

Let H be a Hopf algebra over a field k and let H A → A, h a → h.a, be an action of H on a commutative local Noetherian kalgebra (A, m). We say that this action is linearizable if there exists a minimal system x₁, …, x_n of generators of the maximal ideal m such that h.x_i ε kx₁ + …+ kx_n for all h ε H and i = 1, …, n. In the paper we prove that the actions from a certain class are linearizable (see Theorem 4), and we indicate some consequences of this fact.     

Krylov subspace methods for eigenvalues with special properties and their analysis for normal matrices

In this paper we propose a general approach by which eigenvalues with a special property of a given matrix A can be obtained. In this approach we first determine a scalar function ψ: C → C whose modulus is maximized by the eigenvalues that have the special property. Next, we compute the generalized power iterations uinj + 1 = ψ(A)u_j, j = 0, 1,…, where u₀ is an arbitrary initial vector. Finally, we apply known Krylov subspace methods, such as the Arnoldi and Lanczos methods, to the vector u_n for some sufficiently large n. We can also apply the simultaneous iteration method to the subspace span{x⁽ⁿ⁾₁,…,x⁽ⁿ⁾_k} with some sufficiently large n, where x^(j+1)_m = ψ(A)x^(j)_m, j = 0, 1,…, m = 1,…, k. In all cases the resulting Ritz pairs are approximations to the eigenpairs of A with the special property. We provide a rather thorough convergence analysis of the approach involving all three methods as n → ∞ for the case in which A is a normal matrix. We also discuss the connections and similarities of our approach with the existing methods and approaches in the literature.     

Linear operators preserving unitarily invariant norms of matrices

Let F_{m × n} be the set of all m × n matrices over the field F = C or R Denote by U_n(F) the group of all n × n unitary or orthogonal matrices according as F = C or F-R. A norm N() on Fm ×n, is unitarily invariant if N(UAV) = N(A): for all A ∈ F _m×n U ∈ U_m(F). and V ∈ U_n(F). We characterize those linear operators TF_{m × n} → F_{m × n}which satisfy N (T(A)) = N(A)for all A ∈ F_{m × n}

for a given unitarily invariant norm N(). It is shown that the problem is equivalent to characterizing those operators which preserve certain subsets in F_{m × n} To develop the theory we prove some results concerning unitary operators on F_{m × n} which are of independent interest.     

Flawless 0-sequences and Hilbert functions of Cohen-Macaulay integral domains

Let A = A₀A₁ be a commutative graded ring such that (i) A₀ = k a field, (ii) A = k[A₁] and (iii) dim_k A₁ < ∞. It is well known that the formal power series ∑^∞_{n = 0} (dim_kA_n)λ ⁿ is of the form (h₀ + h₁λ + + h_sλ^s)/(1 − λ)^dimA with each h_iε . We are interested in the sequence (h₀, h₁,…,h_s), called the h-vector of A, when A is a Cohen–Macaulay integral domain. In this paper, after summarizing fundamental results (Section 1), we study h-vectors of certain Gorenstein domains (Section 2) and find some examples of h-vectors arising from integrally closed level domains (Sections 3 and 4).     

Linear recurrence relations in the algebra of matrices and applications

We study here some linear recurrence relations in the algebra of square matrices. With the aid of the Cayley–Hamilton Theorem, we derive some explicit formulas for Aⁿ (nr) and e^tA for every r×r matrix A, in terms of the coefficients of its characteristic polynomial and matrices A^j, where 0jr−1.     

MEROMORPHIC FUNCTIONS SHARING TWO FINITE SETS

Let S1 ＝ {∞} and S2 ＝ {w: Ps(w)＝ 0}, Ps(w) being a uniqueness polynomial under some restricted conditions. Then, for any given nonconstant meromorphic function f, there exist at most finitely many nonconstant meromorphic functions g such that f-1(Si) ＝ g-1(Si)(i ＝ 1,2), where f-1(Si) and g-1(Si) denote the pull-backs of Si considered as a divisor, namely, the inverse images of Si counted with multiplicities, by f and g respectively.     

Matrices with zero trace as commutators of nilpotents

Let C ∈ F^n×n have minimum polynomial m(x). Suppose C is of zero trace and m(x) splits over F. Then, except when n = 2 and m(x) = (x - c)² or when n = 3 and m(x) = x - c)² with c ≠ 0, there exist nilpotents A, B ∈ F^n×n such that C = AB - BA.     

Equality of higher numerical ranges of matrices and a conjecture of Kippenhahn on Hermitian pencils

Let M_n be the algebra of all n × n complex matrices. For 1 k n, the kth numerical range of A M_n is defined by W_k(A) = (1/k)∑_j^k=1x_j^*Ax_j : x₁, …, x_k is an orthonormal set in ⁿ]. It is known that tr A/n = W_n(A) W_n−1(A) W₁(A). We study the condition on A under which W_m(A) = W_k(A) for some given 1 m < k n. It turns out that this study is closely related to a conjecture of Kippenhahn on Hermitian pencils. A new class of counterexamples to the conjecture is constructed, based on the theory of the numerical range.     

Resolutions of Stanley-Reisner rings and Alexander duality

Associated to any simplicial complex Δ on n vertices is a square-free monomial ideal I_Δ in the polynomial ring A = k[x₁, …, x_n], and its quotient k[Δ] = A/I_Δ known as the Stanley-Reisner ring. This note considers a simplicial complex Δ^* which is in a sense a canonical Alexander dual to Δ, previously considered in [1, 5]. Using Alexander duality and a result of Hochster computing the Betti numbers dim_k Tor_i^A (k[Δ],k), it is shown (Proposition 1) that these Betti numbers are computable from the homology of links of faces in Δ^*. As corollaries, we prove that I_Δ has a linear resolution as A-module if and only if Δ^* is Cohen-Macaulay over k, and show how to compute the Betti numbers dim_k Tor_i^A (k[Δ],k) in some cases where Δ^* is wellbehaved (shellable, Cohen-Macaulay, or Buchsbaum). Some other applications of the notion of shellability are also discussed.     

Balanced partitions of two sets of points in the plane

We prove the following theorem. Let m≥2 and q≥1 be integers and let S and T be two disjoint sets of points in the plane such that no three points of ST are on the same line, |S|=2q and |T|=mq. Then ST can be partitioned into q disjoint subsets P₁,P₂,…,P_q satisfying the following two conditions: (i) conv(P_i)∩conv(P_j)=φ for all 1≤i<j≤q, where conv(P_i) denotes the convex hull of P_i; and (ii) |P_i∩S|=2 and |P_i∩T|=m for all 1≤i≤q.     

The kernels of irreducible representations of the general linear group

A derivation for the kernel of the irreducible representation T_(λ) of the general linear group GL_n(C) is given. This is then applied to the problem of determining necessary and sufficient conditions under which T_(λ)(A) = T_(λ)(B), where A and B are linear transformations, not necessarily invertible. Finally, conditions are obtained under which normality of T_(λ)(A) implies normality of A.     

Possible values for 2

From GCH and P^{m(κ)-hypermeasurable} (1 <m<gw), we construct a model satisfying 2_ⁿ = _a(n) and 2_^ω = _ω+m for a monotone a:ω→ω satisfying a(n)>n.     

On the Differential Polynomial of a Graph

We introduce the differential polynomial of a graph. The differential polynomial of a graph G of order n is the polynomial B(G; x) :=∑?(G)k=-nB_k(G) x~(n+k), where B_k(G) denotes the number of vertex subsets of G with differential equal to k. We state some properties of B(G;x) and its coefficients.In particular, we compute the differential polynomial for complete, empty, path, cycle, wheel and double star graphs. We also establish some relationships between B(G; x) and the differential polynomials of graphs which result by removing, adding, and subdividing an edge from G.     

On generalized Boolean functions II

This cycle of papers is based on the concept of generalized Bolean functions introduced by the author in the first article of the series. Every generalized Boolean function f:Bⁿ→B can be written in a manner similar to the canonical disjunctive form using some function defined on A×B, where A is a finite subset of B containing 0 and 1. The set of those functions f is denoted by GBF_n[A]. In this paper the following questions are presented: (1) What is the relationship between GBF_n[A₁] and GBF_n[A₂] when A₁A₂. (2) What can be said about GBF_n[A₁∩A₂] and GBF_n[A₁A₂] in comparison with GBF_n[A₁]∩GBF_n[A₂] and GBF_n[A₁]GBF_n[A₂], respectively.     

A matrix problem associated with computing Cohen-Macaulay type

Suppose A∈M_n×m(F), B∈M_n×t(F) for some field F. Define Г(AB) to be the set of n×n diagonal matrices D such that the column space of DA is contained in the column space of B. In this paper we determine dim Г(AB). For matrices AB of the same rank we provide an algorithm for computing dim Г(AB).     

Strong Laws of Large Numbers for Double Sums of Banach Space Valued Random Elements

For a double array {V_(m,n), m ≥ 1, n ≥ 1} of independent, mean 0 random elements in a real separable Rademacher type p(1 ≤ p ≤ 2) Banach space and an increasing double array {b_(m,n), m ≥1, n ≥ 1} of positive constants, the limit law ■ and in L_p as m∨n→∞ is shown to hold if ■ This strong law of large numbers provides a complete characterization of Rademacher type p Banach spaces. Results of this form are also established when 0 p ≤ 1 where no independence or mean 0 conditions are placed on the random elements and without any geometric conditions placed on the underlying Banach space.     

Efficient robust estimation of parameter in the random censorship model

For an open set Θ of ^k, let \s{P_θ: θ Θ\s} be a parametric family of probabilities modeling the distribution of i.i.d. random variables X₁,…, X_n. Suppose X_i's are subject to right censoring and one is only able to observe the pairs (min(X_i, Y_i), [X_i Y_i]), i = 1,…, n, where [A] denotes the indicator function of the event A, Y₁,…, Y_n are independent of X₁,…, X_n and i.i.d. with unknown distribution Q₀. This paper investigates estimation of the value θ that gives a fitted member of the parametric family when the distributions of X₁ and Y₁ are subject to contamination. The constructed estimators are adaptive under the semi-parametric model and robust against small contaminations: they achieve a lower bound for the local asymptotic minimax risk over Hellinger neighborhoods, in the Hájel—Le Cam sense. The work relies on Beran (1981). The construction employs some results on product-limit estimators.