Convergence proofs and error estimates for an iterative method for conformal mapping

Summary The iterative method as introduced in [8] and [9] for the determination of the conformal mapping of the unit disc onto a domainG is here described explicitly in terms of the operatorK, which assigns to a periodic functionu its periodic conjugate functionK u. It is shown that whenever the boundary curve ofG is parametrized by a function with Lipschitz continuous derivative then the method converges locally in the Sobolev spaceW of 2-periodic absolutely continuous functions with square integrable derivative. If is in a Hölder classC ²⁺, the order of convergence is at least 1+. If is inC ^l+1+ withl1, 0<<1, then the iteration converges inC ^l+. For analytic boundary curves the convergence takes place in a space of analytic functions.For the numerical implementation of the method the operatorK can be approximated by Wittich's method, which can be applied very effectively using fast Fourier transform. The Sobolev norm of the numerical error can be estimated in terms of the numberN of grid points. It isO(N ^1–l–) if is inC ^l+1+, andO (exp (–N/2)) if is an analytic curve. The number in the latter formula is bounded by logR, whereR is the radius of the largest circle into which can be extended analytically such that'(z)0 for |z|<R. The results of some test calculations are reported.     

Round-off error analysis of iterations for large linear systems

Summary We deal with the rounding error analysis of successive approximation iterations for the solution of large linear systemsA _x =b. We prove that Jacobi, Richardson, Gauss-Seidel and SOR iterations arenumerically stable wheneverA=A ^*>0 andA has PropertyA. This means that the computed resultx _k approximates the exact solution with relative error of order A·A ^–1 where is the relative computer precision. However with the exception of Gauss-Seidel iteration the residual vector Ax _k –b is of order A² A ^–1 and hence the remaining three iterations arenot well-behaved.This work was partly done during the author's visit at Carnegie-Mellon University and it was supported in part by the Office of Naval Research under Contract N00014-76-C-0370; NR 044-422 and by the National Science Foundation under Grant MCS75-222-55     

Higher order methods for a singularly perturbed problem

Summary A finite element method using piecewise polynomials of degree k is used to approximate the problem u+u=f, >0 a small parameter. A very irregular mesh is used. On this mesh error estimates of order0(h ^k+1) are obtained uniformly in ,h the maximum stepsize, fork2. The condition number of the system of linear equations one has to solve in order to get the approximation is estimated. Extension of the results to more complicated problems is briefly indicated. Finally, a numerical example is given.Work performed while visiting the IBM Thomas J. Watson Research Center, Yorktown Heights, N.Y.     

Small potential corrections for the discrete eigenvalues of the Sturm-Liouville problem

Summary The inverse Sturm-Liouville problem is the problem of finding a good approximation of a potential functionq such that the eigenvalue problem (^*)–y +qy=y holds on (0, ) fory(0)=y()=0 and a set of given eigenvalues . Since this problem has to be solved numerically by discretization and since the higher discrete eigenvalues strongly deviate from the corresponding Sturm-Liouville eigenvalues , asymptotic corrections for the 's serve to get better estimates forq. Let _k (1kn) be the first eigenvalues of (^*), let _k be the corresponding discrete eigenvalues obtained by the finite element method for (^*) and let _{k k} for the special caseq=0. Then, starting from an asymptotic correction technique proposed by Paine, de Hoog and Anderssen, new estimates for the errors of the corrected discrete eigenvalues are obtained and confirm and improve the knownO(kh ²)(h:=/(n+1)) behaviour. The estimates are based on new Sobolev inequalities and on Fourier analysis and it is shown that for 4+c ₂ k(n+1)/2, wherec ₁ andc ₂ are constants depending onq which tend to 0 for vanishingq.     

Additive functions

1<q<2 L:= _n=1 1/q ⁿ=1/q–1. [0,1] _n()=1, A _n:= _i=1 ^n–1 _i(x)/q^i+1/n x _n(x)=0, _n>. , = _{n=1 n}(x)/qⁿ. F: [0,L]R , F(x)= _{n=1 n}(x)a_n, _n=1 ¦a _n¦<. [0,L]. q(1,2), . , q(1, 2), . .     

On Gauss quadrature formulas with equal weights

Summary It is well known that the Chebyshev weight function (1–x ²)^–1/2 is the only weight function (up to a linear transformation) for which then point Gauss quadrature formula has equal weights for alln. In this paper we describe all weight functions for which thenm point Gauss quadrature formula has equal weights for alln, wherem is fixed.     

The probability of the triangle inequality in probabilistic metric squares

The notion of a random semi-metric space provides an alternate approach to the study of probabilistic metric spaces from the standpoint of random variables instead of distribution functions and permits a new investigation of the triangle inequality. Starting with a probability space (, , P) and an abstract setS, each pair of points,p, q, inS is assigned a random variableX _pq with the interpretation thatX _pq () is the distance betweenp andq at the instant . The probability of the eventJ _pqr = { :X _pr ()X _pq ()+X _qr ()} is studied under distribution function conditions imposed by Menger Spaces (K. Menger, Statistical Metrics, Proc. Nat. Acad. Sci., U.S.A., 28 (1942), 535–537; B. Schweizer and A. Sklar, Statistical Metric Spaces, Pacific J. Math.10 (1960), 313–334). It turns out that for > 0 there are 3 non-negative, identically-distributed random variablesX, Y andZ for whichP(X < Y + Z) < . This and other results show that distribution function triangle inequalities are very weak. Conditional probabilities are introduced to give necessary and sufficient conditions forP(J _pqr ) = 1.     

Approximation of a nondifferentiable nonlinear problem related to MHD equilibria

Summary We present a simple method, based on a variant of the implicit function theorem, which leads to the existence of (a part of) a nontrivial solution branch of the nonlinear eigenvalue problem –u=u ⁺ in ,u=–1 on , where is a two-dimensional domain with boundary . The advantage of this method is that we can apply it for analysing the approximation of the above problem by a finite element method; the error analysis of the discrete problem appears immediately. We give also an iteration scheme which allows to solve the approximate problem.     

An explicit solution of the linear filtering problem for certain classes of nonrational spectral densities

We present an explicit solution of the problem of optimal linear filtering: the recovery of the useful signal(s) at the instantt+, (>0,<0, or=0) from known values of the received signal(s)=(s)+(s) in the past, i.e., at the instantt–s, s0. In doing so we assume the random processes(s) and /gr(s) are stationary and jointly stationary, while the stationary process of noise (s) with zero mean is assumed to be mutually correlated and jointly stationary with the process(s) under the assumption that there exists a common spectral densityf() for these processes.Translated fromTeoriya Sluchaínykh Protsessov, Vol. 14, pp. 83–91, 1986.     

1-Homogeneous Graphs with Cocktail Party μ-Graphs

Let be a graph with diameter d 2. Recall is 1-homogeneous (in the sense of Nomura) whenever for every edge xy of the distance partition{{z V() | (z, y) = i, (x, z) = j} | 0 i, j d}is equitable and its parameters do not depend on the edge xy. Let be 1-homogeneous. Then is distance-regular and also locally strongly regular with parameters (v,k,,), where v = k, k = a ₁, (v – k – 1) = k(k – 1 – ) and c ₂ + 1, since a -graph is a regular graph with valency . If c ₂ = + 1 and c ₂ 1, then is a Terwilliger graph, i.e., all the -graphs of are complete. In [11] we classified the Terwilliger 1-homogeneous graphs with c ₂ 2 and obtained that there are only three such examples. In this article we consider the case c ₂ = + 2 3, i.e., the case when the -graphs of are the Cocktail Party graphs, and obtain that either = 0, = 2 or is one of the following graphs: (i) a Johnson graph J(2m, m) with m 2, (ii) a folded Johnson graph J¯(4m, 2m) with m 3, (iii) a halved m-cube with m 4, (iv) a folded halved (2m)-cube with m 5, (v) a Cocktail Party graph K _{m × 2} with m 3, (vi) the Schläfli graph, (vii) the Gosset graph.     

Equivalence Bimodule Between Non-Commutative Tori

The non-commutative torus C ^*(ⁿ,) is realized as the C*-algebra of sections of a locally trivial C*-algebra bundle over S with fibres isomorphic to C ^*n/S, ₁) for a totally skew multiplier ₁ on ⁿ/S. D. Poguntke [9] proved that A is stably isomorphic to C(S) C(^*( Zⁿ/S, ₁) C(S) A M_kl( C) for a simple non-commutative torus A and an integer kl. It is well-known that a stable isomorphism of two separable C*-algebras is equivalent to the existence of equivalence bimodule between them. We construct an A-C(S) A-equivalence bimodule.     

On certain imbedding and extension theorems for a weighted class of functions

n- , S n–1 m (0m<n), () S . () , ().     

A Hubert space characterization among function spaces

— [0,1] ,E — - e=1 [0,1]. I — _E =1, E=L ₂ x _e =xL ₂ x E.

---

---

This work was prepared when the second author was a visiting professor of the CNR at the University of Firenze. He was supported by the Soros International Fund.     

Onp-Helson sets in R n

p- . E R ⁿ -, f () _p(R ⁿ)., ER ⁿ 2nq ₀, E— - q ₀(q ₀-1). : q₀>2 n1 E R ⁿ 2nq ₀, p- p<₀. , f-[-, ]ⁿ, f A _p(R ⁿ) , _p([-, ]ⁿ) (1 << ).     

Oscillation and non-oscillation theorems for superlinear Emden–Fowler equations of the fourth order

We study the oscillatory behavior of solutions of the fourth-order Emden–Fowler equation: (E) y^(iv)+q(t)|y|sgny=0, where >1 and q(t) is a positive continuous function on [t₀,), t₀>0. Our main results Theorem 2 – if (q(t)t^(3+5)/2)0, then equation (E) has oscillatory solutions; Theorem 3 – if lim_tq(t)t^4+(-1)=0, >0, then every solution y(t) of equation (E) is either non-oscillatory or satisfies limsup_tt^-+i|y⁽ⁱ⁾(t)|= for < and i=0,1,2,3,4. These results complement those given by Kura for equation (E) when q(t)<0 and provide analogues to the results of the second-order equation, y+q(t)|y|sgny=0,>1. Mathematics Subject Classification (2000) 34C10, 34C15     

A graded scale of parametric families of distributions,and parameter estimates based on the sample mean

Summary Denote by _k a class of familiesP={P} of distributions on the line R¹ depending on a general scalar parameter , being an interval of R¹, and such that the moments µ₁()=xdP ,...,µ_2k ()=x ^2k dP are finite, ₁ (), ..., _k (), _k+1 () ..., _k () exist and are continuous, with ₁ () 0, and _j +1 ()= ₁ () _j () +[₂() -₁()²] _j ()/ ₁ (), J=2, ..., k. Let ₁=¯x=x ₁ + ... +x _n/n, ₂=x ₁ ² + ... +x _n ²/n, ..., _k =(x ₁ ^k + ... +x _n ^k/n denote the sample moments constructed for a sample x₁, ..., x_n from a population with distribution Pg. We prove that the estimator of the parameter by the method of moments determined from the equation ₁= ₁() and depending on the observations x₁, ..., x_n only via the sample mean ¯x is asymptotically admissible (and optimal) in the class _k of the estimators determined by the estimator equations of the form ₀ () + ₁ () ₁ + ... + _k () _k =0 if and only ifP _k .The asymptotic admissibility (respectively, optimality) means that the variance of the limit, as n (normal) distribution of an estimator normalized in a standard way is less than the same characteristic for any estimator in the class under consideration for at least one 9 (respectively, for every ).The scales arise of classes _{1 2}... of parametric families and of classes _{1 2} ... of estimators related so that the asymptotic admissibility of an estimator by the method of moments in the class _k is equivalent to the membership of the familyP in the class _k .The intersection consists only of the families of distributions with densities of the form h(x) exp {C₀() + C₁() x } when for the latter the problem of moments is definite, that is, there is no other family with the same moments ₁ (), ₂ (), ...Such scales in the problem of estimating the location parameter were predicted by Linnik about 20 years ago and were constructed by the author in [1] (see also [2, 3]) in exact, not asymptotic, formulation.Translated from Problemy Ustoichivosti Stokhasticheskikh Modelei, pp. 41–47, 1981.     

Discrete Spectrum of Differential Operators in Spectral Gaps under Nonnegative Perturbations of High Order

Let A be a self-adjoint elliptic second-order differential operator, let (, ) be an inner gap in the spectrum of A, and let B(t) = A + tW ^* W, where W is a differential operator of higher order. Conditions are obtained under which the spectrum of the operator B(t) in the gap (, ) is either discrete, or does not accumulate to the right-hand boundary of the spectral gap, or is finite. The quantity N(, A, W, ), (, ), > 0 (the number of eigenvalues of the operator B(t) passing the point (, ) as t increases from 0 to ) is considered. Estimates of N(, A, W, ) are obtained. For the perturbation W ^* W of a special form, the asymptotics of N(, A, W, ) as + is given. Bibliography: 5 titles.     

Isolated points in duals of certain locally compact groups

Summary LetG be a separable locally compact group with dual space. consists of all equivalence classes of irreducible unitary representations ofG, and is endowed with the Fell-topology. We study the topological properties in of the square-integrable representations ofG. [ is square-integrable provided there is a coordinate functiong((g)v, v),gG, for which is inL ²(G) w.r.t. left Haar measure onG.]SupposeG contains an open normal subgroupN of the formeKN ⁿ e whereK is compact. (All groups with a compact invariant neighborhood of the identity, [IN] groups, satisfy this condition.) In this case we show that if is square-integrable then {} is an open point of.Finally, our techniques are used to prove this result for arbitrary (non connected) nilpotent Lie groups.     

Limiting behavior of certain mixtures of distributions

U — [0, 1] Y — . X=[1–U ^1/v /Y], U Y.