Rate of convergence of numerical approximations to homoclinic bifurcation points

For x=f (x, ), x Rⁿ, R, having a hyperbolic or semihyperbolicequilibrium p(), we study the numerical approximation of parametervalues ^* at which there is an orbit homoclinic to p(). We approximate^* by solving a finite-interval boundary value problem on J=[T_–,T₊], T_–<0<T₊, with boundary conditions that sayx(T_–) and x(T₊) are in approximations to appropriate invariantmanifolds of p(). A phase condition is also necessary to makethe solution unique. Using a lemma of Xiao-Biao Lin, we improve,for certain phase conditions, existing estimates on the rateof convergence of the computed homoclinic bifurcation parametervalue , to the true value ^*. The estimates we obtain agree withthe rates of convergence observed in numerical experiments.Unfortunately, the phase condition most commonly used in numericalwork is not covered by our results.     

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.     

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     

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

An elliptic boundary-value problem on a domain with prescribedDirichlet data on _I is approximated using a finite-elementspace of approximation power h^K in the L² norm. It is shownthat the total flux across _I can be approximated with an errorof O(h^K) when is a curved domain in Rⁿ (n = 2 or 3) and isoparametricelements are used. When is a polyhedron, an O(h^2K–2)approximation is given. We use these results to study the finite-elementapproximation of elliptic equations when the prescribed boundarydata on _I is the total flux. Present address: School of Mathematical and Physical Sciences,University of Sussex, Brighton, Sussex BN1 9QH.     

Finite-Element Approximation of a Plasma Equilibrium Problem

The plasma problem studied is: given R⁺ find (, d, u) R ?R ? H¹() such that Let ₁ < ₂ be the first two eigenvalues of the associatedlinear eigenvalue problem: find $$\left(\lambda ,\phi \right)\in\mathrm{R;}\times {\hbox{ H }}_{0}^{1}\left(\Omega \right)$$such that For ₀(0,₂) it is well known that there exists a unique solution(₀, d₀, u₀) to the above problem. We show that the standard continuous piecewise linear Galerkinfinite-element approximatinon $$\left({\lambda }_{0},{\hbox{d }}_{0}^{k},{u}_{0}^{h}\right)$$, for ₀(0,₂), converges atthe optimal rate in the H¹, L², and L norms as h, the mesh length,tends to 0. In addition, we show that dist (, ^h)Ch² ln 1/h,where $${\Gamma }^{\left(h\right)}=\left\{x\in \Omega :{u}_{0}^{\left(h\right)}\left(x\right)=0\right\}$$.Finally we consider a more practical approximation involvingnumerical integration.     

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.     

The Convergence of Operator Approximations at Turning Points

Present address: Department of Mathematics, University of Reading, Reading RG6 2AX. We consider the convergence of solution curves of approximationsto parameter-dependent operator equations of the form G(, x)= 0. Provided G_x(, x) remains non-singular this problem is cateredfor by a simple extension to standard theory. In this paper,however, attention is concentrated on solution curves throughcertain singular points (⁰, x⁰), and the main result is thatconvergence depends on consistency and stability results forthe linear eigenvalue problem G_x(⁰, x⁰)0 = 0.     

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.     

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 .     

Frame duality properties for projective unitary representations

Let be a projective unitary representation of a countable groupG on a separable Hilbert space H. If the set B of Bessel vectorsfor is dense in H, then for any vector x H the analysis operator_x makes sense as a densely defined operator from B to ²(G)-space.Two vectors x and y are called -orthogonal if the range spacesof _x and _y are orthogonal, and they are -weakly equivalent ifthe closures of the ranges of _x and _y are the same. These propertiesare characterized in terms of the commutant of the representation.It is proved that a natural geometric invariant (the orthogonalityindex) of the representation agrees with the cyclic multiplicityof the commutant of (G). These results are then applied to Gaborsystems. A sample result is an alternate proof of the knowntheorem that a Gabor sequence is complete in L²( ^d) ifand only if the corresponding adjoint Gabor sequence is ²-linearlyindependent. Some other applications are also discussed.     

On Mason's conjecture concerning interpolation by polynomials in z and z-1 on an annulus

A polynomial of degree n in z^–1 and n – 1 in z isdefined by an interpolation projection from the space A(N) of functions f analytic in thecircular annulus ^–1<|z| and continuous on its boundaries|z|=^–1, . The points of interpolation are chosen to coincidewith the n roots of zⁿ=–n and the n roots of zⁿ=^–n.We prove Mason's conjecture that the corresponding Lebesguefunction attains its maximal value on the inner circle. We alsoestimate the bound of the Lebesgue constant . It is proved that the following estimate for theoperator norm holds: where _n, is the Lebesgue constant of Gronwall for equally spacedinterpolation on a circle by a polynomial of degree n.     

Recurrence and asymptotics for orthonormal rational functions on an interval

^{Let µ be a positive bounded Borel measure on a subsetI of the real line and = {₁, ..., _n} a sequence of arbitrary ‘complex’poles outside I. Suppose {₁, ..., _n} is the sequence of rationalfunctions with poles in orthonormal on I with respect to µ. First, we are concernedwith reducing the number of different coefficients in the three-termrecurrence relation satisfied by these orthonormal rationalfunctions. Next, we consider the case in which I = [–1, 1] and µ satisfies the Erdos–Turán conditionµ' > 0 a.e. on I (where µ' is the Radon–Nikodymderivative of the measure µ with respect to the Lebesguemeasure) to discuss the convergence of _n+1(x)/_n(x) as n tendsto infinity and to derive asymptotic formulas for the recurrencecoefficients in the three-term recurrence relation. Finally,we give a strong convergence result for _n(x) under the morerestrictive condition that µ satisfies the Szeg condition(1 – x2)–1/2 log µ'(x) L1([– 1, 1]).}     

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||=+.     

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}.     

Central L-Values and Toric Periods for GL(2)

Let be a cusp form on GL(2) over a number field F and let Ebe a quadratic extension of F. Denote by _E the base change of to E and by a unitary character of A^x_E/ E^x. We use the relativetrace formula to give an explicit formula for L(1/2, _E ) interms of period integrals of Gross–Prasad test vectors.We give an application of this formula to equidistribution ofgeodesics on a hyperbolic 3-fold.     

A Krylov subspace algorithm for multiquadric interpolation in many dimensions

We consider the version of multiquadric interpolation wherethe interpolation conditions are the equations s(x_i) = f_i, i= 1,2,..., n, and where the interpolant has the form s(x) =_j=1ⁿ _j (||x – x_j ||² + c²)^1/2 + x R^d, subject to theconstraint _j=1ⁿ _j = 0. The points x_i R^d, the right-hand sidesf_i, i = 1,2,...,n, and the constant c are data. The equationsand the constraint define the parameters _j, j = 1,2,...,n, and. The resultant approximation s f is useful in many applications,but the calculation of the parameters by direct methods requiresO (n³) operations, and n may be large. Therefore iterative proceduresfor this calculation have been studied at Cambridge since 1993,the main task of each iteration being the computation of s(x_i),i = 1,2,...,n, for trial values of the required parameters.These procedures are based on approximations to Lagrange functions,and often they perform very well. For example, ten iterationsusually provide enough accuracy in the case d = 2 and c = 0,for general positions of the data points, but the efficiencydeteriorates if d and c are increased. Convergence can be guaranteedby the inclusion of a Krylov subspace technique that employsthe native semi-norm of multiquadric functions. An algorithmof this kind is specified, its convergence is proved, and carefulattention is given to the choice of the operator that definesthe Krylov subspace, which is analogous to pre-conditioningin the conjugate gradient method. Finally, some numerical resultsare presented and discussed, for values of d and n from theintervals [2,40] and [200,10 000], respectively.     

Near-best Lp Approximations by Real and Complex Chebyshev Series

The expansion of a real or complex function in a series of Chebyshevpolynomials of the first and second kinds is discussed in thecontext of near-best approximation. The discussion covers realand complex approximation on the real interval [–1, 1]as a special example of the complex elliptical contour , as well as complex approximationon an elliptical domain, an ellipse exterior, and an ellipticalannulus (including special cases in which part of the boundarycollapses into a "crack"). Two distinct types of function spacesare considered, namely appropriately weighted L_p measure spacesand analytic function spaces, and resulting approximations areshown in all cases to be near-best in the L_p norm within a relativedistance asymptotic to (4^-2 log n)^2p-1–1 for all p (1p ), where relates to the order of approximation.     

The Isolated Points of G{rho} and the w*-Strongly Exposed Points of P{rho}(G)0

Let G be a separable locally compact group and let be its dualspace with Fell's topology. It is well known that the set P(G)of continuous positive-definite functions on G can be identifiedwith the set of positive linear functionals on the group C^*-algebraC^*(G). We show that if is discrete in , then there exists anonzero positive-definite function associated with such that is a w^*-strongly exposed point of P(G)₀, where P(G)₀={f P(G):f(e)1. Conversely, if some nonzero positive-definite function associatedwith is a w^*-strongly exposed point of P(G)₀, then is isolatedin . Consequently, G is compact if and only if, for every ,there exists a nonzero positive-definite function associatedwith that is a w^*-strongly exposed point of P(G)₀. If, in addition,G is unimodular and , then is isolated in if and only if somenonzero positive-definite function associated with is a w^*-stronglyexposed point of P(G)₀, where is the left regular representationof G and is the reduced dual space of G. We prove that if B(G)has the Radon–Nikodym property, then the set of isolatedpoints of (so square-integrable if G is unimodular) is densein . It is also proved that if G is a separable SIN-group, thenG is amenable if and only if there exists a closed point in. In particular, for a countable discrete non-amenable groupG (for example the free group F₂ on two generators), there isno closed point in its reduced dual space .     

Some Remarks on the Cone of Completely Positive Maps between Von Neumann Algebras

The close relationship between the notions of positive formsand representations for a C*-algebra A is one of the most basicfacts in the subject. In particular the weak containment ofrepresentations is well understood in terms of positive forms:given a representation of A in a Hilbert space H and a positiveform on A, its associated representation is weakly containedin (that is, ker ker ) if and only if belongs to the weak*closure of the cone of all finite sums of coefficients of .Among the results on the subject, let us recall the followingones. Suppose that A is concretely represented in H. Then everypositive form on A is the weak* limit of forms of the typex ^k_{i=1 i}, x_i with the _i in H; moreover if A is a von Neumannsubalgebra of (H) and is normal, there exists a sequence (_i)_{i 1} in H such that (x) = _{i 1 i}, x_i for all x.     

On hearing the shape of a bounded domain with Robin boundary conditions

The asymptotic expansions of the trace of the heat kernel (t)= [sum ]_j=1 exp(-t_j) for small positive t, where {_j} _j=1 arethe eigenvalues of the negative Laplacian -_n = -[sum ]ⁿ_k=1 (/x^k)²in Rⁿ (n = 2 or 3), are studied for a general multiply connectedbounded domain which is surrounded by simply connected boundeddomains _i with smooth boundaries _i (i = 1,...,m), where smoothfunctions Y_i (i = 1,...,m) are assuming the Robin boundary conditions(n_i + Y_i) = 0 on _i. Here /n_i denote differentiations along theinward-pointing normals to _i (i = 1,...,m). Some applicationsof an ideal gas enclosed in the multiply connected bounded containerwith Neumann or Robin boundary conditions are given.