Wavelet bases in

Some years ago, compactly supported divergence-free wavelets were constructed which also gave rise to a stable (biorthogonal) wavelet splitting of . These bases have successfully been used both in the analysis and numerical treatment of the Stokes and Navier-Stokes equations. In this paper, we construct stable wavelet bases for the stream function spaces . Moreover, -free vector wavelets are constructed and analysed. The relationship between and are expressed in terms of these wavelets. We obtain discrete (orthogonal) Hodge decompositions.

Our construction works independently of the space dimension, but in terms of general assumptions on the underlying wavelet systems in that are used as building blocks. We give concrete examples of such bases for tensor product and certain more general domains . As an application, we obtain wavelet multilevel preconditioners in and .

     

Analysis of multilevel methods for eddy current problems

In papers by Arnold, Falk, and Winther, and by Hiptmair, novel multigrid methods for discrete -elliptic boundary value problems have been proposed. Such problems frequently occur in computational electromagnetism, particularly in the context of eddy current simulation.

This paper focuses on the analysis of those nodal multilevel decompositions of the spaces of edge finite elements that form the foundation of the multigrid methods. It provides a significant extension of the existing theory to the case of locally vanishing coefficients and nonconvex domains. In particular, asymptotically uniform convergence of the multigrid method with respect to the number of refinement levels can be established under assumptions that are satisfied in realistic settings for eddy current problems.

The principal idea is to use approximate Helmholtz-decompositions of the function space into an -regular subspace and gradients. The main results of standard multilevel theory for -elliptic problems can then be applied to both subspaces. This yields preliminary decompositions still outside the edge element spaces. Judicious alterations can cure this.

     

Canonical construction of finite elements

The mixed variational formulation of many elliptic boundary value problems involves vector valued function spaces, like, in three dimensions, and . Thus finite element subspaces of these function spaces are indispensable for effective finite element discretization schemes. Given a simplicial triangulation of the computational domain , among others, Raviart, Thomas and Nédélec have found suitable conforming finite elements for and . At first glance, it is hard to detect a common guiding principle behind these approaches. We take a fresh look at the construction of the finite spaces, viewing them from the angle of differential forms. This is motivated by the well-known relationships between differential forms and differential operators: , and can all be regarded as special incarnations of the exterior derivative of a differential form. Moreover, in the realm of differential forms most concepts are basically dimension-independent. Thus, we arrive at a fairly canonical procedure to construct conforming finite element subspaces of function spaces related to differential forms. In any dimension we can give a simple characterization of the local polynomial spaces and degrees of freedom underlying the definition of the finite element spaces. With unprecedented ease we can recover the familiar - and -conforming finite elements, and establish the unisolvence of degrees of freedom. In addition, the use of differential forms makes it possible to establish crucial algebraic properties of the canonical interpolation operators and representation theorems in a single sweep for all kinds of spaces.

     

A priori error estimates for Galerkin approximations to porous medium and fast diffusion equations

Galerkin approximations to solutions of a Cauchy-Dirichlet problem governed by the generalized porous medium equation

on bounded convex domains are considered. The range of the parameter includes the fast diffusion case . Using an Euler finite difference approximation in time, the semi-discrete solution is shown to converge to the exact solution in norm with an error controlled by for and for . For the fully discrete problem, a global convergence rate of in norm is shown for the range . For , a rate of is shown in norm.

     

The monic integer transfinite diameter

We study the problem of finding nonconstant monic integer polynomials, normalized by their degree, with small supremum on an interval . The monic integer transfinite diameter is defined as the infimum of all such supremums. We show that if has length , then .

We make three general conjectures relating to the value of for intervals of length less than . We also conjecture a value for where . We give some partial results, as well as computational evidence, to support these conjectures.

We define functions and , which measure properties of the lengths of intervals with on either side of . Upper and lower bounds are given for these functions.

We also consider the problem of determining when is a Farey interval. We prove that a conjecture of Borwein, Pinner and Pritsker concerning this value is true for an infinite family of Farey intervals.

     

Computation of Stark-Tamagawa units

Let be a totally real number field and let denote an odd prime number. We design an algorithm which computes strong numerical evidence for the validity of the ``Equivariant Tamagawa Number Conjecture' for the -equivariant motive , where is a cyclic extension of degree and group . This conjecture is a very deep refinement of the classical analytic class number formula. In the course of the algorithm, we compute a set of special units which must be considered as a generalization of the (conjecturally existing) Stark units associated to first order vanishing Dirichlet -functions.

     

The convergence of harmonic Ritz values, harmonic Ritz vectors, and refined harmonic Ritz vectors

This paper concerns a harmonic projection method for computing an approximation to an eigenpair of a large matrix . Given a target point and a subspace that contains an approximation to , the harmonic projection method returns an approximation to . Three convergence results are established as the deviation of from approaches zero. First, the harmonic Ritz value converges to if a certain Rayleigh quotient matrix is uniformly nonsingular. Second, the harmonic Ritz vector converges to if the Rayleigh quotient matrix is uniformly nonsingular and remains well separated from the other harmonic Ritz values. Third, better error bounds for the convergence of are derived when converges. However, we show that the harmonic projection method can fail to find the desired eigenvalue --in other words, the method can miss if it is very close to . To this end, we propose to compute the Rayleigh quotient of with respect to and take it as a new approximate eigenvalue. is shown to converge to once tends to , no matter how is close to . Finally, we show that if the Rayleigh quotient matrix is uniformly nonsingular, then the refined harmonic Ritz vector, or more generally the refined eigenvector approximation introduced by the author, converges. We construct examples to illustrate our theory.

     

New quadratic polynomials with high densities of prime values

Hardy and Littlewood's Conjecture F implies that the asymptotic density of prime values of the polynomials , is related to the discriminant of via a quantity The larger is, the higher the asymptotic density of prime values for any quadratic polynomial of discriminant . A technique of Bach allows one to estimate accurately for any , given the class number of the imaginary quadratic order with discriminant , and for any 0$"> given the class number and regulator of the real quadratic order with discriminant . The Manitoba Scalable Sieve Unit (MSSU) has shown us how to rapidly generate many discriminants for which is potentially large, and new methods for evaluating class numbers and regulators of quadratic orders allow us to compute accurate estimates of efficiently, even for values of with as many as decimal digits. Using these methods, we were able to find a number of discriminants for which, under the assumption of the Extended Riemann Hypothesis, is larger than any previously known examples.

     

A new approximation technique for div-curl systems

In this paper, we describe an approximation technique for div-curl systems based in where is a domain in . We formulate this problem as a general variational problem with different test and trial spaces. The analysis requires the verification of an appropriate inf-sup condition. This results in a very weak formulation where the solution space is and the data reside in various negative norm spaces. Subsequently, we consider finite element approximations based on this weak formulation. The main approach of this paper involves the development of ``stable pairs" of discrete test and trial spaces. With this approach, we enlarge the test space so that the discrete inf-sup condition holds and we use a negative-norm least-squares formulation to reduce to a uniquely solvable linear system. This leads to optimal order estimates for problems with minimal regularity which is important since it is possible to construct magnetostatic field problems whose solutions have low Sobolev regularity (e.g., with ). The resulting algebraic equations are symmetric, positive definite and well conditioned. A second approach using a smaller test space which adds terms to the form for stabilization will also be mentioned. Some numerical results are also presented.

     

All first-order averaging techniques for a posteriori finite element error control on unstructured grids are efficient and reliable

All first-order averaging or gradient-recovery operators for lowest-order finite element methods are shown to allow for an efficient a posteriori error estimation in an isotropic, elliptic model problem in a bounded Lipschitz domain in . Given a piecewise constant discrete flux (that is the gradient of a discrete displacement) as an approximation to the unknown exact flux (that is the gradient of the exact displacement), recent results verify efficiency and reliability of

in the sense that is a lower and upper bound of the flux error up to multiplicative constants and higher-order terms. The averaging space consists of piecewise polynomial and globally continuous finite element functions in components with carefully designed boundary conditions. The minimal value is frequently replaced by some averaging operator applied within a simple post-processing to . The result provides a reliable error bound with .

This paper establishes and so equivalence of and . This implies efficiency of for a large class of patchwise averaging techniques which includes the ZZ-gradient-recovery technique. The bound established for tetrahedral finite elements appears striking in that the shape of the elements does not enter: The equivalence is robust with respect to anisotropic meshes. The main arguments in the proof are Ascoli's lemma, a strengthened Cauchy inequality, and elementary calculations with mass matrices.

     

-error estimates for ``shifted' surface spline interpolation on Sobolev space

The accuracy of interpolation by a radial basis function is usually very satisfactory provided that the approximant is reasonably smooth. However, for functions which have smoothness below a certain order associated with the basis function , no approximation power has yet been established. Hence, the purpose of this study is to discuss the -approximation order ( ) of interpolation to functions in the Sobolev space with \max(0,d/2-d/p)$">. We are particularly interested in using the ``shifted' surface spline, which actually includes the cases of the multiquadric and the surface spline. Moreover, we show that the accuracy of the interpolation method can be at least doubled when additional smoothness requirements and boundary conditions are met.

     

A new bound for the smallest

Let denote the number of primes and let denote the usual integral logarithm of . We prove that there are at least integer values of in the vicinity of with . This improves earlier bounds of Skewes, Lehman, and te Riele. We also plot more than 10000 values of in four different regions, including the regions discovered by Lehman, te Riele, and the authors of this paper, and a more distant region in the vicinity of , where appears to exceed by more than . The plots strongly suggest, although upper bounds derived to date for are not sufficient for a proof, that exceeds for at least integers in the vicinity of . If it is possible to improve our bound for by finding a sign change before , our first plot clearly delineates the potential candidates. Finally, we compute the logarithmic density of and find that as departs from the region in the vicinity of , the density is , and that it varies from this by no more than over the next integers. This should be compared to Rubinstein and Sarnak.

     

Multigrid methods for the computation of singular solutions and stress intensity factors I: Corner singularities

We consider the Poisson equation with homogeneous Dirichlet boundary condition on a two-dimensional polygonal domain with re-entrant angles. A multigrid method for the computation of singular solutions and stress intensity factors using piecewise linear functions is analyzed. When , the rate of convergence to the singular solution in the energy norm is shown to be , and the rate of convergence to the stress intensity factors is shown to be , where is the largest re-entrant angle of the domain and can be arbitrarily small. The cost of the algorithm is . When , the algorithm can be modified so that the convergence rate to the stress intensity factors is . In this case the maximum error of the multigrid solution over the vertices of the triangulation is shown to be .

     

Analysis of PSLQ, an integer relation finding algorithm

Let be either the real, complex, or quaternion number system and let be the corresponding integers. Let be a vector in . The vector has an integer relation if there exists a vector , , such that . In this paper we define the parameterized integer relation construction algorithm PSLQ, where the parameter can be freely chosen in a certain interval. Beginning with an arbitrary vector , iterations of PSLQ will produce lower bounds on the norm of any possible relation for . Thus PSLQ can be used to prove that there are no relations for of norm less than a given size. Let be the smallest norm of any relation for . For the real and complex case and each fixed parameter in a certain interval, we prove that PSLQ constructs a relation in less than iterations.

     

Merging the Bramble-Pasciak-Steinbach and the Crouzeix-Thomée criterion for -projection onto finite element spaces

Suppose is a finite-dimensional linear space based on a triangulation of a domain , and let denote the -projection onto . Provided the mass matrix of each element and the surrounding mesh-sizes obey the inequalities due to Bramble, Pasciak, and Steinbach or that neighboring element-sizes obey the global growth-condition due to Crouzeix and Thomée, is -stable: For all we have with a constant that is independent of, e.g., the dimension of .

This paper provides a more flexible version of the Bramble-Pasciak- Steinbach criterion for -stability on an abstract level. In its general version, (i) the criterion is applicable to all kind of finite element spaces and yields, in particular, -stability for nonconforming schemes on arbitrary (shape-regular) meshes; (ii) it is weaker than (i.e., implied by) either the Bramble-Pasciak-Steinbach or the Crouzeix-Thomée criterion for regular triangulations into triangles; (iii) it guarantees -stability of a priori for a class of adaptively-refined triangulations into right isosceles triangles.

     

On quadratic fields with large 3-rank

Davenport and Heilbronn defined a bijection between classes of binary cubic forms and classes of cubic fields, which has been used to tabulate the latter. We give a simpler proof of their theorem then analyze and improve the table-building algorithm. It computes the multiplicities of the general cubic discriminants (real or imaginary) up to in time and space , or more generally in time and space for a freely chosen positive . A variant computes the -ranks of all quadratic fields of discriminant up to with the same time complexity, but using only units of storage. As an application we obtain the first real quadratic fields with , and prove that is the smallest imaginary quadratic field with -rank equal to .

     

Detecting perfect powers by factoring into coprimes

This paper presents an algorithm that, given an integer , finds the largest integer such that is a th power. A previous algorithm by the first author took time where ; more precisely, time ; conjecturally, time . The new algorithm takes time . It relies on relatively complicated subroutines--specifically, on the first author's fast algorithm to factor integers into coprimes--but it allows a proof of the bound without much background; the previous proof of relied on transcendental number theory.

The computation of is the first step, and occasionally the bottleneck, in many number-theoretic algorithms: the Agrawal-Kayal-Saxena primality test, for example, and the number-field sieve for integer factorization.

     

Continuous-time Kreiss resolvent condition on infinite-dimensional spaces

Given the infinitesimal generator of a -semigroup on the Banach space which satisfies the Kreiss resolvent condition, i.e., there exists an such that for all complex with positive real part, we show that for general Banach spaces this condition does not give any information on the growth of the associated -semigroup. For Hilbert spaces the situation is less dramatic. In particular, we show that the semigroup can grow at most like . Furthermore, we show that for every there exists an infinitesimal generator satisfying the Kreiss resolvent condition, but whose semigroup grows at least like . As a consequence, we find that for with the standard Euclidian norm the estimate cannot be replaced by a lower power of or .

     

New techniques for bounds on the total number of prime factors of an odd perfect number

Let denote the sum of the positive divisors of . We say that is perfect if . Currently there are no known odd perfect numbers. It is known that if an odd perfect number exists, then it must be of the form , where are distinct primes and . Define the total number of prime factors of as . Sayers showed that . This was later extended by Iannucci and Sorli to show that . This was extended by the author to show that . Using an idea of Carl Pomerance this paper extends these results. The current new bound is .

     

On the stability of the

We prove the stability in of the projection onto a family of finite element spaces of conforming piecewise linear functions satisfying certain local mesh conditions. We give explicit formulae to check these conditions for a given finite element mesh in any number of spatial dimensions. In particular, stability of the projection in holds for locally quasiuniform geometrically refined meshes as long as the volume of neighboring elements does not change too drastically.