1// Copyright 2020 Google LLC 2// 3// This source code is licensed under the BSD-style license found in the 4// LICENSE file in the root directory of this source tree. 5 6$assert NR % 4 == 0 7$assert ACTIVATION != "MINMAX" or ARCH in ["ARM", "X86", "RELAXED"] 8$assert not FMA or ARCH == "RELAXED" 9$ABC = "0123456789ABCDEFGHIJKLMNOPQRSTUVWXYZ" 10#include <assert.h> 11 12#include <wasm_simd128.h> 13 14#include <xnnpack/gemm.h> 15 16 17$assert ACTIVATION in ["LINEAR", "RELU", "MINMAX"] 18$if ACTIVATION == "MINMAX": 19$ WASM_F32X4_MIN={"ARM": "wasm_f32x4_min", "X86": "wasm_f32x4_pmin", "RELAXED": "__builtin_wasm_relaxed_min_f32x4"}[ARCH] 20$ WASM_F32X4_MAX={"ARM": "wasm_f32x4_max", "X86": "wasm_f32x4_pmax", "RELAXED": "__builtin_wasm_relaxed_max_f32x4"}[ARCH] 21$ACTIVATION_SUFFIX = {"LINEAR": ""}.get(ACTIVATION, "_" + ACTIVATION.lower()) 22$ISA = "wasmsimd" if not FMA and (ACTIVATION in ["LINEAR", "RELU"] or ARCH != "RELAXED") else "wasmrelaxedsimd" 23$ARCH_SUFFIX = "" if not FMA and (ACTIVATION in ["LINEAR", "RELU"] or ARCH == "RELAXED") else "_" + ("fma" if FMA else ARCH.lower()) 24$PARAMS = {"LINEAR": "xnn_f32_default_params", "RELU": "xnn_f32_relu_params", "MINMAX": "xnn_f32_minmax_params"}[ACTIVATION] 25void xnn_f32_gemm${"inc" if INC else ""}${ACTIVATION_SUFFIX}_ukernel_${MR}x${NR}s4__${ISA}${ARCH_SUFFIX}( 26 size_t mr, 27 size_t nc, 28 size_t kc, 29 const float*restrict a, 30 size_t a_stride, 31 const float*restrict w, 32 float*restrict c, 33 size_t cm_stride, 34 size_t cn_stride, 35 $if INC: 36 const float*restrict acc, 37 const union ${PARAMS} params[restrict XNN_MIN_ELEMENTS(1)]) 38{ 39 assert(mr != 0); 40 assert(mr <= ${MR}); 41 assert(nc != 0); 42 assert(kc != 0); 43 assert(kc % sizeof(float) == 0); 44 assert(a != NULL); 45 assert(w != NULL); 46 assert(c != NULL); 47 $if INC: 48 assert(acc != NULL); 49 50 const float* a0 = a; 51 float* c0 = c; 52 $for M in range(1, MR): 53 const float* a${M} = (const float*) ((uintptr_t) a${M-1} + a_stride); 54 float* c${M} = (float*) ((uintptr_t) c${M-1} + cm_stride); 55 $if M % 2 == 0: 56 if XNN_UNPREDICTABLE(mr <= ${M}) { 57 a${M} = a${M-1}; 58 c${M} = c${M-1}; 59 } 60 $elif M + 1 == MR: 61 if XNN_UNPREDICTABLE(mr != ${M+1}) { 62 a${M} = a${M-1}; 63 c${M} = c${M-1}; 64 } 65 $else: 66 if XNN_UNPREDICTABLE(mr < ${M+1}) { 67 a${M} = a${M-1}; 68 c${M} = c${M-1}; 69 } 70 71 $if ACTIVATION == "MINMAX": 72 const v128_t vmin = wasm_v128_load64_splat(params->wasmsimd.min); 73 const v128_t vmax = wasm_v128_load64_splat(params->wasmsimd.max); 74 do { 75 $if INC: 76 $for M in range(MR): 77 $for N in range(0, NR, 4): 78 v128_t vacc${M}x${ABC[N:N+4]} = wasm_v128_load(acc + ${M*NR+N}); 79 acc += ${MR*NR}; 80 $else: 81 $for N in range(0, NR, 4): 82 v128_t vacc0x${ABC[N:N+4]} = wasm_v128_load(w + ${N}); 83 $for M in range(1, MR): 84 $for N in range(0, NR, 4): 85 v128_t vacc${M}x${ABC[N:N+4]} = vacc0x${ABC[N:N+4]}; 86 w += ${NR}; 87 88 size_t k = kc; 89 while (k >= 4 * sizeof(float)) { 90 $for M in range(MR): 91 v128_t va${M} = wasm_v128_load(a${M}); 92 a${M} += 4; 93 94 $for L in range(4): 95 96 $for N in range(0, NR, 4): 97 const v128_t vb${ABC[N:N+4]}c${L} = wasm_v128_load(w + ${L * NR + N}); 98 99 $for N in range(0, NR, 4): 100 $for M in range(MR): 101 $if FMA: 102 vacc${M}x${ABC[N:N+4]} = __builtin_wasm_fma_f32x4(vacc${M}x${ABC[N:N+4]}, va${M}, vb${ABC[N:N+4]}c${L}); 103 $else: 104 vacc${M}x${ABC[N:N+4]} = wasm_f32x4_add(vacc${M}x${ABC[N:N+4]}, wasm_f32x4_mul(va${M}, vb${ABC[N:N+4]}c${L})); 105 106 $if L + 1 != 4: 107 $for M in range(MR): 108 va${M} = wasm_v32x4_shuffle(va${M}, va${M}, 1, 2, 3, 0); 109 110 w += ${4 * NR}; 111 k -= 4 * sizeof(float); 112 } 113 if XNN_UNLIKELY(k != 0) { 114 $for M in range(MR): 115 v128_t va${M} = wasm_v128_load(a${M}); 116 a${M} = (const float*) ((uintptr_t) a${M} + k); 117 118 const v128_t vzero = wasm_f32x4_const_splat(0.0f); 119 $for L in range(4): 120 121 $for N in range(0, NR, 4): 122 const v128_t vb${ABC[N:N+4]}c${L} = wasm_v128_load(w + ${L * NR + N}); 123 124 $for N in range(0, NR, 4): 125 $for M in range(MR): 126 $if FMA: 127 vacc${M}x${ABC[N:N+4]} = __builtin_wasm_fma_f32x4(vacc${M}x${ABC[N:N+4]}, wasm_v128_andnot(va${M}, wasm_f32x4_eq(vb${ABC[N:N+4]}c${L}, vzero)), vb${ABC[N:N+4]}c${L}); 128 $else: 129 vacc${M}x${ABC[N:N+4]} = wasm_f32x4_add(vacc${M}x${ABC[N:N+4]}, wasm_f32x4_mul(wasm_v128_andnot(va${M}, wasm_f32x4_eq(vb${ABC[N:N+4]}c${L}, vzero)), vb${ABC[N:N+4]}c${L})); 130 131 $if L + 1 != 4: 132 $for M in range(MR): 133 va${M} = wasm_v32x4_shuffle(va${M}, va${M}, 1, 2, 3, 0); 134 135 w += ${4 * NR}; 136 } 137 138 $if ACTIVATION == "MINMAX": 139 $for N in range(0, NR, 4): 140 $for M in range(MR): 141 vacc${M}x${ABC[N:N+4]} = ${WASM_F32X4_MAX}(vmin, vacc${M}x${ABC[N:N+4]}); 142 143 $for N in range(0, NR, 4): 144 $for M in range(MR): 145 vacc${M}x${ABC[N:N+4]} = ${WASM_F32X4_MIN}(vmax, vacc${M}x${ABC[N:N+4]}); 146 $elif ACTIVATION == "RELU": 147 const v128_t vzero = wasm_i32x4_const_splat(0); 148 $for N in range(0, NR, 4): 149 $for M in range(MR): 150 vacc${M}x${ABC[N:N+4]} = wasm_i32x4_max(vacc${M}x${ABC[N:N+4]}, vzero); 151 152 if XNN_LIKELY(nc >= ${NR}) { 153 $for M in reversed(range(MR)): 154 wasm_v128_store(c${M}, vacc${M}x${ABC[0:4]}); 155 $for N in range(4, NR, 4): 156 wasm_v128_store(c${M} + ${N}, vacc${M}x${ABC[N:N+4]}); 157 c${M} = (float*) ((uintptr_t) c${M} + cn_stride); 158 159 $for M in reversed(range(MR)): 160 a${M} = (const float*) ((uintptr_t) a${M} - kc); 161 162 nc -= ${NR}; 163 } else { 164 $for LOG2N in reversed(range(NR.bit_length())): 165 $if NR != 1 << LOG2N: 166 if (nc & ${1 << LOG2N}) { 167 $if LOG2N >= 2: 168 $for M in reversed(range(MR)): 169 wasm_v128_store(c${M}, vacc${M}x${ABC[0:4]}); 170 $for N in range(4, 1 << LOG2N, 4): 171 wasm_v128_store(c${M} + ${N}, vacc${M}x${ABC[N:N+4]}); 172 173 $for M in reversed(range(MR)): 174 $for N in range(0, 1 << (LOG2N - 1), 4): 175 vacc${M}x${ABC[N:N+4]} = vacc${M}x${ABC[N + (1 << LOG2N):N + (1 << LOG2N)+4]}; 176 177 $for M in reversed(range(MR)): 178 c${M} += ${1 << LOG2N}; 179 $elif LOG2N == 1: 180 $for M in reversed(range(MR)): 181 *((double*) c${M}) = wasm_f64x2_extract_lane(vacc${M}x${ABC[0:4]}, 0); 182 183 $for M in reversed(range(MR)): 184 vacc${M}x${ABC[0:4]} = wasm_v32x4_shuffle(vacc${M}x${ABC[0:4]}, vacc${M}x${ABC[0:4]}, 2, 3, 2, 3); 185 186 $for M in reversed(range(MR)): 187 c${M} += 2; 188 $elif LOG2N == 0: 189 $for M in reversed(range(MR)): 190 *c${M} = wasm_f32x4_extract_lane(vacc${M}x${ABC[0:4]}, 0); 191 } 192 193 nc = 0; 194 } 195 } while (nc != 0); 196} 197