A note on contractivity in the numerical solution of initial value problems

This paper concerns the stability analysis of numerical methods for approximating the solutions to (stiff) initial value problems. Our analysis includes the case of (nonlinear) systems of differential equations that are essentially more general than the classical test equationU=U, with a complex constant.We explore the relation between two stability concepts, viz. the concepts of contractivity and weak contractivity.General Runge-Kutta methods, one-stage Rosenbrock methods and a notable rational Runge-Kutta method are analysed in some detail.     

Stability analysis of numerical methods for delay differential equations

Summary This paper deals with the stability analysis of step-by-step methods for the numerical solution of delay differential equations. We focus on the behaviour of such methods when they are applied to the linear testproblemU(t)=U(t)+U(t–) with >0 and , complex. A general theorem is presented which can be used to obtain complete characterizations of the stability regions of these methods.     

Ultrafilter translations

We develop a method for extending results about ultrafilters into a more general setting. In this paper we shall be mainly concerned with applications to cardinality logics. For example, assumingV=L, Gödel's Axiom of Constructibility, we prove that if > then the logic with the quantifier there exist many is (,)-compact if and only if either is weakly compact or is singular of cofinality<. As a corollary, for every infinite cardinals and , there exists a (,)-compact non-(,)-compact logic if and only if either < orcf<cf or < is weakly compact.Counterexamples are given showing that the above statements may fail, ifV=L is not assumed.However, without special assumptions, analogous results are obtained for the stronger notion of [,]-compactness.     

The dependence of critical parameter bounds on the monotonicity of a Newton sequence

Summary A new method is proposed for the inclusion of the critical parameter ^* of some convex operator equationu=Tu (appearing e.g. in thermal explosion theory). It is based on the fact that for a fixed Newton's method starting with a suitable subsolution is not monotonically if and only if >^*. Several numerical examples arising from nonlinear boundary value problems illustrate the efficiency of the method.     

A new interpolation procedure for adapting Runge-Kutta methods to delay differential equations

This paper deals with adapting Runge-Kutta methods to differential equations with a lagging argument. A new interpolation procedure is introduced which leads to numerical processes that satisfy an important asymptotic stability condition related to the class of testproblemsU(t)=U(t)+U(t–) with , C, Re()<–||, and >0. Ifc _i denotes theith abscissa of a given Runge-Kutta method, then in thenth stept _n–1t _n :=t _n–1+h of the numerical process our interpolation procedure computes an approximation toU(t _n–1+c _i h–) from approximations that have already been generated by the process at pointst _j–1+c _i h(j=1,2,3,...). For two of these new processes and a standard process we shall consider the convergence behaviour in an actual application to a given, stiff problem.     

Nonnegative Hall Polynomials

The number of subgroups of type and cotype in a finite abelian p-group of type is a polynomialg with integral coefficients. We prove g has nonnegative coefficients for all partitions and if and only if no two parts of differ by more than one. Necessity follows from a few simple facts about Hall-Littlewood symmetric functions; sufficiency relies on properties of certain order-preserving surjections that associate to each subgroup a vector dominated componentwise by . The nonzero components of (H) are the parts of , the type of H; if no two parts of differ by more than one, the nonzero components of – (H) are the parts of , the cotype of H. In fact, we provide an order-theoretic characterization of those isomorphism types of finite abelian p-groups all of whose Hall polynomials have nonnegative coefficients.     

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.     

On the stability properties of Brown's multistep multiderivative methods

Summary Brown introducedk-step methods usingl derivatives. We investigate for whichk andl the methods are stable or unstable. It is seen that to anyl the method becomes unstable fork large enough. All methods withk2(l+1) are stable. Fork=1,2,..., 18 there exists a _k such that the methods are stable for anyl _k and unstable for anyl < _k . The _k are given.     

The dynamics of perturbations of the contracting Lorenz attractor

We show here that by modifying the eigenvalues ₂ < ₃ < 0 < ₁ of the geometric Lorenz attractor, replacing the usualexpanding condition ₃+₁ > 0 by acontracting condition ₃+₁ < 0, we can obtain vector fields exhibiting transitive non-hyperbolic attractors which are persistent in the following measure theoretical sense: They correspond to a positive Lebesgue measure set in a twoparameter space. Actually, there is a codimension-two submanifold in the space of all vector fields, whose elements are full density points for the set of vector fields that exhibit a contracting Lorenz-like attractor in generic two parameter families through them. On the other hand, for an open and dense set of perturbations, the attractor breaks into one or at most two attracting periodic orbits, the singularity, a hyperbolic set and a set of wandering orbits linking these objects.     

Conjugation for polynomial mappings

We consider Keller's functions, namely polynomial functionsf:C ⁿ C ⁿ with detf(x)=1 at allx C ⁿ. Keller conjectured that they are all bijective and have polynomial inverses. The problem is still open.Without loss of generality assumef(0)=0 andf'(0)=I. We study the existence of certain mappingsh , > 1, defined by power series in a ball with center at the origin, such thath(0)=I andh (f(x))=h (x). So eachh conjugates f to its linear part I in a ball where it is injective.We conjecture that for Keller's functionsf of the homogeneous formf(x)=x +g(x),g(sx)=s ^dg(x),g(x)ⁿ=0,xC ⁿ,sC the conjugationh for f is anentire function.     

Solving a nonlinear spectral problem for a matrix

This paper examines the solving of the eigenvalue problem for a matrix M () with a nonlinear occurrence of the spectral parameter. Two methods are suggested for replacing the equation dat M()=0 by a scalar equationf()=0. Here the functionf() is not written formally, but a rule for computingf() at a fixed point of the domain in which the desired roots lie is indicated. Müller's method is used to solve the equationf()=0. The eigenvalue found is refined by Newton's method based on the normalized expansion of matrix M(), and the linearly independent vectors corresponding to it are computed. An ALGOL program and test examples are presented.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 58, pp. 54–66, 1976.     

On a fixed-point theorem on locally convex linear topological spaces

Iff is a self mapping on a closed convex subsetK of a separated quasicomplete locally convex linear topological spaceE such that (i)E is strictly convex, (ii)f (K) is contained in a compact subset ofK and (iii)f satisfies a contraction condition, then it is shown that for eachxK, the sequence of {U ⁿ (x)} _n =1 of iterates, whereU KK is defined byU (y)=f(y)+(1-) y, yK, converges to a fixed point off.     

Spectral parameter asymptotics of the Weyl solutions of Sturm-Liouville equations

In the paper one investigates the dependence of Weyl's solution ,)=c(,)+n()s(,) of the Sturm-Liouville equation y+q()y=²y on the spectral parameter . Under the condition that the potential q is bounded from below and q()exp(c₀+c[in1 ¦¦), it is proved for {ie217-01} for any positive values and A. If q()>1 and {ie217-02} for all >0, then in the semiplane >0 the Weyl solution (, ) is obtained from the Weyl solution (,x) is obtained from the Weyl solution e^ix with zero potential, with the aid of a generalization of B. Ya Levin's transformation operators.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova Akademii Nauk SSSR, Vol. 170, pp. 184–206, 1989.I express my sincere gratitude to L. A. Pastur and I. V. Ostrovskii for valuable advice and discussions.     

Isometries of symmetric operator spaces associated with AFD factors of typeII and symmetric vector-valued spaces

If is a surjective isometry of the separable symmetric operator spaceE(M, ) associated with the approximately finite-dimensional semifinite factorM and if · _E(M,) is not proportional to · _{L 2}, then there exist a unitary operatorUM and a Jordan automorphismJ ofM such that(x)=UJ(x) for allxME(M, ). We characterize also surjective isometries of vector-valued symmetric spacesF((0, 1), E(M, )).Research supported by the Australian Research Council     

The spectral bound of Schrödinger operators

Let V: R ^N [0, ] be a measurable function, and >0 be a parameter. We consider the behaviour of the spectral bound of the operator 1/2–V as a function of . In particular, we give a formula for the limiting value as , in terms of the integrals of V over subsets of R ^N on which the Laplacian with Dirichlet boundary conditions has prescribed values. We also consider the question whether this limiting value is attained for finite .     

Nonlinear eigenvalue approximation

Summary For each in some domainD in the complex plane, letF() be a linear, compact operator on a Banach spaceX and letF be holomorphic in . Assuming that there is a so thatI–F() is not one-to-one, we examine two local methods for approximating the nonlinear eigenvalue . In the Newton method the smallest eigenvalue of the operator pencil [I–F(),F()] is used as increment. We show that under suitable hypotheses the sequence of Newton iterates is locally, quadratically convergent. Second, suppose 0 is an eigenvalue of the operator pencil [I–F(),I] with algebraic multiplicitym. For fixed leth() denote the arithmetic mean of them eigenvalues of the pencil [I–F(),I] which are closest to 0. Thenh is holomorphic in a neighborhood of andh()=0. Under suitable hypotheses the classical Muller's method applied toh converges locally with order approximately 1.84.     

Cardinal Invariants of λ-Topologies

We study the basic cardinal-valued invariants of C (X) such as weight, density, network weight, i-weight, and tightness, where C (X) is the space of all continuous real functions on X in the -topology.     

Reinforcement problems for elliptic equations and variational inequalities

Summary Consider the Dirichlet problem for an elliptic equation in a domain , with coefficients having discontinuity on a surface . Suppose divides into _{1 2}(₂ the inner core), the thickness of ₁ is of order of magnitude , and the modulus of ellipticity in ₁ is of order magnitude ₁. The asymptotic behavior of the solution is studied as 0, ₁ 0, provided lim (/₁) exists. Other questions of this type are studied both for elliptic equations and for elliptic variational inequalities.The second author is partially supported by National Science Foundation Grant 7406375 A01. The third author is partially supported by National Science Foundation Grant MC575-21416 A01.     

A Proof of the Jungnickel-Tonchev Conjecture on Quasi-Multiple Quasi-Symmetric Designs

A proof of the following conjecture of Jungnickel and Tonchev on quasi-multiple quasi-symmetric designs is given: Let D be a design whose parameter set (v,b,r,k,) equals (v,sv,sk,k, s) for some positive integer s and for some integers v,k, that satisfy (v-1) = k(k-1) (that is, these integers satisfy the parametric feasibility conditions for a symmetric (v,k,)-design). Further assume that D is a quasi-symmetric design, that is D has at most two block intersection numbers. If (k, (s-1)) = 1, then the only way D can be constructed is by taking multiple copies of a symmetric (v,k, )-design.     

Une méthode numérique pour le calcul des points de retournement. Application à un problème aux limites non-linéaire

Summary A finite element method (P₁) with numerical integration for approximating the boundary value problem –u=e ^u is considered. It is shown that the discrete problem has a solution branch (with turning point) which converges uniformely to a solution branch of the continuous problem. Error estimates are given; for example it is found that , >0, where ₀ and _h ⁰ are critical values of the parameter for continuous and discrete problems.