Brownian charges around loops

Summary Let (X _t ⁿ ) be a Poisson sequence of independent Brownian motions in ^d ,d3; Let be a compact oriented submanifold of d, of dimensiond–2 and volume ; let t be the sum of the windings of (X _s ⁿ , 0st) around ; then t/t converges in law towards a Cauchy variable of parameter /2. A similar result is valid when the winding is replaced by the integral of a harmonic 1-form in ^d .     

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     

Generalizations of the projection method with applications to SOR theory for hermitian positive semidefinite linear systems

Summary Ann×n complex matrixB is calledparacontracting if B₂1 and 0x[N(I-B)]Bx₂<x₂. We show that a productB=B _k B _k–1 ...B ₁ ofk paracontracting matrices is semiconvergent and give upper bounds on the subdominant eigenvalue ofB in terms of the subdominant singular values of theB _i 's and in terms of the angles between certain subspaces. Our results here extend earlier results due to Halperin and due to Smith, Solomon and Wagner. We also determine necessary and sufficient conditions forn numbers in the interval [0, 1] to form the spectrum of a product of two orthogonal projections and hence characterize the subdominant eigenvalue of such a product. In the final part of the paper we apply the upper bounds mentioned earlier to provide an estimate on the subdominant eigenvalue of the SOR iteration matrix associated with ann×n hermitian positive semidefinite matrixA none of whose diagonal entries vanish.The work of this author was supported in part by NSF Research Grant No. MCS-8400879     

Black holes on the plane drawn by a Wiener process

Summary We say that the discD()R ², of radius , located around the origin isp-covered in timeT by a Wiener processW(·) if for anyzD() there exists a 0tT such thatW(t) is a point of the disc of radiusp, located aroundz. The supremum of those 's (0) is studied for which,D() isp-covered inT.     

Elementarym-harmonic cardinal B-splines

We generalize the notion of B-spline to the thin plate splines and to otherd-dimensional polyharmonic splines as defined in [Duchon, [3]]; for regular nets, we give the main properties of these B-splines: Fourier transform, decay when x , stability, integration property, links between B-splines of different orders or of different dimensions and in particular link with the polynomial B-splines, approximation using B-splines... We show that, in some sense, B-splines may be considered as a regularized form of the Dirac distribution.     

Multilevel large deviations and interacting diffusions

Summary Let ( ^N ) be a sequence of random variables with values in a topological space which satisfy the large deviation principle. For eachM and eachN, let ^{M, N} denote the empirical measure associated withM independent copies of ^N . As a main result, we show that ( ^{M, N} ) also satisfies the large deviation principle asM,N. We derive several representations of the associated rate function. These results are then applied to empirical measure processes ^{M, N} (t) =M ^–1 _i=1 ^N _i ^N (t) 0tT, where ( ₁ ^N ,..., _M ^N (t)) is a system of weakly interacting diffusions with noise intensity 1/N. This is a continuation of our previous work on the McKean-Vlasov limit and related hierarchical models ([4], [5]).Research partially supported by a Natural Science and Engineering Research Council of Canada operating grant     

Distoriton of area and conditioned Brownian motion

Summary In a simply connected planar domainD the expected lifetime of conditioned Brownian motion may be viewed as a function on the set of hyperbolic geodesics for the domain. We show that each hyperbolic geodesic induces a decomposition ofD into disjoint subregions and that the subregions are obtained in a natural way using Euclidean geometric quantities relating toD. The lifetime associated with on each _j is then shown to be bounded by the product of the diameter of the smallest ball containing _j and the diameter of the largest ball in _j . Because this quantity is never larger than, and in general is much smaller than, the area of the largest ball in _j it leads to finite lifetime estimates in a variety of domains of infinite area.Research of the first author was supported in part by NSF Grant DMS-9100811Research of the second author was supported in part by NSF Grant DMS-9105407     

On monotone extensions of boundary data

Summary A functionf C (), is called monotone on if for anyx, y the relation x – y ₊ ^s impliesf(x)f(y). Given a domain with a continuous boundary and given any monotone functionf on we are concerned with the existence and regularity ofmonotone extensions i.e., of functionsF which are monotone on all of and agree withf on . In particular, we show that there is no linear mapping that is capable of producing a monotone extension to arbitrarily given monotone boundary data. Three nonlinear methods for constructing monotone extensions are then presented. Two of these constructions, however, have the common drawback that regardless of how smooth the boundary data may be, the resulting extensions will, in general, only be Lipschitz continuous. This leads us to consider a third and more involved monotonicity preserving extension scheme to prove that, when is the unit square [0, 1]² in ², strictly monotone analytic boundary data admit a monotone analytic extension.Research supported by NSF Grant 8922154Research supported by DARPA: AFOSR #90-0323     

Large deviations for the invariant measure of a reaction-diffusion equation with non-Gaussian perturbations

Summary In this paper we establish a large deviations principle for the invariant measure of the non-Gaussian stochastic partial differential equation (SPDE) _t v =v +f(x,v )+(x,v ) . Here is a strongly-elliptic second-order operator with constant coefficients, h:=DH _xx-h, and the space variablex takes values on the unit circleS ¹. The functionsf and are of sufficient regularity to ensure existence and uniqueness of a solution of the stochastic PDE, and in particular we require that 0<mM wherem andM are some finite positive constants. The perturbationW is a Brownian sheet. It is well-known that under some simple assumptions, the solutionv ² is aC ^k (S ¹)-valued Markov process for each 0<1/2, whereC (S ¹) is the Banach space of real-valued continuous functions onS ¹ which are Hölder-continuous of exponent . We prove, under some further natural assumptions onf and which imply that the zero element ofC (S ¹) is a globally exponentially stable critical point of the unperturbed equation _t ⁰ = ⁰ +f(x,⁰), that has a unique stationary distributionv ^K, on (C (S ¹), (C ^K (S ¹))) when the perturbation parameter is small enough. Some further calculations show that as tends to zero,v ^K, tends tov ^K,0, the point mass centered on the zero element ofC (S ¹). The main goal of this paper is to show that in factv ^K, is governed by a large deviations principle (LDP). Our starting point in establishing the LDP forv ^K, is the LDP for the process , which has been shown in an earlier paper. Our methods of deriving the LDP forv ^K, based on the LDP for are slightly non-standard compared to the corresponding proofs for finite-dimensional stochastic differential equations, since the state spaceC (S ¹) is inherently infinite-dimensional.This work was performed while the author was with the Department of Mathematics, University of Maryland, College Park, MD 20742, USA     

Orthogonal least squares fitting with linear manifolds

Summary The problem is considered of fitting a linear manifold of dimensions with 1sn–1 to a given set of points in ⁿ such that the sum of orthogonal squared distances attains a minimum.     

The central limit theorem for Chébli-Trimèche hypergroups

Let (S _nn>-1) be a random walk on a hypergroup ( ₊ , *), i.e., a Markov chain with transition kernelN(x, A) = _x * (A), where is a fixed probability measure on ₊ such that the second moment exists. Then depending on the growth of the hypergroup two situations can occur: when ( ₊ , *) is of exponential growth then it is shown thatS _n is asymptotically normal. In the case of polynomial growth {more precisely, if the densityA of the Haar measure of ( ₊ , *) satisfies lim[A()/A()]=}, the normalized variablesS _n/[n Var()/(+1)]^1/2 converge to a Rayleigh distribution with parameter .     

Positive Solutions for a Singular Nonlinear Problem on a Bounded Domain in R 2

For a bounded regular Jordan domain in R ², we introduce and study a new class of functions K() related on its Green function G. We exploit the properties of this class to prove the existence and the uniqueness of a positive solution for the singular nonlinear elliptic equation u+(x,u)=0, in D(), with u=0 on and uC(), where is a nonnegative Borel measurable function in ×(0,) that belongs to a convex cone which contains, in particular, all functions (x,t)=q(x)t ^–,>0 with nonnegative functions qK(). Some estimates on the solution are also given.     

Block SSOR preconditionings for high order 3d fe systems

For solving 3D high order hierarchical FE systems the block SSOR preconditioned CG algorithms based on new stripwise block two-color orderings of degrees of freedom and providing for efficient concurrent/vector implementation are suggested. As demonstrated by numerical results for the 3D Navier equations approximated using hierarchical orderp, 2 p 5, FE's the convergence rate of such BSSOR-CG algorithms is only slightly dependent onp and mesh nonunformity.     

Harmonic analysis on tori

We consider measurable subsets {ofR}ⁿ with 0<m()<, and we assume that has a spectral set . (In the special case when is also assumed open, may be obtained as the joint spectrum of a family of commuting self-adjoint operators {H _k: 1kn} in L ² () such that each H _k is an extension of i(/x _k) on C _c (), k=1, ..., n.)It is known that is a fundamental domain for a lattice if is itself a lattice. In this paper, we consider a class of examples where is not assumed to be a lattice. Instead is assumed to have a certain inhomogeneous form, and we prove a necessary and sufficient condition for to be a fundamental domain for some lattice in {ofR}ⁿ. We are thus able to decide the question, fundamental domain or not, by considering only properties of the spectrum . Our criterion is obtained as a corollary to a theorem concerning partitions of sets which have a spectrum of inhomogeneous form.Work supported in part by the NSF.Work supported in part by the NSRC, Denmark.     

Non-Uniformly Elliptic Equations with Natural Growth in the Gradient

We prove the existence of bounded solutions for a class of nonlinear elliptic problems of type–div(a(x,u,Du))=H(x,u,Du)+f, uW ^1,p ₀()L (),where a(x,,)b(||)|| ^p , b is a continuous monotone decreasing function and |H(x,,)| k()|| ^p , k is a continuous monotone increasing function.     

Circularity of Finite Groups without Fixed Points

Let be a fixed point free group given by the presentation where and are relative prime numbers, t = /s and s = gcd( – 1,), and is the order of modulo . We prove that if (1) = 2, and (2) is embeddable into the multiplicative group of some skew field, then is circular. This means that there is some additive group N on which acts fixed point freely, and |((a)+b)((c)+d)| 2 whenever a,b,c,d N, a0c, are such that (a)+b(c)+d.     

Left linear theories — A generalization of module theory

In this paper we introduce left linear theories of exponentN (a set) on the setL as mapsL ×L ^N (l, ) l · L such that for alll L and , L ^N the relation (l · ) =l( · ) holds, where · L ^N is given by ( · )(i) = (i),i N. We assume thatL has a unit, that is an element L ^N withl · =l, for alll L, and · = , for all L ^N . Next, left (resp. right)L-modules andL-M-bimodules and their homomorphisms are defined and lead to categoriesL-Mod, Mod-L, andL-M-Mod. These categories are algebraic categories and their free objects are described explicitly. Finally, Hom(X, Y) andX Y are introduced and their properties are investigated.Herrn Professor Dr. D. Pumplün zum 60. Geburtstag gewidmet     

Theorem on a convergence condition in the spaces

For a given -function (u), a condition on a -function (u) is found such that it is necessary and sufficient for the following to hold: if f_n(x) f(x) and f _n (x)M (n=1, 2, ...) where M>0 is an absolute constant, then f _n (x)–f(x)0(n). An analogous condition for convergence in Orlicz spaces is obtained as a corollary.Translated from Matematicheskie Zametki, Vol. 21, No. 5, pp. 615–626, May, 1977.The author thanks V. A. Skvortsov for his constant attention and guidance on this paper.     

On semiregular conditional distributions

Let (,A,P) denote some probability space and some sub--algebra ofA. It is shown that there exists a semiregular versionQ (A),A, , of the conditional distributionP(A|), AA, i.e., Q (A), (AA fixed) is andAQ (A),AA ( fixed), is a probability charge satisfyingQ (N)=0, , for allP-zero setsN, if and only ifL ₁(,P|) has a lifting, which exists for any sub--algebra ofA ifL ₁(,A P) is separable. Separability ofL ₁(,A,P) implies also the existence of a strongly semiregular versionQ (A),A, , ofP(A|), A , i.e., Q (A), (AA fixed), is -measurable andAQ (A),A ( fixed), is a probability charge. Furthermore,P can be written as P ₁+(1–)P ₂, 01, whereP ₁ are probability measures onA such thatP ₁(A|),AA, has a semiregular version vanishing for anyP-zero setN andP ₂ is singular with respect to any probability measure onA of the type ofP ₁. In the case 0<<1 the probability measuresP _j ,j=1, 2, are uniquely determined. The decomposition can be carried over to the case, where the additional condition thatQ (N)=0 for all and anyP-zero setN is valid, is omitted respectively semiregularity is replaced by (i) strong semiregularity, or (ii) classical regularity. In the last mentioned case (ii) the decomposition is multiplicative.