Zech's logarithm
Template:Short description Template:Use dmy dates Zech logarithms are used to implement addition in finite fields when elements are represented as powers of a generator <math>\alpha</math>.
Zech logarithms are named after Julius Zech,<ref name="Zech_1849"/><ref name="Zech_1863"/><ref name="Zech_1892"/><ref name="Zech_1910"/> and are also called Jacobi logarithms,<ref name="Lidl_1997"/> after Carl G. J. Jacobi who used them for number theoretic investigations.<ref name="Jacoby_1846"/>
DefinitionEdit
Given a primitive element <math>\alpha</math> of a finite field, the Zech logarithm relative to the base <math>\alpha</math> is defined by the equation
- <math>\alpha^{Z_\alpha(n)} = 1 + \alpha^n,</math>
which is often rewritten as
- <math>Z_\alpha(n) = \log_\alpha(1 + \alpha^n).</math>
The choice of base <math>\alpha</math> is usually dropped from the notation when it is clear from the context.
To be more precise, <math>Z_\alpha</math> is a function on the integers modulo the multiplicative order of <math>\alpha</math>, and takes values in the same set. In order to describe every element, it is convenient to formally add a new symbol <math>-\infty</math>, along with the definitions
- <math>\alpha^{-\infty} = 0</math>
- <math>n + (-\infty) = -\infty</math>
- <math>Z_\alpha(-\infty) = 0</math>
- <math>Z_\alpha(e) = -\infty</math>
where <math>e</math> is an integer satisfying <math>\alpha^e = -1</math>, that is <math>e=0</math> for a field of characteristic 2, and <math>e = \frac{q-1}{2}</math> for a field of odd characteristic with <math>q</math> elements.
Using the Zech logarithm, finite field arithmetic can be done in the exponential representation:
- <math>\alpha^m + \alpha^n = \alpha^m \cdot (1 + \alpha^{n-m}) = \alpha^m \cdot \alpha^{Z(n-m)} = \alpha^{m + Z(n-m)} </math>
- <math>-\alpha^n = (-1) \cdot \alpha^n = \alpha^e \cdot \alpha^n = \alpha^{e+n}</math>
- <math>\alpha^m - \alpha^n = \alpha^m + (-\alpha^n) = \alpha^{m + Z(e+n-m)} </math>
- <math>\alpha^m \cdot \alpha^n = \alpha^{m+n}</math>
- <math>\left( \alpha^m \right)^{-1} = \alpha^{-m}</math>
- <math>\alpha^m / \alpha^n = \alpha^m \cdot \left( \alpha^n \right)^{-1} = \alpha^{m - n}</math>
These formulas remain true with our conventions with the symbol <math>-\infty</math>, with the caveat that subtraction of <math>-\infty</math> is undefined. In particular, the addition and subtraction formulas need to treat <math>m = -\infty</math> as a special case.
This can be extended to arithmetic of the projective line by introducing another symbol <math>+\infty</math> satisfying <math>\alpha^{+\infty} = \infty</math> and other rules as appropriate.
For fields of characteristic 2,
- <math>Z_\alpha(n) = m \iff Z_\alpha(m) = n</math>.
UsesEdit
For sufficiently small finite fields, a table of Zech logarithms allows an especially efficient implementation of all finite field arithmetic in terms of a small number of integer addition/subtractions and table look-ups.
The utility of this method diminishes for large fields where one cannot efficiently store the table. This method is also inefficient when doing very few operations in the finite field, because one spends more time computing the table than one does in actual calculation.
ExamplesEdit
Let Template:Nowrap be a root of the primitive polynomial Template:Nowrap. The traditional representation of elements of this field is as polynomials in α of degree 2 or less.
A table of Zech logarithms for this field are Template:Nowrap, Template:Nowrap, Template:Nowrap, Template:Nowrap, Template:Nowrap, Template:Nowrap, Template:Nowrap, and Template:Nowrap. The multiplicative order of α is 7, so the exponential representation works with integers modulo 7.
Since α is a root of Template:Nowrap then that means Template:Nowrap, or if we recall that since all coefficients are in GF(2), subtraction is the same as addition, we obtain Template:Nowrap.
The conversion from exponential to polynomial representations is given by
- <math>\alpha^3 = \alpha^2 + 1</math> (as shown above)
- <math>\alpha^4 = \alpha^3 \alpha = (\alpha^2 + 1)\alpha = \alpha^3 + \alpha = \alpha^2 + \alpha + 1</math>
- <math>\alpha^5 = \alpha^4 \alpha = (\alpha^2 + \alpha + 1)\alpha = \alpha^3 + \alpha^2 + \alpha = \alpha^2 + 1 + \alpha^2 + \alpha = \alpha + 1</math>
- <math>\alpha^6 = \alpha^5 \alpha = (\alpha + 1)\alpha = \alpha^2 + \alpha</math>
Using Zech logarithms to compute αTemplate:Hairsp6 + αTemplate:Hairsp3:
- <math>\alpha^6 + \alpha^3 = \alpha^{6 + Z(-3)} = \alpha^{6 + Z(4)} = \alpha^{6 + 6} = \alpha^{12} = \alpha^5,</math>
or, more efficiently,
- <math>\alpha^6 + \alpha^3 = \alpha^{3 + Z(3)} = \alpha^{3 + 2} = \alpha^5,</math>
and verifying it in the polynomial representation:
- <math>\alpha^6 + \alpha^3 = (\alpha^2 + \alpha) + (\alpha^2 + 1) = \alpha + 1 = \alpha^5.</math>
See alsoEdit
- Gaussian logarithm
- Irish logarithm, a similar technique derived empirically by Percy Ludgate
- Finite field arithmetic
- Logarithm table
ReferencesEdit
Further readingEdit
- Template:Cite book
- Template:Cite journal
- Template:Cite journal [1] [2] [3]
- {{#invoke:citation/CS1|citation
|CitationClass=web }}