Calculating the logarithm rapidly and accurately of points on the standard normal curve--a normal curve with mean of zero and standard deviation of one--is required for applications that involve repeated determinations of the logarithm. This paper points out the relative lack of discussion of the subject in the literature, surveys what there is of its history, presents the algorithm, describes three of its implementations along with empirical results, and provides the means of selecting one of the algorithms, based on requirements of speed and accuracy.
In a brief but insightful review of the literature calculating the logarithm of the standard normal distribution, Linhart identifies a problem with accuracy in the tails of the curve and presents three candidate solutions to remedy the problem. One uses continued fractions, another uses rational Chebyshev approximations, and the third uses a compiler resident error function.
Although the paper focuses on C implementations of the algorithm, a major feature is the analysis of the algorithm independent of its implementation; it makes a contribution to the literature by noting the negative effect of points at the extremes. One of the goals of the paper is to isolate ways to moderate this effect. A valuable feature of the paper is its presentation of testing results regarding error ranges.
The paper provides a valuable resource for applied statisticians and for programmers interested in the implementation of statistical algorithms.