Thelevel-index (LI) representation of numbers, and itsalgorithms forarithmetic operations, were introduced byCharles Clenshaw andFrank Olver in 1984.[1]
The symmetric form of the LI system and its arithmetic operations were presented by Clenshaw and Peter Turner in 1987.[2]
Michael Anuta, Daniel Lozier, Nicolas Schabanel and Turner developed the algorithm forsymmetric level-index (SLI) arithmetic, and a parallel implementation of it. There has been extensive work on developing the SLI arithmetic algorithms and extending them tocomplex andvector arithmetic operations.
The idea of the level-index system is to represent a non-negativereal numberX as
where, and the process of exponentiation is performedℓ times, with.ℓ andf are thelevel andindex ofX respectively.x =ℓ +f is the LI image ofX. For example,
so its LI image is
The symmetric form is used to allow negative exponents, if the magnitude ofX is less than 1. One takessgn(log(X)) orsgn(|X| − |X|−1) and stores it (after substituting +1 for 0 for the reciprocal sign; since forX = 1 = e0 the LI image isx = 1.0 and uniquely definesX = 1, we can do away without a third state and use only one bit for the two states −1 and +1[clarification needed]) as the reciprocal signrX. Mathematically, this is equivalent to taking thereciprocal (multiplicative inverse) of a small-magnitude number, and then finding the SLI image for the reciprocal. Using one bit for the reciprocal sign enables the representation of extremely small numbers.
Asign bit may also be used to allow negative numbers. One takessgn(X) and stores it (after substituting +1 for 0 for the sign; since forX = 0 the LI image isx = 0.0 and uniquely definesX = 0, we can do away without a third state and use only one bit for the two states −1 and +1[clarification needed]) as the signsX. Mathematically, this is equivalent to taking the inverse (additive inverse) of a negative number, and then finding the SLI image for the inverse. Using one bit for the sign enables the representation of negative numbers.
The mapping function is called thegeneralized logarithm function. It is defined as
and it maps onto itself monotonically, thus being invertible on this interval. The inverse, thegeneralized exponential function, is defined by
The density of valuesX represented byx has no discontinuities as we go from levelℓ toℓ + 1 (a very desirable property) since
The generalized logarithm function is closely related to theiterated logarithm used in computer science analysis of algorithms.
Formally, we can define the SLI representation for an arbitrary realX (not 0 or 1) as
wheresX is the sign (additive inversion or not) ofX, andrX is the reciprocal sign (multiplicative inversion or not) as in the following equations:
whereas forX = 0 or 1, we have
For example,
and its SLI representation is