A note on the design of hp-version interior penalty discontinuous Galerkin finite element methods for degenerate problems

^** Email: emmanuil.georgoulis{at}mcs.le.ac.uk^*** Email: al{at}maths.strath.ac.uk We consider a variant of the hp-version interior penalty discontinuousGalerkin finite element method (IP-DGFEM) for second-order problemsof degenerate type. We do not assume uniform ellipticity ofthe diffusion tensor. Moreover, diffusion tensors of arbitraryform are covered in the theory presented. A new, refined recipefor the choice of the discontinuity-penalization parameter (thatis present in the formulation of the IP-DGFEM) is given. Makinguse of the recently introduced augmented Sobolev space framework,we prove an hp-optimal error bound in the energy norm and anh-optimal and slightly p-suboptimal (by only half an order ofp) bound in the L² norm (the latter, for the symmetric versionof the IP-DGFEM), provided that the solution belongs to an augmentedSobolev space.     

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.     

A Fourth-order Tridiagonal Finite Difference Method for General Non-linear Two-point Boundary Value Problems with Mixed Boundary Conditions

We present a new fourth-order finite difference method for thegeneral second-order non-linear differential equation y^N = f(x,y, y') subject to mixed two-point boundary conditions. An interestingfeature of our method is that each discretization of the differentialequation at an interior grid point is based on just three evaluationsof f. We establish, under appropriate conditions, O(h⁴)-convergenceof the finite difference scheme. In the case of linear differentialequations, our finite difference scheme leads to tridiagonallinear systems. Numerical examples are considered to demonstratecomputationally the fourth order of the method.     

The Norm of Powers of the Indefinite Integral Operator on (0, 1)

In this paper we find the norm of powers of the indefinite integraloperator V, acting on L²(0, 1). This answers a question raisedby Halmos, and supplements some recent results of Manakov in[9]. Using results of Stepanov in [13], we show that the operatornorm of Vⁿ is asymptotically equal to the Hilbert–Schmidtnorm as n . 1991 Mathematics Subject Classification 26D15,47B38.     

A Petrov-Galerkin method with quadrature for elliptic boundary value problems

We propose and analyse a fully discrete Petrov–Galerkinmethod with quadrature, for solving second-order, variable coefficient,elliptic boundary value problems on rectangular domains. Inour scheme, the trial space consists of C² splines of degreer 3, the test space consists of C⁰ splines of degree r –2, and we use composite (r – 1)-point Gauss quadrature.We show existence and uniqueness of the approximate solutionand establish optimal order error bounds in H², H¹ and L² norms.     

Rates of Recurrence for Zq and Rq Extensions of Subshifts of Finite Type

In this note we give new asymptotic formulae for certain countingfunctions associated to the periodic behaviour of Z^q and R^qextensions of subshifts of finite type. In the case of the Z^qextensions, these strengthen previous estimates of Marcus andTuncel [9]. For both types of extension, our results complementthe central limit type results of Lalley [6]. Our proof requiresthe application of ideas from thermodynamic formalism. Whilstdeveloping this approach, in Section 2, we take the opportunityto present a counter-example to a related conjecture of Coelho-Filho[2].     

On Fusion in Unipotent Blocks

The context of this note is as follows. One considers a connectedreductive group G and a Frobenius endomorphism F: G G definingG over a finite field of order q. One denotes by G^F the associated(finite) group of fixed points. Let l be a prime not dividing q. We are interested in the l-blocksof the finite group G^F. Such a block is called unipotent ifthere is a unipotent character (see, for instance, [6, Definition12.1]) among its representations in characteristic zero. Roughlyspeaking, it is believed that the study of arbitrary blocksof G^F might be reduced to unipotent blocks (see [2, Théorème2.3], [5, Remark 3.6]). In view of certain conjectures aboutblocks (see, for instance, [9]), it would be interesting tofurther reduce the study of unipotent blocks to the study ofprincipal blocks (blocks containing the trivial character).Our Theorem 7 is a step in that direction: we show that thelocal structure of any unipotent block of G^F is very close tothat of a principal block of a group of related type (notionof ‘control of fusion’, see [13, 49]). 1991 MathematicsSubject Classification 20Cxx.     

Non-variational approximation of discrete eigenvalues of self-adjoint operators

^** Email: L.Boulton{at}ma.hw.ac.uk We establish sufficient conditions for approximation of discreteeigenvalues of self-adjoint operators in the second-order projectionmethod suggested recently in Levitin & Shargorodsky (2004,Spectral pollution and second order relative spectra for self-adjointoperators. IMA J. Numer. Anal., 24, 393–416). We findfairly explicit estimates for the eigenvalue error and studyin detail two concrete model examples. Our results show thatsecond-order projection strategies not only are universallypollution free but also achieve approximation under naturalconditions on the discretising basis.     

Simplices of Maximal Volume or Minimal Total Edge Length in Hyperbolic Space

In this paper, we are mainly concerned with n-dimensional simplicesin hyperbolic space Hⁿ. We will also consider simplices withideal vertices, and we suggest that the reader keeps the Poincaréunit ball model of hyperbolic space in mind, in which the sphereat infinity Hⁿ() corresponds to the bounding sphere of radius1. It is known that all hyperbolic simplices (even the idealones) have finite volume. However, explicit calculation of theirvolume is generally a very difficult problem (see, for example,[1] or [16]). Our first theorem states that, amongst all simplicesin a closed geodesic ball, the simplex of maximal volume isregular. We call a simplex regular if every permutation of itsvertices can be realized by an isometry of Hⁿ. A correspondingresult for simplices in the sphere has been proved by Böröczky[4].     

A Note on a Maximal Function Over Arbitrary Sets of Directions

Let M_S be the universal maximal operator over unit vectors ofarbitrary directions. This operator is not bounded in L²(R²).We consider a sequence of operators over sets of finite equidistributeddirections converging to M_S. We provide a new proof of N. Katz'sbound for such operators. As a corollary, we deduce that M_Sis bounded from some subsets of L² to L². These subsets arecomposed of positive functions whose Fourier transforms havea logarithmic decay or which are supported on a disc. 1991 MathematicsSubject Classification 42B25.     

On the computation of the generalized inverse of a polynomial matrix

The main purpose of this paper is to present a quicker and less memory-expensive algorithm for the generalized inversionof polynomial matrices than those presented earlier (Karampetakis,1997a Computation of the generalized inverse of a polynomialmatrix and applications. Linear Algebr. Appl. 252, 35–60 and Karampetakis, 1997b Generalized inverses of two variable polynomial matrices and applications. Circuit Syst. & Signal Process. 16, 439–453). Received 24 January, 1999. ⁺ karampetakis@ccf.auth.gr     

RIGIDITY OF PERIODIC DIFFEOMORPHISMS OF HOMOTOPY K3 SURFACES

In this paper, we show that homotopy K3 surfaces do not admita periodic diffeomorphism of odd prime order 3 acting triviallyon cohomology. This gives a negative answer for period 3 toProblem 4.124 in Kirby's problem list. In addition, we givean obstruction in terms of the rationality and the sign of thespin numbers to the non-existence of a periodic diffeomorphismof odd prime order acting trivially on cohomology of homotopyK3 surfaces. The main strategy is to calculate the Seiberg–Witteninvariant for the trivial spin^c structure in the presence ofsuch a Z_p-symmetry in two ways: (1) the new interpretation ofthe Seiberg–Witten invariants of Furuta and Fang, and(2) the theorem of Morgan and Szabó on the Seiberg–Witteninvariant of homotopy K3 surfaces for the trivial Spin^c structure.As a consequence, we derive a contradiction for any periodicdiffeomorphism of prime order 3 acting trivially on cohomologyof homotopy K3 surfaces.     

Intertwining Maps from Certain Group Algebras

In [17, 18, 19], we began to investigate the continuity propertiesof homomorphisms from (non-abelian) group algebras. Alreadyin [19], we worked with general intertwining maps [3, 12]. Thesemaps not only provide a unified approach to both homomorphismsand derivations, but also have some significance in their ownright in connection with the cohomology comparison problem [4]. The present paper is a continuation of [17, 18, 19]; this timewe focus on groups which are connected or factorizable in thesense of [26]. In [26], G. A. Willis showed that if G is a connectedor factorizable, locally compact group, then every derivationfrom L¹(G) into a Banach L¹(G)-module is automatically continuous.For general intertwining maps from L¹(G), this conclusion isfalse: if G is connected and, for some nN, has an infinite numberof inequivalent, n-dimensional, irreducible unitary representations,then there is a discontinuous homomorphism from L¹(G into aBanach algebra by [18, Theorem 2.2] (provided that the continuumhypothesis is assumed). Hence, for an arbitrary intertwiningmap from L¹(G), the best we can reasonably hope for is a resultasserting the continuity of on a ‘large’, preferablydense subspace of L¹(G). Even if the target space of is a Banachmodule (which implies that the continuity ideal I() of is closed),it is not a priori evident that is automatically continuous:the proofs of the automatic continuity theorems in [26] relyon the fact that we can always confine ourselves to restrictionsto L¹(G) of derivations from M(G) [25, Lemmas 3.1 and 3.4].It is not clear if this strategy still works for an arbitraryintertwining map from L¹(G) into a Banach L¹(G)-module.     

On the division of polynomial matrices

The main purpose of this paper is to determine two new algorithmsfor the division of the polynomial matrix B(s) R[s]^pxq by A(s) R[s]^pxp (a) based on the Laurent matrix expansion at s = =of the inverse of A(s), i.e. A(s)^–1, and (b) in a waysimilar to the one presented by Gantmacher (1959).     

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.     

Semidiscrete finite element Galerkin approximations to the equations of motion arising in the Oldroyd model

In this paper, a semidiscrete finite element Galerkin methodfor the equations of motion arising in the 2D Oldroyd modelof viscoelastic fluids with zero forcing function is analysed.Some new a priori bounds for the exact solutions are derivedunder realistically assumed conditions on the data. Moreover,the long-time behaviour of the solution is established. By introducinga Stokes–Volterra projection, optimal error bounds forthe velocity in the L(L²) as well as in the L(H¹)-norms andfor the pressure in the L(L²)-norm are derived which are validuniformly in time t > 0.     

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.     

Beardon's Diophantine Equations and Non-Free Mobius Groups

Let µ be a real number. The Möbius group G_µis the matrix group generated by It is known that G_µ is free if |µ| 2 (see [1])or if µ is transcendental (see [3, 8]). Moreover, thereis a set of irrational algebraic numbers µ which is densein (–2, 2) and for which G_µ is non-free [2, p. 528].We may assume that µ > 0, and in this paper we considerrational µ in (0, 2). The following problem is difficult. Let G^nf denote the set of all rational numbers µ in (0,2) for which G_µ is non-free. In 1969 Lyndon and Ullman[8] proved that G^nf contains the elements of the forms p/(p²+ 1) and 1/(p + 1), where p = 1, 2, ..., and that if µ₀ G^nf, then µ₀/p G^nf for p = 1, 2, .... In 1993 Beardon[2] studied problem (P) by means of the words of the form A^rB^s A^t and A^r B^s A^t B^u A^v, and he obtained a sufficient conditionfor solvability of (P), included implicitly in [2, pp. 530–531],by means of the following Diophantine equations: 1991 Mathematics SubjectClassification 20E05, 20H20, 11D09.     

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.     

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.