A Generalised Skolem-Mahler-lech Theorem for Affine Varieties

The Skolem–Mahler–Lech theorem states that if f(n)is a sequence given by a linear recurrence over a field of characteristic0, then the set of m such that f(m) is equal to 0 is the unionof a finite number of arithmetic progressions in m 0 and afinite set. We prove that if X is a subvariety of an affinevariety Y over a field of characteristic 0 and q is a pointin Y, and is an automorphism of Y, then the set of m such that^m(q) lies in X is a union of a finite number of complete doubly-infinitearithmetic progressions and a finite set. We show that thisis a generalisation of the Skolem–Mahler–Lech theorem.     

A Mean Value Property of Poly-Temperatures on a Strip Domain

We consider the iterates of the heat operator on Rⁿ⁺¹={(X, t); X=(x₁, x₂, ..., x_n)Rⁿ, tR}. Let Rⁿ⁺¹ be a domain,and let m1 be an integer. A lower semi-continuous and locallyintegrable function u on is called a poly-supertemperatureof degree m if (–H)^mu0 on (in the sense of distribution). If u and –u are both poly-supertemperatures of degreem, then u is called a poly-temperature of degree m. Since His hypoelliptic, every poly-temperature belongs to C(), andhence (–H)^m u(X, t)=0 (X, t). For the case m=1, we simply call the functions the supertemperatureand the temperature. In this paper, we characterise a poly-temperature and a poly-supertemperatureon a strip D={(X, t);XRⁿ, 0<t<T} by an integral mean on a hyperplane. To state our result precisely,we define a mean A[·, ·]. This plays an essentialrole in our argument.     

Martingales in the Ok Corral

In the model of the OK Corral formulated by Williams and McIlroy[2]: ‘Two lines of gunmen face each other, there beinginitially m on one side, n on the other. Each person involvedis a hopeless shot, but keeps firing at the enemy until eitherhe himself is killed or there is no one left on the other side.’They are interested in the number S of survivors when the shootingceases, and the surprising result, for which they give bothnumerical and heuristic evidence, is that when m=n, [formula] as m (where E denotes expectation). It is the occurrence of this curious power of m, rather thanany application to real gunfights, which makes the Williams–McIlroyprocess of interest. The purpose of this paper is to give anessentially elementary proof of the fact that if m and n arenot too different, then S/m^3/4 has a simple asymptotic distributionthat leads at once to (1.1). 1991 Mathematics Subject Classification60F05.     

Limits Along Parallel Lines and the Classical Fine Topology

The fine topology on Rⁿ (n2) is the coarsest topology for whichall superharmonic functions on Rⁿ are continuous. We refer toDoob [11, 1.XI] for its basic properties and its relationshipto the notion of thinness. This paper presents several theoremsrelating the fine topology to limits of functions along parallellines. (Results of this nature for the minimal fine topologyhave been given by Doob – see [10, Theorem 3.1] or [11,1.XII.23] – and the second author [15].) In particular,we will establish improvements and generalizations of resultsof Lusin and Privalov [18], Evans [12], Rudin [20], Bagemihland Seidel [6], Schneider [21], Berman [7], and Armitage andNelson [4], and will also solve a problem posed by the latterauthors. An early version of our first result is due to Evans [12, p.234], who proved that, if u is a superharmonic function on R³,then there is a set ER²x{0}, of two-dimensional measure 0, suchthat u(x, y,·) is continuous on R whenever (x, y, 0)E.We denote a typical point of Rⁿ by X=(X' x), where X'R^n–1and xR. Let :RⁿR^n–1x{0} denote the projection map givenby (X', x) = (X', 0). For any function f:Rⁿ[–, +] andpoint X we define the vertical and fine cluster sets of f atX respectively by C_V(f;X)={l[–, +]: there is a sequence (t_m) of numbersin R\{x} such that t_mx and f(X', t_m)l}| and C_F(f;X)={l[–, +]: for each neighbourhood N of l in [–,+], the set f^–1(N) is non-thin at X}. Sets which are open in the fine topology will be called finelyopen, and functions which are continuous with respect to thefine topology will be called finely continuous. Corollary 1(ii)below is an improvement of Evans' result.     

Shape-preserving interpolation by G1 and G2 PH quintic splines

The interpolation of a planar sequence of points p₀, ..., p_Nby shape-preserving G¹ or G² PH quintic splines with specifiedend conditions is considered. The shape-preservation propertyis secured by adjusting ‘tension’ parameters thatarise upon relaxing parametric continuity to geometric continuity.In the G² case, the PH spline construction is based on applyingNewton–Raphson iterations to a global system of equations,commencing with a suitable initialization strategy—thisgeneralizes the construction described previously in NumericalAlgorithms 27, 35–60 (2001). As a simpler and cheaperalternative, a shape-preserving G¹ PH quintic spline schemeis also introduced. Although the order of continuity is lower,this has the advantage of allowing construction through purelylocal equations.     

Error expansion for a first-order upwind difference scheme applied to a model convection-diffusion problem

A singularly perturbed convection–diffusion problem isconsidered. The problem is discretized using a simple first-orderupwind difference scheme on general meshes. We derive an expansionof the error of the scheme that enables uniform error boundswith respect to the perturbation parameter in the discrete maximumnorm for both a defect correction method and the Richardsonextrapolation technique. This generalizes and simplifies resultsobtained in earlier publications by Fröhner et al.(2001,Numer. Algorithms, 26, 281–299) and by Natividad &Stynes (2003, Appl. Numer. Math., 45, 315–329). Numericalexperiments complement our theoretical results.     

Time-optimal control problem for parabolic equations with control constraints and infinite number of variables

A time optimal control problem for parabolic equations withan infinite number of variables is considered. A time optimalcontrol problem is replaced by an equivalent one with a performanceindex in the form of integral form. Constraints on controlsare assumed. To obtain the optimality conditions for the Neumannproblem, the generalization of the Dubovitskii–Milyutintheorem given by Walczak (1984, Acta Universitatis LodziensisFolia Mathematica, 187–196; 1984, J. Optim. Theor. Appl.,42, 561–582) was applied.     

Semi-discrete central-upwind schemes with reduced dissipation for Hamilton-Jacobi equations

We introduce a new family of Godunov-type semi-discrete centralschemes for multidimensional Hamilton–Jacobi equations.These schemes are a less dissipative generalization of the central-upwindschemes that have been recently proposed in Kurganov, Noelleand Petrova (2001, SIAM J. Sci. Comput., 23, pp. 707–740).We provide the details of the new family of methods in one,two, and three space dimensions, and then verify their expectedlow-dissipative property in a variety of examples.     

An Euler-Type Version of the Local Steiner Formula for Convex Bodies

The purpose of this note is to establish a new version of thelocal Steiner formula and to give an application to convex bodiesof constant width. This variant of the Steiner formula generalizesresults of Hann [3] and Hug [6], who use much less elementarytechniques than the methods of this paper. In fact, Hann askedfor a simpler proof of these results [4, Problem 2, p. 900].We remark that our formula can be considered as a Euclideananalogue of a spherical result proved in [2, p. 46], and thatour method can also be applied in hyperbolic space. For some remarks on related formulas in certain two-dimensionalMinkowski spaces, see Hann [5, p. 363]. For further information about the notions used below, we referto Schneider's book [9]. Let Kⁿ be the set of all convex bodiesin Euclidean space Rⁿ, that is, the set of all compact, convex,non-empty subsets of Rⁿ. Let S^n–1 be the unit sphere.For KKⁿ, let NorK be the set of all support elements of K, thatis, the pairs (x, u)RⁿxS^n–1 such that x is a boundarypoint of K and u is an outer unit normal vector of K at thepoint x. The support measures (or generalized curvature measures)of K, denoted by ₀(K.), ..., _n–1(K.), are the unique Borelmeasures on RⁿxS^n–1 that are concentrated on NorK andsatisfy [formula] for all integrable functions f:RⁿR; here denotes the Lebesguemeasure on Rⁿ. Equation (1), which is a consequence and a slightgeneralization of Theorem 4.2.1 in Schneider [9], is calledthe local Steiner formula. Our main result is the following.1991 Mathematics Subject Classification 52A20, 52A38, 52A55.     

Stability of uniformly Morse-Smale gradient-like numerical methods for flows

In Garay (1996, Numer. Math., 72, 449–479) and Li (1997b,SIAM J. Math. Anal., 28, 381–388), it was shown that thequalitative properties of a Morse–Smale gradient-likeflow are preserved by its numerical approximations. In thispaper, we show that the qualitative properties of a family ofuniformly Morse–Smale gradient-like numerical methodsare preserved by the approximated flow. The techniques usedin the study of the structural stability theorem for diffeomorphismsare the main tools for this work.     

Oscillation Properties of a Semi-discrete Approximation to the Beam Equation with a Second Order Term

The solution of the equation w(x)u_tt+[p(x)u_xx]_xx–[p(x)u_x]_x=0, 0< x < L, t > 0, where it is assumed that w, p,and q are positive on the interval [0, L], is approximated bythe method of straight lines. The resulting approximation isa linear system of differential equations with coefficient matrixS. The matrix S is studied under a variety of boundary conditionswhich result in a conservative system. In all cases the matrixS is shown to be similar to an oscillation matrix.     

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 .     

Euler Characteristics of Real Varieties

In this paper we show how to associate to any real projectivealgebraic variety Z RP^n–1 a real polynomial F₁:Rⁿ,0 R, 0 with an algebraically isolated singularity, having theproperty that (Z) = (1 – deg (grad F₁), where deg (gradF₁ is the local real degree of the gradient grad F₁:Rⁿ, 0 Rⁿ,0. This degree can be computed algebraically by the method ofEisenbud and Levine, and Khimshiashvili [5]. The variety Z neednot be smooth. This leads to an expression for the Euler characteristic ofany compact algebraic subset of Rⁿ, and the link of a quasihomogeneousmapping f: Rⁿ, 0 Rⁿ, 0 again in terms of the local degree ofa gradient with algebraically isolated singularity. Similar expressions for the Euler characteristic of an arbitraryalgebraic subset of Rⁿ and the link of any polynomial map aregiven in terms of the degrees of algebraically finite gradientmaps. These maps do involve ‘sufficiently small’constants, but the degrees involved ar (theoretically, at least)algebraically computable.     

Kazhdan-Lusztig Cells, q-Schur Algebras and James' Conjecture

We consider the Dipper–James q-Schur algebra S_q(n, r)_k,defined over a field k and with parameter q 0. An understandingof the representation theory of this algebra is of considerableinterest in the representation theory of finite groups of Lietype and quantum groups; see, for example, [6] and [11]. Itis known that S_q(n, r)_k is a semisimple algebra if q is nota root of unity. Much more interesting is the case when S_q(n,r)_k is not semisimple. Then we have a corresponding decompositionmatrix which records the multiplicities of the simple modulesin certain ‘standard modules’ (or ‘Weyl modules’).A major unsolved problem is the explicit determination of thesedecomposition matrices.     

Optimal error estimates for the hp-version interior penalty discontinuous Galerkin finite element method

We consider the hp-version interior penalty discontinuous Galerkinfinite-element method (hp-DGFEM) for second-order linear reaction–diffusionequations. To the best of our knowledge, the sharpest knownerror bounds for the hp-DGFEM are due to Rivière et al.(1999,Comput. Geosci., 3, 337–360) and Houston et al.(2002,SIAM J. Numer. Anal., 99, 2133–2163). These are optimalwith respect to the meshsize h but suboptimal with respect tothe polynomial degree p by half an order of p. We present improvederror bounds in the energy norm, by introducing a new functionspace framework. More specifically, assuming that the solutionsbelong element-wise to an augmented Sobolev space, we deducefully hp-optimal error bounds.     

Hypersurface Singularities with 2-Dimensional Critical Locus

Consider an analytic germ f:(C^m, 0)(C, 0) (m3) whose criticallocus is a 2-dimensional complete intersection with an isolatedsingularity (icis). We prove that the homotopy type of the Milnorfiber of f is a bouquet of spheres, provided that the extendedcodimension of the germ f is finite. This result generalizesthe cases when the dimension of the critical locus is zero [8],respectively one [12]. Notice that if the critical locus isnot an icis, then the Milnor fiber, in general, is not homotopicallyequivalent to a wedge of spheres. For example, the Milnor fiberof the germ f:(C⁴, 0)(C, 0), defined by f(x₁, x₂, x₃, x₄) =x₁x₂x₃x₄ has the homotopy type of S¹xS¹xS¹. On the other hand,the finiteness of the extended codimension seems to be the rightgeneralization of the isolated singularity condition; see forexample [9–12, 17, 18]. In the last few years different types of ‘bouquet theorems’have appeared. Some of them deal with germs f:(X, x)(C, 0) wheref defines an isolated singularity. In some cases, similarlyto the Milnor case [8], F has the homotopy type of a bouquetof (dim X–1)-spheres, for example when X is an icis [2],or X is a complete intersection [5]. Moreover, in [13] Siersmaproved that F has a bouquet decomposition FF₀Sⁿ...Sⁿ (whereF₀ is the complex link of (X, x)), provided that both (X, x)and f have an isolated singularity. Actually, Siersma conjecturedand Tibr proved [16] a more general bouquet theorem for thecase when (X, x) is a stratified space and f defines an isolatedsingularity (in the sense of the stratified spaces). In thiscase F_iF_i, where the F_i are repeated suspensions of complexlinks of strata of X. (If (X, x) has the ‘Milnor property’,then the result has been proved by Lê; for details see[6].) In our situation, the space-germ (X, x) is smooth, but f hasbig singular locus. Surprisingly, for dim Sing f^–1(0)2,the Milnor fiber is again a bouquet (actually, a bouquet ofspheres, maybe of different dimensions). This result is in thespirit of Siersma's paper [12], where dim Sing f^–1(0)= 1. In that case, there is only a rather small topologicalobstruction for the Milnor fiber to be homotopically equivalentto a bouquet of spheres (as explained in Corollary 2.4). Inthe present paper, we attack the dim Sing f^–1(0) = 2 case.In our investigation some results of Zaharia are crucial [17,18].     

Binary Quadratic Forms with Large Discriminants and Sums of Two Squareful Numbers II

Let F = (F₁, ..., F_m) be an m-tuple of primitive positive binaryquadratic forms and let U_F(x) be the number of integers notexceeding x that can be represented simultaneously by all theforms F_j, j = 1, ... , m. Sharp upper and lower bounds for U_F(x)are given uniformly in the discriminants of the quadratic forms. As an application, a problem of Erds is considered. Let V(x)be the number of integers not exceeding x that are representableas a sum of two squareful numbers. Then V(x) = x(log x)^–+o(1)with = 1 – 2^–1/3 = 0.206....     

On the Representation of Primes by Cubic Polynomials in Two Variables

In an earlier paper (see Proc. London Math. Soc. (3) 84 (2002)257–288) we showed that an irreducible integral binarycubic form f(x, y) attains infinitely many prime values, providingthat it has no fixed prime divisor. We now extend this resultby showing that f(m, n) still attains infinitely many primevalues if m and n are restricted by arbitrary congruence conditions,providing that there is still no fixed prime divisor. Two immediate consequences for the solvability of diagonal cubicDiophantine equations are given. 2000 Mathematics Subject Classification11N32 (primary), 11N36, 11R44 (secondary).     

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.     

Multilevel flexible specification of the production function in health economics

^** Email: grassetti{at}stat.unipd.it^*** Email: e.gori{at}dss.uniud.it^**** Email: simona.minotti{at}unicatt.it Previous studies on hospitals' efficiency often refer to quiterestrictive functional forms for the technology (Aigner et al.,1977, J. Econom., 6, 21–37). In this paper, referringto a study about some hospitals in Lombardy, we formulate convenientcorrectives to a statistical model based on the translogarithmicfunction—the most widely used flexible functional form(Christensen et al., 1973, Rev. Econ. Stat., 55, 28–45).More specifically, in order to take into consideration the hierarchicalstructure of the data (as in Gori et al., 2002, Stat. Appl.,14, 247–275), we propose a multilevel model, ignoringfor the moment the one-side error specification, typical ofstochastic frontier analysis (Aigner et al., 1977, J. Econom.,6, 21–37). Given this simplification, however, we areeasily able to take into account some typical econometric problemsas, e.g. heteroscedasticity. The estimated production functioncan be used to identify the technical inefficiency of hospitals(as already seen in previous works), but also to draw some economicconsiderations about scale elasticity, scale efficiency andoptimal resource allocation of the productive units. We willshow, in fact, that for the translogarithmic specification itis possible to obtain the elasticity of the output (regardingan input) at hospital level as a weighted sum of elasticitiesat ward level. Analogous results can be achieved for scale elasticity,which measures how output changes in response to simultaneousinputs variation. In addition, referring to scale efficiencyand to optimal resource allocation, we will consider the resultsof Ray (1998, J. Prod. Anal., 11, 183–194) to our context.The interpretation of the results is surely an interesting administrativeinstrument for decision makers in order to analyse the productiveconditions of each hospital and its single wards and also todecide the preferable interventions.