*
* XPRESS compression routines.
*
+ * See the comments in xpress-decomp.c about the XPRESS format.
+ */
+
+/*
* Copyright (C) 2012 Eric Biggers
*
- * wimlib - Library for working with WIM files
+ * This file is part of wimlib, a library for working with WIM files.
*
- * This library is free software; you can redistribute it and/or modify it under
- * the terms of the GNU Lesser General Public License as published by the Free
- * Software Foundation; either version 2.1 of the License, or (at your option) any
- * later version.
+ * wimlib is free software; you can redistribute it and/or modify it under the
+ * terms of the GNU General Public License as published by the Free
+ * Software Foundation; either version 3 of the License, or (at your option)
+ * any later version.
*
- * This library is distributed in the hope that it will be useful, but WITHOUT ANY
- * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A
- * PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details.
+ * wimlib is distributed in the hope that it will be useful, but WITHOUT ANY
+ * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
+ * A PARTICULAR PURPOSE. See the GNU General Public License for more
+ * details.
*
- * You should have received a copy of the GNU Lesser General Public License along
- * with this library; if not, write to the Free Software Foundation, Inc., 59
- * Temple Place, Suite 330, Boston, MA 02111-1307 USA
+ * You should have received a copy of the GNU General Public License
+ * along with wimlib; if not, see http://www.gnu.org/licenses/.
*/
-/* See the comments in xpress-decomp.c about the XPRESS format. */
-
-
#include "xpress.h"
#include "comp.h"
-#include "huffman.h"
#include <stdlib.h>
#include <string.h>
-static inline u32 bsr32(u32 n)
-{
-#if defined(__x86__) || defined(__x86_64__)
- asm("bsrl %0, %0;"
- : "=r"(n)
- : "0" (n));
- return n;
-#else
- u32 pow = 0;
- while ((n >>= 1) != 0)
- pow++;
- return pow;
-#endif
-}
-
-
-/*
+/*
* Writes @match, which is a match given in the intermediate representation for
* XPRESS matches, to the output stream @ostream.
*
* @codewords and @lens provide the Huffman code that is being used.
*/
-static int xpress_write_match(struct output_bitstream *ostream, u32 match,
+static int xpress_write_match(struct output_bitstream *ostream, u32 match,
const u16 codewords[], const u8 lens[])
{
- uint main_sym;
- uint huff_sym;
- uint offset_bsr;
- uint match_len;
- uint match_offset;
+ u32 adjusted_match_len = match & 0xffff;
+ u32 match_offset = match >> 16;
+ u32 len_hdr = min(adjusted_match_len, 0xf);
+ u32 offset_bsr = bsr32(match_offset);
+ u32 sym = len_hdr | (offset_bsr << 4) | XPRESS_NUM_CHARS;
int ret;
- u8 byte1;
- main_sym = (match & 0xff);
- huff_sym = main_sym + XPRESS_NUM_CHARS;
- ret = bitstream_put_bits(ostream, codewords[huff_sym], lens[huff_sym]);
+ ret = bitstream_put_bits(ostream, codewords[sym], lens[sym]);
if (ret != 0)
return ret;
- offset_bsr = main_sym >> 4;
-
- match_len = (match >> 8) & 0xff;
- match_offset = (match >> 16);
-
-
- match_len -= XPRESS_MIN_MATCH;
- if (match_len >= 0xf) {
- byte1 = (u8)(match_len - 0xf);
+ if (adjusted_match_len >= 0xf) {
+ u8 byte1 = min(adjusted_match_len - 0xf, 0xff);
ret = bitstream_put_byte(ostream, byte1);
if (ret != 0)
return ret;
if (byte1 == 0xff) {
- ret = bitstream_put_two_bytes(ostream, match_len);
+ ret = bitstream_put_two_bytes(ostream, adjusted_match_len);
if (ret != 0)
return ret;
}
}
- return bitstream_put_bits(ostream, match_offset ^ (1 << offset_bsr),
- offset_bsr);
+ return bitstream_put_bits(ostream,
+ match_offset ^ (1 << offset_bsr), offset_bsr);
}
-static int xpress_write_compressed_literals(struct output_bitstream *ostream,
- const u32 match_tab[],
- uint num_matches,
- const u16 codewords[],
+static int xpress_write_compressed_literals(struct output_bitstream *ostream,
+ const u32 match_tab[],
+ unsigned num_matches,
+ const u16 codewords[],
const u8 lens[])
{
- uint i;
- u32 match;
- int ret;
+ for (unsigned i = 0; i < num_matches; i++) {
+ int ret;
+ u32 match = match_tab[i];
- for (i = 0; i < num_matches; i++) {
- match = match_tab[i];
- if (match >= XPRESS_NUM_CHARS) {
- /* match */
- ret = xpress_write_match(ostream, match, codewords,
+ if (match >= XPRESS_NUM_CHARS) /* match */
+ ret = xpress_write_match(ostream, match, codewords,
lens);
- if (ret != 0)
- return ret;
- } else {
- /* literal byte */
- ret = bitstream_put_bits(ostream, codewords[match],
+ else /* literal byte */
+ ret = bitstream_put_bits(ostream, codewords[match],
lens[match]);
- if (ret != 0)
- return ret;
- }
+ if (ret != 0)
+ return ret;
}
- return bitstream_put_bits(ostream, codewords[256], lens[256]);
+ return bitstream_put_bits(ostream, codewords[XPRESS_END_OF_DATA],
+ lens[XPRESS_END_OF_DATA]);
}
static u32 xpress_record_literal(u8 literal, void *__freq_tab)
return literal;
}
-static u32 xpress_record_match(uint match_offset, uint match_len,
- void *__freq_tab, void *ignore)
+static u32 xpress_record_match(unsigned match_offset, unsigned match_len,
+ void *freq_tab, void *ignore)
{
- u32 *freq_tab = __freq_tab;
- u32 len_hdr;
- u32 offset_bsr;
- u32 match;
-
- wimlib_assert(match_len >= XPRESS_MIN_MATCH &&
- match_len <= XPRESS_MAX_MATCH);
- wimlib_assert(match_offset > 0);
+ wimlib_assert(match_len >= XPRESS_MIN_MATCH &&
+ match_len <= XPRESS_MAX_MATCH);
+ wimlib_assert(match_offset >= XPRESS_MIN_OFFSET &&
+ match_offset <= XPRESS_MAX_OFFSET);
- len_hdr = min(match_len - XPRESS_MIN_MATCH, 15);
- offset_bsr = bsr32(match_offset);
- match = (offset_bsr << 4) | len_hdr;
- freq_tab[match + XPRESS_NUM_CHARS]++;
- match |= match_len << 8;
- match |= match_offset << 16;
- return match;
+ /*
+ * The intermediate representation of XPRESS matches is as follows:
+ *
+ * bits description
+ * ---- -----------------------------------------------------------
+ *
+ * 16-31 match offset (XPRESS_MIN_OFFSET < x < XPRESS_MAX_OFFSET)
+ *
+ * 0-15 adjusted match length (0 <= x <= XPRESS_MAX_MATCH - XPRESS_MIN_MATCH)
+ *
+ * Literals are simply represented as themselves and can be
+ * distinguished from matches by the fact that only literals will have
+ * the upper three bytes completely clear. */
+
+ u32 adjusted_match_len = match_len - XPRESS_MIN_MATCH;
+ u32 len_hdr = min(adjusted_match_len, 0xf);
+ u32 offset_bsr = bsr32(match_offset);
+ u32 sym = len_hdr | (offset_bsr << 4) | XPRESS_NUM_CHARS;
+ ((u32*)freq_tab)[sym]++;
+ return adjusted_match_len | (match_offset << 16);
}
static const struct lz_params xpress_lz_params = {
- .min_match = 3,
+ .min_match = XPRESS_MIN_MATCH,
.max_match = XPRESS_MAX_MATCH,
.good_match = 16,
.nice_match = 32,
.too_far = 4096,
};
-/*
+/*
* Performs XPRESS compression on a block of data.
*
* @__uncompressed_data: Pointer to the data to be compressed.
* @compressed_data and @compressed_len_ret will contain the compressed data and
* its length. A return value of nonzero means that compressing the data did
* not reduce its size, and @compressed_data will not contain the full
- * compressed data.
+ * compressed data.
*/
-int xpress_compress(const void *__uncompressed_data, uint uncompressed_len,
- void *__compressed_data, uint *compressed_len_ret)
+int xpress_compress(const void *__uncompressed_data, unsigned uncompressed_len,
+ void *__compressed_data, unsigned *compressed_len_ret)
{
const u8 *uncompressed_data = __uncompressed_data;
u8 *compressed_data = __compressed_data;
u32 match_tab[uncompressed_len];
u32 freq_tab[XPRESS_NUM_SYMBOLS];
u16 codewords[XPRESS_NUM_SYMBOLS];
- u8 lens[XPRESS_NUM_SYMBOLS];
- uint num_matches;
- uint compressed_len;
- uint i;
+ u8 lens[XPRESS_NUM_SYMBOLS];
+ unsigned num_matches;
+ unsigned compressed_len;
+ unsigned i;
int ret;
- XPRESS_DEBUG("uncompressed_len = %u\n", uncompressed_len);
-
- if (uncompressed_len < 300)
+ /* XPRESS requires 256 bytes of overhead for the Huffman tables, so it's
+ * impossible cannot compress 256 bytes or less of data to less than the
+ * input size. */
+ if (uncompressed_len <= XPRESS_NUM_SYMBOLS / 2)
return 1;
ZERO_ARRAY(freq_tab);
+ num_matches = lz_analyze_block(uncompressed_data, uncompressed_len,
+ match_tab, xpress_record_match,
+ xpress_record_literal, freq_tab,
+ NULL, freq_tab,
+ &xpress_lz_params);
- num_matches = lz_analyze_block(uncompressed_data, uncompressed_len,
- match_tab, xpress_record_match,
- xpress_record_literal, freq_tab,
- NULL, freq_tab,
- &xpress_lz_params);
-
- XPRESS_DEBUG("using %u matches\n", num_matches);
-
- freq_tab[256]++;
+ freq_tab[XPRESS_END_OF_DATA]++;
make_canonical_huffman_code(XPRESS_NUM_SYMBOLS, XPRESS_MAX_CODEWORD_LEN,
- freq_tab, lens, codewords);
+ freq_tab, lens, codewords);
/* IMPORTANT NOTE:
*
* bitstream_put_bits() will output 2 bytes at a time in little-endian
* order, which is the order that is needed for the compressed literals.
* However, the bytes in the lengths table are in order, so they need to
- * be written one at a time without using bitstream_put_bits().
+ * be written one at a time without using bitstream_put_bits().
*
* Because of this, init_output_bitstream() is not called until after
* the lengths table is output.
for (i = 0; i < XPRESS_NUM_SYMBOLS; i += 2)
*compressed_data++ = (lens[i] & 0xf) | (lens[i + 1] << 4);
- init_output_bitstream(&ostream, compressed_data, uncompressed_len -
- XPRESS_NUM_SYMBOLS / 2 - 1);
+ init_output_bitstream(&ostream, compressed_data,
+ uncompressed_len - XPRESS_NUM_SYMBOLS / 2 - 1);
- ret = xpress_write_compressed_literals(&ostream, match_tab, num_matches,
- codewords, lens);
+ ret = xpress_write_compressed_literals(&ostream, match_tab,
+ num_matches, codewords, lens);
if (ret != 0)
return ret;
+ /* Flush any bits that are buffered. */
ret = flush_output_bitstream(&ostream);
if (ret != 0)
return ret;
+ /* Assert that there are no output bytes between the ostream.output
+ * pointer and the ostream.next_bit_output pointer. This can only
+ * happen if bytes had been written at the ostream.output pointer before
+ * the last bit word was written to the stream. But, this does not
+ * occur since xpress_write_match() always finishes by writing some bits
+ * (a Huffman symbol), and the bitstream was just flushed. */
+ wimlib_assert(ostream.output - ostream.next_bit_output == 2);
+
+ /* The length of the compressed data is supposed to be the value of the
+ * ostream.output pointer before flushing, which is now the
+ * output.next_bit_output pointer after flushing.
+ *
+ * There will be an extra 2 bytes at the ostream.bit_output pointer,
+ * which is zeroed out. (These 2 bytes may be either the last bytes in
+ * the compressed data, in which case they are actually unnecessary, or
+ * they may precede a number of bytes embedded into the bitstream.) */
+ if (ostream.bit_output >
+ (const u8*)__compressed_data + uncompressed_len - 3)
+ return 1;
+ *(u16*)ostream.bit_output = cpu_to_le16(0);
+ compressed_len = ostream.next_bit_output - (const u8*)__compressed_data;
- compressed_len = ostream.output - (u8*)__compressed_data;
-
- XPRESS_DEBUG("Compressed %u => %u bytes\n",
- uncompressed_len, compressed_len);
+ wimlib_assert(compressed_len <= uncompressed_len - 1);
*compressed_len_ret = compressed_len;
#ifdef ENABLE_VERIFY_COMPRESSION
/* Verify that we really get the same thing back when decompressing. */
- XPRESS_DEBUG("Verifying the compressed data.\n");
u8 buf[uncompressed_len];
- ret = xpress_decompress(__compressed_data, compressed_len, buf,
+ ret = xpress_decompress(__compressed_data, compressed_len, buf,
uncompressed_len);
if (ret != 0) {
- fprintf(stderr, "xpress_compress(): Failed to decompress data "
- "we compressed!\n");
+ ERROR("xpress_compress(): Failed to decompress data we "
+ "compressed");
abort();
}
for (i = 0; i < uncompressed_len; i++) {
if (buf[i] != uncompressed_data[i]) {
- fprintf(stderr, "xpress_compress(): Data we compressed "
- "didn't decompress to the original data "
- "(difference at byte %u of %u)\n",
- i + 1, uncompressed_len);
+ ERROR("xpress_compress(): Data we compressed didn't "
+ "decompress to the original data (difference at "
+ "byte %u of %u)", i + 1, uncompressed_len);
abort();
}
}
- XPRESS_DEBUG("Compression verified to be correct.\n");
#endif
-
return 0;
-
}