1 /* 2 * CDDL HEADER START 3 * 4 * The contents of this file are subject to the terms of the 5 * Common Development and Distribution License (the "License"). 6 * You may not use this file except in compliance with the License. 7 * 8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE 9 * or https://opensource.org/licenses/CDDL-1.0. 10 * See the License for the specific language governing permissions 11 * and limitations under the License. 12 * 13 * When distributing Covered Code, include this CDDL HEADER in each 14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE. 15 * If applicable, add the following below this CDDL HEADER, with the 16 * fields enclosed by brackets "[]" replaced with your own identifying 17 * information: Portions Copyright [yyyy] [name of copyright owner] 18 * 19 * CDDL HEADER END 20 */ 21 /* 22 * Copyright (C) 2016 Romain Dolbeau. All rights reserved. 23 * Copyright (C) 2016 Gvozden Nešković. All rights reserved. 24 */ 25 26 #include <sys/isa_defs.h> 27 28 #if defined(__x86_64) && defined(HAVE_AVX512F) 29 30 #include <sys/types.h> 31 #include <sys/simd.h> 32 #include <sys/debug.h> 33 34 #ifdef __linux__ 35 #define __asm __asm__ __volatile__ 36 #endif 37 38 #define _REG_CNT(_0, _1, _2, _3, _4, _5, _6, _7, N, ...) N 39 #define REG_CNT(r...) _REG_CNT(r, 8, 7, 6, 5, 4, 3, 2, 1) 40 41 #define VR0_(REG, ...) "zmm"#REG 42 #define VR1_(_1, REG, ...) "zmm"#REG 43 #define VR2_(_1, _2, REG, ...) "zmm"#REG 44 #define VR3_(_1, _2, _3, REG, ...) "zmm"#REG 45 #define VR4_(_1, _2, _3, _4, REG, ...) "zmm"#REG 46 #define VR5_(_1, _2, _3, _4, _5, REG, ...) "zmm"#REG 47 #define VR6_(_1, _2, _3, _4, _5, _6, REG, ...) "zmm"#REG 48 #define VR7_(_1, _2, _3, _4, _5, _6, _7, REG, ...) "zmm"#REG 49 50 #define VR0(r...) VR0_(r) 51 #define VR1(r...) VR1_(r) 52 #define VR2(r...) VR2_(r, 1) 53 #define VR3(r...) VR3_(r, 1, 2) 54 #define VR4(r...) VR4_(r, 1, 2) 55 #define VR5(r...) VR5_(r, 1, 2, 3) 56 #define VR6(r...) VR6_(r, 1, 2, 3, 4) 57 #define VR7(r...) VR7_(r, 1, 2, 3, 4, 5) 58 59 #define VRy0_(REG, ...) "ymm"#REG 60 #define VRy1_(_1, REG, ...) "ymm"#REG 61 #define VRy2_(_1, _2, REG, ...) "ymm"#REG 62 #define VRy3_(_1, _2, _3, REG, ...) "ymm"#REG 63 #define VRy4_(_1, _2, _3, _4, REG, ...) "ymm"#REG 64 #define VRy5_(_1, _2, _3, _4, _5, REG, ...) "ymm"#REG 65 #define VRy6_(_1, _2, _3, _4, _5, _6, REG, ...) "ymm"#REG 66 #define VRy7_(_1, _2, _3, _4, _5, _6, _7, REG, ...) "ymm"#REG 67 68 #define VRy0(r...) VRy0_(r) 69 #define VRy1(r...) VRy1_(r) 70 #define VRy2(r...) VRy2_(r, 1) 71 #define VRy3(r...) VRy3_(r, 1, 2) 72 #define VRy4(r...) VRy4_(r, 1, 2) 73 #define VRy5(r...) VRy5_(r, 1, 2, 3) 74 #define VRy6(r...) VRy6_(r, 1, 2, 3, 4) 75 #define VRy7(r...) VRy7_(r, 1, 2, 3, 4, 5) 76 77 #define R_01(REG1, REG2, ...) REG1, REG2 78 #define _R_23(_0, _1, REG2, REG3, ...) REG2, REG3 79 #define R_23(REG...) _R_23(REG, 1, 2, 3) 80 81 #define ELEM_SIZE 64 82 83 typedef struct v { 84 uint8_t b[ELEM_SIZE] __attribute__((aligned(ELEM_SIZE))); 85 } v_t; 86 87 88 #define XOR_ACC(src, r...) \ 89 { \ 90 switch (REG_CNT(r)) { \ 91 case 4: \ 92 __asm( \ 93 "vpxorq 0x00(%[SRC]), %%" VR0(r)", %%" VR0(r) "\n" \ 94 "vpxorq 0x40(%[SRC]), %%" VR1(r)", %%" VR1(r) "\n" \ 95 "vpxorq 0x80(%[SRC]), %%" VR2(r)", %%" VR2(r) "\n" \ 96 "vpxorq 0xc0(%[SRC]), %%" VR3(r)", %%" VR3(r) "\n" \ 97 : : [SRC] "r" (src)); \ 98 break; \ 99 } \ 100 } 101 102 #define XOR(r...) \ 103 { \ 104 switch (REG_CNT(r)) { \ 105 case 8: \ 106 __asm( \ 107 "vpxorq %" VR0(r) ", %" VR4(r)", %" VR4(r) "\n" \ 108 "vpxorq %" VR1(r) ", %" VR5(r)", %" VR5(r) "\n" \ 109 "vpxorq %" VR2(r) ", %" VR6(r)", %" VR6(r) "\n" \ 110 "vpxorq %" VR3(r) ", %" VR7(r)", %" VR7(r)); \ 111 break; \ 112 case 4: \ 113 __asm( \ 114 "vpxorq %" VR0(r) ", %" VR2(r)", %" VR2(r) "\n" \ 115 "vpxorq %" VR1(r) ", %" VR3(r)", %" VR3(r)); \ 116 break; \ 117 } \ 118 } 119 120 121 #define ZERO(r...) XOR(r, r) 122 123 124 #define COPY(r...) \ 125 { \ 126 switch (REG_CNT(r)) { \ 127 case 8: \ 128 __asm( \ 129 "vmovdqa64 %" VR0(r) ", %" VR4(r) "\n" \ 130 "vmovdqa64 %" VR1(r) ", %" VR5(r) "\n" \ 131 "vmovdqa64 %" VR2(r) ", %" VR6(r) "\n" \ 132 "vmovdqa64 %" VR3(r) ", %" VR7(r)); \ 133 break; \ 134 case 4: \ 135 __asm( \ 136 "vmovdqa64 %" VR0(r) ", %" VR2(r) "\n" \ 137 "vmovdqa64 %" VR1(r) ", %" VR3(r)); \ 138 break; \ 139 } \ 140 } 141 142 #define LOAD(src, r...) \ 143 { \ 144 switch (REG_CNT(r)) { \ 145 case 4: \ 146 __asm( \ 147 "vmovdqa64 0x00(%[SRC]), %%" VR0(r) "\n" \ 148 "vmovdqa64 0x40(%[SRC]), %%" VR1(r) "\n" \ 149 "vmovdqa64 0x80(%[SRC]), %%" VR2(r) "\n" \ 150 "vmovdqa64 0xc0(%[SRC]), %%" VR3(r) "\n" \ 151 : : [SRC] "r" (src)); \ 152 break; \ 153 } \ 154 } 155 156 #define STORE(dst, r...) \ 157 { \ 158 switch (REG_CNT(r)) { \ 159 case 4: \ 160 __asm( \ 161 "vmovdqa64 %%" VR0(r) ", 0x00(%[DST])\n" \ 162 "vmovdqa64 %%" VR1(r) ", 0x40(%[DST])\n" \ 163 "vmovdqa64 %%" VR2(r) ", 0x80(%[DST])\n" \ 164 "vmovdqa64 %%" VR3(r) ", 0xc0(%[DST])\n" \ 165 : : [DST] "r" (dst)); \ 166 break; \ 167 } \ 168 } 169 170 #define MUL2_SETUP() \ 171 { \ 172 __asm("vmovq %0, %%xmm31" :: "r"(0x1d1d1d1d1d1d1d1d)); \ 173 __asm("vpbroadcastq %xmm31, %zmm31"); \ 174 __asm("vmovq %0, %%xmm30" :: "r"(0x8080808080808080)); \ 175 __asm("vpbroadcastq %xmm30, %zmm30"); \ 176 __asm("vmovq %0, %%xmm29" :: "r"(0xfefefefefefefefe)); \ 177 __asm("vpbroadcastq %xmm29, %zmm29"); \ 178 } 179 180 #define _MUL2(r...) \ 181 { \ 182 switch (REG_CNT(r)) { \ 183 case 2: \ 184 __asm( \ 185 "vpandq %" VR0(r)", %zmm30, %zmm26\n" \ 186 "vpandq %" VR1(r)", %zmm30, %zmm25\n" \ 187 "vpsrlq $7, %zmm26, %zmm28\n" \ 188 "vpsrlq $7, %zmm25, %zmm27\n" \ 189 "vpsllq $1, %zmm26, %zmm26\n" \ 190 "vpsllq $1, %zmm25, %zmm25\n" \ 191 "vpsubq %zmm28, %zmm26, %zmm26\n" \ 192 "vpsubq %zmm27, %zmm25, %zmm25\n" \ 193 "vpsllq $1, %" VR0(r)", %" VR0(r) "\n" \ 194 "vpsllq $1, %" VR1(r)", %" VR1(r) "\n" \ 195 "vpandq %zmm26, %zmm31, %zmm26\n" \ 196 "vpandq %zmm25, %zmm31, %zmm25\n" \ 197 "vpternlogd $0x6c,%zmm29, %zmm26, %" VR0(r) "\n" \ 198 "vpternlogd $0x6c,%zmm29, %zmm25, %" VR1(r)); \ 199 break; \ 200 default: \ 201 VERIFY(0); \ 202 } \ 203 } 204 205 #define MUL2(r...) \ 206 { \ 207 switch (REG_CNT(r)) { \ 208 case 4: \ 209 _MUL2(R_01(r)); \ 210 _MUL2(R_23(r)); \ 211 break; \ 212 case 2: \ 213 _MUL2(r); \ 214 break; \ 215 } \ 216 } 217 218 #define MUL4(r...) \ 219 { \ 220 MUL2(r); \ 221 MUL2(r); \ 222 } 223 224 225 /* General multiplication by adding powers of two */ 226 227 #define _mul_x2_in 21, 22 228 #define _mul_x2_acc 23, 24 229 230 #define _MUL_PARAM(x, in, acc) \ 231 { \ 232 if (x & 0x01) { COPY(in, acc); } else { ZERO(acc); } \ 233 if (x & 0xfe) { MUL2(in); } \ 234 if (x & 0x02) { XOR(in, acc); } \ 235 if (x & 0xfc) { MUL2(in); } \ 236 if (x & 0x04) { XOR(in, acc); } \ 237 if (x & 0xf8) { MUL2(in); } \ 238 if (x & 0x08) { XOR(in, acc); } \ 239 if (x & 0xf0) { MUL2(in); } \ 240 if (x & 0x10) { XOR(in, acc); } \ 241 if (x & 0xe0) { MUL2(in); } \ 242 if (x & 0x20) { XOR(in, acc); } \ 243 if (x & 0xc0) { MUL2(in); } \ 244 if (x & 0x40) { XOR(in, acc); } \ 245 if (x & 0x80) { MUL2(in); XOR(in, acc); } \ 246 } 247 248 #define MUL_x2_DEFINE(x) \ 249 static void \ 250 mul_x2_ ## x(void) { _MUL_PARAM(x, _mul_x2_in, _mul_x2_acc); } 251 252 253 MUL_x2_DEFINE(0); MUL_x2_DEFINE(1); MUL_x2_DEFINE(2); MUL_x2_DEFINE(3); 254 MUL_x2_DEFINE(4); MUL_x2_DEFINE(5); MUL_x2_DEFINE(6); MUL_x2_DEFINE(7); 255 MUL_x2_DEFINE(8); MUL_x2_DEFINE(9); MUL_x2_DEFINE(10); MUL_x2_DEFINE(11); 256 MUL_x2_DEFINE(12); MUL_x2_DEFINE(13); MUL_x2_DEFINE(14); MUL_x2_DEFINE(15); 257 MUL_x2_DEFINE(16); MUL_x2_DEFINE(17); MUL_x2_DEFINE(18); MUL_x2_DEFINE(19); 258 MUL_x2_DEFINE(20); MUL_x2_DEFINE(21); MUL_x2_DEFINE(22); MUL_x2_DEFINE(23); 259 MUL_x2_DEFINE(24); MUL_x2_DEFINE(25); MUL_x2_DEFINE(26); MUL_x2_DEFINE(27); 260 MUL_x2_DEFINE(28); MUL_x2_DEFINE(29); MUL_x2_DEFINE(30); MUL_x2_DEFINE(31); 261 MUL_x2_DEFINE(32); MUL_x2_DEFINE(33); MUL_x2_DEFINE(34); MUL_x2_DEFINE(35); 262 MUL_x2_DEFINE(36); MUL_x2_DEFINE(37); MUL_x2_DEFINE(38); MUL_x2_DEFINE(39); 263 MUL_x2_DEFINE(40); MUL_x2_DEFINE(41); MUL_x2_DEFINE(42); MUL_x2_DEFINE(43); 264 MUL_x2_DEFINE(44); MUL_x2_DEFINE(45); MUL_x2_DEFINE(46); MUL_x2_DEFINE(47); 265 MUL_x2_DEFINE(48); MUL_x2_DEFINE(49); MUL_x2_DEFINE(50); MUL_x2_DEFINE(51); 266 MUL_x2_DEFINE(52); MUL_x2_DEFINE(53); MUL_x2_DEFINE(54); MUL_x2_DEFINE(55); 267 MUL_x2_DEFINE(56); MUL_x2_DEFINE(57); MUL_x2_DEFINE(58); MUL_x2_DEFINE(59); 268 MUL_x2_DEFINE(60); MUL_x2_DEFINE(61); MUL_x2_DEFINE(62); MUL_x2_DEFINE(63); 269 MUL_x2_DEFINE(64); MUL_x2_DEFINE(65); MUL_x2_DEFINE(66); MUL_x2_DEFINE(67); 270 MUL_x2_DEFINE(68); MUL_x2_DEFINE(69); MUL_x2_DEFINE(70); MUL_x2_DEFINE(71); 271 MUL_x2_DEFINE(72); MUL_x2_DEFINE(73); MUL_x2_DEFINE(74); MUL_x2_DEFINE(75); 272 MUL_x2_DEFINE(76); MUL_x2_DEFINE(77); MUL_x2_DEFINE(78); MUL_x2_DEFINE(79); 273 MUL_x2_DEFINE(80); MUL_x2_DEFINE(81); MUL_x2_DEFINE(82); MUL_x2_DEFINE(83); 274 MUL_x2_DEFINE(84); MUL_x2_DEFINE(85); MUL_x2_DEFINE(86); MUL_x2_DEFINE(87); 275 MUL_x2_DEFINE(88); MUL_x2_DEFINE(89); MUL_x2_DEFINE(90); MUL_x2_DEFINE(91); 276 MUL_x2_DEFINE(92); MUL_x2_DEFINE(93); MUL_x2_DEFINE(94); MUL_x2_DEFINE(95); 277 MUL_x2_DEFINE(96); MUL_x2_DEFINE(97); MUL_x2_DEFINE(98); MUL_x2_DEFINE(99); 278 MUL_x2_DEFINE(100); MUL_x2_DEFINE(101); MUL_x2_DEFINE(102); MUL_x2_DEFINE(103); 279 MUL_x2_DEFINE(104); MUL_x2_DEFINE(105); MUL_x2_DEFINE(106); MUL_x2_DEFINE(107); 280 MUL_x2_DEFINE(108); MUL_x2_DEFINE(109); MUL_x2_DEFINE(110); MUL_x2_DEFINE(111); 281 MUL_x2_DEFINE(112); MUL_x2_DEFINE(113); MUL_x2_DEFINE(114); MUL_x2_DEFINE(115); 282 MUL_x2_DEFINE(116); MUL_x2_DEFINE(117); MUL_x2_DEFINE(118); MUL_x2_DEFINE(119); 283 MUL_x2_DEFINE(120); MUL_x2_DEFINE(121); MUL_x2_DEFINE(122); MUL_x2_DEFINE(123); 284 MUL_x2_DEFINE(124); MUL_x2_DEFINE(125); MUL_x2_DEFINE(126); MUL_x2_DEFINE(127); 285 MUL_x2_DEFINE(128); MUL_x2_DEFINE(129); MUL_x2_DEFINE(130); MUL_x2_DEFINE(131); 286 MUL_x2_DEFINE(132); MUL_x2_DEFINE(133); MUL_x2_DEFINE(134); MUL_x2_DEFINE(135); 287 MUL_x2_DEFINE(136); MUL_x2_DEFINE(137); MUL_x2_DEFINE(138); MUL_x2_DEFINE(139); 288 MUL_x2_DEFINE(140); MUL_x2_DEFINE(141); MUL_x2_DEFINE(142); MUL_x2_DEFINE(143); 289 MUL_x2_DEFINE(144); MUL_x2_DEFINE(145); MUL_x2_DEFINE(146); MUL_x2_DEFINE(147); 290 MUL_x2_DEFINE(148); MUL_x2_DEFINE(149); MUL_x2_DEFINE(150); MUL_x2_DEFINE(151); 291 MUL_x2_DEFINE(152); MUL_x2_DEFINE(153); MUL_x2_DEFINE(154); MUL_x2_DEFINE(155); 292 MUL_x2_DEFINE(156); MUL_x2_DEFINE(157); MUL_x2_DEFINE(158); MUL_x2_DEFINE(159); 293 MUL_x2_DEFINE(160); MUL_x2_DEFINE(161); MUL_x2_DEFINE(162); MUL_x2_DEFINE(163); 294 MUL_x2_DEFINE(164); MUL_x2_DEFINE(165); MUL_x2_DEFINE(166); MUL_x2_DEFINE(167); 295 MUL_x2_DEFINE(168); MUL_x2_DEFINE(169); MUL_x2_DEFINE(170); MUL_x2_DEFINE(171); 296 MUL_x2_DEFINE(172); MUL_x2_DEFINE(173); MUL_x2_DEFINE(174); MUL_x2_DEFINE(175); 297 MUL_x2_DEFINE(176); MUL_x2_DEFINE(177); MUL_x2_DEFINE(178); MUL_x2_DEFINE(179); 298 MUL_x2_DEFINE(180); MUL_x2_DEFINE(181); MUL_x2_DEFINE(182); MUL_x2_DEFINE(183); 299 MUL_x2_DEFINE(184); MUL_x2_DEFINE(185); MUL_x2_DEFINE(186); MUL_x2_DEFINE(187); 300 MUL_x2_DEFINE(188); MUL_x2_DEFINE(189); MUL_x2_DEFINE(190); MUL_x2_DEFINE(191); 301 MUL_x2_DEFINE(192); MUL_x2_DEFINE(193); MUL_x2_DEFINE(194); MUL_x2_DEFINE(195); 302 MUL_x2_DEFINE(196); MUL_x2_DEFINE(197); MUL_x2_DEFINE(198); MUL_x2_DEFINE(199); 303 MUL_x2_DEFINE(200); MUL_x2_DEFINE(201); MUL_x2_DEFINE(202); MUL_x2_DEFINE(203); 304 MUL_x2_DEFINE(204); MUL_x2_DEFINE(205); MUL_x2_DEFINE(206); MUL_x2_DEFINE(207); 305 MUL_x2_DEFINE(208); MUL_x2_DEFINE(209); MUL_x2_DEFINE(210); MUL_x2_DEFINE(211); 306 MUL_x2_DEFINE(212); MUL_x2_DEFINE(213); MUL_x2_DEFINE(214); MUL_x2_DEFINE(215); 307 MUL_x2_DEFINE(216); MUL_x2_DEFINE(217); MUL_x2_DEFINE(218); MUL_x2_DEFINE(219); 308 MUL_x2_DEFINE(220); MUL_x2_DEFINE(221); MUL_x2_DEFINE(222); MUL_x2_DEFINE(223); 309 MUL_x2_DEFINE(224); MUL_x2_DEFINE(225); MUL_x2_DEFINE(226); MUL_x2_DEFINE(227); 310 MUL_x2_DEFINE(228); MUL_x2_DEFINE(229); MUL_x2_DEFINE(230); MUL_x2_DEFINE(231); 311 MUL_x2_DEFINE(232); MUL_x2_DEFINE(233); MUL_x2_DEFINE(234); MUL_x2_DEFINE(235); 312 MUL_x2_DEFINE(236); MUL_x2_DEFINE(237); MUL_x2_DEFINE(238); MUL_x2_DEFINE(239); 313 MUL_x2_DEFINE(240); MUL_x2_DEFINE(241); MUL_x2_DEFINE(242); MUL_x2_DEFINE(243); 314 MUL_x2_DEFINE(244); MUL_x2_DEFINE(245); MUL_x2_DEFINE(246); MUL_x2_DEFINE(247); 315 MUL_x2_DEFINE(248); MUL_x2_DEFINE(249); MUL_x2_DEFINE(250); MUL_x2_DEFINE(251); 316 MUL_x2_DEFINE(252); MUL_x2_DEFINE(253); MUL_x2_DEFINE(254); MUL_x2_DEFINE(255); 317 318 319 typedef void (*mul_fn_ptr_t)(void); 320 321 static const mul_fn_ptr_t __attribute__((aligned(256))) 322 gf_x2_mul_fns[256] = { 323 mul_x2_0, mul_x2_1, mul_x2_2, mul_x2_3, mul_x2_4, mul_x2_5, 324 mul_x2_6, mul_x2_7, mul_x2_8, mul_x2_9, mul_x2_10, mul_x2_11, 325 mul_x2_12, mul_x2_13, mul_x2_14, mul_x2_15, mul_x2_16, mul_x2_17, 326 mul_x2_18, mul_x2_19, mul_x2_20, mul_x2_21, mul_x2_22, mul_x2_23, 327 mul_x2_24, mul_x2_25, mul_x2_26, mul_x2_27, mul_x2_28, mul_x2_29, 328 mul_x2_30, mul_x2_31, mul_x2_32, mul_x2_33, mul_x2_34, mul_x2_35, 329 mul_x2_36, mul_x2_37, mul_x2_38, mul_x2_39, mul_x2_40, mul_x2_41, 330 mul_x2_42, mul_x2_43, mul_x2_44, mul_x2_45, mul_x2_46, mul_x2_47, 331 mul_x2_48, mul_x2_49, mul_x2_50, mul_x2_51, mul_x2_52, mul_x2_53, 332 mul_x2_54, mul_x2_55, mul_x2_56, mul_x2_57, mul_x2_58, mul_x2_59, 333 mul_x2_60, mul_x2_61, mul_x2_62, mul_x2_63, mul_x2_64, mul_x2_65, 334 mul_x2_66, mul_x2_67, mul_x2_68, mul_x2_69, mul_x2_70, mul_x2_71, 335 mul_x2_72, mul_x2_73, mul_x2_74, mul_x2_75, mul_x2_76, mul_x2_77, 336 mul_x2_78, mul_x2_79, mul_x2_80, mul_x2_81, mul_x2_82, mul_x2_83, 337 mul_x2_84, mul_x2_85, mul_x2_86, mul_x2_87, mul_x2_88, mul_x2_89, 338 mul_x2_90, mul_x2_91, mul_x2_92, mul_x2_93, mul_x2_94, mul_x2_95, 339 mul_x2_96, mul_x2_97, mul_x2_98, mul_x2_99, mul_x2_100, mul_x2_101, 340 mul_x2_102, mul_x2_103, mul_x2_104, mul_x2_105, mul_x2_106, mul_x2_107, 341 mul_x2_108, mul_x2_109, mul_x2_110, mul_x2_111, mul_x2_112, mul_x2_113, 342 mul_x2_114, mul_x2_115, mul_x2_116, mul_x2_117, mul_x2_118, mul_x2_119, 343 mul_x2_120, mul_x2_121, mul_x2_122, mul_x2_123, mul_x2_124, mul_x2_125, 344 mul_x2_126, mul_x2_127, mul_x2_128, mul_x2_129, mul_x2_130, mul_x2_131, 345 mul_x2_132, mul_x2_133, mul_x2_134, mul_x2_135, mul_x2_136, mul_x2_137, 346 mul_x2_138, mul_x2_139, mul_x2_140, mul_x2_141, mul_x2_142, mul_x2_143, 347 mul_x2_144, mul_x2_145, mul_x2_146, mul_x2_147, mul_x2_148, mul_x2_149, 348 mul_x2_150, mul_x2_151, mul_x2_152, mul_x2_153, mul_x2_154, mul_x2_155, 349 mul_x2_156, mul_x2_157, mul_x2_158, mul_x2_159, mul_x2_160, mul_x2_161, 350 mul_x2_162, mul_x2_163, mul_x2_164, mul_x2_165, mul_x2_166, mul_x2_167, 351 mul_x2_168, mul_x2_169, mul_x2_170, mul_x2_171, mul_x2_172, mul_x2_173, 352 mul_x2_174, mul_x2_175, mul_x2_176, mul_x2_177, mul_x2_178, mul_x2_179, 353 mul_x2_180, mul_x2_181, mul_x2_182, mul_x2_183, mul_x2_184, mul_x2_185, 354 mul_x2_186, mul_x2_187, mul_x2_188, mul_x2_189, mul_x2_190, mul_x2_191, 355 mul_x2_192, mul_x2_193, mul_x2_194, mul_x2_195, mul_x2_196, mul_x2_197, 356 mul_x2_198, mul_x2_199, mul_x2_200, mul_x2_201, mul_x2_202, mul_x2_203, 357 mul_x2_204, mul_x2_205, mul_x2_206, mul_x2_207, mul_x2_208, mul_x2_209, 358 mul_x2_210, mul_x2_211, mul_x2_212, mul_x2_213, mul_x2_214, mul_x2_215, 359 mul_x2_216, mul_x2_217, mul_x2_218, mul_x2_219, mul_x2_220, mul_x2_221, 360 mul_x2_222, mul_x2_223, mul_x2_224, mul_x2_225, mul_x2_226, mul_x2_227, 361 mul_x2_228, mul_x2_229, mul_x2_230, mul_x2_231, mul_x2_232, mul_x2_233, 362 mul_x2_234, mul_x2_235, mul_x2_236, mul_x2_237, mul_x2_238, mul_x2_239, 363 mul_x2_240, mul_x2_241, mul_x2_242, mul_x2_243, mul_x2_244, mul_x2_245, 364 mul_x2_246, mul_x2_247, mul_x2_248, mul_x2_249, mul_x2_250, mul_x2_251, 365 mul_x2_252, mul_x2_253, mul_x2_254, mul_x2_255 366 }; 367 368 #define MUL(c, r...) \ 369 { \ 370 switch (REG_CNT(r)) { \ 371 case 4: \ 372 COPY(R_01(r), _mul_x2_in); \ 373 gf_x2_mul_fns[c](); \ 374 COPY(_mul_x2_acc, R_01(r)); \ 375 COPY(R_23(r), _mul_x2_in); \ 376 gf_x2_mul_fns[c](); \ 377 COPY(_mul_x2_acc, R_23(r)); \ 378 break; \ 379 default: \ 380 VERIFY(0); \ 381 } \ 382 } 383 384 385 #define raidz_math_begin() kfpu_begin() 386 #define raidz_math_end() kfpu_end() 387 388 389 #define SYN_STRIDE 4 390 391 #define ZERO_STRIDE 4 392 #define ZERO_DEFINE() {} 393 #define ZERO_D 0, 1, 2, 3 394 395 #define COPY_STRIDE 4 396 #define COPY_DEFINE() {} 397 #define COPY_D 0, 1, 2, 3 398 399 #define ADD_STRIDE 4 400 #define ADD_DEFINE() {} 401 #define ADD_D 0, 1, 2, 3 402 403 #define MUL_STRIDE 4 404 #define MUL_DEFINE() MUL2_SETUP() 405 #define MUL_D 0, 1, 2, 3 406 407 #define GEN_P_STRIDE 4 408 #define GEN_P_DEFINE() {} 409 #define GEN_P_P 0, 1, 2, 3 410 411 #define GEN_PQ_STRIDE 4 412 #define GEN_PQ_DEFINE() {} 413 #define GEN_PQ_D 0, 1, 2, 3 414 #define GEN_PQ_C 4, 5, 6, 7 415 416 #define GEN_PQR_STRIDE 4 417 #define GEN_PQR_DEFINE() {} 418 #define GEN_PQR_D 0, 1, 2, 3 419 #define GEN_PQR_C 4, 5, 6, 7 420 421 #define SYN_Q_DEFINE() {} 422 #define SYN_Q_D 0, 1, 2, 3 423 #define SYN_Q_X 4, 5, 6, 7 424 425 #define SYN_R_DEFINE() {} 426 #define SYN_R_D 0, 1, 2, 3 427 #define SYN_R_X 4, 5, 6, 7 428 429 #define SYN_PQ_DEFINE() {} 430 #define SYN_PQ_D 0, 1, 2, 3 431 #define SYN_PQ_X 4, 5, 6, 7 432 433 #define REC_PQ_STRIDE 4 434 #define REC_PQ_DEFINE() MUL2_SETUP() 435 #define REC_PQ_X 0, 1, 2, 3 436 #define REC_PQ_Y 4, 5, 6, 7 437 #define REC_PQ_T 8, 9, 10, 11 438 439 #define SYN_PR_DEFINE() {} 440 #define SYN_PR_D 0, 1, 2, 3 441 #define SYN_PR_X 4, 5, 6, 7 442 443 #define REC_PR_STRIDE 4 444 #define REC_PR_DEFINE() MUL2_SETUP() 445 #define REC_PR_X 0, 1, 2, 3 446 #define REC_PR_Y 4, 5, 6, 7 447 #define REC_PR_T 8, 9, 10, 11 448 449 #define SYN_QR_DEFINE() {} 450 #define SYN_QR_D 0, 1, 2, 3 451 #define SYN_QR_X 4, 5, 6, 7 452 453 #define REC_QR_STRIDE 4 454 #define REC_QR_DEFINE() MUL2_SETUP() 455 #define REC_QR_X 0, 1, 2, 3 456 #define REC_QR_Y 4, 5, 6, 7 457 #define REC_QR_T 8, 9, 10, 11 458 459 #define SYN_PQR_DEFINE() {} 460 #define SYN_PQR_D 0, 1, 2, 3 461 #define SYN_PQR_X 4, 5, 6, 7 462 463 #define REC_PQR_STRIDE 4 464 #define REC_PQR_DEFINE() MUL2_SETUP() 465 #define REC_PQR_X 0, 1, 2, 3 466 #define REC_PQR_Y 4, 5, 6, 7 467 #define REC_PQR_Z 8, 9, 10, 11 468 #define REC_PQR_XS 12, 13, 14, 15 469 #define REC_PQR_YS 16, 17, 18, 19 470 471 472 #include <sys/vdev_raidz_impl.h> 473 #include "vdev_raidz_math_impl.h" 474 475 DEFINE_GEN_METHODS(avx512f); 476 DEFINE_REC_METHODS(avx512f); 477 478 static boolean_t 479 raidz_will_avx512f_work(void) 480 { 481 return (kfpu_allowed() && zfs_avx_available() && 482 zfs_avx2_available() && zfs_avx512f_available()); 483 } 484 485 const raidz_impl_ops_t vdev_raidz_avx512f_impl = { 486 .init = NULL, 487 .fini = NULL, 488 .gen = RAIDZ_GEN_METHODS(avx512f), 489 .rec = RAIDZ_REC_METHODS(avx512f), 490 .is_supported = &raidz_will_avx512f_work, 491 .name = "avx512f" 492 }; 493 494 #endif /* defined(__x86_64) && defined(HAVE_AVX512F) */ 495