The Fast Numerical Solution of Very Large Elliptic Difference Schemes

Let L_kv_k = g_k be a system of difference equations discretizingan elliptic boundary value problem. Assume the system to be"very large", that means that the number of unknowns exceedsthe capacity of storage. We present a method for solving theproblem with much less storage requirement. For two-dimensionalproblems the size of the needed storage decreases from O(h^–2)to (or even O(h^–5/4)). The computational work increasesonly by a factor about six. The technique can be generalizedto nonlinear problems. The algorithm is also useful for computerswith a small number of parallel processors.     

On Positive Multipeak Solutions of a Nonlinear Elliptic Problem

In this paper we continue our investigation in [5, 7, 8] onmultipeak solutions to the problem –²u+u=Q(x)|u|^q–2u, xR^N, uH¹(R^N) (1.1) where = ^N_i=1²/x²_i is the Laplace operator in R^N, 2 < q < for N = 1, 2, 2 < q < 2N/(N–2) for N3, and Q(x)is a bounded positive continuous function on R^N satisfying thefollowing conditions. (Q₁) Q has a strict local minimum at some point x₀R^N, that is,for some > 0 Q(x)>Q(x₀) for all 0 < |x–x₀| < . (Q₂) There are constants C, > 0 such that |Q(x)–Q(y)|C|x–y| for all |x–x₀| , |y–y₀| . Our aim here is to show that corresponding to each strict localminimum point x₀ of Q(x) in R^N, and for each positive integerk, (1.1) has a positive solution with k-peaks concentratingnear x₀, provided is sufficiently small, that is, a solutionwith k-maximum points converging to x₀, while vanishing as 0 everywhere else in R^N.     

A Note on the WKB Method for Difference Equations

We consider the asymptotic solution of the second-order differenceequation y_n + 1 –2y_n + y_n–1 + Q_ny_n = 0, where Q_n= N^–Q(n/N), 0 < < 2, Q(s) being a differentiablefunction of s, and N a large parameter such that Q(n/N) variesby order unity as n varies by order N. A discrete WKB methodis proposed, the form of the asymptotic expansion being similarto that used in the conventional WKB method. A particular Q(s)is studied, for which results of the discrete WKB method arein agreement with the results from the approach due to Bremmer(1951).     

Asymptotics for the Heat Content of a Planar Region with a Fractal Polygonal Boundary

Let k 3 be an integer. For 0<s<1, let D_s R² be the setthat is constructed iteratively as follows. Take a regular openk-gon with sides of unit length, attach regular open k-gonswith sides of length s to the middles of the edges, and so on.At each stage of the iteration the k-gons that are added area factor s smaller than the previous generation and are attachedto the outer edges of the family grown so far. The set D_s isdefined to be the interior of the closure of the union of allthe k-gons. It is easy to see that there must exist some s_k> 0 such that no k-gons overlap if and only if 0 < s s_k. We derive an explicit formula for s_k. The set D_s is open, bounded, connected and has a fractal polygonalboundary. Let denote the heat content of D_s at time t when D_s initially has temperature 0and D_s is kept at temperature 1. We derive the complete short-timeexpansion of up to terms that are exponentially small in 1/t. It turns out that there arethree regimes, corresponding to 0<s<1/(k–1), s=1/(k–1),and 1/(k–1)<s s_k. For s 1/(k–1) the expansionhas the form where p_s is a log (1/s²)-periodic function, d_s=log (k–1)/log(1/s) is a similarity dimension, A_s and B are constants relatedto the edges and vertices, respectively, of D_s, and r_s is anerror exponent. For s=1/(k–1), the t^1/2-term carries anadditional log t. 1991 Mathematics Subject Classification: 11D25,11G05, 14G05.     

On the Local and Superlinear Convergence of Quasi-Newton Methods

This paper presents a local convergence analysis for severalwell-known quasi-Newton methods when used, without line searches,in an iteration of the form to solve for x^* such that Fx^* = 0. The basic idea behind theproofs is that under certain reasonable conditions on x_o, Fand x_o, the errors in the sequence of approximations {H_k} toF'(x^*)^–1 can be shown to be of bounded deterioration inthat these errors, while not ensured to decrease, can increaseonly in a controlled way. Despite the fact that H_k is not shownto approach F'(x^*)^–1, the methods considered, includingthose based on the single-rank Broyden and double-rank Davidon-Fletcher-Powellformulae, generate locally Q-superlinearly convergent sequences{x_k}.     

Shapiro's Cyclic Sum

Shapiro's cyclic sum is defined by , If K is the cone in Rⁿ of points withnon-negative coordinates, it is shown that the minimum of Ein K is a fixed point of T², where T is the non-linear operatordefined by (Tx)_i = x_n–i+1/(x_n–i+2 + x_n–i+3)²for i = 1,2,...,n. It is conjectured that Tx = S^kx, where Sis the shift operator in Rⁿ, and a proof is given under someadditional hypotheses. One of the consequences is a simple proofthat at the minimum point, a_i(x) = a_n–i+1–k(x) fori = 1,2,...,n.     

The streamline-diffusion method for nonconforming Qrot1 elements on rectangular tensor-product meshes

When the streamline–diffusion finite element method isapplied to convection–diffusion problems using nonconformingtrial spaces, it has previously been observed that stabilityand convergence problems may occur. It has consequently beenproposed that certain jump terms should be added to the bilinearform to obtain the same stability and convergence behaviouras in the conforming case. The analysis in this paper showsthat for the Q^rot₁ 1 element on rectangular shape-regular tensor-productmeshes, no jump terms are needed to stabilize the method. Inthis case moreover, for smooth solutions we derive in the streamline–diffusionnorm convergence of order h^3/2 (uniformly in the diffusion coefficientof the problem), where h is the mesh diameter. (This estimateis already known for the conforming case.) Our analysis alsoshows that similar stability and convergence results fail tohold true for analogous piecewise linear nonconforming elements.     

Long-Time Behavior of Solutions of the Fast Diffusion Equations with Critical Absorption Terms

This paper is devoted to the long-time behavior of solutionsto the Cauchy problem of the porous medium equation u_t = (u^m)– u^p in Rⁿ x (0,) with (1 – 2/n)₊ < m < 1and the critical exponent p = m + 2/n. For the strictly positiveinitial data u(x,0) = O(1 + |x|)^–k with n + mn(2 –n + nm)/(2[2 – m + mn(1 – m)]) k < 2/(1 –m), we prove that the solution of the above Cauchy problem convergesto a fundamental solution of u_t = (u^m) with an additional logarithmicanomalous decay exponent in time as t .     

Ganea's Conjecture on Lusternik-Schnirelmann Category

A series of complexes Q_p indexed by all primes p is constructedwith catQ_p=2 and catQ_pxSⁿ=2 for either n2 or n=1 and p=2. Thisdisproves Ganea's conjecture on Lusternik–Schnirelmann(LS) category. 1991 Mathematics Subject Classification 55M30.     

Stability in the Sense of Bounded Average Power

An R^m-valued sequence (x_k): = (x_k : k = 1, 2, ...), e.g. generatedrecursively by x_k = f_k (x_k–_k, U_k), is called ‘averagepth power bounded’ if (1/K) is bounded uniformly in K= 1, 2,.... (The case p = 2 may correspond to ‘power’in the physical sense.) This is a notion of stability. Givenestimates of the form: f_k (x, u) < a x + ¶ _k conditionsare obtained on the coefficient sequence (a_k) and the inputestimates e_k:=¶_k (u_k) which ensure this form of stabilityfor the output (x_k). In particular, a condition (utilized inan application to adaptive control) is obtained which imposes(i) a bound b on (a_k) and a ‘sparsity measure’ m(K) on #{kK: a_k>} as K ( >1) (ii) average pth power boundednesson (e_k), and (iii) a growth condition on (e_k) related to b andm (•). This condition is sharp.     

Necessary and Sufficient Conditions for Exponential Stability and Ultimate Boundedness of Systems Governed by Stochastic Partial Differential Equations

Consider the following infinite dimensional stochastic evolutionequation over some Hilbert space H with norm |·|: It is proved that under certain mild assumptions, the strongsolution X_t(x₀)VHV*, t 0, is mean square exponentially stableif and only if there exists a Lyapunov functional (·,·):HxR₊R¹ which satisfies the following conditions: (i)c₁|x|²–k₁e^–µ₁t(x,t)c₂|x|²+k₂+k₂e^–µ₂t; (ii) L(x,t)–c₃(x,t)+k₃e^–µ₃t, xV, t0; where L is the infinitesimal generator of the Markov processX_t and c_i, k_i, µ_i, i = 1, 2, 3, are positive constants.As a by-product, the characterization of exponential ultimateboundedness of the strong solution is established as the nulldecay rates (that is, µ_i = 0) are considered.     

A multigrid-type method for thin plate spline interpolation on a circle

We consider the problem of thin plate spline interpolation ton equally spaced points on a circle, where the number of datapoints is sufficiently large for work of O(n³ to be unacceptable.We develop an iterative multigrid-type method, each iterationcomprising ngrid stages, and n being an integer multiple of2^ngrid–1. We let the first grid, V₁ be the full set ofdata points, V say, and each subsequent (coarser) grid, V_k,k=2, 3,...,ngrid, contain exactly half of the data points ofthe preceding (finer) grid, these data points being equallyspaced. At each stage of the iteration, we correct our current approximationto the thin plate spline interpolant by an estimate of the interpolantto the current residuals on V_k, where the correction is constructedfrom Lagrange functions of interpolation on small local subsetsof p data points in V_k. When the coarsest grid is reached, however,then the interpolation problem is solved exactly on its q=n/2^ngrid–1points. The iterative process continues until the maximum residualdoes not exceed a specified tolerance. Each iteration has the effect of premultiplying the vector ofresiduals by an n x n matrix R, and thus convergence will dependupon the spectral radius, (R), of this matrix. We investigatethe dependence of the spectral radius on the values of n, p,and q. In all the cases we have considered, we find (R) <<1, and thus rapid convergence is assured.     

Determination of a Convex Body from Minkowski Sums of its Projections

For a convex body K in R^d and 1 K d – 1, let P_K (K)be the Minkowski sum (average) of all orthogonal projectionsof K onto k-dimensional subspaces of R^d. It is Known that theoperator P_k is injective if kd/2, k=3 for all d, and if k =2, d 14. It is shown that P_2k (K) determines a convex body K among allcentrally symmetric convex bodies and P_2k+1(K) determines aconvex body K among all bodies of constant width. Correspondingstability results are also given. Furthermore, it is shown thatany convex body K is determined by the two sets P_k (K) and P_k'(K) if 1 < k < k'. Concerning the range of P_k , 1 k d–2, it is shown that its closure (in the Hausdorff-metric)does not contain any polytopes other than singletons.     

Hypersurfaces in a Unit Sphere Sn+1(1) with Constant Scalar Curvature

The paper considers n-dimensional hypersurfaces with constantscalar curvature of a unit sphere S^n–1(1). The hypersurfaceS^k(c₁)xS^n–k(c₂) in a unit sphere Sⁿ⁺¹(1) is characterized,and it is shown that there exist many compact hypersurfaceswith constant scalar curvature in a unit sphere Sⁿ⁺¹(1) whichare not congruent to each other in it. In particular, it isproved that if M is an n-dimensional (n > 3) complete locallyconformally flat hypersurface with constant scalar curvaturen(n–1)r in a unit sphere Sⁿ⁺¹(1), then r > 1–2/n,and (1) when r (n–2)/(n–1), if then M is isometric to S¹xS^n–1(c),where S is the squared norm of the second fundamental form ofM; (2) there are no complete hypersurfaces in Sⁿ⁺¹(1) with constantscalar curvature n(n–1)r and with two distinct principalcurvatures, one of which is simple, such that r = (n–2)/(n–1)and      

A non-conforming finite element method with anisotropic mesh grading for the Stokes problem in domains with edges

The solution of the Stokes problem in three-dimensional domainswith edges has anisotropic singular behaviour which is treatednumerically by using anisotropic finite element meshes. Thevelocity is approximated by Crouzeix–Raviart (nonconformingP₁ ) elements and the pressure by piecewise constants. Thismethod is stable for general meshes (without minimal or maximalangle condition). Denoting by N_e the number of elements in themesh, the interpolation and consistency errors are of the optimalorder h N_e^–1/3 which is proved for tensor product meshes.As a by-product, we analyse also nonconforming prismatic elementswith P₁ [oplus ] span {x₃²} as the local space for the velocitywhere x₃ is the direction of the edge.     

Multipliers of Qs spaces and the corona theorem

For s>0, let Q_s be the space of all analytic functions onthe unit disc such that |f'(z)|²(1–|z|²)^s dA(z) isan s-Carleson measure. Here we prove that the corona theoremholds for the algebra of pointwise multipliers of Q_s.     

The Direct Limits Of The Banach-mazur Compacta

Let 1 p . For each n-dimensional Banach space E = (E, || ·||), we define a norm || · ||_p on E x R as follows: [formula] It is shown that the correspondence (E, || · ||) (Ex R, || · ||_p) defines a topological embedding of oneBanach–Mazur compactum, BM(n), into another, BM(n 1),and hence we obtain a tower of Banach–Mazur compacta:BM(1) BM(2) BM(3) ···. Let BM_p be thedirect limit of this tower. We prove that BM_p is homeomorphicto Q = dir lim Qⁿ, where Q = [0, 1] is the Hilbert cube. 1991Mathematics Subject Classification 46B04, 46B20, 52A21, 57N20,54H15.     

Exotic smooth structures on Formula

Motivated by Stipsicz and Szabó's exotic 4-manifoldswith b₂⁺ = 3 and b₂^– = 8, we construct a family of simplyconnected smooth 4-manifolds with b₂⁺ = 3 and b₂^– = 8.As a corollary, we conclude that the topological 4-manifold     

(k, l)-Algebraic Stability of Runge-Kutta Methods

For ordinary differential equations satisfying a one-sided Lipschitzcondition with Lipschitz constant v, the solutions satisfy with l=hv, so that, in the case of Runge-Kutta methods, estimatesof the form ||y_n||²k(l)||y_n–1||² are desirable. Burrage(1986) has investigated the behaviour of the error-boundingfunction k for positive l for the family of s-stage Gauss methodsof order 2s, and has shown that k(l)=exp 2l+O(l³) (l0) for s3.In this paper, we extend the analysis of k to any irreduciblealgebraically stable Runge-Kutta method, and obtain resultsabout the maximum order of k as an approximation to exp 2l.As a particular example, we investigate the function k for allalgebraically stable methods of order 2s–1.     

Serre's Theorem on the Cohomology Algebra of a p-Group

The purpose of this note is to give a proof of a theorem ofSerre, which states that if G is a p-group which is not elementaryabelian, then there exist an integer m and non-zero elementsx₁, ..., x_m H¹ (G, Z/p) such that with ß the Bockstein homomorphism. Denote by m_G thesmallest integer m satisfying the above property. The theoremwas originally proved by Serre [5], without any bound on m_G.Later, in [2], Kroll showed that m_G p^k – 1, with k =dim_Z/pH¹ (G, Z/p). Serre, in [6], also showed that m_G (p^k –1)/(p – 1). In [3], using the Evens norm map, Okuyamaand Sasaki gave a proof with a slight improvement on Serre'sbound; it follows from their proof (see, for example, [1, Theorem4.7.3]) that m_G (p + 1)p^k–2. However, m_G can be sharpenedfurther, as we see below. For convenience, write H^*(G, Z/p) = H^*(G). For every x_i H¹(G),set 1991 Mathematics SubjectClassification 20J06.