The PSLQ Algorithm for Empirical Data
The celebrated integer relation finding algorithm PSLQ has been successfully used in many applications. However, the PSLQ was only analyzed theoretically for exact input. When the input data are irrational numbers, they must be approximate ones due to finite precision in computer. That is, when the algorithm takes empirical data (inexact data with error bounded) instead of exact real numbers as its input, how do we ensure theoretically the output of the algorithm to be an exact integer relation? In this paper, we investigate the PSLQ algorithm for empirical data as its input. First, we give a termination condition for this case. Secondly we analyze a perturbation on the hyperplane matrix constructed from the input data and hence disclose a relationship between the accuracy of the input data and the output quality (an upper bound on the absolute value of the inner product of the exact data and the computed integer relation). Further, we also analyze the computational complexity for PSLQ with empirical data. Examples on transcendental numbers and algebraic numbers show the meaningfulness of our error control strategies.
READ FULL TEXT