Let {x 1 , x 2 , · · · , x n } be a vector of real numbers. An integer relation algorithm is a computational scheme to find the n integers a k , if they exist, such that a 1 x 1 + a 2 x 2 + · · · + a n x n = 0. In the past few years, integer relation algorithms have been utilized to discover new results in mathematics and physics. Existing programs for this purpose require very large amounts of computer time, due in part to the requirement for multiprecision arithmetic, yet are poorly suited for parallel processing. This paper presents a new integer relation algorithm designed for parallel computer systems, but as a bonus it also gives superior results on single processor systems. Single-and multi-level implementations of this algorithm are described, together with performance results on a parallel computer system. Several applications of these programs are discussed, including some new results in number theory, quantum field theory and chaos theory.
IntroductionLet x = (x 1 , x 2 , · · · , x n ) be a vector of real numbers. x is said to possess an integer relation if there exist integers a i , not all zero, such that a 1 x 1 + a 2 x 2 + · · · + a n x n = 0. By an integer relation algorithm, we mean a practical computational scheme that can recover (provided the computer implementation has sufficient numeric precision) the vector of integers a i , if it exists, or can produce bounds within which no integer relation exists.The problem of finding integer relations among a set of real numbers was first studied by Euclid, who gave an iterative scheme which, when applied to two real numbers, either terminates, yielding an exact relation, or produces an infinite sequence of approximate relations. The generalization of this problem for n > 2 was attempted by Euler, Jacobi, Poincaré, Minkowski, Perron, Brun, Bernstein, among others. The first integer relation algorithm with the required properties mentioned above was discovered in 1977 by Ferguson and Forcade [18]. Since then, a number of other integer relation algorithms have been discovered, including the "HJLS" algorithm [19] (which is based on the LLL algorithm), and the "PSLQ" algorithm.