#include "crypto_core_wforcesntrup1013.h"
#include "crypto_decode_int16.h"
#include "params.h"
#include <immintrin.h>

#define int16 int16_t


/* out = in if bottom bits of in have weight w */
/* otherwise out = (1,1,...,1,0,0,...,0) */
int PQCLEAN_SNTRUP1013_AVX2_crypto_core_wforcesntrup1013(unsigned char *out, const unsigned char *in) {
    int16 weight;
    int16 mask;
    __m256i maskvec;
    int i;

    crypto_core_weight((unsigned char *) &weight, in);
    PQCLEAN_SNTRUP1013_AVX2_crypto_decode_int16(&weight, (unsigned char *) &weight);

    mask = (int16) ((weight - w) | (w - weight));
    mask >>= 15;
    maskvec = _mm256_set1_epi16((short) ~mask);

    i = w - 32;
    for (;;) {
        do {
            __m256i x = _mm256_loadu_si256((__m256i *) in);
            x ^= _mm256_set1_epi8(1);
            x &= maskvec;
            x ^= _mm256_set1_epi8(1);
            _mm256_storeu_si256((__m256i *) out, x);
            in += 32;
            out += 32;
            i -= 32;
        } while (i >= 0);
        if (i <= -32) {
            break;
        }
        in += i;
        out += i;
    }

    i = p - w - 32;
    for (;;) {
        do {
            __m256i x = _mm256_loadu_si256((__m256i *) in);
            x &= maskvec;
            _mm256_storeu_si256((__m256i *) out, x);
            in += 32;
            out += 32;
            i -= 32;
        } while (i >= 0);
        if (i <= -32) {
            break;
        }
        in += i;
        out += i;
    }

    return 0;
}