Field of the Invention
This invention relates generally to data processing systems and, more particularly, to a system for decoding and correcting errors in data using an error correction code.
Background of the Invention
Data stored on magnetic media, such as a magnetic disks, are typically stored in encoded form, so that errors in the stored data can possibly be corrected. The errors may occur, for example, because of inter-symbol interference, a defect in the disk, or noise. As the density of the data stored on the disk increases, more errors are likely, and the system is required to correct greater numbers of errors, which include greater numbers of burst errors. A burst error is typically defined as a contiguous number of symbols in which the first symbol and the last symbol are erroneous. The speed with which the system corrects the errors, including the burst errors, is important to the overall speed with which the system processes the data.
Prior to recording, multiple-bit data symbols are encoded using an error correction code (ECC). When the data symbols are retrieved from the disk and demodulated, the ECC is employed to, as the name implies, correct the erroneous data.
Specifically, before a string of k data symbols is written to a disk, it is mathematically encoded using an (n, k) ECC to form n-k ECC symbols. The ECC symbols are then appended to the data string to form an n-symbol error correction code word, which is then written to, or stored, on the disk. When the data are read from the disk, the code words containing the data symbols and ECC symbols are retrieved and mathematically decoded. During decoding, errors in the data are detected and, if possible, corrected through manipulation of the ECC symbols for a detailed description of decoding see, Peterson and Weldon, Error Correction Codes, 2nd Ed. MIT Press, 1972!.
To correct multiple errors in strings of data symbols, the system typically uses an ECC that efficiently and effectively utilizes the various mathematical properties of sets of symbols known as Galois fields. Galois fields are represented "GF (PM)" , where "P"is a prime number and "M" can be thought of as the number of digits, base "P", in each element or symbol in the field. P usually has the value 2 in digital computer and disk drive applications and, therefore, M is the number of bits in each symbol. The ECC's commonly used with the Galois Fields are Reed Solomon codes or BCH codes.
There are essentially four major steps in decoding a corrupted code word of a high rate Reed-Solomon code or a BCH code. The system first determines error syndromes that are based on the results of a manipulation of the ECC symbols. Next, using the error syndromes, the system determines an error locator polynomial, which is a polynomial that has the same degree as the number of errors. The system then finds the roots of the error locator polynomial and from each root determines the location of an associated error in the code word. Finally, the system finds error values for the error locations.
The steps of determining the syndromes and finding the error locations are the most time consuming in the error correction process. This invention relates to the step of finding the error locations.
Codes are generally interleaved to provide maximum burst error protection. Otherwise, there are too many errors for the (non-interleaved) code to handle or it takes too long for the error correcting system to determine the error locations using the conventional trial and error method.
Interleaving segments the data string such that adjacent symbols are encoded by the ECC into different code words. Thus, the burst errors become independent errors spread out between the interleaved code words, and can be readily and quickly corrected.
Interleaving the codes, however, requires separately encoding each code segment and separately decoding each segment. Further, interleaving requires , for each interleaved code, extra ECC symbols for error detection. Alternatively , error detection could be performed using an error detection code that overlays the ECC , such as a parity check to guard against incorrect correction of the data. The interleaved codes must, however, correct the error detection symbols in order to ensure that the error detection code can perform its error detection function. The interleaved codes are thus either inefficient for requiring, for each interleaved code, ECC symbols for error detection or time consuming for requiring additional operations that are associated with an overlaid error detection code.
Another problem with interleaved codes is that the error correction fails if there are too many errors in one of the data segments, even if the other data segments are essentially error-free.
What is needed is a system that reduces the time it takes to determine the locations and lengths of burst errors using a non-interleaved code, such that the errors can be corrected in real time.
Summary of the Invention
The invention is a multiple-solid-burst error correcting system for determining the number and locations of "solid burst" errors in a high-rate Reed Solomon or BCH code. A solid burst is one that includes a succession of erroneous symbols. This differs from the conventional definition of a burst error, which may include one or more non-erroneous symbols in between first and last erroneous symbols. Bursts that include non-erroneous symbols are treated as multiple solid bursts by the multiple-solid-burst error correcting system. An erroneous symbol that has only adjacent error-free symbols is a solid burst that consists of one symbol.
The system determines the greatest common divisor of the error locator polynomial .sigma.(x), which has roots .alpha..sup.-i.sbsp.k that correspond to error locations i.sub.k, and a mapping error locator polynomial .sigma.(.alpha.*x) that maps the error locations i.sub.k, to locations i.sub.k+ 1. The roots that are common to both the error locator polynomial and the mapping error locator polynomial, that is, the roots that are included in the greatest common divisor, d(x), correspond to adjacent error locations that are contained in the solid bursts. The roots of the non-common factors, p.sub.1 (x), of the error locator polynomial correspond to the first locations of the respective solid bursts and the roots of the non-common factors p.sub.2 (x) of the mapping error locator polynomial correspond to one location beyond the end locations of the solid bursts. The system determines the roots of only the non-common factors to determine the locations that are included in the solid bursts. The system then labels all locations between associated first and last locations as part of a solid burst.
If the error locator polynomial has degree "e" and the greatest common divisor has degree "d", the polynomials p.sub.1 (X) and p.sub.2 (x) each have degree e-d=p. The system finds the roots of the lower-degree p.sub.1 (x), and uses p.sub.2 (x), to test for the end locations of the bursts. Accordingly, the system finds the roots of one degree-p polynomial rather than the roots of a degree-e polynomial, as is required by prior systems.
More specifically, the system determines the greatest common divisor of the error locator polynomial
where "* " represents multiplication over GF(2.sup.M), and a mapping error locator polynomial
where the addition of exponents is modulo the code length n. The factors that are common to both polynomials are part of the greatest common divisor, and correspond to adjacent error locations that are part of the solid burst errors, as discussed above.
The factors that are included only in eqn. 1 correspond to error locations i.sub.k that are, respectively, the first error locations in the solid bursts. Further, the roots .alpha..sup.-(i.sbsp.j.sup.+1) of the actors that are included only in eqn. 2 correspond respectively to error locations that are one location beyond the last locations of the solid bursts. These i.sub.j+ 1 locations do not correspond to the roots of .sigma.(x), and thus, do not contain erroneous symbols.
The system uses a fast, iterative mechanism for determining the greatest common divisor, d(x), of eqns. 1 and 2, as discussed below. Once the system determines the divisor, the system next determines the locations of the first and last symbols of the respective solid bursts by finding the roots of polynomial p.sub.1 (x) where: ##EQU1## and testing to determine if adjacent locations are roots of p.sub.2 (x), where: ##EQU2## The two polynomials p.sub.1 (x) and p.sub.2 (x), respectively, contain the non-common factors of .sigma.(x) and .sigma.(.alpha.*x), and both are degree e-d=p, where d is the degree of the greatest common divisor, which is equal to the number of common factors.
The degree p of p.sub.1 (x) and p.sub.2 (x) is the number of solid bursts. There is generally a relatively small number of solid bursts, and thus, the degree of p.sub.1 (x) and p.sub.2 (x) is relatively low. Accordingly, known fast techniques for finding the roots of the lower-degree polynomial p.sub.1 (x) can be used to determine the first locations of the solid bursts and adjacent error locations are tested as roots of p.sub.2 (x). Thus, the system need not use the more time consuming methods that would be required if it had, also, to find the d common roots.
Brief Description of the Drawings
For a fuller understanding of the nature and objects of the invention, reference should be had to the following detailed description taken in connection with the accompanying drawings, in which:
FIG. 1 is a functional block diagram of a decoder in the multiple-solid-burst error correcting system;
FIG. 2 is a flow chart of the operations of a GCD processor included in the system of FIG. 1;
FIG. 3 is a chart that illustrates solid bursts and a greatest common divisor;
FIG. 4 is a more detailed flow chart of the operations of an error locator processor of FIG. 1 in determining the locations contained within respective solid bursts; and
FIGS. 5A and 5B a are more detailed flow chart of the operations of the GCD processor in determining a greatest common divisor.
Description of the Preferred Embodiment
The mathematical operations of addition, subtraction, multiplication and division discussed herein are Galois Field operations over the applicable Galois Field.
Referring now to FIG. 1, a decoder 10 includes a syndrome generator 12 that operates in a conventional manner to produce a set of error syndromes for a received code word. If the syndromes are all zero, the syndrome generator 12 determines that the code word is error-free. Otherwise, the syndrome generator sends the set of syndromes to an error locator polynomial generator 14, which produces in a conventional manner from the syndromes a degree "e" error locator polynomial, .sigma.(x):
where e is the number of errors in the code word.
The roots, .sigma..sup.-i.sbsp.k, of the error locator polynomial, and thus each factor of .sigma.(x), correspond to error locations i.sub.k in the code word.
Referring also to FIG. 2, the error locator polynomial generator 14 determines if e.ltoreq..beta., where .beta. is a predetermined number of errors that can be corrected using a direct solution, that is, can be corrected in hardware {step 100}. If so, the error locator polynomial generator 14 sends the error locator polynomial to an error locator processor 16, which selects and implements an appropriate direct solution process for determining the e roots of the error locator polynomial {step 102}. For example, if e =3, the error locator processor 16 implements a process for determining the roots of a degree-three polynomial.
If e>.beta., the error locator polynomial generator 14 sends the error locator polynomial to a GCD processor 18, which produces a mapping error locator polynomial, .sigma.(.alpha.*x), that maps the roots, .alpha..sup.-i .sbsp.k of the error locator polynomial .sigma.(x) to .alpha..sup.-(i.sbsp.k.sup.+1), where the addition of the exponents is modulo the code length n {step 104}. The mapping error locator polynomial .sigma.(.alpha.*x) thus has roots that correspond to code word locations that are shifted one symbol from, or are adjacent to, the error locations determined by .sigma.(x), as depicted in FIG. 3.
The GCD processor 18 next determines the degree "d" greatest common divisor, d(x), of .sigma.(x) and .sigma.(.alpha.*x).{step 106}. The d common roots indicate error locations that are part of the solid bursts, since they correspond to adjacent error locations. The e-d=p roots, .alpha..sup.-i .sbsp.k, of the error locator polynomial .sigma.(x) that are not common with those of the mapping error locator polynomial .sigma.(.alpha.*x) correspond, respectively, to the first locations of the solid bursts. The p roots .alpha..sup.-(i.sbsp.j.sup.+1) of the mapping error locator polynomial .sigma.(.alpha.*x) that are not common with those of the error locator polynomial .sigma.(x) correspond, respectively, to locations that are one location beyond the last locations of the solid bursts. Thus, the system determines the first and last locations of each solid burst from the non-common roots, and labels as erroneous all of the symbols that are in between the associated first and last symbols.
To find the error locations, the system determines the first locations of the respective solid bursts by finding the roots .alpha..sup.-i .sbsp.k of the degree-p polynomial p.sub.1 (x) which consists of the non-common factors of .sigma.(x). Next, the system, for each root of p.sub.1 (x), tests adjacent locations .alpha..sup.-(i.sbsp.k.sup.+1), .alpha..sup.-(i.sbsp.k.sup.+2) . . . as roots of p.sub.2 (x) to find the end location of the associated solid burst. The system can do this much faster than it can find the roots of a higher degree polynomial by a similar trial and error technique. Further, the system need not find the roots of the d common factors, since these roots correspond to the locations between the various associated first and last locations.
The GCD processor 18 produces the mapping error locator polynomial .sigma.(.alpha.*x) from the error locator polynomial:
by multiplying the coefficients .sigma..sub.j by .alpha..sup.j.
The resulting polynomial is the mapping error locator polynomial, which maps .alpha..sup.-i.sbsp.k to .alpha..sup.-(i.sbsp.k.sup.+1).
The GCD processor 18 uses a novel, iterative operation based on the Euclidean Algorithm to determine a greatest common divisor, d(x), as discussed below with reference to FIG. 5. The iterative operation determines not only d(x), but also the two polynomials p.sub.1 (x) and p.sub.2 (x) where ##EQU3##
Once the GCD processor 18 determines d(x), p.sub.1 (x) and p.sub.2 (x), the processor determines if the degree, p, of p.sub.1 (X), which corresponds to the number of solid bursts, is less than or equal to .beta., to determine if the roots of p.sub.1 (x) can be found using a direct solution {step 108}. If so, the GCD processor sends the polynomial p.sub.1 (x) to the error location processor 16, which directly and quickly determines the p roots using known techniques for finding the roots of a degree-p polynomial {step 110}.
If p>.beta., it means that the system cannot use a direct solution technique to find the error locations and the system must use other, necessarily slower, techniques to find the error locations {step 112}.
As discussed above, every location in a solid burst contains an erroneous symbol. An error location that has only adjacent error-free symbols is a solid burst error of length 1. The degree, d, of the greatest common divisor, d(x), is thus equal to
Where b.sub.m is the length of the m.sup.th solid burst error. Accordingly, if there are one or more long burst errors, the degree of p.sub.1 (x), which is e-d=p, is relatively small. Recall that e is he degree of .sigma.(x), or the total number of erroneous symbols in the codeword.
Once the roots .alpha..sup.-i.sbsp.k of the degree-p polynomial p.sub.1 (x) have been determined, the error locator processor 16 determines which error locations are contained in the p solid bursts. The processor labels the locations that correspond to the roots .alpha..sup.-i.sbsp.k of p.sub.1 (x), respectively, as the first locations of the solid bursts and for each burst tests the consecutive roots .alpha..sup.-(i.sbsp.k+1), .alpha..sup.-(i.sbsp.k+2), . . . until a root .alpha..sup.31 (i.sbsp.j.sup.+1) of p.sub.2 (x) is found. The system then labels the location i.sub.j as the last location of the solid burst and the locations in between i.sub.k and i.sub.j as part of the same solid burst. The system repeats this for all the roots of p.sub.1 (x). {step 112}. In contrast, prior error correction systems must determine all e roots of the error locator polynomial in order to determine the error locations.
Referring now to FIG. 4, the error locator processor 16 selects a first root .alpha..sup.-i.sbsp.kt of p.sub.1 (x), where t=1 is the number of the solid burst. The processor sets "b," the burst error length, to b=1 and determines if .alpha..sup.-(i kt.sup.+b is a root of p.sub.2 (x) {step 204}. If so, the error locator processor determines that the first burst error has a length of b=1, and thus, an end location that corresponds to .alpha..sup.-(i.sbsp.kt.sup.+b-1) or to .alpha..sup.-i.sbsp.kt {step 206}. Otherwise, the processor sets b=b+1 and again determines if .alpha..sup.-(i.sbsp.kt .sup.+b) is a root of p.sub.2 (x) {step 208}. The processor repeats these operations until it finds a root of p.sub.2 (x). The processor then determines that the location, i.sub.k.sbsb.t+b- 1 , which corresponds to .alpha..sup.-(i.sbsp.k.sbsb.t.sup.+b-1) is the end location of the t.sup.th solid burst error. The processor next labels the b locations i.sub.k.sbsb.t to i.sub.k.sbsb.t.sup.+b- 1as part of the t.sup.th solid burst error {step 206}.
Setting t=t+1, the error location processor 14 determines if t>p {steps 210, 212}. If not, the processor repeats steps 202-208 to determine the first and last locations of the next solid burst error, and labels those locations and the locations between them as part of the solid burst. The processor then repeats steps 202-212 until t exceeds p, which is the number of solid bursts errors in the code word {steps 212-213}.
Referring now to FIG. 5, the operations of the GCD processor 18 in determining the greatest common divisor of .sigma.(x) and .sigma.(.alpha.*x) are discussed. The GCD processor 18 iteratively implements the Euclidean Algorithm, to relatively quickly determine the greatest common divisor, d(x).
The Euclidean Division Algorithm states that for any two polynomials s(x) and t(x), with the degree of t(x) less than or equal to that of s(x):
where q(x) is the quotient of ##EQU4## and r(x) is the remainder, with the degree of r(x) less than that of t(x).
If d(x) is a common divisor of s(x) and t(x), d(x) also divides s(x)-q(x)t(x) , that is, d(x) divides the remainder r(x). Thus, d(x) can be recursively determined from the known s(x) and t(x).
Let s.sub.o (x)=s(x) and t.sub.o (x)=t(x) and q.sub.i (x) be defined as the quotient of s.sub.i (x) and t.sub.i (x). Equation 3 then becomes:
for i>0. For i>0, the degree of s.sub.i (x) is greater than the degree of t.sub.i (x) and the degree of t.sub.i (x) is greater than the degree of t.sub.i+1 (x). To make equation 4 recursive, let
and
where the degree of t.sub.i+1 (x) is greater than the degree of t.sub.i+2 (x). The degree of t.sub.i (x) thus monotonically decreases as i increases.
When the remainder t.sub.n+1 (x)=0, it implies that s.sub.n+1 (x)=t.sub.n (x) is one of the greatest common divisors of s(x) and t(x). In Galois Fields, there are a plurality of greatest common divisors. However, each greatest common divisor differs from another greatest common divisor by a constant non-zero factor, c, over the Galois Field. Thus all the greatest common divisors are of the form c.sub.j *t.sub.n (x).
By convention, a polynomial of the form w(x)=p.sub.k x.sup.k +p.sub.k-1 x.sup.k-1. . . +b has degree k; a polynomial w(x)=px+b has degree 1, and a polynomial w(x)=b has degree 0. Also by convention, a polynomial w(x)=0 has degree -.infin..
The iteration of equation 6 continues until the degree of t.sub.n+1 (x) equals either 0 or -.infin.. If t.sub.n+1 (x) has degree 0, s.sub.0 (x) and t.sub.0 (x) are relatively prime, i.e., they have no non-trivial common factors of degree 1 or higher. If .sigma.(x)=s.sub.0 (x) and .sigma.(.alpha.*x)=t.sub.0 (x)=t.sub.0 (x), .sigma.(x) and .sigma.(.alpha.*x) are relatively prime, and every solid burst error has length 1. There are thus too many errors for the ECC to correct using a direct solution, since the degree p of p.sub.1 (x) is equal to e and e>.beta..
If the degree of t.sub.n+1 (x) is -.infin., i.e., t.sub.n+1 (x)=0, all of the greatest common divisors of s(x) and t(x) have the form c.sub.j *t.sub.n (x). There are thus two lower-degree polynomials p.sub.1 (x) and p.sub.2 (x) that consist of the non-common factors of .sigma.(x) and .sigma.(.alpha.*x).
To recursively solve for the greatest common divisor, equations 5 and 6 are used, where .sigma.(x)=s.sub.0 (x) and .sigma.(.alpha.*x)=t.sub.0 (x) i=0,1. . .
In matrix form these equations are ##EQU5## for i.gtoreq.0. The two-by-two matrix of eqn 7 is non-singular for any q.sub.i (x). Thus the inverse of the matrix can be determined: ##EQU6##
The product of the two 2-by-2 matrices is by definition another 2-by-2 matrix. Accordingly, any s.sub.i (x), t.sub.i (x) pair, i.gtoreq.0, is related to s(x)=s.sub.0 (x) and t(x)=t.sub.0 (x) by a 2-by-2 matrix. This 2-by-2 matrix can be obtained iteratively using the Euclidean Algorithm. ##EQU7##
The determinant of ##EQU8## is -1 for all i. Accordingly, for even i's, .vertline.A.sub.i .vertline.=+1and for odd i's .vertline.A.sub.i .vertline.=-1.
Then, by induction ##EQU9## As discussed above, A.sub.i (x) is updated to A.sub.i+1 (x) in accordance with the Euclidean Algorithm: ##EQU10## and to update, the GCD processor 18 renames ##EQU11##
Then the GCD processor continues updating until either t.sub.i+1 (x)=0, that is, has degree -.infin., or t.sub.i+1 (x) has degree 0.
The inverse of A.sub.i (x) is ##EQU12## where j=i and c.sub.j equals -1 for odd j's and +1 for even j's. In GF(2.sup.M), c.sub.j =1 for all values of j. ##EQU13##
If t.sub.i+1 (x) has degree -.infin., that is, if t.sub.i+1 (x)=0, equation 10 becomes: ##EQU14##
The greatest common divisors are thus of the form c.sub.j * s.sub.i+1 (x), which is the same as c.sub.j *t.sub.i (x). {steps 300+310 }. Since .sigma.(x)=s.sub.0 (x) and .sigma.(.alpha.*x)=t.sub.0 (x), the term a.sub.22.sbsb.i +1 (x) consists of all the factors of .sigma.(x) that are not common with .sigma.(.alpha.*x), and a.sub.21.sbsb.i+1 (x) consists of all the factors of .sigma.(.alpha.*x) that are not common with .sigma.(x). Thus, p.sub.1 (x) =a.sub.22.sbsb.i+1 (x) and p.sub.2 (x) =a.sub.21.sbsb.i+1 (x). Further, a.sub.22.sbsb.i+1 (x) consists of the factors with roots that correspond, respectively, to the first locations of the solid bursts and a.sub.21.sbsb.i+1 (x) consists of the factors with roots that correspond, respectively, to locations that are one location beyond the last locations of the solid bursts.
As discussed above, the error locator processor associates a last location with a first location and labels these locations and all the locations in between them as erroneous and part of the same solid burst error. The system can thus quickly correct up to .beta. solid burst errors of any length, since all of the locations included in the burst errors contain erroneous symbols.
Once the system has determined the error locations, the system then uses conventional methods to determine the error values. The system thereafter uses the error values associated with the error locations included in the solid bursts, to correct the erroneous symbols.