// Auto-generated file. Do not edit! // Template: src/qs8-gemm/MRx16c4-avx512vnni.c.in // Generator: tools/xngen // // Copyright 2023 Google LLC // // This source code is licensed under the BSD-style license found in the // LICENSE file in the root directory of this source tree. #include #include #include "xnnpack/common.h" #include "xnnpack/gemm.h" #include "xnnpack/intrinsics-polyfill.h" #include "xnnpack/math.h" #include "xnnpack/unaligned.h" void xnn_qd8_f32_qc8w_gemm_minmax_ukernel_10x16c4__avx512vnni( size_t mr, size_t nc, size_t kc, const int8_t* restrict a, size_t a_stride, const void* restrict w, float* restrict c, size_t cm_stride, size_t cn_stride, const union xnn_f32_minmax_params params[restrict XNN_MIN_ELEMENTS(1)], const struct xnn_qd8_quantization_params quantization_params[restrict XNN_MIN_ELEMENTS(1)]) XNN_OOB_READS { assert(mr != 0); assert(mr <= 10); assert(nc != 0); assert(kc != 0); assert(kc % sizeof(int8_t) == 0); assert(a != NULL); assert(w != NULL); assert(c != NULL); kc = round_up_po2(kc, 4 * sizeof(int8_t)); const int8_t* a0 = a; float* c0 = c; const int8_t* a1 = (const int8_t*) ((uintptr_t) a0 + a_stride); float* c1 = (float*) ((uintptr_t) c0 + cm_stride); if XNN_UNPREDICTABLE(mr < 2) { a1 = a0; c1 = c0; } const int8_t* a2 = (const int8_t*) ((uintptr_t) a1 + a_stride); float* c2 = (float*) ((uintptr_t) c1 + cm_stride); if XNN_UNPREDICTABLE(mr <= 2) { a2 = a1; c2 = c1; } const int8_t* a3 = (const int8_t*) ((uintptr_t) a2 + a_stride); float* c3 = (float*) ((uintptr_t) c2 + cm_stride); if XNN_UNPREDICTABLE(mr < 4) { a3 = a2; c3 = c2; } const int8_t* a4 = (const int8_t*) ((uintptr_t) a3 + a_stride); float* c4 = (float*) ((uintptr_t) c3 + cm_stride); if XNN_UNPREDICTABLE(mr <= 4) { a4 = a3; c4 = c3; } const int8_t* a5 = (const int8_t*) ((uintptr_t) a4 + a_stride); float* c5 = (float*) ((uintptr_t) c4 + cm_stride); if XNN_UNPREDICTABLE(mr < 6) { a5 = a4; c5 = c4; } const int8_t* a6 = (const int8_t*) ((uintptr_t) a5 + a_stride); float* c6 = (float*) ((uintptr_t) c5 + cm_stride); if XNN_UNPREDICTABLE(mr <= 6) { a6 = a5; c6 = c5; } const int8_t* a7 = (const int8_t*) ((uintptr_t) a6 + a_stride); float* c7 = (float*) ((uintptr_t) c6 + cm_stride); if XNN_UNPREDICTABLE(mr < 8) { a7 = a6; c7 = c6; } const int8_t* a8 = (const int8_t*) ((uintptr_t) a7 + a_stride); float* c8 = (float*) ((uintptr_t) c7 + cm_stride); if XNN_UNPREDICTABLE(mr <= 8) { a8 = a7; c8 = c7; } const int8_t* a9 = (const int8_t*) ((uintptr_t) a8 + a_stride); float* c9 = (float*) ((uintptr_t) c8 + cm_stride); if XNN_UNPREDICTABLE(mr != 10) { a9 = a8; c9 = c8; } const __m512i vinput_zero_point0 = _mm512_set1_epi32((int) quantization_params[0].zero_point); const __m512i vinput_zero_point1 = _mm512_set1_epi32((int) quantization_params[1].zero_point); const __m512i vinput_zero_point2 = _mm512_set1_epi32((int) quantization_params[2].zero_point); const __m512i vinput_zero_point3 = _mm512_set1_epi32((int) quantization_params[3].zero_point); const __m512i vinput_zero_point4 = _mm512_set1_epi32((int) quantization_params[4].zero_point); const __m512i vinput_zero_point5 = _mm512_set1_epi32((int) quantization_params[5].zero_point); const __m512i vinput_zero_point6 = _mm512_set1_epi32((int) quantization_params[6].zero_point); const __m512i vinput_zero_point7 = _mm512_set1_epi32((int) quantization_params[7].zero_point); const __m512i vinput_zero_point8 = _mm512_set1_epi32((int) quantization_params[8].zero_point); const __m512i vinput_zero_point9 = _mm512_set1_epi32((int) quantization_params[9].zero_point); const __m512 voutput_min = _mm512_set1_ps(params->scalar.min); const __m512 voutput_max = _mm512_set1_ps(params->scalar.max); // XNN_FORCE_REALIZATION(voutput_min); // XNN_FORCE_REALIZATION(voutput_max); do { const __m512i vksum0123456789ABCDEF = _mm512_load_epi32(w); __m512i vacc0x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point0); __m512i vacc1x0x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc1x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point1); __m512i vacc1x1x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc2x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point2); __m512i vacc1x2x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc3x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point3); __m512i vacc1x3x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc4x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point4); __m512i vacc1x4x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc5x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point5); __m512i vacc1x5x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc6x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point6); __m512i vacc1x6x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc7x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point7); __m512i vacc1x7x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc8x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point8); __m512i vacc1x8x0123456789ABCDEF = _mm512_setzero_epi32(); __m512i vacc9x0123456789ABCDEF = _mm512_mullo_epi32(vksum0123456789ABCDEF, vinput_zero_point9); __m512i vacc1x9x0123456789ABCDEF = _mm512_setzero_epi32(); w = (const int32_t*) w + 16; size_t k = kc; while (k >= 8 * sizeof(int8_t)) { const __m512i va0x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a0)); const __m512i va0x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a0 + 4)); a0 += 8; const __m512i va1x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a1)); const __m512i va1x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a1 + 4)); a1 += 8; const __m512i va2x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a2)); const __m512i va2x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a2 + 4)); a2 += 8; const __m512i va3x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a3)); const __m512i va3x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a3 + 4)); a3 += 8; const __m512i va4x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a4)); const __m512i va4x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a4 + 4)); a4 += 8; const __m512i va5x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a5)); const __m512i va5x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a5 + 4)); a5 += 8; const __m512i va6x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a6)); const __m512i va6x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a6 + 4)); a6 += 8; const __m512i va7x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a7)); const __m512i va7x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a7 + 4)); a7 += 8; const __m512i va8x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a8)); const __m512i va8x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a8 + 4)); a8 += 8; const __m512i va9x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a9)); const __m512i va9x4567 = _mm512_set1_epi32((int) unaligned_load_u32(a9 + 4)); a9 += 8; const __m512i vb0123456789ABCDEFx0123 = _mm512_load_si512(w); const __m512i vb0123456789ABCDEFx4567 = _mm512_load_si512((const int8_t*) w + 64); vacc0x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc0x0123456789ABCDEF, va0x0123, vb0123456789ABCDEFx0123); vacc1x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x0123456789ABCDEF, va1x0123, vb0123456789ABCDEFx0123); vacc2x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc2x0123456789ABCDEF, va2x0123, vb0123456789ABCDEFx0123); vacc3x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc3x0123456789ABCDEF, va3x0123, vb0123456789ABCDEFx0123); vacc4x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc4x0123456789ABCDEF, va4x0123, vb0123456789ABCDEFx0123); vacc5x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc5x0123456789ABCDEF, va5x0123, vb0123456789ABCDEFx0123); vacc6x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc6x0123456789ABCDEF, va6x0123, vb0123456789ABCDEFx0123); vacc7x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc7x0123456789ABCDEF, va7x0123, vb0123456789ABCDEFx0123); vacc8x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc8x0123456789ABCDEF, va8x0123, vb0123456789ABCDEFx0123); vacc9x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc9x0123456789ABCDEF, va9x0123, vb0123456789ABCDEFx0123); vacc1x0x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x0x0123456789ABCDEF, va0x4567, vb0123456789ABCDEFx4567); vacc1x1x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x1x0123456789ABCDEF, va1x4567, vb0123456789ABCDEFx4567); vacc1x2x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x2x0123456789ABCDEF, va2x4567, vb0123456789ABCDEFx4567); vacc1x3x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x3x0123456789ABCDEF, va3x4567, vb0123456789ABCDEFx4567); vacc1x4x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x4x0123456789ABCDEF, va4x4567, vb0123456789ABCDEFx4567); vacc1x5x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x5x0123456789ABCDEF, va5x4567, vb0123456789ABCDEFx4567); vacc1x6x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x6x0123456789ABCDEF, va6x4567, vb0123456789ABCDEFx4567); vacc1x7x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x7x0123456789ABCDEF, va7x4567, vb0123456789ABCDEFx4567); vacc1x8x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x8x0123456789ABCDEF, va8x4567, vb0123456789ABCDEFx4567); vacc1x9x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x9x0123456789ABCDEF, va9x4567, vb0123456789ABCDEFx4567); w = (const int8_t*) w + 128; k -= 8 * sizeof(int8_t); } vacc0x0123456789ABCDEF = _mm512_add_epi32(vacc0x0123456789ABCDEF, vacc1x0x0123456789ABCDEF); vacc1x0123456789ABCDEF = _mm512_add_epi32(vacc1x0123456789ABCDEF, vacc1x1x0123456789ABCDEF); vacc2x0123456789ABCDEF = _mm512_add_epi32(vacc2x0123456789ABCDEF, vacc1x2x0123456789ABCDEF); vacc3x0123456789ABCDEF = _mm512_add_epi32(vacc3x0123456789ABCDEF, vacc1x3x0123456789ABCDEF); vacc4x0123456789ABCDEF = _mm512_add_epi32(vacc4x0123456789ABCDEF, vacc1x4x0123456789ABCDEF); vacc5x0123456789ABCDEF = _mm512_add_epi32(vacc5x0123456789ABCDEF, vacc1x5x0123456789ABCDEF); vacc6x0123456789ABCDEF = _mm512_add_epi32(vacc6x0123456789ABCDEF, vacc1x6x0123456789ABCDEF); vacc7x0123456789ABCDEF = _mm512_add_epi32(vacc7x0123456789ABCDEF, vacc1x7x0123456789ABCDEF); vacc8x0123456789ABCDEF = _mm512_add_epi32(vacc8x0123456789ABCDEF, vacc1x8x0123456789ABCDEF); vacc9x0123456789ABCDEF = _mm512_add_epi32(vacc9x0123456789ABCDEF, vacc1x9x0123456789ABCDEF); if (k != 0) { const __m512i va0x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a0)); a0 += 4; const __m512i va1x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a1)); a1 += 4; const __m512i va2x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a2)); a2 += 4; const __m512i va3x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a3)); a3 += 4; const __m512i va4x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a4)); a4 += 4; const __m512i va5x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a5)); a5 += 4; const __m512i va6x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a6)); a6 += 4; const __m512i va7x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a7)); a7 += 4; const __m512i va8x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a8)); a8 += 4; const __m512i va9x0123 = _mm512_set1_epi32((int) unaligned_load_u32(a9)); a9 += 4; const __m512i vb0123456789ABCDEF = _mm512_load_si512(w); vacc0x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc0x0123456789ABCDEF, va0x0123, vb0123456789ABCDEF); vacc1x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc1x0123456789ABCDEF, va1x0123, vb0123456789ABCDEF); vacc2x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc2x0123456789ABCDEF, va2x0123, vb0123456789ABCDEF); vacc3x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc3x0123456789ABCDEF, va3x0123, vb0123456789ABCDEF); vacc4x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc4x0123456789ABCDEF, va4x0123, vb0123456789ABCDEF); vacc5x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc5x0123456789ABCDEF, va5x0123, vb0123456789ABCDEF); vacc6x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc6x0123456789ABCDEF, va6x0123, vb0123456789ABCDEF); vacc7x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc7x0123456789ABCDEF, va7x0123, vb0123456789ABCDEF); vacc8x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc8x0123456789ABCDEF, va8x0123, vb0123456789ABCDEF); vacc9x0123456789ABCDEF = _mm512_dpbusd_epi32(vacc9x0123456789ABCDEF, va9x0123, vb0123456789ABCDEF); w = (const int8_t*) w + 64; k -= 4 * sizeof(int8_t); } __m512 vscaled0x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc0x0123456789ABCDEF); __m512 vscaled1x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc1x0123456789ABCDEF); __m512 vscaled2x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc2x0123456789ABCDEF); __m512 vscaled3x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc3x0123456789ABCDEF); __m512 vscaled4x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc4x0123456789ABCDEF); __m512 vscaled5x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc5x0123456789ABCDEF); __m512 vscaled6x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc6x0123456789ABCDEF); __m512 vscaled7x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc7x0123456789ABCDEF); __m512 vscaled8x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc8x0123456789ABCDEF); __m512 vscaled9x0123456789ABCDEF = _mm512_cvtepi32_ps(vacc9x0123456789ABCDEF); vscaled0x0123456789ABCDEF = _mm512_mul_ps(vscaled0x0123456789ABCDEF, _mm512_set1_ps(quantization_params[0].inv_scale)); vscaled1x0123456789ABCDEF = _mm512_mul_ps(vscaled1x0123456789ABCDEF, _mm512_set1_ps(quantization_params[1].inv_scale)); vscaled2x0123456789ABCDEF = _mm512_mul_ps(vscaled2x0123456789ABCDEF, _mm512_set1_ps(quantization_params[2].inv_scale)); vscaled3x0123456789ABCDEF = _mm512_mul_ps(vscaled3x0123456789ABCDEF, _mm512_set1_ps(quantization_params[3].inv_scale)); vscaled4x0123456789ABCDEF = _mm512_mul_ps(vscaled4x0123456789ABCDEF, _mm512_set1_ps(quantization_params[4].inv_scale)); vscaled5x0123456789ABCDEF = _mm512_mul_ps(vscaled5x0123456789ABCDEF, _mm512_set1_ps(quantization_params[5].inv_scale)); vscaled6x0123456789ABCDEF = _mm512_mul_ps(vscaled6x0123456789ABCDEF, _mm512_set1_ps(quantization_params[6].inv_scale)); vscaled7x0123456789ABCDEF = _mm512_mul_ps(vscaled7x0123456789ABCDEF, _mm512_set1_ps(quantization_params[7].inv_scale)); vscaled8x0123456789ABCDEF = _mm512_mul_ps(vscaled8x0123456789ABCDEF, _mm512_set1_ps(quantization_params[8].inv_scale)); vscaled9x0123456789ABCDEF = _mm512_mul_ps(vscaled9x0123456789ABCDEF, _mm512_set1_ps(quantization_params[9].inv_scale)); const __m512 vfilter_output_scale0123456789ABCDEF = _mm512_load_ps((const float*) w); const __m512 vbias0123456789ABCDEF = _mm512_load_ps((const float*) w + 16); w = (const float*) w + 32; vscaled0x0123456789ABCDEF = _mm512_fmadd_ps(vscaled0x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled1x0123456789ABCDEF = _mm512_fmadd_ps(vscaled1x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled2x0123456789ABCDEF = _mm512_fmadd_ps(vscaled2x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled3x0123456789ABCDEF = _mm512_fmadd_ps(vscaled3x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled4x0123456789ABCDEF = _mm512_fmadd_ps(vscaled4x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled5x0123456789ABCDEF = _mm512_fmadd_ps(vscaled5x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled6x0123456789ABCDEF = _mm512_fmadd_ps(vscaled6x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled7x0123456789ABCDEF = _mm512_fmadd_ps(vscaled7x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled8x0123456789ABCDEF = _mm512_fmadd_ps(vscaled8x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled9x0123456789ABCDEF = _mm512_fmadd_ps(vscaled9x0123456789ABCDEF, vfilter_output_scale0123456789ABCDEF, vbias0123456789ABCDEF); vscaled0x0123456789ABCDEF = _mm512_max_ps(vscaled0x0123456789ABCDEF, voutput_min); vscaled1x0123456789ABCDEF = _mm512_max_ps(vscaled1x0123456789ABCDEF, voutput_min); vscaled2x0123456789ABCDEF = _mm512_max_ps(vscaled2x0123456789ABCDEF, voutput_min); vscaled3x0123456789ABCDEF = _mm512_max_ps(vscaled3x0123456789ABCDEF, voutput_min); vscaled4x0123456789ABCDEF = _mm512_max_ps(vscaled4x0123456789ABCDEF, voutput_min); vscaled5x0123456789ABCDEF = _mm512_max_ps(vscaled5x0123456789ABCDEF, voutput_min); vscaled6x0123456789ABCDEF = _mm512_max_ps(vscaled6x0123456789ABCDEF, voutput_min); vscaled7x0123456789ABCDEF = _mm512_max_ps(vscaled7x0123456789ABCDEF, voutput_min); vscaled8x0123456789ABCDEF = _mm512_max_ps(vscaled8x0123456789ABCDEF, voutput_min); vscaled9x0123456789ABCDEF = _mm512_max_ps(vscaled9x0123456789ABCDEF, voutput_min); vscaled0x0123456789ABCDEF = _mm512_min_ps(vscaled0x0123456789ABCDEF, voutput_max); vscaled1x0123456789ABCDEF = _mm512_min_ps(vscaled1x0123456789ABCDEF, voutput_max); vscaled2x0123456789ABCDEF = _mm512_min_ps(vscaled2x0123456789ABCDEF, voutput_max); vscaled3x0123456789ABCDEF = _mm512_min_ps(vscaled3x0123456789ABCDEF, voutput_max); vscaled4x0123456789ABCDEF = _mm512_min_ps(vscaled4x0123456789ABCDEF, voutput_max); vscaled5x0123456789ABCDEF = _mm512_min_ps(vscaled5x0123456789ABCDEF, voutput_max); vscaled6x0123456789ABCDEF = _mm512_min_ps(vscaled6x0123456789ABCDEF, voutput_max); vscaled7x0123456789ABCDEF = _mm512_min_ps(vscaled7x0123456789ABCDEF, voutput_max); vscaled8x0123456789ABCDEF = _mm512_min_ps(vscaled8x0123456789ABCDEF, voutput_max); vscaled9x0123456789ABCDEF = _mm512_min_ps(vscaled9x0123456789ABCDEF, voutput_max); if XNN_LIKELY(nc >= 16) { _mm512_storeu_ps(c0, vscaled0x0123456789ABCDEF); _mm512_storeu_ps(c1, vscaled1x0123456789ABCDEF); _mm512_storeu_ps(c2, vscaled2x0123456789ABCDEF); _mm512_storeu_ps(c3, vscaled3x0123456789ABCDEF); _mm512_storeu_ps(c4, vscaled4x0123456789ABCDEF); _mm512_storeu_ps(c5, vscaled5x0123456789ABCDEF); _mm512_storeu_ps(c6, vscaled6x0123456789ABCDEF); _mm512_storeu_ps(c7, vscaled7x0123456789ABCDEF); _mm512_storeu_ps(c8, vscaled8x0123456789ABCDEF); _mm512_storeu_ps(c9, vscaled9x0123456789ABCDEF); a0 = (const int8_t*) ((uintptr_t) a0 - kc); a1 = (const int8_t*) ((uintptr_t) a1 - kc); a2 = (const int8_t*) ((uintptr_t) a2 - kc); a3 = (const int8_t*) ((uintptr_t) a3 - kc); a4 = (const int8_t*) ((uintptr_t) a4 - kc); a5 = (const int8_t*) ((uintptr_t) a5 - kc); a6 = (const int8_t*) ((uintptr_t) a6 - kc); a7 = (const int8_t*) ((uintptr_t) a7 - kc); a8 = (const int8_t*) ((uintptr_t) a8 - kc); a9 = (const int8_t*) ((uintptr_t) a9 - kc); c0 = (float*) ((uintptr_t) c0 + cn_stride); c1 = (float*) ((uintptr_t) c1 + cn_stride); c2 = (float*) ((uintptr_t) c2 + cn_stride); c3 = (float*) ((uintptr_t) c3 + cn_stride); c4 = (float*) ((uintptr_t) c4 + cn_stride); c5 = (float*) ((uintptr_t) c5 + cn_stride); c6 = (float*) ((uintptr_t) c6 + cn_stride); c7 = (float*) ((uintptr_t) c7 + cn_stride); c8 = (float*) ((uintptr_t) c8 + cn_stride); c9 = (float*) ((uintptr_t) c9 + cn_stride); nc -= 16; } else { // Prepare mask for valid 32-bit elements (depends on nc). const __mmask16 vmask = _cvtu32_mask16((UINT32_C(1) << nc) - 1); _mm512_mask_storeu_ps(c0, vmask, vscaled0x0123456789ABCDEF); _mm512_mask_storeu_ps(c1, vmask, vscaled1x0123456789ABCDEF); _mm512_mask_storeu_ps(c2, vmask, vscaled2x0123456789ABCDEF); _mm512_mask_storeu_ps(c3, vmask, vscaled3x0123456789ABCDEF); _mm512_mask_storeu_ps(c4, vmask, vscaled4x0123456789ABCDEF); _mm512_mask_storeu_ps(c5, vmask, vscaled5x0123456789ABCDEF); _mm512_mask_storeu_ps(c6, vmask, vscaled6x0123456789ABCDEF); _mm512_mask_storeu_ps(c7, vmask, vscaled7x0123456789ABCDEF); _mm512_mask_storeu_ps(c8, vmask, vscaled8x0123456789ABCDEF); _mm512_mask_storeu_ps(c9, vmask, vscaled9x0123456789ABCDEF); nc = 0; } } while (nc != 0); }