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.     

High-accuracy P-stable Methods for y' = f(t, y)

We obtain a one-parameter family of sixth-order P-stable methodsfor the numerical integration of periodic or near-periodic differentialequations that are defined by initial-value problems of theform: y^" = f(t, y), y(t₀)= y₀, y^'(t₀)= y^’₀. Our P-stablemethods are symmetric and involve three function evaluationsper step (periteration, in case f(t, y) is non-linear in y).For non-linear problems, starting values for the solution ofthe implicit equations by modified Newton's method are suggestedand illustrated by an example.     

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.     

Symmetric Groups as Products of Abelian Subgroups

A proof is given that the full symmetric group over any infiniteset is the product of finitely many Abelian subgroups. In fact,289 subgroups suffice. Sharp bounds are also obtained on theminimal number k, such that the finite symmetric group S_n isthe product of k Abelian subgroups. Using this, S_n is provedto be the product of 72n^1/2(log n)^3/2 cyclic subgroups. 2000Mathematics Subject Classification 20B30, 20D40.     

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.     

The Solubility of Diagonal Cubic Diophantine Equations

This paper proves conditional existence results for non-trivialsolutions of the equation where the coefficients a_i and the unknowns X_i are taken to berational integers. No such results were previously known for n6. The proofs useelementary facts about the 3-descent procedure for ellipticcurves of the form E_A: X³ + Y³ = AZ³. Thus, when n=4, and the a_i are each prime, and are all congruentto 2 modulo 3, it is shown that (*) will have non-trivial solutions,providing that the Selmer conjecture holds for the curves E_A.One may replace the Selmer conjecture by an appropriate formof the Generalized Riemann Hypothesis, when n=5 and the a_i areagain taken to be primes, all congruent to 8 modulo 9. Finally,when n=5, one may require only that the a_i be square-free andcoprime to 3, providing one assumes both the Selmer conjectureand a special case of Schinzel's conjecture (on the representationof primes by cubic polynomials). 1991 Mathematics Subject Classification:11D25, 11G05, 14G05.     

Stationary Critical Points of the Heat Flow in Spaces of Constant Curvature

The paper considers stationary critical points of the heat flowin sphere S^N and in hyperbolic space H^N, and proves severalresults corresponding to those in Euclidean space R^N which havebeen proved by Magnanini and Sakaguchi. To be precise, it isshown that a solution u of the heat equation has a stationarycritical point, if and only if u satisfies some balance lawwith respect to the point for any time. In Cauchy problems forthe heat equation, it is shown that the solution u has a stationarycritical point if and only if the initial data satisfies thebalance law with respect to the point. Furthermore, one point,say x₀, is fixed and initial-boundary value problems are consideredfor the heat equation on bounded domains containing x₀. It isshown that for any initial data satisfying the balance law withrespect to x₀ (or being centrosymmetric with respect to x₀)the corresponding solution always has x₀ as a stationary criticalpoint, if and only if the domain is a geodesic ball centredat x₀ (or is centrosymmetric with respect to x₀, respectively).     

On Artin's Conjecture, I: Systems of Diagonal Forms

As a special case of a well-known conjecture of Artin, it isexpected that a system of R additive forms of degree k, say [formula] with integer coefficients a_ij, has a non-trivial solution inQ_p for all primes p whenever [formula] Here we adopt the convention that a solution of (1) is non-trivialif not all the x_i are 0. To date, this has been verified onlywhen R=1, by Davenport and Lewis [4], and for odd k when R=2,by Davenport and Lewis [7]. For larger values of R, and in particularwhen k is even, more severe conditions on N are required toassure the existence of p-adic solutions of (1) for all primesp. In another important contribution, Davenport and Lewis [6]showed that the conditions [formula] are sufficient. There have been a number of refinements of theseresults. Schmidt [13] obtained N>>R²k³ log k, and Low,Pitman and Wolff [10] improved the work of Davenport and Lewisby showing the weaker constraints [formula] to be sufficient for p-adic solubility of (1). A noticeable feature of these results is that for even k, onealways encounters a factor k³ log k, in spite of the expectedk² in (2). In this paper we show that one can reach the expectedorder of magnitude k². 1991 Mathematics Subject Classification11D72, 11D79.     

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     

Canonical curves in P3

Let Hilb^6t–3(P³) be the Hilbert scheme of closed 1-dimensionalsubschemes of degree 6 and arithmetic genus 4 in P³. Let H bethe component of Hilb^6t–3(P³) whose generic point correspondsto a canonical curve, that is, a complete intersection of aquadric and a cubic surface in P³. Let F be the vector spaceof linear forms in the variables z₁, z₂, z₃, z₄. Denote by F_dthe vector space of homogeneous forms of degree d. Set X = (f₂,f₃)where f₂ P(F₂) is a quadric surface, and f₃ P(F₃/f₂ ·F) is a cubic modulo f₂. Wehave a rational map, : X ... Hdefined by (f₂,f₃) f₂ f₃. It fails to be regular along thelocus where f₂ and f₃ acquire a common linear component. Ourmain result gives an explicit resolution of the indeterminaciesof as well as of the singularities of H. 2000 Mathematical Subject Classification: 14C05, 14N05, 14N10,14N15.     

Arbitrary-norm hyperplane separation by variable neighbourhood search

^** Email: alejandro.karam{at}hec.ca^*** Email: gilles.caporossi{at}gerad.ca^**** Email: pierre.hansen{at}gerad.ca We consider the problem of separating two sets of points ina Euclidean space with a hyperplane that minimizes the sum ofp-norm distances to the plane of points lying on the ‘wrong’side of the plane. A variable neighbourhood search heuristicis used to determine the plane coefficients. For a set of exampleswith L₁-norm, L₂-norm and L-norm, for which the exact solutioncan be computed, we show that our algorithm finds it in mostcases and gets good approximations in the others. The use ofour heuristic solutions for problems in these norms can dramaticallyaccelerate exact algorithms. Our method can be applied on verylarge instances that are intractable by exact algorithms. Sincethe proposed approach works for truly arbitrary norms (otherthan the traditional 1, 2 and ), we can explore for the firsttime the effects of the choice of p on the generalization propertiesof p-norm hyperplane separation.     

A Mod Two Analogue of a Conjecture of Cooke

The mod two cohomology of the three connective covering of S³has the form F₂[X_2n] E(Sq¹X_2n) where x_2n is in degree 2n and n = 2. If F denotes the homotopytheoretic fibre of the map S³ B²S¹ of degree 2, then the mod2 cohomology of F is also of the same form for n = 1. Notice(cf. Section 7 of the present paper) that the existence of spaceswhose cohomology has this form for high values of n would immediatelyprovide Arf invariant elements in the stable stem. Hence, itis worthwhile to determine for what values of n the above algebracan be realized as the mod2 cohomology of some space. The purposeof this paper is to construct a further example of a space withsuch a cohomology algebra for n = 4 and to show that no othervalues of n are admissible. More precisely, we prove the following.     

On the Discreteness and Convergence in n-Dimensional Mobius Groups

Throughout this paper, we adopt the same notations as in [1,6, 8] such as the Möbius group M(Rⁿ), the Clifford algebraC_n–1, the Clifford matrix group SL(2, _n), the Cliffordnorm of ||A||=(|a|²+|b|²+|c|²+|d|²) (1) and the Clifford metric of SL(2, _n) or of the Möbius groupM(Rⁿ) d(A₁,A₂)=||A₁–A₂||(|a₁–a₂|²+|b₁–b₂|²+|c₁–c₂|²+|d₁–d₂|²)(2) where |·| is the norm of a Clifford number and represents f_i M(), i = 1,2, and so on. In addition, we adopt some notions in [6, 12]:the elementary group, the uniformly bounded torsion, and soon. For example, the definition of the uniformly bounded torsionis as follows.     

A note on least-squares mixed finite elements in relation to standard and mixed finite elements

^** Email: brandts{at}science.uva.nl The least-squares mixed finite-element method for second-orderelliptic problems yields an approximation u_h V_h H₀¹() of thepotential u together with an approximation p_{h h} H(div ; )of the vector field p = – Au. Comparing u_h with the standardfinite-element approximation of u in V_h, and p_h with the mixedfinite-element approximation of p, it turns out that they arehigher-order perturbations of each other. In other words, theyare ‘superclose’. Refined a priori bounds and superconvergenceresults can now be proved. Also, the local mass conservationerror is of higher order than could be concluded from the standarda priori analysis.     

Harmonic morphisms on heaven spaces

We prove that any (real or complex) analytic horizontally conformalsubmersion from a three-dimensional conformal manifold (M³,c_M) to a two-dimensional conformal manifold (N², c_N) can be,locally, ‘extended’ to a unique harmonic morphismfrom the (eaven)-space (H⁴, g) of (M³, c_N) to (N², c_N). Moreover,any positive harmonic morphism with two-dimensional fibres from(H⁴, g) is obtained in this way.     

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.     

Integral Representations of Derivatives and Integrals with Respect to the Order of the Bessel Functions Jv(t), Iv(t), the Anger Function Jv(t) and the Integral Bessel Function Jiv(t)

^*To whom Correspondence should be addressed. On sabbatical leave in the University of Alberta, Department of Chemical Engineering, Edmonton, Alberta, Canada, T6G 2G6. Integral representations of integrals and derivatives with respectto the order of the Bessel functions J_v(t) I_v(t), the integralBessel function Ji_v(t) and the Anger function J_v(t) are presented.The Laplace transform technique is applied to derive them. Theintegral representations permit the evaluation of a number oftrigonometric integrals.     

On the Sloan iteration applied to integral equations of the first kind

When the piecewise constant collocation method is used to solvean integral equation of the first kind with logarithmic kernel,the convergence rate is O(h) in the L₂ norm. In this note weshow that O(h³) or O(h⁵) convergence in any Sobolev norm (andthus, for example, in L) may be obtained by a simple cheap postprocessingof the original collocation solution. The construction of thepostprocessor is based on writing the first kind equation asa second kind equation, and applying the Sloan iteration tothe latter equation. The theoretical convergence rates are verifiedin a numerical example.     

The Cohomology of the Sylow 2-Subgroup of J2

The Hall-Janko-Wales group J₂ is one of the twenty-six sporadicfinite simple groups. The cohomology of its Sylow 2-subgroupS_J is computed, an important step in calculating the mod 2 cohomologyof J₂. The spectral sequence corresponding to the central extensionfor S_J is described and shown to collapse at the eighth page.The group S_J contains two subgroups (the central product of a dihedral and a quaternionic group)and 2²⁺⁴ (the Sylow 2-subgroup of the matrix group PSL₃(₄))which detect the cohomology of S_J. The cohomology relationsfor the subgroup 2²⁺⁴ are computed.     

Higher string topology on general spaces

In this paper, I give a generalized analogue of the string topologyresults of Chas and Sullivan, and of Cohen and Jones. For afinite simplicial complex X and k 1, I construct a spectrumMaps(S^k, X)^S(X), which is obtained by taking a generalizationof the Spivak bundle on X (which however is not a stable spherebundle unless X is a Poincaré space), pulling back toMaps(S^k, X) and quotienting out the section at infinity. I showthat the corresponding chain complex is naturally homotopy equivalentto an algebra over the (k + 1)-dimensional unframed little diskoperad C_{k + 1}. I also prove a conjecture of Kontsevich, whichstates that the Quillen cohomology of a based C_k-algebra (inthe category of chain complexes) is equivalent to a shift ofits Hochschild cohomology, as well as prove that the operadC_*C_k is Koszul-dual to itself up to a shift in the derived category.This gives one a natural notion of (derived) Koszul dual C_*C_k-algebras.I show that the cochain complex of X and the chain complex of^k X are Koszul dual to each other as C_*C_k-algebras, and thatthe chain complex of Maps(S^k, X)^S(X) is naturally equivalentto their (equivalent) Hochschild cohomology in the categoryof C_* C_k-algebras. 2000 Mathematics Subject Classification 55P48(primary), 16E40, 55N45, 18D50 (secondary).