Small into isomorphisms onL ∞ spaces

IfT is an isomorphism ofL (A, ) intoL (B, ) which satisfies the condition T T ^–11+, where (A, ) is a -finite measure space, thenT/T is close to an isometry with an error less than 4.     

An algorithm for best approximate solutions ofAx=b with a smooth strictly convex norm

Summary In this paper, overdetermined systems ofm linear equations inn unknowns are considered. With ^m equipped with a smooth strictly convex norm, ·, an iterative algorithm for finding the best approximate solution of the linear system which minimizes the ·-error is given. The convergence of the algorithm is established and numerical results are presented for the case when · is anl _p norm, 1<p<.Portions of this paper are taken from the author's Ph.D. thesis at Michigan State University     

Numerical stability of the Chebyshev method for the solution of large linear systems

Summary This paper contains the rounding error analysis for the Chebyshev method for the solution of large linear systemsAx+g=0 whereA=A ^* is positive definite. We prove that the Chebyshev method in floating point arithmetic is numerically stable, which means that the computed sequence {x _k} approximates the solution such that x _k – is of order AA ^–1 where is the relative computer precision.We also point out that in general the Chebyshev method is not well-behaved, which means that the computed residualsr _k=Ax _k+g are of order A²A ^–1.This work was supported in part by the Office of Naval Research under Contract N0014-67-0314-0010, NR 044-422 and by the National Science Foundation under Grant GJ32111     

Stability inL 1 of some integral operators

A class of Markov operators appearing in biomathematics is investigated. It is proved that these operators are asymptotic stable inL ¹, i.e. lim _n P ⁿ f=0 forfL ¹ and f(x) dx=0.     

Upper bounds for theL p -norms of martingales

Summary Let (f _n ) be a martingale. We establish a relationship between exponential bounds for the probabilities of the typeP(|f _n |>·T(f _n )) and the size of the constantC _p appearing in the inequality f ^* _p C _p T ^*(f) _p , for some quasi-linear operators acting on martingales.This research was supported in part by NSF Grant, no. DMS-8902418On leave from Academy of Physical Education, Warsaw, Poland     

The functional equation of multiplicative derivation is superstable on standard operator algebras

LetX be a real or complex infinite dimensional Banach space andA a standard operator algebra onX. Denote byB(X) the algebra of all bounded linear operators onX. Let : _{+ +} be a function with the property lim _t (t)t ^–1=0. Assume that a mappingD:A B(X) satisfies D(AB)–AD(B)–D(A)B<(A B) for all operatorsA, B D (no linearity or continuity ofD is assumed). ThenD is of the formD(A)=AT–TA for someTB(X).This work was supported by the Research Council of Slovenia     

Convergence of random products of compact contractions in Hilbert space

Let {T₁, ..., T_N} be a finite set of linear contraction mappings of a Hilbert space H into itself, and let r be a mapping from the natural numbers N to {1, ..., N}. One can form S_n=T_r(n)...T_r(1) which could be described as a random product of the T_i's. Roughly, the S_n converge strongly in the mean, but additional side conditions are necessary to ensure uniform, strong or weak convergence. We examine contractions with three such conditions. (W): x_n1, Tx_n1 implies (I-T)x_n0 weakly, (S): x_n1, Tx_n1 implies (I-T)x_n0 strongly, and (K): there exists a constant K>0 such that for all x, (I-T)x²K(x²–Tx²).We have three main results in the event that the T_i's are compact contractions. First, if r assumes each value infinitely often, then S_n converges uniformly to the projection Q on the subspace _i= ₁ ^N [x|T_ix=x]. Secondly we prove that for such compact contractions, the three conditions (W), (S), and (K) are equivalent. Finally if S=S(T₁, ..., T_N) denotes the algebraic semigroup generated by the T_i's, then there exists a fixed positive constant K such that each element in S satisfies (K) with that K.     

The measure of a translated ball in uniformly convex spaces

Let (E, ¦·¦) be a uniformly convex Banach space with the modulus of uniform convexity of power type. Let be the convolution of the distribution of a random series inE with independent one-dimensional components and an arbitrary probability measure onE. Under some assumptions about the components and the smoothness of the norm we show that there exists a constant such that |{·<t}–{·+r<t}|r ^q , whereq depends on the properties of the norm. We specify it in the case ofL spaces, >1.     

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     

Resolvent of the Toeplitz operator may increase arbitrarily fast

It is proved that for every sequence of points _n from the unit circle, _n1, and for an arbitrary sequence of positive numbers A_n, A_n, there exists a continuous real function u, such that for the Toeplitz operator T (acting in the Hardy space H²) with the symbol =e ^iu we have the estimates (T–_nI)^–1>A_n, n.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova, AN SSSR, Vol. 157, pp. 175–177, 1987.     

Attracting properties for one dimensional flows of a general barotropic viscous fluid. Periodic flows

Summary We consider the motion of a barotropic compressible fluid in a one dimensional bounded region with impermeable boundary, see equation (1.1). Here, u(t, q) denotes the velocity and v(t, q) the specific volume. The quantity log v(t, q) measures the displacement of v(t, q) with respect to the equilibrium v 1. For the sake of brevity we denote here different norms by the simbol . We show that there is a positive constant r₀=r₀(), a small ball B₁ (r) (with radius R₁ (r), ), and a large ball B(r) (with radius R(r), ) such that the following holds, for each r [0, r₀ [(i) If f(t) < r for all t 0, and if (u(0), log v(0))R(r) (i.e. (u(0), log v(0)) B(r)) then, for sufficiently large values of t, (u(t), log v(t))R₁ (r); (ii) The solutions starting at time t=0 from the large ball B(r) have all the same asymptotic behaviour (see (1.11)); (iii) If f is T-periodic then there is a (unique) T-periodic solution (u(t), log v(t)) inside the small ball B₁ (r). This periodic solution atracts all solutions which intersect the large ball B(r). Periodic solutions had been previously studied only for very specific pressure laws, namely p(v)-log v and p(v)-v^–1.     

An extremal property of outer functions

Our main result is the following: iff (z) is in the space H², and F(z) is its outer part, then F⁽ⁿ⁾_H2F⁽ⁿ⁾_H2(n=1,2,...), the left side being finite if the right side is finite. Under certain essential restrictions, this. inequality was proved by B. I. Korenblyum and V. S. Korolevich [1].Translated from Matematicheskie Zametki, Vol. 10, No. 1, pp. 53–56, July, 1971.     

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     

On the condition number of linear least squares problems in a weighted Frobenius norm

LetA be anm × n, m n full rank real matrix andb a real vector of sizem. We give in this paper an explicit formula for the condition number of the linear least squares problem (LLSP) defined by min Ax–b₂,x ⁿ . Let and be two positive real numbers, we choose the weighted Frobenius norm [A, b] _F on the data and the usual Euclidean norm on the solution. A straightforward generalization of the backward error of [9] to this norm is also provided. This allows us to carry out a first order estimate of the forward error for the LLSP with this norm. This enables us to perform a complete backward error analysis in the chosen norms.Finally, some numerical results are presented in the last section on matrices from the collection of [5]. Three algorithms have been tested: the QR factorization, the Normal Equations (NE), the Semi-Normal Equations (SNE).     

A Note on the Optimal L 2-Estimate of the Finite Volume Element Method

In this note, the optimal L ²-error estimate of the finite volume element method (FVE) for elliptic boundary value problem is discussed. It is shown that u–u _{h 0}Ch ²|ln h|^1/2f_1,1 and u–u _{h 0}Ch ²f_1,p , p>1, where u is the solution of the variational problem of the second order elliptic partial differential equation, u _h is the solution of the FVE scheme for solving the problem, and f is the given function in the right-hand side of the equation.     

On some variational problems for 2-dimensional Hermitian metrics

Summary Let (M, J, g) be a compact complex 2-dimensional Hermitian manifold with the Kähler form , and the torsion 1-form defined by d = . In this note we obtain the Euler-Lagrange equations for the variational functionals defined by ² and d², whereg runs in the space of all the Hermitian metrics onM. In the first case, the extremals are precisely the Kähler metrics [Gd]. In the second case, we also write down a formula for the second variation.Communicated by J. Szenthe     

A superfast algorithm for multi-dimensional Padé systems

For a vector ofk+1 matrix power series, a superfast algorithm is given for the computation of multi-dimensional Padé systems. The algorithm provides a method for obtaining matrix Padé, matrix Hermite Padé and matrix simultaneous Padé approximants. When the matrix power series is normal or perfect, the algorithm is shown to calculate multi-dimensional matrix Padé systems of type (n ₀,...,n _k ) inO(n · log²n) block-matrix operations, where n=n ₀+...+n _k . Whenk=1 and the power series is scalar, this is the same complexity as that of other superfast algorithms for computing Padé systems. Whenk>1, the fastest methods presently compute these matrix Padé approximants with a complexity ofO(n²). The algorithm succeeds also in the non-normal and non-perfect case, but with a possibility of an increase in the cost complexity.Supported in part by NSERC grant No. A8035.Partially supported by NSERC operating grant No. 6194.     

On the Monotone Error Rule for Parameter Choice in Iterative and Continuous Regularization Methods

We consider in Hilbert spaces linear ill-posed problems Ax = y with noisy data y satisfying y –y. Regularized approximations x _r to the minimum-norm solution x of Ax = y are constructed by continuous regularization methods or by iterative methods. For the choice of the regularization parameter r (the stopping index n in iterative methods) the following monotone error rule (ME rule) is used: we choose r = r _ME (n = n _ME) as the largest r-value with the guaranteed monotonical decrease of the error x _r – x for r [0, r _ME] (x _n – x <#60; x _n–1 – x for n = 1, 2, ..., n _ME). Main attention is paid to iterative methods of gradient type and to nonstationary implicit iteration methods. As shown, the ME rule leads for many methods to order optimal error bounds. Comparisons with other rules for the choice of the stopping index are made and numerical examples are given.This revised version was published online in October 2005 with corrections to the Cover Date.     

Comparison theorems for nonlinear evolution equations and their application to the investigation of the dependence of the solutions of the equation ut=δϕ(u) ON ϕ and on the boundary conditions

One obtains estimates of the form, whereu. are generalized solutions of the equationsdu/dt=Au, du/dt=Bu whileA, B are non-linear,m-dissipative operators in a Banach space, and there exists an operatorP:D(A)D(B), such thatPw · W+BPw –Aw, uniformly on some setw. These results are applied to the investigation of the dependence of the solutions of the Cauchy, Dirichlet problems and of the problem with the boundary condition –du/dn=(u) for the equation u₁=(u) on the continuous nondecreasing functions and.Translated from Problemy Matematicheskogo Analiza, No. 9, pp. 183–198, 1984.The author is sincerely grateful to O. A. Ladyzhenskaya and N. N. Ural'tseva for their interest in this paper and for useful discussions.     

Least squares algorithms for finding solutions of overdetermined linear equations which minimize error in an abstract norm

In this paper overdetermined linear equations inn unknowns are considered. Ifm is the number of equations, let _m be equipped with a smooth strictly convex norm, ·. Algorithms for finding best-fit solutions of the system which minimize the ·-error are given. The algorithms are iterative and in particular apply to the important case where _m is given thel ^p -norms, (1<p<). The algorithms consist of obtaining a least square solution i.e. carrying out an orthogonal projection at each stage of the iteration and solving a non-linear equation in a single real variable. The convergence of the algorithms are proved in the paper.