for (i = 0; i < m; i++)

    {

           alpha_to[i] = mask;

           index_of[alpha_to[i]] = i;

           if (p[i] != 0) alpha_to[m] ^= mask;

           mask <<= 1;

    } index_of[alpha_to[m]] = m; mask >>= 1;

    for (i = m+1; i < n; i++)

    {

           if (alpha_to[i-1] >= mask)

                 alpha_to[i] = alpha_to[m] ^ ((alpha_to[i-1]^mask)<<1);

           else

                 alpha_to[i] = alpha_to[i-1]<<1;

           index_of[alpha_to[i]] = i;

    } index_of[0] = -1;

}

Сама же функция умножения выглядит тривиально, укладываясь всего в пяток строк. В большинстве программных реализаций кодера/декодера Рида-Соломона, которые мне только доводилось видеть, операция умножения даже не выносится в отдельную процедуру, а реализуется непосредственно по месту вызова.

Листинг 4. Функция быстрого табличного умножения полиномов в полях Галуа