A continuity property of multivariate Lagrange interpolation

Let be a sequence of interpolation schemes in of degree (i.e. for each one has unique interpolation by a polynomial of total degree and total order . Suppose that the points of tend to as and the Lagrange-Hermite interpolants, , satisfy for all monomials with . Theorem: for all functions of class in a neighborhood of . (Here denotes the Taylor series of at 0 to order .) Specific examples are given to show the optimality of this result.

     

Optimal information for approximating periodic analytic functions

Let be a strip in the complex plane. For fixed integer let denote the class of -periodic functions , which are analytic in and satisfy in . Denote by the subset of functions from that are real-valued on the real axis. Given a function , we try to recover at a fixed point by an algorithm on the basis of the information

where , are the Fourier coefficients of . We find the intrinsic error of recovery

Furthermore the -dimensional optimal information error, optimal sampling error and -widths of in , the space of continuous functions on , are determined. The optimal sampling error turns out to be strictly greater than the optimal information error. Finally the same problems are investigated for the class , consisting of all -periodic functions, which are analytic in with -integrable boundary values. In the case sampling fails to yield optimal information as well in odd as in even dimensions.

     

Further tabulation of the Erdös-Selfridge function

For a positive integer , the Erdös-Selfridge function is the least integer such that all prime factors of exceed . This paper describes a rapid method of tabulating using VLSI based sieving hardware. We investigate the number of admissible residues for each modulus in the underlying sieving problem and relate this number to the size of . A table of values of for is provided.

     

Computational experiences on the distances of polynomials to irreducible polynomials

In this paper we deal with a problem of Turán concerning the `distance' of polynomials to irreducible polynomials. Using computational methods we prove that for any monic polynomial of degree there exists a monic polynomial with deg() = deg() such that is irreducible over and the `distance' of and is .

     

Bounds for Multiplicative Cosets over Fields of Prime Order

Let be a positive integer and suppose that is an odd prime with . Suppose that and consider the polynomial . If this polynomial has any roots in , where the coset representatives for are taken to be all integers with , then these roots will form a coset of the multiplicative subgroup of consisting of the th roots of unity mod . Let be a coset of in , and define . In the paper ``Numbers Having Small th Roots mod ' (Mathematics of Computation, Vol. 61, No. 203 (1993),pp. 393-413), Robinson gives upper bounds for of the form , where is the Euler phi-function. This paper gives lower bounds that are of the same form, and seeks to sharpen the constants in the upper bounds of Robinson. The upper bounds of Robinson are proven to be optimal when is a power of or when

     

On a Problem of Byrnes concerning Polynomials with Restricted Coefficients

We consider a question of Byrnes concerning the minimal degree of a polynomial with all coefficients in which has a zero of a given order at . For , we prove his conjecture that the monic polynomial of this type of minimal degree is given by , but we disprove this for . We prove that a polynomial of this type must have , which is in sharp contrast with the situation when one allows coefficients in . The proofs use simple number theoretic ideas and depend ultimately on the fact that .

     

On the rapid computation of various polylogarithmic constants

We give algorithms for the computation of the -th digit of certain transcendental numbers in various bases. These algorithms can be easily implemented (multiple precision arithmetic is not needed), require virtually no memory, and feature run times that scale nearly linearly with the order of the digit desired. They make it feasible to compute, for example, the billionth binary digit of or on a modest work station in a few hours run time. We demonstrate this technique by computing the ten billionth hexadecimal digit of , the billionth hexadecimal digits of and , and the ten billionth decimal digit of . These calculations rest on the observation that very special types of identities exist for certain numbers like , , and . These are essentially polylogarithmic ladders in an integer base. A number of these identities that we derive in this work appear to be new, for example the critical identity for :

     

On Wendt's determinant

Wendt's determinant of order is the circulant determinant whose -th entry is the binomial coefficient , for . We give a formula for , when is even not divisible by 6, in terms of the discriminant of a polynomial , with rational coefficients, associated to . In particular, when where is a prime , this yields a factorization of involving a Fermat quotient, a power of and the 6-th power of an integer.

     

On some inequalities for the incomplete gamma function

Let be a positive real number. We determine all real numbers and such that the inequalities

are valid for all . And, we determine all real numbers and such that

hold for all .

     

On searching for solutions of the Diophantine equation

We propose a new search algorithm to solve the equation for a fixed value of . By parametrizing min, this algorithm obtains and (if they exist) by solving a quadratic equation derived from divisors of . By using several efficient number-theoretic sieves, the new algorithm is much faster on average than previous straightforward algorithms. We performed a computer search for 51 values of below 1000 (except ) for which no solution has previously been found. We found eight new integer solutions for and in the range of .

     

On the asymptotic spectrum of Hermitian block Toeplitz matrices with Toeplitz blocks

We study the asymptotic behaviour of the eigenvalues of Hermitian block Toeplitz matrices , with Toeplitz blocks. Such matrices are generated by the Fourier coefficients of an integrable bivariate function , and we study their eigenvalues for large and , relating their behaviour to some properties of as a function; in particular we show that, for any fixed , the first eigenvalues of tend to , while the last tend to , so extending to the block case a well-known result due to Szegö. In the case the 's are positive-definite, we study the asymptotic spectrum of , where is a block Toeplitz preconditioner for the conjugate gradient method, applied to solve the system , obtaining strict estimates, when and are fixed, and exact limit values, when and tend to infinity, for both the condition number and the conjugate gradient convergence factor of the previous matrices. Extensions to the case of a deeper nesting level of the block structure are also discussed.

     

Computing Canonical Heights With Little (or No) Factorization

Let be an elliptic curve with discriminant , and let . The standard method for computing the canonical height is as a sum of local heights . There are well-known series for computing the archimedean height , and the non-archimedean heights are easily computed as soon as all prime factors of have been determined. However, for curves with large coefficients it may be difficult or impossible to factor . In this note we give a method for computing the non-archimedean contribution to which is quite practical and requires little or no factorization. We also give some numerical examples illustrating the algorithm.

     

New estimates for Ritz vectors

The following estimate for the Rayleigh-Ritz method is proved:

Here is a bounded self-adjoint operator in a real Hilbert/euclidian space, one of its eigenpairs, a trial subspace for the Rayleigh-Ritz method, and a Ritz pair. This inequality makes it possible to analyze the fine structure of the error of the Rayleigh-Ritz method, in particular, it shows that if an eigenvector is close to the trial subspace with accuracy and a Ritz vector is an approximation to another eigenvector, with a different eigenvalue. Generalizations of the estimate to the cases of eigenspaces and invariant subspaces are suggested, and estimates of approximation of eigenspaces and invariant subspaces are proved.

     

Equivalent formulae for the supremum and stability of weighted pseudoinverses

During recent decades, there have been a great number of research articles studying interior-point methods for solving problems in mathematical programming and constrained optimization. Stewart and O'Leary obtained an upper bound for scaled pseudoinverses of a matrix where is a set of diagonal positive definite matrices. We improved their results to obtain the supremum of scaled pseudoinverses and derived the stability property of scaled pseudoinverses. Forsgren further generalized these results to derive the supremum of weighted pseudoinverses where is a set of diagonally dominant positive semidefinite matrices, by using a signature decomposition of weighting matrices and by applying the Binet-Cauchy formula and Cramer's rule for determinants. The results are also extended to equality constrained linear least squares problems. In this paper we extend Forsgren's results to a general complex matrix to establish several equivalent formulae for , where is a set of diagonally dominant positive semidefinite matrices, or a set of weighting matrices arising from solving equality constrained least squares problems. We also discuss the stability property of these weighted pseudoinverses.

     

A search for Wieferich and Wilson primes

An odd prime is called a Wieferich prime if

alternatively, a Wilson prime if

To date, the only known Wieferich primes are and , while the only known Wilson primes are , and . We report that there exist no new Wieferich primes , and no new Wilson primes . It is elementary that both defining congruences above hold merely (mod ), and it is sometimes estimated on heuristic grounds that the ``probability" that is Wieferich (independently: that is Wilson) is about . We provide some statistical data relevant to occurrences of small values of the pertinent Fermat and Wilson quotients (mod ).

     

Pure product polynomials and the Prouhet-Tarry-Escott problem

An -factor pure product is a polynomial which can be expressed in the form for some natural numbers . We define the norm of a polynomial to be the sum of the absolute values of the coefficients. It is known that every -factor pure product has norm at least . We describe three algorithms for determining the least norm an -factor pure product can have. We report results of our computations using one of these algorithms which include the result that every -factor pure product has norm strictly greater than if is , , , or .

     

Density of Carmichael numbers with three prime factors

We get an upper bound of on the number of Carmichael numbers with exactly three prime factors.

     

A unified approach to evaluation algorithms for multivariate polynomials

We present a unified framework for most of the known and a few new evaluation algorithms for multivariate polynomials expressed in a wide variety of bases including the Bernstein-Bézier, multinomial (or Taylor), Lagrange and Newton bases. This unification is achieved by considering evaluation algorithms for multivariate polynomials expressed in terms of L-bases, a class of bases that include the Bernstein-Bézier, multinomial, and a rich subclass of Lagrange and Newton bases. All of the known evaluation algorithms can be generated either by considering up recursive evaluation algorithms for L-bases or by examining change of basis algorithms for L-bases. For polynomials of degree in variables, the class of up recursive evaluation algorithms includes a parallel up recurrence algorithm with computational complexity , a nested multiplication algorithm with computational complexity and a ladder recurrence algorithm with computational complexity . These algorithms also generate a new generalization of the Aitken-Neville algorithm for evaluation of multivariate polynomials expressed in terms of Lagrange L-bases. The second class of algorithms, based on certain change of basis algorithms between L-bases, include a nested multiplication algorithm with computational complexity , a divided difference algorithm, a forward difference algorithm, and a Lagrange evaluation algorithm with computational complexities , and per point respectively for the evaluation of multivariate polynomials at several points.

     

Numerical solution of isospectral flows

In this paper we are concerned with the problem of solving numerically isospectral flows. These flows are characterized by the differential equation

where is a symmetric matrix, is a skew-symmetric matrix function of and is the Lie bracket operator. We show that standard Runge-Kutta schemes fail in recovering the main qualitative feature of these flows, that is isospectrality, since they cannot recover arbitrary cubic conservation laws. This failure motivates us to introduce an alternative approach and establish a framework for generation of isospectral methods of arbitrarily high order.

     

Modular forms which behave like theta series

In this paper, we determine all modular forms of weights , , for the full modular group which behave like theta series, i.e., which have in their Fourier expansions, the constant term and all other Fourier coefficients are non-negative rational integers. In fact, we give convex regions in (resp. in ) for the cases (resp. for the cases ). Corresponding to each lattice point in these regions, we get a modular form with the above property. As an application, we determine the possible exceptions of quadratic forms in the respective dimensions.