/* ******************************************************************

   Huffman decoder, part of New Generation Entropy library

   Copyright (C) 2013-2016, Yann Collet.

   BSD 2-Clause License (http://www.opensource.org/licenses/bsd-license.php)

   Redistribution and use in source and binary forms, with or without

   modification, are permitted provided that the following conditions are

   met:

       * Redistributions of source code must retain the above copyright

   notice, this list of conditions and the following disclaimer.

       * Redistributions in binary form must reproduce the above

   copyright notice, this list of conditions and the following disclaimer

   in the documentation and/or other materials provided with the

   distribution.

   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS

   "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT

   LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR

   A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT

   OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,

   SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT

   LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,

   DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY

   THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT

   (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE

   OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

    You can contact the author at :

    - FSE+HUF source repository : https://github.com/Cyan4973/FiniteStateEntropy

    - Public forum : https://groups.google.com/forum/#!forum/lz4c

****************************************************************** */

/* **************************************************************

*  Dependencies

****************************************************************/

#include <string.h>     /* memcpy, memset */

#include "bitstream.h"  /* BIT_* */

#include "compiler.h"

#include "fse.h"        /* header compression */

#define HUF_STATIC_LINKING_ONLY

#include "huf.h"

#include "error_private.h"

/* **************************************************************

*  Error Management

****************************************************************/

#define HUF_isError ERR_isError

#define HUF_STATIC_ASSERT(c) { enum { HUF_static_assert = 1/(int)(!!(c)) }; }   /* use only *after* variable declarations */

#define CHECK_F(f) { size_t const err_ = (f); if (HUF_isError(err_)) return err_; }

/* **************************************************************

*  Byte alignment for workSpace management

****************************************************************/

#define HUF_ALIGN(x, a)         HUF_ALIGN_MASK((x), (a) - 1)

#define HUF_ALIGN_MASK(x, mask) (((x) + (mask)) & ~(mask))

/*-***************************/

/*  generic DTableDesc       */

/*-***************************/

typedef struct { BYTE maxTableLog; BYTE tableType; BYTE tableLog; BYTE reserved; } DTableDesc;

static DTableDesc HUF_getDTableDesc(const HUF_DTable* table)

{

    DTableDesc dtd;

    memcpy(&dtd, table, sizeof(dtd));

    return dtd;

}

/*-***************************/

/*  single-symbol decoding   */

/*-***************************/

typedef struct { BYTE byte; BYTE nbBits; } HUF_DEltX2;   /* single-symbol decoding */

size_t HUF_readDTableX2_wksp(HUF_DTable* DTable, const void* src, size_t srcSize, void* workSpace, size_t wkspSize)

{

    U32 tableLog = 0;

    U32 nbSymbols = 0;

    size_t iSize;

    void* const dtPtr = DTable + 1;

    HUF_DEltX2* const dt = (HUF_DEltX2*)dtPtr;

    U32* rankVal;

    BYTE* huffWeight;

    size_t spaceUsed32 = 0;

    rankVal = (U32 *)workSpace + spaceUsed32;

    spaceUsed32 += HUF_TABLELOG_ABSOLUTEMAX + 1;

    huffWeight = (BYTE *)((U32 *)workSpace + spaceUsed32);

    spaceUsed32 += HUF_ALIGN(HUF_SYMBOLVALUE_MAX + 1, sizeof(U32)) >> 2;

    if ((spaceUsed32 << 2) > wkspSize) return ERROR(tableLog_tooLarge);

    HUF_STATIC_ASSERT(sizeof(DTableDesc) == sizeof(HUF_DTable));

    /* memset(huffWeight, 0, sizeof(huffWeight)); */   /* is not necessary, even though some analyzer complain ... */

    iSize = HUF_readStats(huffWeight, HUF_SYMBOLVALUE_MAX + 1, rankVal, &nbSymbols, &tableLog, src, srcSize);

    if (HUF_isError(iSize)) return iSize;

    /* Table header */

    {   DTableDesc dtd = HUF_getDTableDesc(DTable);

        if (tableLog > (U32)(dtd.maxTableLog+1)) return ERROR(tableLog_tooLarge);   /* DTable too small, Huffman tree cannot fit in */

        dtd.tableType = 0;

        dtd.tableLog = (BYTE)tableLog;

        memcpy(DTable, &dtd, sizeof(dtd));

    }

    /* Calculate starting value for each rank */

    {   U32 n, nextRankStart = 0;

        for (n=1; n<tableLog+1; n++) {

            U32 const current = nextRankStart;

            nextRankStart += (rankVal[n] << (n-1));

            rankVal[n] = current;

    }   }

    /* fill DTable */

    {   U32 n;

        for (n=0; n<nbSymbols; n++) {

            U32 const w = huffWeight[n];

            U32 const length = (1 << w) >> 1;

            U32 u;

            HUF_DEltX2 D;

            D.byte = (BYTE)n; D.nbBits = (BYTE)(tableLog + 1 - w);

            for (u = rankVal[w]; u < rankVal[w] + length; u++)

                dt[u] = D;

            rankVal[w] += length;

    }   }

    return iSize;

}

size_t HUF_readDTableX2(HUF_DTable* DTable, const void* src, size_t srcSize)

{

    U32 workSpace[HUF_DECOMPRESS_WORKSPACE_SIZE_U32];

    return HUF_readDTableX2_wksp(DTable, src, srcSize,

                                 workSpace, sizeof(workSpace));

}

typedef struct { U16 sequence; BYTE nbBits; BYTE length; } HUF_DEltX4;  /* double-symbols decoding */

FORCE_INLINE_TEMPLATE BYTE

HUF_decodeSymbolX2(BIT_DStream_t* Dstream, const HUF_DEltX2* dt, const U32 dtLog)

{

    size_t const val = BIT_lookBitsFast(Dstream, dtLog); /* note : dtLog >= 1 */

    BYTE const c = dt[val].byte;

    BIT_skipBits(Dstream, dt[val].nbBits);

    return c;

}

#define HUF_DECODE_SYMBOLX2_0(ptr, DStreamPtr) \

    *ptr++ = HUF_decodeSymbolX2(DStreamPtr, dt, dtLog)

#define HUF_DECODE_SYMBOLX2_1(ptr, DStreamPtr)  \

    if (MEM_64bits() || (HUF_TABLELOG_MAX<=12)) \

        HUF_DECODE_SYMBOLX2_0(ptr, DStreamPtr)

#define HUF_DECODE_SYMBOLX2_2(ptr, DStreamPtr) \

    if (MEM_64bits()) \

        HUF_DECODE_SYMBOLX2_0(ptr, DStreamPtr)

HINT_INLINE size_t

HUF_decodeStreamX2(BYTE* p, BIT_DStream_t* const bitDPtr, BYTE* const pEnd, const HUF_DEltX2* const dt, const U32 dtLog)

{

    BYTE* const pStart = p;

    /* up to 4 symbols at a time */

    while ((BIT_reloadDStream(bitDPtr) == BIT_DStream_unfinished) & (p < pEnd-3)) {

        HUF_DECODE_SYMBOLX2_2(p, bitDPtr);

        HUF_DECODE_SYMBOLX2_1(p, bitDPtr);

        HUF_DECODE_SYMBOLX2_2(p, bitDPtr);

        HUF_DECODE_SYMBOLX2_0(p, bitDPtr);

    }

    /* [0-3] symbols remaining */

    if (MEM_32bits())

        while ((BIT_reloadDStream(bitDPtr) == BIT_DStream_unfinished) & (p < pEnd))

            HUF_DECODE_SYMBOLX2_0(p, bitDPtr);

    /* no more data to retrieve from bitstream, no need to reload */

    while (p < pEnd)

        HUF_DECODE_SYMBOLX2_0(p, bitDPtr);

    return pEnd-pStart;

}

FORCE_INLINE_TEMPLATE size_t

HUF_decompress1X2_usingDTable_internal_body(

          void* dst,  size_t dstSize,

    const void* cSrc, size_t cSrcSize,

    const HUF_DTable* DTable)

{

    BYTE* op = (BYTE*)dst;

    BYTE* const oend = op + dstSize;

    const void* dtPtr = DTable + 1;

    const HUF_DEltX2* const dt = (const HUF_DEltX2*)dtPtr;

    BIT_DStream_t bitD;

    DTableDesc const dtd = HUF_getDTableDesc(DTable);

    U32 const dtLog = dtd.tableLog;

    CHECK_F( BIT_initDStream(&bitD, cSrc, cSrcSize) );

    HUF_decodeStreamX2(op, &bitD, oend, dt, dtLog);

    if (!BIT_endOfDStream(&bitD)) return ERROR(corruption_detected);

    return dstSize;

}

FORCE_INLINE_TEMPLATE size_t

HUF_decompress4X2_usingDTable_internal_body(

          void* dst,  size_t dstSize,

    const void* cSrc, size_t cSrcSize,

    const HUF_DTable* DTable)

{

    /* Check */

    if (cSrcSize < 10) return ERROR(corruption_detected);  /* strict minimum : jump table + 1 byte per stream */

    {   const BYTE* const istart = (const BYTE*) cSrc;

        BYTE* const ostart = (BYTE*) dst;

        BYTE* const oend = ostart + dstSize;

        const void* const dtPtr = DTable + 1;

        const HUF_DEltX2* const dt = (const HUF_DEltX2*)dtPtr;

        /* Init */

        BIT_DStream_t bitD1;

        BIT_DStream_t bitD2;

        BIT_DStream_t bitD3;

        BIT_DStream_t bitD4;

        size_t const length1 = MEM_readLE16(istart);

        size_t const length2 = MEM_readLE16(istart+2);

        size_t const length3 = MEM_readLE16(istart+4);

        size_t const length4 = cSrcSize - (length1 + length2 + length3 + 6);

        const BYTE* const istart1 = istart + 6;  /* jumpTable */

        const BYTE* const istart2 = istart1 + length1;

        const BYTE* const istart3 = istart2 + length2;

        const BYTE* const istart4 = istart3 + length3;

        const size_t segmentSize = (dstSize+3) / 4;

        BYTE* const opStart2 = ostart + segmentSize;

        BYTE* const opStart3 = opStart2 + segmentSize;

        BYTE* const opStart4 = opStart3 + segmentSize;

        BYTE* op1 = ostart;

        BYTE* op2 = opStart2;

        BYTE* op3 = opStart3;

        BYTE* op4 = opStart4;

        U32 endSignal = BIT_DStream_unfinished;

        DTableDesc const dtd = HUF_getDTableDesc(DTable);

        U32 const dtLog = dtd.tableLog;

        if (length4 > cSrcSize) return ERROR(corruption_detected);   /* overflow */

        CHECK_F( BIT_initDStream(&bitD1, istart1, length1) );

        CHECK_F( BIT_initDStream(&bitD2, istart2, length2) );

        CHECK_F( BIT_initDStream(&bitD3, istart3, length3) );

        CHECK_F( BIT_initDStream(&bitD4, istart4, length4) );

        /* up to 16 symbols per loop (4 symbols per stream) in 64-bit mode */

        endSignal = BIT_reloadDStream(&bitD1) | BIT_reloadDStream(&bitD2) | BIT_reloadDStream(&bitD3) | BIT_reloadDStream(&bitD4);

        while ( (endSignal==BIT_DStream_unfinished) && (op4<(oend-3)) ) {

            HUF_DECODE_SYMBOLX2_2(op1, &bitD1);

            HUF_DECODE_SYMBOLX2_2(op2, &bitD2);

            HUF_DECODE_SYMBOLX2_2(op3, &bitD3);

            HUF_DECODE_SYMBOLX2_2(op4, &bitD4);

            HUF_DECODE_SYMBOLX2_1(op1, &bitD1);

            HUF_DECODE_SYMBOLX2_1(op2, &bitD2);

            HUF_DECODE_SYMBOLX2_1(op3, &bitD3);

            HUF_DECODE_SYMBOLX2_1(op4, &bitD4);

            HUF_DECODE_SYMBOLX2_2(op1, &bitD1);

            HUF_DECODE_SYMBOLX2_2(op2, &bitD2);

            HUF_DECODE_SYMBOLX2_2(op3, &bitD3);

            HUF_DECODE_SYMBOLX2_2(op4, &bitD4);

            HUF_DECODE_SYMBOLX2_0(op1, &bitD1);

            HUF_DECODE_SYMBOLX2_0(op2, &bitD2);

            HUF_DECODE_SYMBOLX2_0(op3, &bitD3);

            HUF_DECODE_SYMBOLX2_0(op4, &bitD4);

            BIT_reloadDStream(&bitD1);

            BIT_reloadDStream(&bitD2);

            BIT_reloadDStream(&bitD3);

            BIT_reloadDStream(&bitD4);

        }

        /* check corruption */

        /* note : should not be necessary : op# advance in lock step, and we control op4.

         *        but curiously, binary generated by gcc 7.2 & 7.3 with -mbmi2 runs faster when >=1 test is present */

        if (op1 > opStart2) return ERROR(corruption_detected);

        if (op2 > opStart3) return ERROR(corruption_detected);

        if (op3 > opStart4) return ERROR(corruption_detected);

        /* note : op4 supposed already verified within main loop */

        /* finish bitStreams one by one */

        HUF_decodeStreamX2(op1, &bitD1, opStart2, dt, dtLog);

        HUF_decodeStreamX2(op2, &bitD2, opStart3, dt, dtLog);

        HUF_decodeStreamX2(op3, &bitD3, opStart4, dt, dtLog);

        HUF_decodeStreamX2(op4, &bitD4, oend,     dt, dtLog);

        /* check */

        { U32 const endCheck = BIT_endOfDStream(&bitD1) & BIT_endOfDStream(&bitD2) & BIT_endOfDStream(&bitD3) & BIT_endOfDStream(&bitD4);

          if (!endCheck) return ERROR(corruption_detected); }

        /* decoded size */

        return dstSize;

    }

}

FORCE_INLINE_TEMPLATE U32

HUF_decodeSymbolX4(void* op, BIT_DStream_t* DStream, const HUF_DEltX4* dt, const U32 dtLog)

{

    size_t const val = BIT_lookBitsFast(DStream, dtLog);   /* note : dtLog >= 1 */

    memcpy(op, dt+val, 2);

    BIT_skipBits(DStream, dt[val].nbBits);

    return dt[val].length;

}

FORCE_INLINE_TEMPLATE U32

HUF_decodeLastSymbolX4(void* op, BIT_DStream_t* DStream, const HUF_DEltX4* dt, const U32 dtLog)

{

    size_t const val = BIT_lookBitsFast(DStream, dtLog);   /* note : dtLog >= 1 */

    memcpy(op, dt+val, 1);

    if (dt[val].length==1) BIT_skipBits(DStream, dt[val].nbBits);

    else {

        if (DStream->bitsConsumed < (sizeof(DStream->bitContainer)*8)) {

            BIT_skipBits(DStream, dt[val].nbBits);

            if (DStream->bitsConsumed > (sizeof(DStream->bitContainer)*8))

                /* ugly hack; works only because it's the last symbol. Note : can't easily extract nbBits from just this symbol */

                DStream->bitsConsumed = (sizeof(DStream->bitContainer)*8);

    }   }

    return 1;

}

#define HUF_DECODE_SYMBOLX4_0(ptr, DStreamPtr) \

    ptr += HUF_decodeSymbolX4(ptr, DStreamPtr, dt, dtLog)

#define HUF_DECODE_SYMBOLX4_1(ptr, DStreamPtr) \

    if (MEM_64bits() || (HUF_TABLELOG_MAX<=12)) \

        ptr += HUF_decodeSymbolX4(ptr, DStreamPtr, dt, dtLog)

#define HUF_DECODE_SYMBOLX4_2(ptr, DStreamPtr) \

    if (MEM_64bits()) \

        ptr += HUF_decodeSymbolX4(ptr, DStreamPtr, dt, dtLog)

HINT_INLINE size_t

HUF_decodeStreamX4(BYTE* p, BIT_DStream_t* bitDPtr, BYTE* const pEnd,

                const HUF_DEltX4* const dt, const U32 dtLog)

{

    BYTE* const pStart = p;

    /* up to 8 symbols at a time */

    while ((BIT_reloadDStream(bitDPtr) == BIT_DStream_unfinished) & (p < pEnd-(sizeof(bitDPtr->bitContainer)-1))) {

        HUF_DECODE_SYMBOLX4_2(p, bitDPtr);

        HUF_DECODE_SYMBOLX4_1(p, bitDPtr);

        HUF_DECODE_SYMBOLX4_2(p, bitDPtr);

        HUF_DECODE_SYMBOLX4_0(p, bitDPtr);

    }

    /* closer to end : up to 2 symbols at a time */

    while ((BIT_reloadDStream(bitDPtr) == BIT_DStream_unfinished) & (p <= pEnd-2))

        HUF_DECODE_SYMBOLX4_0(p, bitDPtr);

    while (p <= pEnd-2)

        HUF_DECODE_SYMBOLX4_0(p, bitDPtr);   /* no need to reload : reached the end of DStream */

    if (p < pEnd)

        p += HUF_decodeLastSymbolX4(p, bitDPtr, dt, dtLog);

    return p-pStart;

}

FORCE_INLINE_TEMPLATE size_t

HUF_decompress1X4_usingDTable_internal_body(

          void* dst,  size_t dstSize,

    const void* cSrc, size_t cSrcSize,

    const HUF_DTable* DTable)

{

    BIT_DStream_t bitD;

    /* Init */

    CHECK_F( BIT_initDStream(&bitD, cSrc, cSrcSize) );

    /* decode */

    {   BYTE* const ostart = (BYTE*) dst;

        BYTE* const oend = ostart + dstSize;

        const void* const dtPtr = DTable+1;   /* force compiler to not use strict-aliasing */

        const HUF_DEltX4* const dt = (const HUF_DEltX4*)dtPtr;

        DTableDesc const dtd = HUF_getDTableDesc(DTable);

        HUF_decodeStreamX4(ostart, &bitD, oend, dt, dtd.tableLog);

    }

    /* check */

    if (!BIT_endOfDStream(&bitD)) return ERROR(corruption_detected);

    /* decoded size */

    return dstSize;

}

FORCE_INLINE_TEMPLATE size_t

HUF_decompress4X4_usingDTable_internal_body(

          void* dst,  size_t dstSize,

    const void* cSrc, size_t cSrcSize,

    const HUF_DTable* DTable)

{

    if (cSrcSize < 10) return ERROR(corruption_detected);   /* strict minimum : jump table + 1 byte per stream */

    {   const BYTE* const istart = (const BYTE*) cSrc;

        BYTE* const ostart = (BYTE*) dst;

        BYTE* const oend = ostart + dstSize;

        const void* const dtPtr = DTable+1;

        const HUF_DEltX4* const dt = (const HUF_DEltX4*)dtPtr;

        /* Init */

        BIT_DStream_t bitD1;

        BIT_DStream_t bitD2;

        BIT_DStream_t bitD3;

        BIT_DStream_t bitD4;

        size_t const length1 = MEM_readLE16(istart);

        size_t const length2 = MEM_readLE16(istart+2);

        size_t const length3 = MEM_readLE16(istart+4);

        size_t const length4 = cSrcSize - (length1 + length2 + length3 + 6);

        const BYTE* const istart1 = istart + 6;  /* jumpTable */

        const BYTE* const istart2 = istart1 + length1;

        const BYTE* const istart3 = istart2 + length2;

        const BYTE* const istart4 = istart3 + length3;

        size_t const segmentSize = (dstSize+3) / 4;

        BYTE* const opStart2 = ostart + segmentSize;

        BYTE* const opStart3 = opStart2 + segmentSize;

        BYTE* const opStart4 = opStart3 + segmentSize;

        BYTE* op1 = ostart;

        BYTE* op2 = opStart2;

        BYTE* op3 = opStart3;

        BYTE* op4 = opStart4;

        U32 endSignal;

        DTableDesc const dtd = HUF_getDTableDesc(DTable);

        U32 const dtLog = dtd.tableLog;

        if (length4 > cSrcSize) return ERROR(corruption_detected);   /* overflow */

        CHECK_F( BIT_initDStream(&bitD1, istart1, length1) );

        CHECK_F( BIT_initDStream(&bitD2, istart2, length2) );

        CHECK_F( BIT_initDStream(&bitD3, istart3, length3) );

        CHECK_F( BIT_initDStream(&bitD4, istart4, length4) );

        /* 16-32 symbols per loop (4-8 symbols per stream) */

        endSignal = BIT_reloadDStream(&bitD1) | BIT_reloadDStream(&bitD2) | BIT_reloadDStream(&bitD3) | BIT_reloadDStream(&bitD4);

        for ( ; (endSignal==BIT_DStream_unfinished) & (op4<(oend-(sizeof(bitD4.bitContainer)-1))) ; ) {

            HUF_DECODE_SYMBOLX4_2(op1, &bitD1);

            HUF_DECODE_SYMBOLX4_2(op2, &bitD2);

            HUF_DECODE_SYMBOLX4_2(op3, &bitD3);

            HUF_DECODE_SYMBOLX4_2(op4, &bitD4);

            HUF_DECODE_SYMBOLX4_1(op1, &bitD1);

            HUF_DECODE_SYMBOLX4_1(op2, &bitD2);

            HUF_DECODE_SYMBOLX4_1(op3, &bitD3);

            HUF_DECODE_SYMBOLX4_1(op4, &bitD4);

            HUF_DECODE_SYMBOLX4_2(op1, &bitD1);

            HUF_DECODE_SYMBOLX4_2(op2, &bitD2);

            HUF_DECODE_SYMBOLX4_2(op3, &bitD3);

            HUF_DECODE_SYMBOLX4_2(op4, &bitD4);

            HUF_DECODE_SYMBOLX4_0(op1, &bitD1);

            HUF_DECODE_SYMBOLX4_0(op2, &bitD2);

            HUF_DECODE_SYMBOLX4_0(op3, &bitD3);

            HUF_DECODE_SYMBOLX4_0(op4, &bitD4);

            endSignal = BIT_reloadDStream(&bitD1) | BIT_reloadDStream(&bitD2) | BIT_reloadDStream(&bitD3) | BIT_reloadDStream(&bitD4);

        }

        /* check corruption */

        if (op1 > opStart2) return ERROR(corruption_detected);

        if (op2 > opStart3) return ERROR(corruption_detected);

        if (op3 > opStart4) return ERROR(corruption_detected);

        /* note : op4 already verified within main loop */

        /* finish bitStreams one by one */

        HUF_decodeStreamX4(op1, &bitD1, opStart2, dt, dtLog);

        HUF_decodeStreamX4(op2, &bitD2, opStart3, dt, dtLog);

        HUF_decodeStreamX4(op3, &bitD3, opStart4, dt, dtLog);

        HUF_decodeStreamX4(op4, &bitD4, oend,     dt, dtLog);

        /* check */

        { U32 const endCheck = BIT_endOfDStream(&bitD1) & BIT_endOfDStream(&bitD2) & BIT_endOfDStream(&bitD3) & BIT_endOfDStream(&bitD4);

          if (!endCheck) return ERROR(corruption_detected); }

        /* decoded size */

        return dstSize;

    }

}

typedef size_t (*HUF_decompress_usingDTable_t)(void *dst, size_t dstSize,

                                               const void *cSrc,

                                               size_t cSrcSize,

                                               const HUF_DTable *DTable);

#if DYNAMIC_BMI2

#define X(fn)                                                               \

                                                                            \

    static size_t fn##_default(                                             \

                  void* dst,  size_t dstSize,                               \

            const void* cSrc, size_t cSrcSize,                              \

            const HUF_DTable* DTable)                                       \

    {                                                                       \

        return fn##_body(dst, dstSize, cSrc, cSrcSize, DTable);             \

    }                                                                       \

                                                                            \

    static TARGET_ATTRIBUTE("bmi2") size_t fn##_bmi2(                       \

                  void* dst,  size_t dstSize,                               \