static float InnerProductSIMD16Ext()

in cpp/src/Spaces/InnerProduct.h [190:217]

• 20 lines of code
• 2 McCabe index (conditional complexity)

static float InnerProductSIMD16Ext(const float *pVect1, const float *pVect2,
                                   const size_t qty) {
  float PORTABLE_ALIGN64 TmpRes[16];

  size_t qty16 = qty / 16;

  const float *pEnd1 = pVect1 + 16 * qty16;

  __m512 sum512 = _mm512_set1_ps(0);

  while (pVect1 < pEnd1) {
    //_mm_prefetch((char*)(pVect2 + 16), _MM_HINT_T0);

    __m512 v1 = _mm512_loadu_ps(pVect1);
    pVect1 += 16;
    __m512 v2 = _mm512_loadu_ps(pVect2);
    pVect2 += 16;
    sum512 = _mm512_add_ps(sum512, _mm512_mul_ps(v1, v2));
  }

  _mm512_store_ps(TmpRes, sum512);
  float sum = TmpRes[0] + TmpRes[1] + TmpRes[2] + TmpRes[3] + TmpRes[4] +
              TmpRes[5] + TmpRes[6] + TmpRes[7] + TmpRes[8] + TmpRes[9] +
              TmpRes[10] + TmpRes[11] + TmpRes[12] + TmpRes[13] + TmpRes[14] +
              TmpRes[15];

  return 1.0f - sum;
}