| 1 |
|
| 2 | |
| 3 | * huff0 huffman codec,
|
| 4 | * part of Finite State Entropy library
|
| 5 | * Copyright (c) Meta Platforms, Inc. and affiliates.
|
| 6 | *
|
| 7 | * You can contact the author at :
|
| 8 | * - Source repository : https://github.com/Cyan4973/FiniteStateEntropy
|
| 9 | *
|
| 10 | * This source code is licensed under both the BSD-style license (found in the
|
| 11 | * LICENSE file in the root directory of this source tree) and the GPLv2 (found
|
| 12 | * in the COPYING file in the root directory of this source tree).
|
| 13 | * You may select, at your option, one of the above-listed licenses.
|
| 14 | ******************************************************************
|
| 15 |
|
| 16 | #ifndef HUF_H_298734234
|
| 17 | #define HUF_H_298734234
|
| 18 |
|
| 19 |
|
| 20 | #include "zstd_deps.h"
|
| 21 | #include "mem.h"
|
| 22 | #define FSE_STATIC_LINKING_ONLY
|
| 23 | #include "fse.h"
|
| 24 |
|
| 25 |
|
| 26 | #define HUF_BLOCKSIZE_MAX (128 * 1024)
|
| 27 | size_t HUF_compressBound(size_t size);
|
| 28 |
|
| 29 |
|
| 30 | unsigned HUF_isError(size_t code);
|
| 31 | const char* HUF_getErrorName(size_t code);
|
| 32 |
|
| 33 |
|
| 34 | #define HUF_WORKSPACE_SIZE ((8 << 10) + 512 )
|
| 35 | #define HUF_WORKSPACE_SIZE_U64 (HUF_WORKSPACE_SIZE / sizeof(U64))
|
| 36 |
|
| 37 |
|
| 38 | #define HUF_TABLELOG_MAX 12
|
| 39 | #define HUF_TABLELOG_DEFAULT 11
|
| 40 | #define HUF_SYMBOLVALUE_MAX 255
|
| 41 |
|
| 42 | #define HUF_TABLELOG_ABSOLUTEMAX 12
|
| 43 | #if (HUF_TABLELOG_MAX > HUF_TABLELOG_ABSOLUTEMAX)
|
| 44 | # error "HUF_TABLELOG_MAX is too large !"
|
| 45 | #endif
|
| 46 |
|
| 47 |
|
| 48 | |
| 49 | * Static allocation
|
| 50 | *****************************************
|
| 51 |
|
| 52 | #define HUF_CTABLEBOUND 129
|
| 53 | #define HUF_BLOCKBOUND(size) (size + (size>>8) + 8)
|
| 54 | #define HUF_COMPRESSBOUND(size) (HUF_CTABLEBOUND + HUF_BLOCKBOUND(size))
|
| 55 |
|
| 56 |
|
| 57 |
|
| 58 | typedef size_t HUF_CElt;
|
| 59 | #define HUF_CTABLE_SIZE_ST(maxSymbolValue) ((maxSymbolValue)+2)
|
| 60 | #define HUF_CTABLE_SIZE(maxSymbolValue) (HUF_CTABLE_SIZE_ST(maxSymbolValue) * sizeof(size_t))
|
| 61 | #define HUF_CREATE_STATIC_CTABLE(name, maxSymbolValue) \
|
| 62 | HUF_CElt name[HUF_CTABLE_SIZE_ST(maxSymbolValue)]
|
| 63 |
|
| 64 |
|
| 65 | typedef U32 HUF_DTable;
|
| 66 | #define HUF_DTABLE_SIZE(maxTableLog) (1 + (1<<(maxTableLog)))
|
| 67 | #define HUF_CREATE_STATIC_DTABLEX1(DTable, maxTableLog) \
|
| 68 | HUF_DTable DTable[HUF_DTABLE_SIZE((maxTableLog)-1)] = { ((U32)((maxTableLog)-1) * 0x01000001) }
|
| 69 | #define HUF_CREATE_STATIC_DTABLEX2(DTable, maxTableLog) \
|
| 70 | HUF_DTable DTable[HUF_DTABLE_SIZE(maxTableLog)] = { ((U32)(maxTableLog) * 0x01000001) }
|
| 71 |
|
| 72 |
|
| 73 | |
| 74 | * Advanced decompression functions
|
| 75 | *****************************************
|
| 76 |
|
| 77 | |
| 78 | * Huffman flags bitset.
|
| 79 | * For all flags, 0 is the default value.
|
| 80 |
|
| 81 | typedef enum {
|
| 82 | |
| 83 | * If compiled with DYNAMIC_BMI2: Set flag only if the CPU supports BMI2 at runtime.
|
| 84 | * Otherwise: Ignored.
|
| 85 |
|
| 86 | HUF_flags_bmi2 = (1 << 0),
|
| 87 | |
| 88 | * If set: Test possible table depths to find the one that produces the smallest header + encoded size.
|
| 89 | * If unset: Use heuristic to find the table depth.
|
| 90 |
|
| 91 | HUF_flags_optimalDepth = (1 << 1),
|
| 92 | |
| 93 | * If set: If the previous table can encode the input, always reuse the previous table.
|
| 94 | * If unset: If the previous table can encode the input, reuse the previous table if it results in a smaller output.
|
| 95 |
|
| 96 | HUF_flags_preferRepeat = (1 << 2),
|
| 97 | |
| 98 | * If set: Sample the input and check if the sample is uncompressible, if it is then don't attempt to compress.
|
| 99 | * If unset: Always histogram the entire input.
|
| 100 |
|
| 101 | HUF_flags_suspectUncompressible = (1 << 3),
|
| 102 | |
| 103 | * If set: Don't use assembly implementations
|
| 104 | * If unset: Allow using assembly implementations
|
| 105 |
|
| 106 | HUF_flags_disableAsm = (1 << 4),
|
| 107 | |
| 108 | * If set: Don't use the fast decoding loop, always use the fallback decoding loop.
|
| 109 | * If unset: Use the fast decoding loop when possible.
|
| 110 |
|
| 111 | HUF_flags_disableFast = (1 << 5)
|
| 112 | } HUF_flags_e;
|
| 113 |
|
| 114 |
|
| 115 | |
| 116 | * HUF detailed API
|
| 117 | * ***************************************
|
| 118 | #define HUF_OPTIMAL_DEPTH_THRESHOLD ZSTD_btultra
|
| 119 |
|
| 120 | |
| 121 | * 1. count symbol occurrence from source[] into table count[] using FSE_count() (exposed within "fse.h")
|
| 122 | * 2. (optional) refine tableLog using HUF_optimalTableLog()
|
| 123 | * 3. build Huffman table from count using HUF_buildCTable()
|
| 124 | * 4. save Huffman table to memory buffer using HUF_writeCTable()
|
| 125 | * 5. encode the data stream using HUF_compress4X_usingCTable()
|
| 126 | *
|
| 127 | * The following API allows targeting specific sub-functions for advanced tasks.
|
| 128 | * For example, it's possible to compress several blocks using the same 'CTable',
|
| 129 | * or to save and regenerate 'CTable' using external methods.
|
| 130 |
|
| 131 | unsigned HUF_minTableLog(unsigned symbolCardinality);
|
| 132 | unsigned HUF_cardinality(const unsigned* count, unsigned maxSymbolValue);
|
| 133 | unsigned HUF_optimalTableLog(unsigned maxTableLog, size_t srcSize, unsigned maxSymbolValue, void* workSpace,
|
| 134 | size_t wkspSize, HUF_CElt* table, const unsigned* count, int flags);
|
| 135 | size_t HUF_writeCTable_wksp(void* dst, size_t maxDstSize, const HUF_CElt* CTable, unsigned maxSymbolValue, unsigned huffLog, void* workspace, size_t workspaceSize);
|
| 136 | size_t HUF_compress4X_usingCTable(void* dst, size_t dstSize, const void* src, size_t srcSize, const HUF_CElt* CTable, int flags);
|
| 137 | size_t HUF_estimateCompressedSize(const HUF_CElt* CTable, const unsigned* count, unsigned maxSymbolValue);
|
| 138 | int HUF_validateCTable(const HUF_CElt* CTable, const unsigned* count, unsigned maxSymbolValue);
|
| 139 |
|
| 140 | typedef enum {
|
| 141 | HUF_repeat_none,
|
| 142 | HUF_repeat_check,
|
| 143 | HUF_repeat_valid
|
| 144 | } HUF_repeat;
|
| 145 |
|
| 146 | |
| 147 | * Same as HUF_compress4X_wksp(), but considers using hufTable if *repeat != HUF_repeat_none.
|
| 148 | * If it uses hufTable it does not modify hufTable or repeat.
|
| 149 | * If it doesn't, it sets *repeat = HUF_repeat_none, and it sets hufTable to the table used.
|
| 150 | * If preferRepeat then the old table will always be used if valid.
|
| 151 | * If suspectUncompressible then some sampling checks will be run to potentially skip huffman coding
|
| 152 | size_t HUF_compress4X_repeat(void* dst, size_t dstSize,
|
| 153 | const void* src, size_t srcSize,
|
| 154 | unsigned maxSymbolValue, unsigned tableLog,
|
| 155 | void* workSpace, size_t wkspSize,
|
| 156 | HUF_CElt* hufTable, HUF_repeat* repeat, int flags);
|
| 157 |
|
| 158 | |
| 159 | * Same as HUF_buildCTable(), but using externally allocated scratch buffer.
|
| 160 | * `workSpace` must be aligned on 4-bytes boundaries, and its size must be >= HUF_CTABLE_WORKSPACE_SIZE.
|
| 161 |
|
| 162 | #define HUF_CTABLE_WORKSPACE_SIZE_U32 ((4 * (HUF_SYMBOLVALUE_MAX + 1)) + 192)
|
| 163 | #define HUF_CTABLE_WORKSPACE_SIZE (HUF_CTABLE_WORKSPACE_SIZE_U32 * sizeof(unsigned))
|
| 164 | size_t HUF_buildCTable_wksp (HUF_CElt* tree,
|
| 165 | const unsigned* count, U32 maxSymbolValue, U32 maxNbBits,
|
| 166 | void* workSpace, size_t wkspSize);
|
| 167 |
|
| 168 | |
| 169 | * Read compact Huffman tree, saved by HUF_writeCTable().
|
| 170 | * `huffWeight` is destination buffer.
|
| 171 | * @return : size read from `src` , or an error Code .
|
| 172 | * Note : Needed by HUF_readCTable() and HUF_readDTableXn() .
|
| 173 | size_t HUF_readStats(BYTE* huffWeight, size_t hwSize,
|
| 174 | U32* rankStats, U32* nbSymbolsPtr, U32* tableLogPtr,
|
| 175 | const void* src, size_t srcSize);
|
| 176 |
|
| 177 | |
| 178 | * Same as HUF_readStats() but takes an external workspace which must be
|
| 179 | * 4-byte aligned and its size must be >= HUF_READ_STATS_WORKSPACE_SIZE.
|
| 180 | * If the CPU has BMI2 support, pass bmi2=1, otherwise pass bmi2=0.
|
| 181 |
|
| 182 | #define HUF_READ_STATS_WORKSPACE_SIZE_U32 FSE_DECOMPRESS_WKSP_SIZE_U32(6, HUF_TABLELOG_MAX-1)
|
| 183 | #define HUF_READ_STATS_WORKSPACE_SIZE (HUF_READ_STATS_WORKSPACE_SIZE_U32 * sizeof(unsigned))
|
| 184 | size_t HUF_readStats_wksp(BYTE* huffWeight, size_t hwSize,
|
| 185 | U32* rankStats, U32* nbSymbolsPtr, U32* tableLogPtr,
|
| 186 | const void* src, size_t srcSize,
|
| 187 | void* workspace, size_t wkspSize,
|
| 188 | int flags);
|
| 189 |
|
| 190 | |
| 191 | * Loading a CTable saved with HUF_writeCTable()
|
| 192 | size_t HUF_readCTable (HUF_CElt* CTable, unsigned* maxSymbolValuePtr, const void* src, size_t srcSize, unsigned *hasZeroWeights);
|
| 193 |
|
| 194 | |
| 195 | * Read nbBits from CTable symbolTable, for symbol `symbolValue` presumed <= HUF_SYMBOLVALUE_MAX
|
| 196 | * Note 1 : If symbolValue > HUF_readCTableHeader(symbolTable).maxSymbolValue, returns 0
|
| 197 | * Note 2 : is not inlined, as HUF_CElt definition is private
|
| 198 |
|
| 199 | U32 HUF_getNbBitsFromCTable(const HUF_CElt* symbolTable, U32 symbolValue);
|
| 200 |
|
| 201 | typedef struct {
|
| 202 | BYTE tableLog;
|
| 203 | BYTE maxSymbolValue;
|
| 204 | BYTE unused[sizeof(size_t) - 2];
|
| 205 | } HUF_CTableHeader;
|
| 206 |
|
| 207 | |
| 208 | * @returns The header from the CTable specifying the tableLog and the maxSymbolValue.
|
| 209 |
|
| 210 | HUF_CTableHeader HUF_readCTableHeader(HUF_CElt const* ctable);
|
| 211 |
|
| 212 | |
| 213 | * HUF_decompress() does the following:
|
| 214 | * 1. select the decompression algorithm (X1, X2) based on pre-computed heuristics
|
| 215 | * 2. build Huffman table from save, using HUF_readDTableX?()
|
| 216 | * 3. decode 1 or 4 segments in parallel using HUF_decompress?X?_usingDTable()
|
| 217 |
|
| 218 |
|
| 219 | |
| 220 | * Tells which decoder is likely to decode faster,
|
| 221 | * based on a set of pre-computed metrics.
|
| 222 | * @return : 0==HUF_decompress4X1, 1==HUF_decompress4X2 .
|
| 223 | * Assumption : 0 < dstSize <= 128 KB
|
| 224 | U32 HUF_selectDecoder (size_t dstSize, size_t cSrcSize);
|
| 225 |
|
| 226 | |
| 227 | * The minimum workspace size for the `workSpace` used in
|
| 228 | * HUF_readDTableX1_wksp() and HUF_readDTableX2_wksp().
|
| 229 | *
|
| 230 | * The space used depends on HUF_TABLELOG_MAX, ranging from ~1500 bytes when
|
| 231 | * HUF_TABLE_LOG_MAX=12 to ~1850 bytes when HUF_TABLE_LOG_MAX=15.
|
| 232 | * Buffer overflow errors may potentially occur if code modifications result in
|
| 233 | * a required workspace size greater than that specified in the following
|
| 234 | * macro.
|
| 235 |
|
| 236 | #define HUF_DECOMPRESS_WORKSPACE_SIZE ((2 << 10) + (1 << 9))
|
| 237 | #define HUF_DECOMPRESS_WORKSPACE_SIZE_U32 (HUF_DECOMPRESS_WORKSPACE_SIZE / sizeof(U32))
|
| 238 |
|
| 239 |
|
| 240 | |
| 241 |
|
| 242 | |
| 243 |
|
| 244 | size_t HUF_compress1X_usingCTable(void* dst, size_t dstSize, const void* src, size_t srcSize, const HUF_CElt* CTable, int flags);
|
| 245 | |
| 246 | * Same as HUF_compress1X_wksp(), but considers using hufTable if *repeat != HUF_repeat_none.
|
| 247 | * If it uses hufTable it does not modify hufTable or repeat.
|
| 248 | * If it doesn't, it sets *repeat = HUF_repeat_none, and it sets hufTable to the table used.
|
| 249 | * If preferRepeat then the old table will always be used if valid.
|
| 250 | * If suspectUncompressible then some sampling checks will be run to potentially skip huffman coding
|
| 251 | size_t HUF_compress1X_repeat(void* dst, size_t dstSize,
|
| 252 | const void* src, size_t srcSize,
|
| 253 | unsigned maxSymbolValue, unsigned tableLog,
|
| 254 | void* workSpace, size_t wkspSize,
|
| 255 | HUF_CElt* hufTable, HUF_repeat* repeat, int flags);
|
| 256 |
|
| 257 | size_t HUF_decompress1X_DCtx_wksp(HUF_DTable* dctx, void* dst, size_t dstSize, const void* cSrc, size_t cSrcSize, void* workSpace, size_t wkspSize, int flags);
|
| 258 | #ifndef HUF_FORCE_DECOMPRESS_X1
|
| 259 | size_t HUF_decompress1X2_DCtx_wksp(HUF_DTable* dctx, void* dst, size_t dstSize, const void* cSrc, size_t cSrcSize, void* workSpace, size_t wkspSize, int flags);
|
| 260 | #endif
|
| 261 |
|
| 262 | |
| 263 | * If the CPU has BMI2 support, pass bmi2=1, otherwise pass bmi2=0.
|
| 264 |
|
| 265 | size_t HUF_decompress1X_usingDTable(void* dst, size_t maxDstSize, const void* cSrc, size_t cSrcSize, const HUF_DTable* DTable, int flags);
|
| 266 | #ifndef HUF_FORCE_DECOMPRESS_X2
|
| 267 | size_t HUF_decompress1X1_DCtx_wksp(HUF_DTable* dctx, void* dst, size_t dstSize, const void* cSrc, size_t cSrcSize, void* workSpace, size_t wkspSize, int flags);
|
| 268 | #endif
|
| 269 | size_t HUF_decompress4X_usingDTable(void* dst, size_t maxDstSize, const void* cSrc, size_t cSrcSize, const HUF_DTable* DTable, int flags);
|
| 270 | size_t HUF_decompress4X_hufOnly_wksp(HUF_DTable* dctx, void* dst, size_t dstSize, const void* cSrc, size_t cSrcSize, void* workSpace, size_t wkspSize, int flags);
|
| 271 | #ifndef HUF_FORCE_DECOMPRESS_X2
|
| 272 | size_t HUF_readDTableX1_wksp(HUF_DTable* DTable, const void* src, size_t srcSize, void* workSpace, size_t wkspSize, int flags);
|
| 273 | #endif
|
| 274 | #ifndef HUF_FORCE_DECOMPRESS_X1
|
| 275 | size_t HUF_readDTableX2_wksp(HUF_DTable* DTable, const void* src, size_t srcSize, void* workSpace, size_t wkspSize, int flags);
|
| 276 | #endif
|
| 277 |
|
| 278 | #endif
|