galois.ReedSolomon.encode(message: ndarray | FieldArray, parity_only: bool = False) → ndarray | FieldArray

Encodes the message \(\mathbf{m}\) into the Reed-Solomon codeword \(\mathbf{c}\).

Parameters
message: ndarray | FieldArray

The message as either a \(k\)-length vector or \((N, k)\) matrix, where \(N\) is the number of messages. For systematic codes, message lengths less than \(k\) may be provided to produce shortened codewords.

parity_only: bool = False

Optionally specify whether to return only the parity symbols. This only applies to systematic codes. The default is False.

Returns

The codeword as either a \(n\)-length vector or \((N, n)\) matrix. The return type matches the message type. If parity_only=True, the parity symbols are returned as either a \(n - k\)-length vector or \((N, n-k)\) matrix.

Notes

The message vector \(\mathbf{m}\) is defined as \(\mathbf{m} = [m_{k-1}, \dots, m_1, m_0] \in \mathrm{GF}(q)^k\), which corresponds to the message polynomial \(m(x) = m_{k-1} x^{k-1} + \dots + m_1 x + m_0\). The codeword vector \(\mathbf{c}\) is defined as \(\mathbf{c} = [c_{n-1}, \dots, c_1, c_0] \in \mathrm{GF}(q)^n\), which corresponds to the codeword polynomial \(c(x) = c_{n-1} x^{n-1} + \dots + c_1 x + c_0\).

The codeword vector is computed from the message vector by \(\mathbf{c} = \mathbf{m}\mathbf{G}\), where \(\mathbf{G}\) is the generator matrix. The equivalent polynomial operation is \(c(x) = m(x)g(x)\). For systematic codes, \(\mathbf{G} = [\mathbf{I}\ |\ \mathbf{P}]\) such that \(\mathbf{c} = [\mathbf{m}\ |\ \mathbf{p}]\). And in polynomial form, \(p(x) = -(m(x) x^{n-k}\ \textrm{mod}\ g(x))\) with \(c(x) = m(x)x^{n-k} + p(x)\). For systematic and non-systematic codes, each codeword is a multiple of the generator polynomial, i.e. \(g(x)\ |\ c(x)\).

For the shortened \(\textrm{RS}(n-s, k-s)\) code (only applicable for systematic codes), pass \(k-s\) symbols into encode() to return the \(n-s\)-symbol codeword.

Examples

Encode a single message using the \(\textrm{RS}(15, 9)\) code.

In [1]: rs = galois.ReedSolomon(15, 9)

In [2]: GF = rs.field

In [3]: m = GF.Random(rs.k); m
Out[3]: GF([ 5,  7,  1, 10,  6, 14,  8,  2, 11], order=2^4)

In [4]: c = rs.encode(m); c
Out[4]: GF([ 5,  7,  1, 10,  6, 14,  8,  2, 11, 10, 15, 11,  1, 14,  1], order=2^4)

Compute the parity symbols only.

In [5]: p = rs.encode(m, parity_only=True); p
Out[5]: GF([10, 15, 11,  1, 14,  1], order=2^4)

Encode a single message using the shortened \(\textrm{RS}(11, 5)\) code.

In [6]: rs = galois.ReedSolomon(15, 9)

In [7]: GF = rs.field

In [8]: m = GF.Random(rs.k - 4); m
Out[8]: GF([2, 9, 4, 5, 9], order=2^4)

In [9]: c = rs.encode(m); c
Out[9]: GF([ 2,  9,  4,  5,  9,  9,  7, 10,  0, 13,  7], order=2^4)

Compute the parity symbols only.

In [10]: p = rs.encode(m, parity_only=True); p
Out[10]: GF([ 9,  7, 10,  0, 13,  7], order=2^4)

Encode a matrix of three messages using the \(\textrm{RS}(15, 9)\) code.

In [11]: rs = galois.ReedSolomon(15, 9)

In [12]: GF = rs.field

In [13]: m = GF.Random((3, rs.k)); m
Out[13]: 
GF([[ 1, 10,  9,  4,  7,  3,  0,  8, 11],
    [ 5,  6, 11,  5, 10,  3,  4,  4,  6],
    [15,  1, 15,  9, 11,  8,  4,  6, 14]], order=2^4)

In [14]: c = rs.encode(m); c
Out[14]: 
GF([[ 1, 10,  9,  4,  7,  3,  0,  8, 11,  5,  8, 13, 15,  9, 12],
    [ 5,  6, 11,  5, 10,  3,  4,  4,  6,  4,  7, 11, 12, 14,  1],
    [15,  1, 15,  9, 11,  8,  4,  6, 14,  8, 10, 12,  9,  9,  8]],
   order=2^4)

Compute the parity symbols only.

In [15]: p = rs.encode(m, parity_only=True); p
Out[15]: 
GF([[ 5,  8, 13, 15,  9, 12],
    [ 4,  7, 11, 12, 14,  1],
    [ 8, 10, 12,  9,  9,  8]], order=2^4)

Encode a matrix of three messages using the shortened \(\textrm{RS}(11, 5)\) code.

In [16]: rs = galois.ReedSolomon(15, 9)

In [17]: GF = rs.field

In [18]: m = GF.Random((3, rs.k - 4)); m
Out[18]: 
GF([[15, 12, 10, 15,  6],
    [12,  1,  1, 11,  4],
    [ 6,  1,  9,  9,  4]], order=2^4)

In [19]: c = rs.encode(m); c
Out[19]: 
GF([[15, 12, 10, 15,  6,  7, 14,  4,  1,  7, 15],
    [12,  1,  1, 11,  4, 12, 15,  7, 14,  4,  5],
    [ 6,  1,  9,  9,  4, 15,  1, 10,  9, 12,  8]], order=2^4)

Compute the parity symbols only.

In [20]: p = rs.encode(m, parity_only=True); p
Out[20]: 
GF([[ 7, 14,  4,  1,  7, 15],
    [12, 15,  7, 14,  4,  5],
    [15,  1, 10,  9, 12,  8]], order=2^4)
Last update: Jul 28, 2022