Convergence of Univariate Quasi-Interpolation Using Multiquadrics

Quasi-interpolants to a function f: RR on an infinite regularmesh of spacing h can be defined by where :RR is a function with fast decay for large argument. In the approach employing the radial-basis-function : RR, thefunction is a finite linear combination of basis functions(|•–jh|) (jZ). Choosing Hardy's multiquadrics (r)=(r²+c²)^?,we show that sufficiently fast-decaying exist that render quasi-interpolationexact for linear polynomials f. Then, approximating f C²(R),we obtain uniform convergence of s to f as (h, c)0, and convergenceof s' to f' as (h, c²/h)0. However, when c stays bounded awayfrom 0 as h0, there are f C(R) for which s does not convergeto f as h0. We also show that, for all which vanish at infinity but arenot integrable over R, there are no finite linear combinations of the given basis functions allowing the construction of admissiblequasi-interpolants. This includes the case of the inverse multiquadncs(r)=(r²+c²)^–?.     

An Order of Convergence for Some Radial Basis Functions

The quasi-interpolant to a function f : RⁿR on an infinite regulargrid of spacing h can be defined by where : RⁿR is a function which decays quickly for large argument.In the case of radial basis functions has the form where : R⁺R is known as a radial basis function and, in general,?_j R (j = 1,...,m) and x^j Rⁿ (j = 1,...,m), though here onlythe particular case x^j Zⁿ (j = 1,..., m) is considered. Thispaper concentrates on the case (r) = r, a generalization oflinear interpolation, although some of the analysis is moregeneral. It is proved that, if n is odd, then there is a function such that the maximum difference between a sufficiently smoothfunction and its quasi-interpolant is bounded by a constantmultiple of hⁿ⁺¹. This is done by first showing that such aquasi-interpolation formula can reproduce polynomials of degreen.     

Boundary integral methods for singularly perturbed boundary value problems

In this paper we consider boundary integral methods appliedto boundary value problems for the positive definite Helmholtz-typeproblem –U + ²U = 0 in a bounded or unbounded domain,with the parameter real and possibly large. Applications arisein the implementation of space–time boundary integralmethods for the heat equation, where is proportional to 1/(t),and t is the time step. The corresponding layer potentials arisingfrom this problem depend nonlinearly on the parameter and havekernels which become highly peaked as , causing standard discretizationschemes to fail. We propose a new collocation method with arobust convergence rate as . Numerical experiments on a modelproblem verify the theoretical results.     

Convergent and spurious solutions of nonlinear elliptic equations

In this paper we investigate finite element approximations ofnonlinear elliptic equations in three dimensions. By applyingand extending the results of Lopez-Marcos and Sanz-Serna, weprove that the finite element approximation on a mesh of sizeh, has a solution U_k which converges to an exact solution ofthe differential equation as h0. This solution is unique withina suitably defined stability ball Bh. For the particular nonlinearequation u + (u + u^p) we show that the size of B_h depends uponh only if p > 5 when it tends to zero as h 0. In this casewe prove the existence of spurious solutions V_h of the Galerkinapproximation which become unbounded in the maximum norm ash0. The stability ball B_h then acts to separate the convergentand the spurious solutions. We present the results of some numericalexperiments to substantiate our claims.     

Fast evaluation of the gamma function for small rational fractions using complete elliptic integrals of the first kind

Much simplified expressions for certain complete elliptic integralsin terms of the beta function are produced. The reverse procedureof expressing gamma functions in terms of complete ellipticintegrals allows us to use the arithmetic-geometric mean iterationto compute gamma functions, using quadratically convergent iterations.Explicit algorithms are given for computing (n/4) and (n/6),where 1m3 and 1n5.     

A Residually Finite Version of Rips's Construction

For each finitely presented group Q a short exact sequence 1 N G Q 1 is produced, such that G is residually finite andN is finitely generated.     

Spurious solutions of numerical methods for initial value problems

It is well known that some numerical methods for initial valueproblems admit spurious limit sets. Here the existence and behaviourof spurious solutions of Runge-Kutta, linear multistep and predictor-correctormethods are studied in the limit as the step-size h0. In particular,it is shown that for ordinary differential equations definedby globally Lipschitz vector fields, spurious fixed points andperiod 2 solutions cannot exist for h arbitrarily small, whilstfor locally Lipschitz vector fields, spurious solutions mayexist for h arbitrarily small, but must become unbounded ash0. The existence of spurious solutions is also studied forvector fields merely assumed to be continuous, and an exampleis given, showing that in this case spurious solutions may remainbounded as h0. It is shown that if spurious fixed points orperiod 2 solutions of continuous problems exist for h arbitrarilysmall, then as h0 spurious solutions either converge to steadysolutions of the underlying differential equation or divergeto infinity. A necessary condition for the bifurcation spurioussolutions from h=0 is derived. To prove the above results forimplicit Runge-Kutta methods, an additional assumption on theiteration scheme used to solve the nonlinear equations definingthe method is needed; an example of a Runge-Kutta method whichgenerates a bounded spurious solution for a smooth problem withh arbitrarily small is given, showing that such an assumptionis necessary. It is also shown that an Adams-Bashforth/Adams-Moultonpredictor-corrector method in PC^m implementation can generatespurious fixed point solutions for any m.     

Metrical Theorems on Prime Values of the Integer Parts of Real Sequences II

Let [ ] denote the integer part. Among other results in [3]we gave a complete solution to the following problem. PROBLEM. Given an increasing sequence a_n R⁺, n = 1, 2, ...,where a_n as n , are there infinitely many primes in the sequence[a_n] for almost all ?     

Total Flux Estimates for a Finite-Element Approximation of Parabolic Equations

An initial-boundary-value problem for a parabolic equation ina domain x (0, T) with prescribed Dirichlet data on is approximatedusing a continuous-time Galerkin finite-element scheme. It isshown that the total flux across ₁= can be approximated withan error of O(h^k) when is a curved domain in Rⁿ (n = 2 or 3)and isoparametric elements having approximation power h^k inthe L² norm are used.     

Extremal Properties of Contraction Semigroups on Hilbert and Banach Spaces

Let f be a unit vector and T = {T(t) = e^tA: t 0} be a (C₀)contraction semigroup generated by A on a complex Hilbert spaceX. If |T(t)f,f| 1 as t then f is an eigenvector of A correspondingto a purely imaginary eigenvalue. If one allows X to be a Banachspace, the same situation can be considered by replacing T(t)f,fby (T(t)f) where is a unit vector in X^* dual to f. If |(T(t)f)| 1, as t , is f an eigenvector of A? The answer is sometimesyes and sometimes no.     

Finite-Element Approximation of Elliptic Equations with a Neumann or Robin Condition on a Curved Boundary

This paper considers a finite-element approximation of a second-orderself adjoint elliptic equation in a region Rⁿ (with n=2 or 3)having a curved boundary on which a Neumann or Robin conditionis prescribed. If the finite-element space defined over , a union of elements, has approximation power h^kin the L² norm, and if the region of integration is approximatedby ^h with dist (, ^h)Ch^k, then it is shown that one retains optimalrates of convergence for the error in the H¹ and L² norms, whetherQ^h is fitted or unfitted , provided that the numerical integration scheme has sufficientaccuracy.     

A Classification of Horizontally Homothetic Submersions from Space Forms of Nonnegative Curvature

In this paper, it is proved that for n 2, any horizontallyhomothetic submersion : Rⁿ⁺¹ (Nⁿ, h) is a Riemannian submersionup to a homothety. It is also shown that if : Sⁿ⁺¹ (Nⁿ, h)is a horizontally homothetic submersion, then n = 2m, (Nⁿ, h)is isometric to CP^m and, up to a homothety, is a standard Hopffibration S^2m+1 CP^m. 2000 Mathematics Subject Classification53C20, 53C12.     

Finite Element Approximation of a Rigid Punch Indenting a Membrane

Optimal order H¹ and L error bounds are obtained for a continuouspiecewise linear finite element approximation of an obstacleproblem, where the obstacle's height as well as the contactzone, _c, are a priori unknown. The problem models the indentationof a membrane by a rigid punch. For R², given ,g R⁺ and an obstacle defined over E we consider the minimization of |v|²_1,+gµover (v, µ) H¹₀() x R subject to v+µ on E. In additionwe show under certain nondegeneracy conditions that dist (_c,^h_c)Ch ln 1/h, where ^h_c is the finite element approximation to_c. Finally we show that the resulting algebraic problem canbe solved using a projected SOR algorithm.     

Contractions et Hyperdistributions a Spectre de Carleson

Let =(_n)_n1 be a log concave sequence such that lim inf_n+n/n^c>0for some c>0 and ((log _n)/n)_n1 is nonincreasing for some<1/2. We show that, if T is a contraction on the Hilbertspace with spectrum a Carleson set, and if ||T^–n||=O(_n)as n tends to + with _n11/(n log _n)=+, then T is unitary. Onthe other hand, if _n11/(n log _n)<+, then there exists a (non-unitary)contraction T on the Hilbert space such that the spectrum ofT is a Carleson set, ||T^–n||=O(_n) as n tends to +, andlim sup_n+||T^–n||=+.     

Excision in Cyclic Type Homology of Frechet Algebras

It is proved that every topologically pure extension of Fréchetalgebras 0 I A A/I 0 such that I is strongly H-unital hasthe excision property in continuous (co)homology of the followingtypes: bar, naive-Hochschild, Hochschild, cyclic, and periodiccyclic. In particular, the property holds for every extensionof Fréchet algebras such that I has a left or right boundedapproximate identity.     

Discontinuous Galerkin finite element approximation of quasilinear elliptic boundary value problems I: the scalar case

^** Email: Paul.Houston{at}mcs.le.ac.uk^*** Email: Janice.Robson{at}comlab.ox.ac.uk^**** Email: Endre.Suli{at}comlab.ox.ac.uk We develop a one-parameter family of hp-version discontinuousGalerkin finite element methods, parameterised by [–1,1], for the numerical solution of quasilinear elliptic equationsin divergence form on a bounded open set ^d, d 2. In particular,we consider the analysis of the family for the equation –·{µ(x, |u|)u} = f(x) subject to mixed Dirichlet–Neumannboundary conditions on . It is assumed that µ is a real-valuedfunction, µ C( x [0, )), and thereexist positive constants m_µ and M_µ such that m_µ(t– s) µ(x, t)t – µ(x, s)s M_µ(t– s) for t s 0 and all x . Using a result from the theory of monotone operators for any valueof [–1, 1], the corresponding method is shown to havea unique solution u_DG in the finite element space. If u C¹() H^k(), k 2, then with discontinuous piecewise polynomials ofdegree p 1, the error between u and u_DG, measured in the brokenH¹()-norm, is (h^s–1/p^k–3/2), where 1 s min {p+ 1, k}.     

Invariant Subspaces of C1.-Contractions with Non-Reductive Unitary Extensions

Let T be a contraction acting on the Hilbert space H such thatlim_n||Tⁿh||0, for every nonzero h;H. It is proved that if theunitary operator attached to T in a canonic way contains thebilateral shift, then T has a non-trivial invariant subspace.Furthermore, if in addition lim_n||T^*nh||0 holds for every nonzeroh H, then T is shown to be reflexive.     

Fitted and Unfitted Finite-Element Methods for Elliptic Equations with Smooth Interfaces

This paper considers the finite-element approximation of theelliptic interface problem: -?(u) + cu = f in Rⁿ (n = 2 or3), with u = 0 on , where is discontinuous across a smoothsurface in the interior of . First we show that, if the meshis isoparametrically fitted to using simplicial elements ofdegree k - 1, with k 2, then the standard Galerkin method achievesthe optimal rate of convergence in the H¹ and L² norms overthe approximations _l⁴ of _l where _{l 2}. Second, since itmay be computationally inconvenient to fit the mesh to , weanalyse a fully practical piecewise linear approximation ofa related penalized problem, as introduced by Babuska (1970),based on a mesh that is independent of . We show that, by choosingthe penalty parameter appropriately, this approximation convergesto u at the optimal rate in the H¹ norm over _l⁴ and in the L²norm over any interior domain _l^* satisfying _l^* _l^** _l⁴ for somedomain _l^**. Present address: School of Mathematical and Physical Sciences,University of Sussex, Brighton BN1 9QH     

Error Bounds for Hermite Interpolation by Quadratic Splines on an {alpha}-Triangulation

For l, an -triangulation F of a planar domain is such that,for every T F, there holds 1 R_T/2r_T , where R_T (resp. r_T)denotes the radius of the circumscribed (resp. inscribed) circleof the triangle T. When T is varying in F the centre of itsinscribed circle is varying in a compact interior to T and itsorthogonal projections on the sides are varying in compact intervalsinterior to these sides. Precise results are given about thesizes of these compacts and are used for the computation oferror constants in the problem of Hermite interpolation by Powell-Sabinquadratic finite elements, bringing to the fore their dependenceon the parameter .     

On orders Solely of Abelian Groups II

Let C'(x) denote the number of integers n x such that thereis no non-abelian group of order n, but there exists a non-cyclicgroup of order n. Here it is shown that C' (x) , x , where denotes Euler's constant.