Littlewood polynomials with high order zeros

Let be the minimal length of a polynomial with coefficients divisible by . Byrnes noted that for each , and asked whether in fact . Boyd showed that for all , but . He further showed that , and that is one of the 5 numbers , or . Here we prove that . Similarly, let be the maximal power of dividing some polynomial of degree with coefficients. Boyd was able to find for . In this paper we determine for .

     

Analysis of some low order quadrilateral Reissner-Mindlin plate elements

Four quadrilateral elements for the Reissner-Mindlin plate model are considered. The elements are the stabilized MITC4 element of Lyly, Stenberg and Vihinen  (1933), the MIN4 element of Tessler and Hughes (1983), the Q4BL element of Zienkiewicz et al. (1993) and the FMIN4 element of Kikuchi and Ishii (1999). For all elements except the Q4BL element, a unifying variational formulation is introduced, and optimal H and L error bounds uniform in the plate thickness are proven. Moreover, we propose a modified Q4BL element and show that it admits the optimal H and L error bounds uniform in the plate thickness. In particular, we study the convergence behavior of all elements regarding the mesh distortion.

     

Integer transfinite diameter and polynomials with small Mahler measure

In this work, we show how suitable generalizations of the integer transfinite diameter of some compact sets in give very good bounds for coefficients of polynomials with small Mahler measure. By this way, we give the list of all monic irreducible primitive polynomials of of degree at most with Mahler measure less than and of degree and with Mahler measure less than .

     

Limiting set of second order spectra

Let be a self-adjoint operator acting on a Hilbert space . A complex number is in the second order spectrum of relative to a finite-dimensional subspace iff the truncation to of is not invertible. This definition was first introduced in Davies, 1998, and according to the results of Levin and Shargorodsky in 2004, these sets provide a method for estimating eigenvalues free from the problems of spectral pollution. In this paper we investigate various aspects related to the issue of approximation using second order spectra. Our main result shows that under fairly mild hypothesis on the uniform limit of these sets, as increases towards , contain the isolated eigenvalues of of finite multiplicity. Therefore, unlike the majority of the standard methods, second order spectra combine nonpollution and approximation at a very high level of generality.

     

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.

     

Analysis of a bilinear finite element for shallow shells. II: Consistency error

We consider a bilinear reduced-strain finite element of the MITC family for a shallow Reissner-Naghdi type shell. We estimate the consistency error of the element in both membrane- and bending-dominated states of deformation. We prove that in the membrane-dominated case, under severe assumptions on the domain, the finite element mesh and the regularity of the solution, an error bound can be obtained if the contribution of transverse shear is neglected. Here is the thickness of the shell, the mesh spacing, and a smoothness parameter. In the bending-dominated case, the uniformly optimal bound is achievable but requires that membrane and transverse shear strains are of order as . In this case we also show that under sufficient regularity assumptions the asymptotic consistency error has the bound .

     

On the optimal convergence rate of universal and nonuniversal algorithms for multivariate integration and approximation

We study the maximal rate of convergence (mrc) of algorithms for (multivariate) integration and approximation of -variate functions from reproducing kernel Hilbert spaces . Here is an arbitrary kernel all of whose partial derivatives up to order satisfy a Hölder-type condition with exponent . Algorithms use function values and we analyze their rate of convergence as tends to infinity. We focus on universal algorithms which depend on , , and but not on the specific kernel , and nonuniversal algorithms which may depend additionally on .

For universal algorithms the mrc is for both integration and approximation, and for nonuniversal algorithms it is for integration and with for approximation. Hence, the mrc for universal algorithms suffers from the curse of dimensionality if is large relative to , whereas the mrc for nonuniversal algorithms does not since it is always at least for integration, and for approximation. This is the price we have to pay for using universal algorithms. On the other hand, if is large relative to , then the mrc for universal and nonuniversal algorithms is approximately the same.

We also consider the case when we have the additional knowledge that the kernel has product structure, . Here are some univariate kernels whose all derivatives up to order satisfy a Hölder-type condition with exponent . Then the mrc for universal algorithms is for both integration and approximation, and for nonuniversal algorithms it is for integration and with for approximation. If or for all , then the mrc is at least , and the curse of dimensionality is not present. Hence, the product form of reproducing kernels breaks the curse of dimensionality even for universal algorithms.

     

Optimal anisotropic meshes for minimizing interpolation errors in -norm

In this paper, we present a new optimal interpolation error estimate in norm ( ) for finite element simplicial meshes in any spatial dimension. A sufficient condition for a mesh to be nearly optimal is that it is quasi-uniform under a new metric defined by a modified Hessian matrix of the function to be interpolated. We also give new functionals for the global moving mesh method and obtain optimal monitor functions from the viewpoint of minimizing interpolation error in the norm. Some numerical examples are also given to support the theoretical estimates.

     

The 192 solutions of the Heun equation

A machine-generated list of local solutions of the Heun equation is given. They are analogous to Kummer's  solutions of the Gauss hypergeometric equation, since the two equations are canonical Fuchsian differential equations on the Riemann sphere with four and three singular points, respectively. Tabulation is facilitated by the identification of the automorphism group of the equation with  singular points as the Coxeter group  . Each of the expressions is labeled by an element of  . Of the ,  are equivalent expressions for the local Heun function  , and it is shown that the resulting order- group of transformations of  is isomorphic to the symmetric group . The isomorphism encodes each transformation as a permutation of an abstract four-element set, not identical to the set of singular points.

     

Modular exponentiation via the explicit Chinese remainder theorem

Fix pairwise coprime positive integers . We propose representing integers modulo , where is any positive integer up to roughly , as vectors . We use this representation to obtain a new result on the parallel complexity of modular exponentiation: there is an algorithm for the Common CRCW PRAM that, given positive integers , , and in binary, of total bit length , computes in time using processors. For comparison, a parallelization of the standard binary algorithm takes superlinear time; Adleman and Kompella gave an expected time algorithm using processors; von zur Gathen gave an NC algorithm for the highly special case that is polynomially smooth.

     

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.

     

Complexity of inverting the Euler function

Given an integer , how hard is it to find the set of all integers such that , where is the Euler totient function? We present a certain basic algorithm which, given the prime number factorization of , in polynomial time ``on average' (that is, ), finds the set of all such solutions . In fact, in the worst case this set of solutions is exponential in , and so cannot be constructed by a polynomial time algorithm. In the opposite direction, we show, under a widely accepted number theoretic conjecture, that the PARTITION PROBLEM, an NP-complete problem, can be reduced in polynomial (in the input size) time to the problem of deciding whether has a solution, for polynomially (in the input size of the PARTITION PROBLEM) many values of (where the prime factorizations of these are given). What this means is that the problem of deciding whether there even exists a solution to , let alone finding any or all such solutions, is very likely to be intractable. Finally, we establish close links between the problem of inverting the Euler function and the integer factorization problem.

     

Uniform error estimates in the finite element method for a singularly perturbed reaction-diffusion problem

Consider the problem with homogeneous Neumann boundary condition in a bounded smooth domain in . The whole range is treated. The Galerkin finite element method is used on a globally quasi-uniform mesh of size ; the mesh is fixed and independent of .

A precise analysis of how the error at each point depends on and is presented. As an application, first order error estimates in , which are uniform with respect to , are given.

     

Computing the arithmetic genus of Hilbert modular fourfolds

The Hilbert modular fourfold determined by the totally real quartic number field is a desingularization of a natural compactification of the quotient space , where PSL acts on by fractional linear transformations via the four embeddings of into . The arithmetic genus, equal to one plus the dimension of the space of Hilbert modular cusp forms of weight , is a birational invariant useful in the classification of these varieties. In this work, we describe an algorithm allowing for the automated computation of the arithmetic genus and give sample results.

     

Practical solution of the Diophantine equation

Let be an odd prime and , positive integers. In this note we prove that the problem of the determination of the integer solutions to the equation can be easily reduced to the resolution of the unit equation over . The solutions of the latter equation are given by Wildanger's algorithm.

     

Order of convergence of second order schemes based on the minmod limiter

Many second order accurate nonoscillatory schemes are based on the minmod limiter, e.g., the Nessyahu-Tadmor scheme. It is well known that the -error of monotone finite difference methods for the linear advection equation is of order for initial data in , . For second or higher order nonoscillatory schemes very little is known because they are nonlinear even for the simple advection equation. In this paper, in the case of a linear advection equation with monotone initial data, it is shown that the order of the -error for a class of second order schemes based on the minmod limiter is of order at least in contrast to the order for any formally first order scheme.

     

A dual finite element complex on the barycentric refinement

Given a two dimensional oriented surface equipped with a simplicial mesh, the standard lowest order finite element spaces provide a complex centered on Raviart-Thomas divergence conforming vector fields. It can be seen as a realization of the simplicial cochain complex. We construct a new complex of finite element spaces on the barycentric refinement of the mesh which can be seen as a realization of the simplicial chain complex on the original (unrefined) mesh, such that the duality is non-degenerate on for each . In particular is a space of -conforming vector fields which is dual to Raviart-Thomas -conforming elements. When interpreted in terms of differential forms, these two complexes provide a finite-dimensional analogue of Hodge duality.

     

On splitting methods for Schrödinger-Poisson and cubic nonlinear Schrödinger equations

We give an error analysis of Strang-type splitting integrators for nonlinear Schrödinger equations. For Schrödinger-Poisson equations with an -regular solution, a first-order error bound in the norm is shown and used to derive a second-order error bound in the norm. For the cubic Schrödinger equation with an -regular solution, first-order convergence in the norm is used to obtain second-order convergence in the norm. Basic tools in the error analysis are Lie-commutator bounds for estimating the local error and -conditional stability for error propagation, where for the Schrödinger-Poisson system and for the cubic Schrödinger equation.

     

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.

     

Efficient computation of root numbers and class numbers of parametrized families of real abelian number fields

Let be a parametrized family of simplest real cyclic cubic, quartic, quintic or sextic number fields of known regulators, e.g., the so-called simplest cubic and quartic fields associated with the polynomials and . We give explicit formulas for powers of the Gaussian sums attached to the characters associated with these simplest number fields. We deduce a method for computing the exact values of these Gaussian sums. These values are then used to efficiently compute class numbers of simplest fields. Finally, such class number computations yield many examples of real cyclotomic fields of prime conductors and class numbers greater than or equal to . However, in accordance with Vandiver's conjecture, we found no example of for which divides .