Chien search explained

In abstract algebra, the Chien search, named after Robert Tienwen Chien, is a fast algorithm for determining roots of polynomials defined over a finite field. Chien search is commonly used to find the roots of error-locator polynomials encountered in decoding Reed-Solomon codes and BCH codes.

Algorithm

The problem is to find the roots of the polynomial (over the finite field): \Lambda(x) = \lambda_0 + \lambda_1 x + \lambda_2 x^2 + \cdots + \lambda_t x^t

The roots may be found using brute force: there are a finite number of, so the polynomial can be evaluated for each element . If the polynomial evaluates to zero, then that element is a root.

For the trivial case, only the coefficient need be tested for zero. Below, the only concern will be for non-zero .

A straightforward evaluation of the polynomial involves general multiplications and additions. A more efficient scheme would use Horner's method for general multiplications and additions. Both of these approaches may evaluate the elements of the finite field in any order.

Chien search improves upon the above by selecting a specific order for the non-zero elements. In particular, the finite field has a (constant) generator element . Chien tests the elements in the generator's order . Consequently, Chien search needs only multiplications by constants and additions. The multiplications by constants are less complex than general multiplications.

The Chien search is based on two observations:

\beta

may be expressed as
i\beta
\alpha
for some

i\beta

, where

\alpha

is a primitive element of

GF(q)

,

i\beta

is the power number of primitive element

\alpha

. Thus the powers

\alphai

for

0\leqi<(q-1)

cover the entire field (excluding the zero element).

\begin \Lambda(\alpha^i) &=& \lambda_0 &+& \lambda_1 (\alpha^i) &+& \lambda_2 (\alpha^i)^2 &+& \cdots &+& \lambda_t (\alpha^i)^t \\ &\triangleq& \gamma_ &+& \gamma_ &+& \gamma_ &+& \cdots &+& \gamma_ \\ \Lambda(\alpha^) &=& \lambda_0 &+& \lambda_1 (\alpha^) &+& \lambda_2 (\alpha^)^2 &+& \cdots &+& \lambda_t (\alpha^)^t \\ &=& \lambda_0 &+& \lambda_1 (\alpha^i)\,\alpha &+& \lambda_2 (\alpha^i)^2\,\alpha^2 &+& \cdots &+& \lambda_t (\alpha^i)^t\,\alpha^t \\ &=& \gamma_ &+& \gamma_\,\alpha &+& \gamma_\,\alpha^2 &+& \cdots &+& \gamma_\,\alpha^t \\ &\triangleq& \gamma_ &+& \gamma_ &+& \gamma_ &+& \cdots &+& \gamma_\end

In other words, we may define each

Λ(\alphai)

as the sum of a set of terms

\{\gammaj,i\mid0\leqj\leqt\}

, from which the next set of coefficients may be derived thus: \gamma_ = \gamma_\,\alpha^j

In this way, we may start at

i=0

with

\gammaj,0=λj

, and iterate through each value of

i

up to

(q-1)

. If at any stage the resultant summation is zero, i.e. \sum_^t \gamma_ = 0,then

Λ(\alphai)=0

also, so

\alphai

is a root. In this way, we check every element in the field.

When implemented in hardware, this approach significantly reduces the complexity, as all multiplications consist of one variable and one constant, rather than two variables as in the brute-force approach.