diff options
author | Mark Adler <madler@alumni.caltech.edu> | 2011-09-09 23:15:17 -0700 |
---|---|---|
committer | Mark Adler <madler@alumni.caltech.edu> | 2011-09-09 23:15:17 -0700 |
commit | e26a448e9673d67dc2866e11a48d24fc352e5f80 (patch) | |
tree | a2dc10debfd34979dde501afbf1cce6a6c3cbcb7 /algorithm.doc | |
parent | 423eb40306489f9c88f7dba32c2f69179166730b (diff) | |
download | zlib-1.0.2.tar.gz zlib-1.0.2.tar.bz2 zlib-1.0.2.zip |
zlib 1.0.2v1.0.2
Diffstat (limited to 'algorithm.doc')
-rw-r--r-- | algorithm.doc | 105 |
1 files changed, 105 insertions, 0 deletions
diff --git a/algorithm.doc b/algorithm.doc new file mode 100644 index 0000000..156058a --- /dev/null +++ b/algorithm.doc | |||
@@ -0,0 +1,105 @@ | |||
1 | 1. Compression algorithm (deflate) | ||
2 | |||
3 | The deflation algorithm used by zlib (also zip and gzip) is a variation of | ||
4 | LZ77 (Lempel-Ziv 1977, see reference below). It finds duplicated strings in | ||
5 | the input data. The second occurrence of a string is replaced by a | ||
6 | pointer to the previous string, in the form of a pair (distance, | ||
7 | length). Distances are limited to 32K bytes, and lengths are limited | ||
8 | to 258 bytes. When a string does not occur anywhere in the previous | ||
9 | 32K bytes, it is emitted as a sequence of literal bytes. (In this | ||
10 | description, 'string' must be taken as an arbitrary sequence of bytes, | ||
11 | and is not restricted to printable characters.) | ||
12 | |||
13 | Literals or match lengths are compressed with one Huffman tree, and | ||
14 | match distances are compressed with another tree. The trees are stored | ||
15 | in a compact form at the start of each block. The blocks can have any | ||
16 | size (except that the compressed data for one block must fit in | ||
17 | available memory). A block is terminated when deflate() determines that | ||
18 | it would be useful to start another block with fresh trees. (This is | ||
19 | somewhat similar to compress.) | ||
20 | |||
21 | Duplicated strings are found using a hash table. All input strings of | ||
22 | length 3 are inserted in the hash table. A hash index is computed for | ||
23 | the next 3 bytes. If the hash chain for this index is not empty, all | ||
24 | strings in the chain are compared with the current input string, and | ||
25 | the longest match is selected. | ||
26 | |||
27 | The hash chains are searched starting with the most recent strings, to | ||
28 | favor small distances and thus take advantage of the Huffman encoding. | ||
29 | The hash chains are singly linked. There are no deletions from the | ||
30 | hash chains, the algorithm simply discards matches that are too old. | ||
31 | |||
32 | To avoid a worst-case situation, very long hash chains are arbitrarily | ||
33 | truncated at a certain length, determined by a runtime option (level | ||
34 | parameter of deflateInit). So deflate() does not always find the longest | ||
35 | possible match but generally finds a match which is long enough. | ||
36 | |||
37 | deflate() also defers the selection of matches with a lazy evaluation | ||
38 | mechanism. After a match of length N has been found, deflate() searches for a | ||
39 | longer match at the next input byte. If a longer match is found, the | ||
40 | previous match is truncated to a length of one (thus producing a single | ||
41 | literal byte) and the longer match is emitted afterwards. Otherwise, | ||
42 | the original match is kept, and the next match search is attempted only | ||
43 | N steps later. | ||
44 | |||
45 | The lazy match evaluation is also subject to a runtime parameter. If | ||
46 | the current match is long enough, deflate() reduces the search for a longer | ||
47 | match, thus speeding up the whole process. If compression ratio is more | ||
48 | important than speed, deflate() attempts a complete second search even if | ||
49 | the first match is already long enough. | ||
50 | |||
51 | The lazy match evaluation is not performed for the fastest compression | ||
52 | modes (level parameter 1 to 3). For these fast modes, new strings | ||
53 | are inserted in the hash table only when no match was found, or | ||
54 | when the match is not too long. This degrades the compression ratio | ||
55 | but saves time since there are both fewer insertions and fewer searches. | ||
56 | |||
57 | |||
58 | 2. Decompression algorithm (inflate) | ||
59 | |||
60 | The real question is given a Huffman tree, how to decode fast. The most | ||
61 | important realization is that shorter codes are much more common than | ||
62 | longer codes, so pay attention to decoding the short codes fast, and let | ||
63 | the long codes take longer to decode. | ||
64 | |||
65 | inflate() sets up a first level table that covers some number of bits of | ||
66 | input less than the length of longest code. It gets that many bits from the | ||
67 | stream, and looks it up in the table. The table will tell if the next | ||
68 | code is that many bits or less and how many, and if it is, it will tell | ||
69 | the value, else it will point to the next level table for which inflate() | ||
70 | grabs more bits and tries to decode a longer code. | ||
71 | |||
72 | How many bits to make the first lookup is a tradeoff between the time it | ||
73 | takes to decode and the time it takes to build the table. If building the | ||
74 | table took no time (and if you had infinite memory), then there would only | ||
75 | be a first level table to cover all the way to the longest code. However, | ||
76 | building the table ends up taking a lot longer for more bits since short | ||
77 | codes are replicated many times in such a table. What inflate() does is | ||
78 | simply to make the number of bits in the first table a variable, and set it | ||
79 | for the maximum speed. | ||
80 | |||
81 | inflate() sends new trees relatively often, so it is possibly set for a | ||
82 | smaller first level table than an application that has only one tree for | ||
83 | all the data. For inflate, which has 286 possible codes for the | ||
84 | literal/length tree, the size of the first table is nine bits. Also the | ||
85 | distance trees have 30 possible values, and the size of the first table is | ||
86 | six bits. Note that for each of those cases, the table ended up one bit | ||
87 | longer than the "average" code length, i.e. the code length of an | ||
88 | approximately flat code which would be a little more than eight bits for | ||
89 | 286 symbols and a little less than five bits for 30 symbols. It would be | ||
90 | interesting to see if optimizing the first level table for other | ||
91 | applications gave values within a bit or two of the flat code size. | ||
92 | |||
93 | |||
94 | Jean-loup Gailly Mark Adler | ||
95 | gzip@prep.ai.mit.edu madler@alumni.caltech.edu | ||
96 | |||
97 | |||
98 | References: | ||
99 | |||
100 | [LZ77] Ziv J., Lempel A., "A Universal Algorithm for Sequential Data | ||
101 | Compression", IEEE Transactions on Information Theory", Vol. 23, No. 3, | ||
102 | pp. 337-343. | ||
103 | |||
104 | "DEFLATE Compressed Data Format Specification" available in | ||
105 | ftp://ds.internic.net/rfc/rfc1951.txt | ||