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.

     

The global decay to discrete shocks for scalar monotone schemes

Given a family of discrete shocks of a monotone scheme, we prove that the discrete shock profile with rational shock speed is asymptotically stable with respect to the topology : if , then as under no restriction conditions of the initial data to the interval . The asymptotic wave profile is uniquely identified from the above family by a mass parameter.

     

Corrigenda and addition to ``Computer verification of the Ankeny-Artin-Chowla conjecture for all primes less than '

An error in the program for verifying the Ankeny-Artin-Chowla (AAC) conjecture is reported. As a result, in the case of primes which are , the AAC conjecture has been verified using a different multiple of the regulator of the quadratic field than was meant. However, since any multiple of this regulator is suitable for this purpose, provided that it is smaller than , the main result that the AAC conjecture is true for all the primes which are , remains valid.

As an addition, we have verified the AAC conjecture for all the primes between and , with the corrected program.

     

A convergence and stability study of the iterated Lubkin transformation and the -algorithm

In this paper we analyze the convergence and stability of the iterated Lubkin transformation and the -algorithm as these are being applied to sequences whose members behave like as , where and are complex scalars and is a nonnegative integer. We study the three different cases in which (i) , , and (logarithmic sequences), (ii) and (linear sequences), and (iii) (factorial sequences). We show that both methods accelerate the convergence of all three types of sequences. We show also that both methods are stable on linear and factorial sequences, and they are unstable on logarithmic sequences. On the basis of this analysis we propose ways of improving accuracy and stability in problematic cases. Finally, we provide a comparison of these results with analogous results corresponding to the Levin -transformation.

     

Lower bounds for the total stopping time of iterates

The total stopping time of a positive integer is the minimal number of iterates of the function needed to reach the value , and is if no iterate of reaches . It is shown that there are infinitely many positive integers having a finite total stopping time such that 6.14316 \log n.$"> The proof involves a search of trees to depth 60, A heuristic argument suggests that for any constant , a search of all trees to sufficient depth could produce a proof that there are infinitely many such that \gamma\log n.$">It would require a very large computation to search trees to a sufficient depth to produce a proof that the expected behavior of a ``random' iterate, which is occurs infinitely often.

     

Newman polynomials with prescribed vanishing and integer sets with distinct subset sums

We study the problem of determining the minimal degree of a polynomial that has all coefficients in and a zero of multiplicity at . We show that a greedy solution is optimal precisely when , mirroring a result of Boyd on polynomials with coefficients. We then examine polynomials of the form , where is a set of positive odd integers with distinct subset sums, and we develop algorithms to determine the minimal degree of such a polynomial. We determine that satisfies inequalities of the form . Last, we consider the related problem of finding a set of positive integers with distinct subset sums and minimal largest element and show that the Conway-Guy sequence yields the optimal solution for , extending some computations of Lunnon.

     

About the sharpness of the stability estimates in the Kreiss matrix theorem

One of the conditions in the Kreiss matrix theorem involves the resolvent of the matrices under consideration. This so-called resolvent condition is known to imply, for all , the upper bounds and . Here is the spectral norm, is the constant occurring in the resolvent condition, and the order of is equal to .

It is a long-standing problem whether these upper bounds can be sharpened, for all fixed 1$">, to bounds in which the right-hand members grow much slower than linearly with and with , respectively. In this paper it is shown that such a sharpening is impossible. The following result is proved: for each 0$">, there are fixed values 0, K>1$"> and a sequence of matrices , satisfying the resolvent condition, such that for .

The result proved in this paper is also relevant to matrices whose -pseudospectra lie at a distance not exceeding from the unit disk for all 0$">.

     

Sums of heights of algebraic numbers

For , we consider the set . The polynomials are in , with only mild restrictions, and is the Weil height of . We show that this set is dense in for some effectively computable limit point .

     

Nontrivial Galois module structure of cyclotomic fields

We say a tame Galois field extension with Galois group has trivial Galois module structure if the rings of integers have the property that is a free -module. The work of Greither, Replogle, Rubin, and Srivastav shows that for each algebraic number field other than the rational numbers there will exist infinitely many primes so that for each there is a tame Galois field extension of degree so that has nontrivial Galois module structure. However, the proof does not directly yield specific primes for a given algebraic number field For any cyclotomic field we find an explicit so that there is a tame degree extension with nontrivial Galois module structure.

     

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.

     

A fast algorithm for testing reducibility of trinomials mod 2 and some new primitive trinomials of degree 3021377

The standard algorithm for testing reducibility of a trinomial of prime degree over requires bits of memory. We describe a new algorithm which requires only bits of memory and significantly fewer memory references and bit-operations than the standard algorithm.

If is a Mersenne prime, then an irreducible trinomial of degree is necessarily primitive. We give primitive trinomials for the Mersenne exponents , , and . The results for extend and correct some computations of Kumada et al. The two results for are primitive trinomials of the highest known degree.

     

Strong tractability of multivariate integration using quasi-Monte Carlo algorithms

We study quasi-Monte Carlo algorithms based on low discrepancy sequences for multivariate integration. We consider the problem of how the minimal number of function evaluations needed to reduce the worst-case error from its initial error by a factor of depends on and the dimension . Strong tractability means that it does not depend on and is bounded by a polynomial in . The least possible value of the power of is called the -exponent of strong tractability. Sloan and Wozniakowski established a necessary and sufficient condition of strong tractability in weighted Sobolev spaces, and showed that the -exponent of strong tractability is between 1 and 2. However, their proof is not constructive.

In this paper we prove in a constructive way that multivariate integration in some weighted Sobolev spaces is strongly tractable with -exponent equal to 1, which is the best possible value under a stronger assumption than Sloan and Wozniakowski's assumption. We show that quasi-Monte Carlo algorithms using Niederreiter's -sequences and Sobol sequences achieve the optimal convergence order for any 0$"> independent of the dimension with a worst case deterministic guarantee (where is the number of function evaluations). This implies that strong tractability with the best -exponent can be achieved in appropriate weighted Sobolev spaces by using Niederreiter's -sequences and Sobol sequences.

     

Newton-Cotes integration for approximating Stieltjes (generalized Euler) constants

In the Laurent expansion

of the Riemann-Hurwitz zeta function, the coefficients are known as Stieltjes, or generalized Euler, constants. [When , (the Riemann zeta function), and .] We present a new approach to high-precision approximation of . Plots of our results reveal much structure in the growth of the generalized Euler constants. Our results when for , and when for (for such as 53/100, 1/2, etc.) suggest that published bounds on the growth of the Stieltjes constants can be much improved, and lead to several conjectures. Defining , we conjecture that is attained: for any given , for some (and similarly that, given and , is within of for infinitely many ). In addition we conjecture that satisfies for 1$">. We also conjecture that , a special case of a more general conjecture relating the values of and for . Finally, it is known that for . Using this to define for all real 0$">, we conjecture that for nonintegral , is precisely times the -th (Weyl) fractional derivative at of the entire function . We also conjecture that , now defined for all real arguments 0$">, is smooth. Our numerical method uses Newton-Cotes integration formulae for very high-degree interpolating polynomials; it differs in implementation from, but compares in error bounding to, Euler-Maclaurin summation based methods.

     

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

     

Enhanced accuracy by post-processing for finite element methods for hyperbolic equations

We consider the enhancement of accuracy, by means of a simple post-processing technique, for finite element approximations to transient hyperbolic equations. The post-processing is a convolution with a kernel whose support has measure of order one in the case of arbitrary unstructured meshes; if the mesh is locally translation invariant, the support of the kernel is a cube whose edges are of size of the order of only. For example, when polynomials of degree are used in the discontinuous Galerkin (DG) method, and the exact solution is globally smooth, the DG method is of order in the -norm, whereas the post-processed approximation is of order ; if the exact solution is in only, in which case no order of convergence is available for the DG method, the post-processed approximation converges with order in , where is a subdomain over which the exact solution is smooth. Numerical results displaying the sharpness of the estimates are presented.

     

Linear quintuple-product identities

In the first part of this paper, series and product representations of four single-variable triple products , , , and four single-variable quintuple products , , , are defined. Reduced forms and reduction formulas for these eight functions are given, along with formulas which connect them. The second part of the paper contains a systematic computer search for linear trinomial identities. The complete set of such families is found to consist of two 2-parameter families, which are proved using the formulas in the first part of the paper.

     

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.

     

On the convergence of certain Gauss-type quadrature formulas for unbounded intervals

We consider the convergence of Gauss-type quadrature formulas for the integral , where is a weight function on the half line . The -point Gauss-type quadrature formulas are constructed such that they are exact in the set of Laurent polynomials }, where is a sequence of integers satisfying and . It is proved that under certain Carleman-type conditions for the weight and when or goes to , then convergence holds for all functions for which is integrable on . Some numerical experiments compare the convergence of these quadrature formulas with the convergence of the classical Gauss quadrature formulas for the half line.

     

Comparison theorems of Kolmogorov type and exact values of -widths on Hardy-Sobolev classes

Let be a strip in complex plane. denotes those -periodic, real-valued functions on which are analytic in the strip and satisfy the condition , . Osipenko and Wilderotter obtained the exact values of the Kolmogorov, linear, Gel'fand, and information -widths of in , , and 2-widths of in , , .

In this paper we continue their work. Firstly, we establish a comparison theorem of Kolmogorov type on , from which we get an inequality of Landau-Kolmogorov type. Secondly, we apply these results to determine the exact values of the Gel'fand -width of in , . Finally, we calculate the exact values of Kolmogorov -width, linear -width, and information -width of in , , .