mirror of
https://github.com/qdrant/qdrant.git
synced 2026-07-25 20:21:08 -05:00
b414a402bb
* # This is a combination of 7 commits. * SameAsStorage default value * fix coderabbit warnings * neon for u8 bq * sse for u8 bq * fix windows build * rename function * add comments * fmt * fix arm build * review remarks
206 lines
7.2 KiB
C
206 lines
7.2 KiB
C
#include <stdlib.h>
|
|
#include <stdint.h>
|
|
#include <immintrin.h>
|
|
|
|
#include "export_macro.h"
|
|
|
|
#define HSUM256_PS(X, R) \
|
|
float R = 0.0f; \
|
|
{ \
|
|
__m128 x128 = _mm_add_ps(_mm256_extractf128_ps(X, 1), _mm256_castps256_ps128(X)); \
|
|
__m128 x64 = _mm_add_ps(x128, _mm_movehl_ps(x128, x128)); \
|
|
__m128 x32 = _mm_add_ss(x64, _mm_shuffle_ps(x64, x64, 0x55)); \
|
|
R = _mm_cvtss_f32(x32); \
|
|
}
|
|
|
|
#define HSUM256_EPI32(X, R) \
|
|
int R = 0; \
|
|
{ \
|
|
__m128i x128 = _mm_add_epi32(_mm256_extractf128_si256(X, 1), _mm256_castsi256_si128(X)); \
|
|
__m128i x64 = _mm_add_epi32(x128, _mm_srli_si128(x128, 8)); \
|
|
__m128i x32 = _mm_add_epi32(x64, _mm_srli_si128(x64, 4)); \
|
|
R = _mm_cvtsi128_si32(x32); \
|
|
}
|
|
|
|
EXPORT float impl_score_dot_avx(
|
|
const uint8_t* query_ptr,
|
|
const uint8_t* vector_ptr,
|
|
uint32_t dim
|
|
) {
|
|
const __m256i* v_ptr = (const __m256i*)vector_ptr;
|
|
const __m256i* q_ptr = (const __m256i*)query_ptr;
|
|
|
|
__m256i mul1 = _mm256_setzero_si256();
|
|
__m256i mask_epu32 = _mm256_set1_epi32(0xFFFF);
|
|
for (uint32_t _i = 0; _i < dim / 32; _i++) {
|
|
__m256i v = _mm256_loadu_si256(v_ptr);
|
|
__m256i q = _mm256_loadu_si256(q_ptr);
|
|
v_ptr++;
|
|
q_ptr++;
|
|
|
|
__m256i s = _mm256_maddubs_epi16(v, q);
|
|
__m256i s_low = _mm256_cvtepi16_epi32(_mm256_castsi256_si128(s));
|
|
__m256i s_high = _mm256_cvtepi16_epi32(_mm256_extractf128_si256(s, 1));
|
|
mul1 = _mm256_add_epi32(mul1, s_low);
|
|
mul1 = _mm256_add_epi32(mul1, s_high);
|
|
}
|
|
|
|
// the vector sizes are assumed to be multiples of 16, check if one last 16-element part remaining
|
|
if (dim % 32 != 0) {
|
|
__m128i v_short = _mm_loadu_si128((const __m128i*)v_ptr);
|
|
__m128i q_short = _mm_loadu_si128((const __m128i*)q_ptr);
|
|
|
|
__m256i v1 = _mm256_cvtepu8_epi16(v_short);
|
|
__m256i q1 = _mm256_cvtepu8_epi16(q_short);
|
|
|
|
__m256i s = _mm256_mullo_epi16(v1, q1);
|
|
mul1 = _mm256_add_epi32(mul1, _mm256_and_si256(s, mask_epu32));
|
|
mul1 = _mm256_add_epi32(mul1, _mm256_srli_epi32(s, 16));
|
|
}
|
|
__m256 mul_ps = _mm256_cvtepi32_ps(mul1);
|
|
HSUM256_PS(mul_ps, mul_scalar);
|
|
return mul_scalar;
|
|
}
|
|
|
|
EXPORT float impl_score_l1_avx(
|
|
const uint8_t* query_ptr,
|
|
const uint8_t* vector_ptr,
|
|
uint32_t dim
|
|
) {
|
|
const __m256i* v_ptr = (const __m256i*)vector_ptr;
|
|
const __m256i* q_ptr = (const __m256i*)query_ptr;
|
|
|
|
uint32_t m = dim - (dim % 32);
|
|
__m256i sum256 = _mm256_setzero_si256();
|
|
|
|
for (uint32_t i = 0; i < m; i += 32) {
|
|
__m256i v = _mm256_loadu_si256(v_ptr);
|
|
__m256i q = _mm256_loadu_si256(q_ptr);
|
|
v_ptr++;
|
|
q_ptr++;
|
|
|
|
// Compute the difference in both directions and take the maximum for abs
|
|
__m256i diff1 = _mm256_subs_epu8(v, q);
|
|
__m256i diff2 = _mm256_subs_epu8(q, v);
|
|
|
|
__m256i abs_diff = _mm256_max_epu8(diff1, diff2);
|
|
|
|
__m256i abs_diff16_lo = _mm256_unpacklo_epi8(abs_diff, _mm256_setzero_si256());
|
|
__m256i abs_diff16_hi = _mm256_unpackhi_epi8(abs_diff, _mm256_setzero_si256());
|
|
|
|
sum256 = _mm256_add_epi16(sum256, abs_diff16_lo);
|
|
sum256 = _mm256_add_epi16(sum256, abs_diff16_hi);
|
|
}
|
|
|
|
// the vector sizes are assumed to be multiples of 16, check if one last 16-element part remaining
|
|
if (m < dim) {
|
|
__m128i v_short = _mm_loadu_si128((const __m128i * ) v_ptr);
|
|
__m128i q_short = _mm_loadu_si128((const __m128i * ) q_ptr);
|
|
|
|
__m128i diff1 = _mm_subs_epu8(v_short, q_short);
|
|
__m128i diff2 = _mm_subs_epu8(q_short, v_short);
|
|
|
|
__m128i abs_diff = _mm_max_epu8(diff1, diff2);
|
|
|
|
__m128i abs_diff16_lo_128 = _mm_unpacklo_epi8(abs_diff, _mm_setzero_si128());
|
|
__m128i abs_diff16_hi_128 = _mm_unpackhi_epi8(abs_diff, _mm_setzero_si128());
|
|
|
|
__m256i abs_diff16_lo = _mm256_cvtepu16_epi32(abs_diff16_lo_128);
|
|
__m256i abs_diff16_hi = _mm256_cvtepu16_epi32(abs_diff16_hi_128);
|
|
|
|
sum256 = _mm256_add_epi16(sum256, abs_diff16_lo);
|
|
sum256 = _mm256_add_epi16(sum256, abs_diff16_hi);
|
|
}
|
|
|
|
__m256i sum_epi32 = _mm256_add_epi32(
|
|
_mm256_unpacklo_epi16(sum256, _mm256_setzero_si256()),
|
|
_mm256_unpackhi_epi16(sum256, _mm256_setzero_si256()));
|
|
|
|
HSUM256_EPI32(sum_epi32, sum);
|
|
|
|
return (float) sum;
|
|
}
|
|
|
|
EXPORT uint32_t impl_xor_popcnt_scalar8_avx_uint128(
|
|
const uint8_t* query_ptr,
|
|
const uint8_t* vector_ptr,
|
|
uint32_t count
|
|
) {
|
|
const uint64_t* v_ptr = (const uint64_t*)vector_ptr;
|
|
const uint64_t* q_ptr = (const uint64_t*)query_ptr;
|
|
|
|
__m256i sum1 = _mm256_set1_epi32(0);
|
|
__m256i sum2 = _mm256_set1_epi32(0);
|
|
for (uint32_t _i = 0; _i < count; _i++) {
|
|
uint64_t v_1 = *v_ptr;
|
|
uint64_t v_2 = *(v_ptr + 1);
|
|
|
|
__m256i popcnt1 = _mm256_set_epi32(
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 0)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 2)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 4)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 6)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 1)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 3)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 5)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 7))
|
|
);
|
|
sum1 = _mm256_add_epi32(sum1, popcnt1);
|
|
|
|
__m256i popcnt2 = _mm256_set_epi32(
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 8)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 10)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 12)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 14)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 9)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 11)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 13)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 15))
|
|
);
|
|
sum2 = _mm256_add_epi32(sum2, popcnt2);
|
|
|
|
v_ptr += 2;
|
|
q_ptr += 16;
|
|
}
|
|
__m256i factor1 = _mm256_set_epi32(1, 2, 4, 8, 1, 2, 4, 8);
|
|
__m256i factor2 = _mm256_set_epi32(16, 32, 64, 128, 16, 32, 64, 128);
|
|
__m256 result1_mm256 = _mm256_cvtepi32_ps(_mm256_mullo_epi32(sum1, factor1));
|
|
__m256 result2_mm256 = _mm256_cvtepi32_ps(_mm256_mullo_epi32(sum2, factor2));
|
|
HSUM256_PS(_mm256_add_ps(result1_mm256, result2_mm256), mul_scalar);
|
|
return (uint32_t)mul_scalar;
|
|
}
|
|
|
|
EXPORT uint32_t impl_xor_popcnt_scalar4_avx_uint128(
|
|
const uint8_t* query_ptr,
|
|
const uint8_t* vector_ptr,
|
|
uint32_t count
|
|
) {
|
|
const uint64_t* v_ptr = (const uint64_t*)vector_ptr;
|
|
const uint64_t* q_ptr = (const uint64_t*)query_ptr;
|
|
|
|
__m256i sum = _mm256_set1_epi32(0);
|
|
for (uint32_t _i = 0; _i < count; _i++) {
|
|
uint64_t v_1 = *v_ptr;
|
|
uint64_t v_2 = *(v_ptr + 1);
|
|
|
|
__m256i popcnt = _mm256_set_epi32(
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 0)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 2)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 4)),
|
|
_mm_popcnt_u64(v_1 ^ *(q_ptr + 6)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 1)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 3)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 5)),
|
|
_mm_popcnt_u64(v_2 ^ *(q_ptr + 7))
|
|
);
|
|
sum = _mm256_add_epi32(sum, popcnt);
|
|
|
|
v_ptr += 2;
|
|
q_ptr += 8;
|
|
}
|
|
__m256i factor = _mm256_set_epi32(1, 2, 4, 8, 1, 2, 4, 8);
|
|
__m256 result_mm256 = _mm256_cvtepi32_ps(_mm256_mullo_epi32(sum, factor));
|
|
HSUM256_PS(result_mm256, mul_scalar);
|
|
return (uint32_t)mul_scalar;
|
|
}
|