The present work relates generally to linear block codes and, more particularly, to single error correcting, double error detecting (SECDED) codes over GF(q).
SECDED codes play a pivotal role in data transfer/communication applications. One example is data storage applications. Advances in memory technology have provided the capability of storing increasingly large amounts of data in individual memory packages (or modules). This makes it possible to store/retrieve from a single memory module data that previously was stored/retrieved by accessing multiple modules, providing improvements in performance, power consumption, and scalability. When accessing a single module, it is advantageous to address effectively errors that have a single cause. For example, if several I/O terminals of a module provide respective I/O paths, a failed transmitter driving one of the terminals may be used multiple times per memory access. The bits transferred via any given terminal in a single burst may be considered to be a data symbol. If the terminal transfers data in 4-bit bursts, for example, then the symbol size is four bits. A sixteen-byte transfer might be accomplished by transferring a 4-bit burst (i.e., one symbol) on each of 32 terminals of the memory module. Within knowable limits, a SECDED code that appends four parity check symbols of four bits each is capable of correcting any error in a single data symbol (and detecting errors in two data symbols). Thus, errors repeatedly caused by a single failed transmitter associated with a single terminal may be corrected.
The aforementioned example SECDED code has a Hamming distance of at least four, in terms of symbols. This is a specific example of a linear [n, k, d]q code over a general alphabet GF(q), with n=36 (data+parity check) symbols, k=32 data symbols, code distance d=4, and symbol alphabet size q=16 (24 with 4 bits/symbol), that is, a [36, 32, 4]16 code.
One approach to construct an [n, k, d]q code is a random search of parity check matrices in GF(q). Considering the aforementioned example of a [36, 32, 4]16 code, there are 15×36 different single symbol errors, while the syndrome from the four symbol parity has 65536 possibilities. It thus appears that a random search would likely yield a code. Nevertheless, testing as many as 15×230 H matrices obtained via random search has failed to yield a [36, 32, 4]16 code.
Another approach is to construct a [35, 32, 3]16 code, and extend it with a single parity check to a [36, 32, 4]16 code. Since GCD(35,16)=1, construction of a BCH code with n=35 and d=3 may be considered. BCH codes provide a minimum distance guarantee that is determined by the generator polynomial. However, investigation of BCH codes with n=35 and q=16 reveals that the most efficient BCH codes with d=3 are [35, 31, 3]16 codes, which cannot be extended to a [36, 32, 4]16 code.
Another approach is the use of an existing bound, such as the Gilbert-Varshamov (GV) bound, that indicates feasible codes. However, for the code distance d=4, the GV bound indicates that only k=27.5 is achievable.
Another approach is constructing a [36, 32, 4]16 code as four simultaneous [36, 32, 4]2 codes, or two simultaneous [36, 32, 4]4 codes. However, a [36, 32, 4]2 code does not exist because it violates the Hamming bound. Although a linear [36, 32, 4]4 code would not violate the Hamming bound, it has been found that [36, 32, 4]4 codes do not exist.
It is desirable in view of the foregoing to provide for the construction of distance d linear codes more efficiently than by random searching.
The present work provides for non-deterministically generating the parity check matrix for a linear [n, k, d]q code over a general alphabet GF(q). The search space is organized as a tree, and the search can terminate in one of two states, one if it has taken an unsuccessful path and run out of search space, another if it has found a suitable parity check matrix. Also provided is a simple decoder for linear SECDED codes, with an efficient parallel implementation. The present work further provides for generating, from an input parity check matrix, another parity check matrix whose corresponding code has the same distance as that of the input matrix, is a subcode of the input matrix, and requires the minimum number of logical AND gates among all subcodes that encode the same set of symbols.
In the retrieval path 19, the parity check symbols at 17 are exclusive-ORed (XORed) at 12 with the parity check symbols that were stored (and are now retrieved) with the data symbols. The result of the XOR operation at 12 is the syndrome, which is decoded by a syndrome decoder 13. The syndrome decoder 13, which uses information from entries in the matrix H, produces information 15 which may be used by an error corrector, together with the retrieved data symbols, to correct single symbol errors and output correspondingly corrected data. The syndrome decoder 13 also produces information 16 which indicates detection of errors in two (and often more) symbols.
The above-described operations performed by the check symbol generator 11, the syndrome decoder 13 and the error corrector 14 are generally known in the art. Thus, in some embodiments, these components operate in accordance with conventional techniques. However, with respect to the parity check matrix H used by the check symbol generator 11 and the syndrome decoder 13, the present work provides improved techniques for non-deterministically generating this matrix H for a linear SECDED code. The present work exploits the known relationship between the code distance of a linear code and the requirement that there must be linear independence of the column vectors of the parity check matrix for the code. More specifically, if H is the parity check matrix of a linear code C, then the code C has distance d if and only if any non-trivial linear combination of d−1 column vectors of H is not zero. (Column vectors may also be referred to herein simply as columns, or vectors.) A parity check matrix H is described in standard form as
H=[A
r×k
I
r×r],
where A is a matrix describing the parity checks, I is the identity matrix concatenated with A to form H, k is the number of data symbols transferred, and r is the number of parity check symbols generated for the k data symbols. The aforementioned property that no non-trivial linear combination of d−1 column vectors of H results in zero is referred to herein as “d−1 linear independence”, and the columns of H are said to be “d−1 linearly independent.”
Considering the aforementioned example using 4-bit symbols, this provides sixteen (24) available values for a symbol. That is, the size of the symbol alphabet, referred to herein as q, is q=16. As shown generally in
The iterative filter 31 is shown in more detail in
where al represents any of the q symbols in the symbol alphabet, and yl represents any column of I. If the selected vector is not of the form defined by the test at 53 (i.e., satisfies d−1 linear independence), then the vector passes the test, and operations proceed to 55. If the selected vector is of the form defined by the test at 53 (i.e., does not satisfy d−1 linear independence), then the vector fails the test, and is eliminated from P at 54, after which operations proceed to 55. As indicated at 55, the above-described operations are repeated until all vectors of P have been tested at 53, and eliminated at 54 as appropriate. After all vectors of P have been tested, the set P0 is formed at 56 as the set of vectors that have not been eliminated from P.
Noting that a column vector is d−1 linearly dependent with columns of I if and only if that vector has at least d−2 zeros, some embodiments form the set P0 by simply eliminating from P all vectors that contain at least d−2 zeros. Consider, for example, the aforementioned scenario where each element of each column vector contains four bits. Each column vector wherein all four bits of at least d−2 elements are zero is eliminated from P, and the result is P0.
At 63, if the iteration index i=k−1, then all columns of A have been populated, and the matrix is complete. Otherwise, another vector is selected from Pi at 64. The test shown at 65 is applied to the selected vector. The test 65 is whether the selected vector is of the form
where al represents any of the q symbols in the symbol alphabet, al represents any of the q symbols in the symbol alphabet, and yl represents any column of I or Ai. If the selected vector is not of the form defined by the test at 65 (i.e., satisfies d−1 linear independence), then the vector passes the test, and operations proceed to 67. If the selected vector is of the form defined by the test at 65 (i.e., does not satisfy d−1 linear independence), then the vector fails the test, and is eliminated from Pi at 66, after which operations proceed to 67. As indicated at 67, the operations described above at 64-66 are repeated until all vectors of Pi have been tested at 65, and eliminated as appropriate at 66. After all vectors of Pi have been tested, the set Pi+1 is formed at 68 as the set of vectors that have not been eliminated from Pi.
Thereafter, a test is performed at 69 to determine whether Pi+1 contains enough vectors to populate the currently unpopulated columns of A. In the test 69, |P1+i| denotes the number of vectors in Pi+1, and q is the size of the symbol alphabet. For each vector x that is yet to be concatenated to A at 62, at least the q−1 non-zero multiples of that vector will be eliminated at 66, so the maximum number of vectors still available for concatenation to A is bounded by |Pi+1|/q−1. There are at this point i+1 columns already populated in A, with k−i−1 columns remaining to be populated. The test 69 determines whether the maximum number of available vectors is less than the number of unpopulated columns of A. If not, then the iteration index is updated at 600, and operations return to 61 to begin the next iteration. Otherwise, the test 69 determines that there are not enough vectors in Pi+1 to complete the construction of A, so the current attempt to generate A fails.
The operations shown in
In contrast to the above-described techniques of the present work, conventional approaches use a random code search when attempting to find an [n, k, d]q code that is not a member of the currently known block codes or their modifications. The smaller the ratio of the total number of check matrices for [n, k, d]q codes to the total number of k×r matrices, the less likely is the success of a random search technique. The present work provides a more systematic search for check matrices and is significantly more efficient than a random search in finding [n, k, d]q codes if these codes are rare. Although the time needed to verify the d−1 linear independence of a column increases with successive iterations in the present work, the number of acceptable columns decreases with successive iterations. This makes it progressively more likely to choose columns that are acceptable as the time needed to verify column acceptability progressively increases. Also, because candidate columns may be identified as not acceptable while the matrix is being built, the present work determines relatively early if a matrix under construction will not work.
As an example, a random brute force search, using fifteen 2.27 GHz, 8 GByte processors for two weeks, has failed to produce a linear [36, 32, 4]16 code. (Relaxing the check symbol requirement to r=5, or the distance requirement to d=3, has been found to enable random brute force searching to produce the respectively associated [37, 32, 4]16 and [36, 32, 3]16 codes in about an hour.) In contrast, the techniques of the present work have produced a linear [36, 32, 4]16 code in under a minute using a single processor. A linear [68, 64, 4]16 code has also been produced according to the present work.
In some embodiments, the syndrome decoder 13 in
For a SECDED code, only errors that have Hamming weights 0, 1, and 2 need to be distinguished. This forms the essence of the decoding technique. The following relationships among s, e, and h(i) may be written:
h(i)|siff∃e∈GF(q):s=eh(i).
If h(i)|s and s≠0, then e=s|h(i) is well defined. Accordingly, as shown at 71-73 in
Some embodiments achieve improved computational efficiency by assuming that H is in standard form, H=[Ar×kIr×r], and that errors in the parity check symbols will not be corrected. As shown at 81-83 in
Given s,h∈GF(q)r, evaluating whether the condition h(i)∈s is satisfied is equivalent to satisfying the following conditions, where i and j are symbol indices of the syndrome s and the column h(i)
hj=0sj=0i∀j
∀i,j:hi,hj≠0 require (hi)−1si=(hj)−1sj
For purposes of evaluating the foregoing conditions, the multiplicative inverses of the possible non-zero symbol values in H (e.g., (hi)−1 and (hj)−1) may be pre-computed. For example, in the aforementioned case of 4-bit symbols, each symbol value of each column may be represented as a single hexadecimal digit. The respective multiplicative inverses for the hexadecimal digits 1, 2, . . . , E, F are: 1 9 E D B 7 6 F 2 C 5 A 4 3 8.
Applying the conditions defined above to the symbols of an example column vector h=[4 F 0 3]T in a four-row (r=4) H matrix, it can be seen that the h(i)|s condition is satisfied for a given syndrome s if
s3=0
4−1s1=Ds1=F−1s2=8s2=3−1s4=Es4
where s1, s2, s3 and s4 are the syndrome symbols.
In some embodiments, pre-computed multiplicative inverses are stored, and the syndrome decoder 13 of
For a systematic [k+r,k,d]q linear code with parity check matrix in standard form H=[|I], and input x=[x1 . . . xk], the check symbol generator 11 computes the parity check symbols as follows
As is well known in the art, the addition operations in GF(2m) computation may be implemented using XOR logic gates, and the multiplication operations may be implemented using XOR and AND logic gates. The weight of an element a∈GF(2m), denoted as w(a), is defined as the sum of the Hamming weights of vector space representations a,a−1. The weight of a matrix is the sum of the weights of its elements. The weight of H is directly proportional to the total number of AND gates required for both encoding and decoding. Example embodiments of the present work minimize the weight of H, thereby minimizing the total number of AND gates in the check symbol generator 11 and the syndrome decoder 13.
Some embodiments produce a check matrix H(k
As shown in
For each column vector h(i),i=1, . . . , kmax of H,
In some embodiments, the matrix H used in the processing of
Although example embodiments of the present work are described above in detail, this does not limit the scope of the present work, which can be practiced in a variety of embodiments.
This application is a divisional application of U.S. Non-Provisional Application No. 14/029,379, Design, Decoding and Optimized Implementation of SECDED Codes over GF(q), by H. Lee Ward et al., filed Sep. 17, 2013 (allowed Feb. 6, 2014), which is a divisional application of U.S. Non-Provisional Application No. 13/116,976, Design, Decoding and Optimized Implementation of SECDED Codes over GF(q), by H. Lee Ward et al., filed May 26, 2011, now U.S. Pat. No. 8,566,684, issued Oct. 22, 2013, both of which are incorporated herein by reference.
This invention was developed under Contract DE-AC04-94AL85000 between Sandia Corporation and the U.S. Department of Energy. The U.S. Government has certain rights in this invention.
Number | Date | Country | |
---|---|---|---|
Parent | 14029379 | Sep 2013 | US |
Child | 14281567 | US | |
Parent | 13116976 | May 2011 | US |
Child | 14029379 | US |