| 1 | 1. Compression algorithm (deflate)
|
|---|
| 2 |
|
|---|
| 3 | The deflation algorithm used by gzip (also zip and zlib) 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 the behavior of LZW-based _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
|
|---|
| 39 | a 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 process of lazy evaluation begins again. Otherwise,
|
|---|
| 42 | the original match is kept, and the next match search is attempted only N
|
|---|
| 43 | 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 | 2.1 Introduction
|
|---|
| 61 |
|
|---|
| 62 | The real question is, given a Huffman tree, how to decode fast. The most
|
|---|
| 63 | important realization is that shorter codes are much more common than
|
|---|
| 64 | longer codes, so pay attention to decoding the short codes fast, and let
|
|---|
| 65 | the long codes take longer to decode.
|
|---|
| 66 |
|
|---|
| 67 | inflate() sets up a first level table that covers some number of bits of
|
|---|
| 68 | input less than the length of longest code. It gets that many bits from the
|
|---|
| 69 | stream, and looks it up in the table. The table will tell if the next
|
|---|
| 70 | code is that many bits or less and how many, and if it is, it will tell
|
|---|
| 71 | the value, else it will point to the next level table for which inflate()
|
|---|
| 72 | grabs more bits and tries to decode a longer code.
|
|---|
| 73 |
|
|---|
| 74 | How many bits to make the first lookup is a tradeoff between the time it
|
|---|
| 75 | takes to decode and the time it takes to build the table. If building the
|
|---|
| 76 | table took no time (and if you had infinite memory), then there would only
|
|---|
| 77 | be a first level table to cover all the way to the longest code. However,
|
|---|
| 78 | building the table ends up taking a lot longer for more bits since short
|
|---|
| 79 | codes are replicated many times in such a table. What inflate() does is
|
|---|
| 80 | simply to make the number of bits in the first table a variable, and set it
|
|---|
| 81 | for the maximum speed.
|
|---|
| 82 |
|
|---|
| 83 | inflate() sends new trees relatively often, so it is possibly set for a
|
|---|
| 84 | smaller first level table than an application that has only one tree for
|
|---|
|
|---|