be/src/util/simd/bits.h

Source
// Licensed to the Apache Software Foundation (ASF) under one
// or more contributor license agreements.  See the NOTICE file
// distributed with this work for additional information
// regarding copyright ownership.  The ASF licenses this file
// to you under the Apache License, Version 2.0 (the
// "License"); you may not use this file except in compliance
// with the License.  You may obtain a copy of the License at
//
//   http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing,
// software distributed under the License is distributed on an
// "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
// KIND, either express or implied.  See the License for the
// specific language governing permissions and limitations
// under the License.

#pragma once

#include <cstdint>
#include <cstring>
#include <type_traits>
#include <vector>

#if defined(__ARM_NEON)
#include <arm_neon.h>
#endif

#include "util/sse_util.hpp"

namespace doris::simd {
consteval auto bits_mask_length() {
#if defined(__ARM_NEON) && defined(__aarch64__)
    return 16;
#else
    return 32;
#endif
}

#if defined(__ARM_NEON) && defined(__aarch64__)
inline uint64_t get_nibble_mask(uint8x16_t values) {
    // It produces 4-bit out of each byte, alternating between the high 4-bits and low 4-bits of the 16-byte vector.
    // Given that the comparison operators give a 16-byte result of 0x00 or 0xff, the result is close to being a PMOVMSKB,
    // the only difference is that every matching bit is repeated 4 times and is a 64-bit integer.
    // https://community.arm.com/arm-community-blogs/b/infrastructure-solutions-blog/posts/porting-x86-vector-bitmask-optimizations-to-arm-neon?CommentId=af187ac6-ae00-4e4d-bbf0-e142187aa92e
    return vget_lane_u64(vreinterpret_u64_u8(vshrn_n_u16(vreinterpretq_u16_u8(values), 4)), 0);
}
/*
Input 16 bytes of data and convert it into a 64-bit integer, where one bit appears 4 times.
Compare with bytes32_mask_to_bits32_mask, a u8 array with a length of 32
  std::vector<uint8_t> vec = {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1,
                                0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0};

bytes32_mask_to_bits32_mask   0100 0000 0000 0000,1101 0000 0000 0011


                            (1101 0000 0000 0011)
bytes16_mask_to_bits64_mask   1111 1111 0000 1111,0000 0000 0000 0000,0000 0000 0000 0000,0000 0000 1111 1111
                            (0100 0000 0000 0000)
                              0000 1111 0000 0000,0000 0000 0000 0000,0000 0000 0000 0000,0000 0000 0000 0000
*/

inline uint64_t bytes16_mask_to_bits64_mask(const uint8_t* data) {
    const uint8x16_t vfilter = vld1q_u8(data);
    return get_nibble_mask(vmvnq_u8(vceqzq_u8(vfilter)));
}
#endif

inline uint32_t bytes32_mask_to_bits32_mask(const uint8_t* data) {
#ifdef __AVX2__
    auto zero32 = _mm256_setzero_si256();
    auto mask = static_cast<uint32_t>(_mm256_movemask_epi8(
            _mm256_cmpgt_epi8(_mm256_loadu_si256(reinterpret_cast<const __m256i*>(data)), zero32)));
#elif defined(__SSE2__)
    auto zero16 = _mm_setzero_si128();
    uint32_t mask =
            (static_cast<uint32_t>(_mm_movemask_epi8(_mm_cmpgt_epi8(
                    _mm_loadu_si128(reinterpret_cast<const __m128i*>(data)), zero16)))) |
            ((static_cast<uint32_t>(_mm_movemask_epi8(_mm_cmpgt_epi8(
                      _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + 16)), zero16)))
              << 16) &
             0xffff0000);
#else
    uint32_t mask = 0;
    for (std::size_t i = 0; i < 32; ++i) {
        mask |= static_cast<uint32_t>(1 == *(data + i)) << i;
    }
#endif
    return mask;
}

inline auto bytes_mask_to_bits_mask(const uint8_t* data) {
#if defined(__ARM_NEON) && defined(__aarch64__)
    return bytes16_mask_to_bits64_mask(data);
#else
    return bytes32_mask_to_bits32_mask(data);
#endif
}

inline constexpr auto bits_mask_all() {
#if defined(__ARM_NEON) && defined(__aarch64__)
    return 0xffff'ffff'ffff'ffffULL;
#else
    return 0xffffffff;
#endif
}

template <typename Func>
void iterate_through_bits_mask(Func func, decltype(bytes_mask_to_bits_mask(nullptr)) mask) {
#if defined(__ARM_NEON) && defined(__aarch64__)
    mask &= 0x8888'8888'8888'8888ULL;
    while (mask) {
        const auto index = __builtin_ctzll(mask) >> 2;
        func(index);
        mask &= mask - 1;
    }

#else
    while (mask) {
        const auto bit_pos = __builtin_ctzll(mask);
        func(bit_pos);
        mask = mask & (mask - 1);
    }
#endif
}

struct FilterCounts {
    size_t input_rows = 0;
    size_t output_rows = 0;
};

// Merge a filter in-place and count nonzero bytes before and after the merge.
// Nullable filters preserve the scalar expression data[i] &= !null_map[i] & filter[i].
template <bool nullable = false>
inline FilterCounts filter_and_count(uint8_t* __restrict data, const uint8_t* __restrict filter,
                                     size_t size, const uint8_t* __restrict null_map = nullptr) {
    FilterCounts counts;
    size_t i = 0;
#if defined(__SSE2__) || defined(__aarch64__)
    const __m128i zero = _mm_setzero_si128();
    __m128i input_sum = zero;
    __m128i output_sum = zero;
    const size_t vector_size = size / 64 * 64;
    for (; i < vector_size; i += 64) {
        __m128i input_zeros = zero;
        __m128i output_zeros = zero;
        for (size_t offset = 0; offset < 64; offset += 16) {
            const auto before =
                    _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + i + offset));
            const auto filter_values =
                    _mm_loadu_si128(reinterpret_cast<const __m128i*>(filter + i + offset));
            auto after = _mm_and_si128(before, filter_values);
            if constexpr (nullable) {
                const auto nulls =
                        _mm_loadu_si128(reinterpret_cast<const __m128i*>(null_map + i + offset));
                const auto valid = _mm_and_si128(_mm_cmpeq_epi8(nulls, zero), _mm_set1_epi8(1));
                after = _mm_and_si128(after, valid);
            }
            // Each byte accumulates at most four zeros. Reduce once per 64 rows
            // without expanding bytes into size_t lanes or emulating movemask on ARM.
            input_zeros = _mm_sub_epi8(input_zeros, _mm_cmpeq_epi8(before, zero));
            output_zeros = _mm_sub_epi8(output_zeros, _mm_cmpeq_epi8(after, zero));
            _mm_storeu_si128(reinterpret_cast<__m128i*>(data + i + offset), after);
        }
        input_sum = _mm_add_epi64(input_sum, _mm_sad_epu8(input_zeros, zero));
        output_sum = _mm_add_epi64(output_sum, _mm_sad_epu8(output_zeros, zero));
    }
    counts.input_rows = vector_size - static_cast<size_t>(_mm_cvtsi128_si64(_mm_add_epi64(
                                              input_sum, _mm_srli_si128(input_sum, 8))));
    counts.output_rows = vector_size - static_cast<size_t>(_mm_cvtsi128_si64(_mm_add_epi64(
                                               output_sum, _mm_srli_si128(output_sum, 8))));
#endif
    for (; i < size; ++i) {
        const uint8_t before = data[i];
        uint8_t after = before & filter[i];
        if constexpr (nullable) {
            after &= !null_map[i];
        }
        counts.input_rows += before != 0;
        counts.output_rows += after != 0;
        data[i] = after;
    }
    return counts;
}

template <typename T>
    requires requires { std::is_unsigned_v<T>; }
inline T count_zero_num(const int8_t* __restrict data, T size) {
    T num = 0;
    const int8_t* end = data + size;
#if defined(__ARM_NEON)
    const int8_t* end64 = data + (size / 64 * 64);

    for (; data < end64; data += 64) {
        auto a0 = vshrq_n_u8(vceqzq_s8(vld1q_s8(data)), 7);
        auto a1 = vshrq_n_u8(vceqzq_s8(vld1q_s8(data + 16)), 7);
        auto a2 = vshrq_n_u8(vceqzq_s8(vld1q_s8(data + 32)), 7);
        auto a3 = vshrq_n_u8(vceqzq_s8(vld1q_s8(data + 48)), 7);

        auto s0 = vaddq_u8(a0, a1);
        auto s1 = vaddq_u8(a2, a3);
        auto s = vaddq_u8(s0, s1);
        num += vaddvq_u8(s);
    }
#elif defined(__SSE2__) && defined(__POPCNT__)
    const __m128i zero16 = _mm_setzero_si128();
    const int8_t* end64 = data + (size / 64 * 64);

    for (; data < end64; data += 64) {
        num += __builtin_popcountll(
                static_cast<uint64_t>(_mm_movemask_epi8(_mm_cmpeq_epi8(
                        _mm_loadu_si128(reinterpret_cast<const __m128i*>(data)), zero16))) |
                (static_cast<uint64_t>(_mm_movemask_epi8(_mm_cmpeq_epi8(
                         _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + 16)), zero16)))
                 << 16U) |
                (static_cast<uint64_t>(_mm_movemask_epi8(_mm_cmpeq_epi8(
                         _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + 32)), zero16)))
                 << 32U) |
                (static_cast<uint64_t>(_mm_movemask_epi8(_mm_cmpeq_epi8(
                         _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + 48)), zero16)))
                 << 48U));
    }
#endif
    for (; data < end; ++data) {
        num += (*data == 0);
    }
    return num;
}

template <typename T>
    requires requires { std::is_unsigned_v<T>; }
inline T count_zero_num(const int8_t* __restrict data, const uint8_t* __restrict null_map, T size) {
    T num = 0;
    const int8_t* end = data + size;
#if defined(__ARM_NEON)
    const int8_t* end64 = data + (size / 64 * 64);

    for (; data < end64; data += 64, null_map += 64) {
        auto a0 = vshrq_n_u8(vceqzq_s8(vld1q_s8(data)), 7);
        auto a1 = vshrq_n_u8(vceqzq_s8(vld1q_s8(data + 16)), 7);
        auto a2 = vshrq_n_u8(vceqzq_s8(vld1q_s8(data + 32)), 7);
        auto a3 = vshrq_n_u8(vceqzq_s8(vld1q_s8(data + 48)), 7);

        auto r0 = vorrq_u8(a0, vld1q_u8(null_map));
        auto r1 = vorrq_u8(a1, vld1q_u8(null_map + 16));
        auto r2 = vorrq_u8(a2, vld1q_u8(null_map + 32));
        auto r3 = vorrq_u8(a3, vld1q_u8(null_map + 48));

        auto s0 = vaddq_u8(r0, r1);
        auto s1 = vaddq_u8(r2, r3);
        auto s = vaddq_u8(s0, s1);
        num += vaddvq_u8(s);
    }
#elif defined(__SSE2__) && defined(__POPCNT__)
    const __m128i zero16 = _mm_setzero_si128();
    const __m128i one16 = _mm_set1_epi8(1);
    const int8_t* end64 = data + (size / 64 * 64);

    for (; data < end64; data += 64, null_map += 64) {
        num += __builtin_popcountll(
                static_cast<uint64_t>(_mm_movemask_epi8(_mm_or_si128(
                        _mm_cmpeq_epi8(_mm_loadu_si128(reinterpret_cast<const __m128i*>(data)),
                                       zero16),
                        _mm_cmpeq_epi8(_mm_loadu_si128(reinterpret_cast<const __m128i*>(null_map)),
                                       one16)))) |
                (static_cast<uint64_t>(_mm_movemask_epi8(_mm_or_si128(
                         _mm_cmpeq_epi8(
                                 _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + 16)),
                                 zero16),
                         _mm_cmpeq_epi8(
                                 _mm_loadu_si128(reinterpret_cast<const __m128i*>(null_map + 16)),
                                 one16))))
                 << 16U) |
                (static_cast<uint64_t>(_mm_movemask_epi8(_mm_or_si128(
                         _mm_cmpeq_epi8(
                                 _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + 32)),
                                 zero16),
                         _mm_cmpeq_epi8(
                                 _mm_loadu_si128(reinterpret_cast<const __m128i*>(null_map + 32)),
                                 one16))))
                 << 32U) |
                (static_cast<uint64_t>(_mm_movemask_epi8(_mm_or_si128(
                        _mm_cmpeq_epi8(_mm_loadu_si128(reinterpret_cast<const __m128i*>(data + 48)),
                                       zero16),
                        _mm_cmpeq_epi8(
                                _mm_loadu_si128(reinterpret_cast<const __m128i*>(null_map + 48)),
                                one16)))))
                        << 48U);
    }
#endif
    for (; data < end; ++data, ++null_map) {
        num += ((*data == 0) | *null_map);
    }
    return num;
}

// TODO: compare with different SIMD implements
template <class T>
static size_t find_byte(const std::vector<T>& vec, size_t start, T byte) {
    if (start >= vec.size()) {
        return start;
    }
    const void* p = std::memchr((const void*)(vec.data() + start), byte, vec.size() - start);
    if (p == nullptr) {
        return vec.size();
    }
    return (T*)p - vec.data();
}

template <class T>
static size_t find_byte(const T* data, size_t start, size_t end, T byte) {
    if (start >= end) {
        return start;
    }
    const void* p = std::memchr((const void*)(data + start), byte, end - start);
    if (p == nullptr) {
        return end;
    }
    return (T*)p - data;
}

inline size_t find_one(const std::vector<uint8_t>& vec, size_t start) {
    return find_byte<uint8_t>(vec, start, 1);
}

inline size_t find_one(const uint8_t* data, size_t start, size_t end) {
    return find_byte<uint8_t>(data, start, end, 1);
}

inline size_t find_zero(const std::vector<uint8_t>& vec, size_t start) {
    return find_byte<uint8_t>(vec, start, 0);
}

inline bool contain_one(const uint8_t* __restrict data, size_t size) {
    size_t i = 0;
#if defined(__AVX2__)
    for (; i + 32 <= size; i += 32) {
        __m256i chunk = _mm256_loadu_si256(reinterpret_cast<const __m256i*>(data + i));
        if (!_mm256_testz_si256(chunk, chunk)) {
            return true;
        }
    }
#elif defined(__SSE2__)
    const __m128i zero = _mm_setzero_si128();
    for (; i + 16 <= size; i += 16) {
        __m128i chunk = _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + i));
        if (_mm_movemask_epi8(_mm_cmpeq_epi8(chunk, zero)) != 0xFFFF) {
            return true;
        }
    }
#elif defined(__ARM_NEON) && defined(__aarch64__)
    for (; i + 16 <= size; i += 16) {
        if (vmaxvq_u8(vld1q_u8(data + i)) != 0) {
            return true;
        }
    }
#endif
    for (; i < size; ++i) {
        if (data[i]) {
            return true;
        }
    }
    return false;
}

inline bool contain_zero(const uint8_t* __restrict data, size_t size) {
    size_t i = 0;
#if defined(__AVX2__)
    const __m256i zero = _mm256_setzero_si256();
    for (; i + 32 <= size; i += 32) {
        __m256i chunk = _mm256_loadu_si256(reinterpret_cast<const __m256i*>(data + i));
        if (_mm256_movemask_epi8(_mm256_cmpeq_epi8(chunk, zero)) != 0) {
            return true;
        }
    }
#elif defined(__SSE2__)
    const __m128i zero = _mm_setzero_si128();
    for (; i + 16 <= size; i += 16) {
        __m128i chunk = _mm_loadu_si128(reinterpret_cast<const __m128i*>(data + i));
        if (_mm_movemask_epi8(_mm_cmpeq_epi8(chunk, zero)) != 0) {
            return true;
        }
    }
#elif defined(__ARM_NEON) && defined(__aarch64__)
    for (; i + 16 <= size; i += 16) {
        if (vminvq_u8(vld1q_u8(data + i)) == 0) {
            return true;
        }
    }
#endif
    for (; i < size; ++i) {
        if (!data[i]) {
            return true;
        }
    }
    return false;
}

} // namespace doris::simd

Coverage Report

Created: 2026-11-06 10:43