Biquadratic reciprocity and a Lucasian primality test

Let be the sequence defined from a given initial value, the seed, , by the recurrence . Then, for a suitable seed , the number (where is odd) is prime iff . In general depends both on and on . We describe a slight modification of this test which determines primality of numbers with a seed which depends only on , provided . In particular, when , odd, we have a test with a single seed depending only on , in contrast with the unmodified test, which, as proved by W. Bosma in Explicit primality criteria for , Math. Comp. 61 (1993), 97-109, needs infinitely many seeds. The proof of validity uses biquadratic reciprocity.

     

Inverse inequalities on non-quasi-uniform meshes and application to the mortar element method

We present a range of mesh-dependent inequalities for piecewise constant and continuous piecewise linear finite element functions defined on locally refined shape-regular (but possibly non-quasi-uniform) meshes. These inequalities involve norms of the form for positive and negative and , where is a function which reflects the local mesh diameter in an appropriate way. The only global parameter involved is , the total number of degrees of freedom in the finite element space, and we avoid estimates involving either the global maximum or minimum mesh diameter. Our inequalities include new variants of inverse inequalities as well as trace and extension theorems. They can be used in several areas of finite element analysis to extend results--previously known only for quasi-uniform meshes--to the locally refined case. Here we describe applications to (i) the theory of nonlinear approximation and (ii) the stability of the mortar element method for locally refined meshes.

     

Bounds for computing the tame kernel

The tame kernel of the of a number field  is the kernel of some explicit map , where the product runs over all finite primes  of  and is the residue class field at . When is a set of primes of , containing the infinite ones, we can consider the -unit group  of . Then has a natural image in . The tame kernel is contained in this image if  contains all finite primes of  up to some bound. This is a theorem due to Bass and Tate. An explicit bound for imaginary quadratic fields was given by Browkin. In this article we give a bound, valid for any number field, that is smaller than Browkin's bound in the imaginary quadratic case and has better asymptotics. A simplified version of this bound says that we only have to include in  all primes with norm up to  , where  is the discriminant of . Using this bound, one can find explicit generators for the tame kernel, and a ``long enough' search would also yield all relations. Unfortunately, we have no explicit formula to describe what ``long enough' means. However, using theorems from Keune, we can show that the tame kernel is computable.

     

Weak coupling of solutions of first-order least-squares method

A theoretical analysis of a first-order least-squares finite element method for second-order self-adjoint elliptic problems is presented. We investigate the coupling effect of the approximate solutions for the primary function and for the flux . We prove that the accuracy of the approximate solution for the primary function is weakly affected by the flux . That is, the bound for is dependent on , but only through the best approximation for multiplied by a factor of meshsize . Similarly, we provide that the bound for is dependent on , but only through the best approximation for multiplied by a factor of the meshsize . This weak coupling is not true for the non-selfadjoint case. We provide the numerical experiment supporting the theorems in this paper.

     

Finite element superconvergence on Shishkin mesh for 2-D convection-diffusion problems

In this work, the bilinear finite element method on a Shishkin mesh for convection-diffusion problems is analyzed in the two-dimensional setting. A superconvergence rate in a discrete -weighted energy norm is established under certain regularity assumptions. This convergence rate is uniformly valid with respect to the singular perturbation parameter . Numerical tests indicate that the rate is sharp for the boundary layer terms. As a by-product, an -uniform convergence of the same order is obtained for the -norm. Furthermore, under the same regularity assumption, an -uniform convergence of order in the norm is proved for some mesh points in the boundary layer region.

     

Lax theorem and finite volume schemes

This work addresses a theory of convergence for finite volume methods applied to linear equations. A non-consistent model problem posed in an abstract Banach space is proved to be convergent. Then various examples show that the functional framework is non-empty. Convergence with a rate of all TVD schemes for linear advection in 1D is an application of the general result. Using duality techniques and assuming enough regularity of the solution, convergence of the upwind finite volume scheme for linear advection on a 2D triangular mesh is proved in , : provided the solution is in , it proves a rate of convergence in .

     

Counting primes in residue classes

We explain how the Meissel-Lehmer-Lagarias-Miller-Odlyzko method for computing can be used for computing efficiently , the number of primes congruent to modulo up to . As an application, we computed the number of prime numbers of the form less than for several values of up to and found a new region where is less than near .

     

Numerical simulation of stochastic evolution equations associated to quantum Markov semigroups

We address the problem of approximating numerically the solutions of stochastic evolution equations on Hilbert spaces , with respect to Brownian motions, arising in the unraveling of backward quantum master equations. In particular, we study the computation of mean values of , where is a linear operator. First, we introduce estimates on the behavior of . Then we characterize the error induced by the substitution of with the solution of a convenient stochastic ordinary differential equation. It allows us to establish the rate of convergence of to , where denotes the explicit Euler method. Finally, we consider an extrapolation method based on the Euler scheme. An application to the quantum harmonic oscillator system is included.

     

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.

     

An estimate for the number of integers without large prime factors

denotes the number of positive integers and free of prime factors y$">. Hildebrand and Tenenbaum provided a good approximation of . However, their method requires the solution to the equation , and therefore it needs a large amount of time for the numerical solution of the above equation for large . Hildebrand also showed approximates for , where and is the unique solution to . Let be defined by for 0$">. We show approximates , and also approximates , where . Using these approximations, we give a simple method which approximates within a factor in the range , where is any positive constant.

     

Elliptic curves with nonsplit mod representations

We calculate explicitly the -invariants of the elliptic curves corresponding to rational points on the modular curve by giving an expression defined over of the -function in terms of the function field generators and of the elliptic curve . As a result we exhibit infinitely many elliptic curves over with nonsplit mod representations.

     

On the total number of prime factors of an odd perfect number

We say is perfect if , where denotes the sum of the positive divisors of . No odd perfect numbers are known, but it is well known that if such a number exists, it must have prime factorization of the form , where , , ..., are distinct primes and . We prove that if or for all , , then . We also prove as our main result that , where . This improves a result of Sayers given in 1986.

     

Computational scales of Sobolev norms with application to preconditioning

This paper provides a framework for developing computationally efficient multilevel preconditioners and representations for Sobolev norms. Specifically, given a Hilbert space and a nested sequence of subspaces , we construct operators which are spectrally equivalent to those of the form . Here , , are positive numbers and is the orthogonal projector onto with . We first present abstract results which show when is spectrally equivalent to a similarly constructed operator defined in terms of an approximation of , for . We show that these results lead to efficient preconditioners for discretizations of differential and pseudo-differential operators of positive and negative order. These results extend to sums of operators. For example, singularly perturbed problems such as can be preconditioned uniformly independently of the parameter . We also show how to precondition an operator which results from Tikhonov regularization of a problem with noisy data. Finally, we describe how the technique provides computationally efficient bounded discrete extensions which have applications to domain decomposition.

     

Real zeros of real odd Dirichlet -functions

Let be a real odd Dirichlet character of modulus , and let be the associated Dirichlet -function. As a consequence of the work of Low and Purdy, it is known that if and , , , then has no positive real zeros. By a simple extension of their ideas and the advantage of thirty years of advances in computational power, we are able to prove that if , then has no positive real zeros.

     

The Euler implicit/explicit scheme for the 2D time-dependent Navier-Stokes equations with smooth or non-smooth initial data

This paper considers the stability and convergence results for the Euler implicit/explicit scheme applied to the spatially discretized two-dimensional (2D) time-dependent Navier-Stokes equations. A Galerkin finite element spatial discretization is assumed, and the temporal treatment is implicit/explict scheme, which is implicit for the linear terms and explicit for the nonlinear term. Here the stability condition depends on the smoothness of the initial data , i.e., the time step condition is in the case of , in the case of and in the case of for mesh size and some positive constant . We provide the -stability of the scheme under the stability condition with and obtain the optimal error estimate of the numerical velocity and the optimal error estimate of the numerical pressure under the stability condition with .

     

Chebyshev's bias for composite numbers with restricted prime divisors

Let denote the number of primes with . Chebyshev's bias is the phenomenon for which ``more often' \pi(x;d,r)$">, than the other way around, where is a quadratic nonresidue mod and is a quadratic residue mod . If for every up to some large number, then one expects that for every . Here denotes the number of integers such that every prime divisor of satisfies . In this paper we develop some tools to deal with this type of problem and apply them to show that, for example, for every .

In the process we express the so-called second order Landau-Ramanujan constant as an infinite series and show that the same type of formula holds for a much larger class of constants.

     

Stabilized finite element method based on the Crank-Nicolson extrapolation scheme for the time-dependent Navier-Stokes equations

This paper provides an error analysis for the Crank-Nicolson extrapolation scheme of time discretization applied to the spatially discrete stabilized finite element approximation of the two-dimensional time-dependent Navier-Stokes problem, where the finite element space pair for the approximation of the velocity and the pressure is constructed by the low-order finite element: the quadrilateral element or the triangle element with mesh size . Error estimates of the numerical solution to the exact solution with are derived.

     

Asymptotics of recurrence coefficients for orthonormal polynomials on the line--Magnus's method revisited

We use Freud equations to obtain the main term in the asymptotic expansion of the recurrence coefficients associated with orthonormal polynomials for weights on the real line where is an even polynomial of fixed degree with nonnegative coefficients or where . Here for some real -1$">.