1 /* 2 * M-profile MVE Operations 3 * 4 * Copyright (c) 2021 Linaro, Ltd. 5 * 6 * This library is free software; you can redistribute it and/or 7 * modify it under the terms of the GNU Lesser General Public 8 * License as published by the Free Software Foundation; either 9 * version 2.1 of the License, or (at your option) any later version. 10 * 11 * This library is distributed in the hope that it will be useful, 12 * but WITHOUT ANY WARRANTY; without even the implied warranty of 13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 14 * Lesser General Public License for more details. 15 * 16 * You should have received a copy of the GNU Lesser General Public 17 * License along with this library; if not, see <http://www.gnu.org/licenses/>. 18 */ 19 20 #include "qemu/osdep.h" 21 #include "qemu/int128.h" 22 #include "cpu.h" 23 #include "internals.h" 24 #include "vec_internal.h" 25 #include "exec/helper-proto.h" 26 #include "exec/cpu_ldst.h" 27 #include "exec/exec-all.h" 28 #include "tcg/tcg.h" 29 30 static uint16_t mve_element_mask(CPUARMState *env) 31 { 32 /* 33 * Return the mask of which elements in the MVE vector should be 34 * updated. This is a combination of multiple things: 35 * (1) by default, we update every lane in the vector 36 * (2) VPT predication stores its state in the VPR register; 37 * (3) low-overhead-branch tail predication will mask out part 38 * the vector on the final iteration of the loop 39 * (4) if EPSR.ECI is set then we must execute only some beats 40 * of the insn 41 * We combine all these into a 16-bit result with the same semantics 42 * as VPR.P0: 0 to mask the lane, 1 if it is active. 43 * 8-bit vector ops will look at all bits of the result; 44 * 16-bit ops will look at bits 0, 2, 4, ...; 45 * 32-bit ops will look at bits 0, 4, 8 and 12. 46 * Compare pseudocode GetCurInstrBeat(), though that only returns 47 * the 4-bit slice of the mask corresponding to a single beat. 48 */ 49 uint16_t mask = FIELD_EX32(env->v7m.vpr, V7M_VPR, P0); 50 51 if (!(env->v7m.vpr & R_V7M_VPR_MASK01_MASK)) { 52 mask |= 0xff; 53 } 54 if (!(env->v7m.vpr & R_V7M_VPR_MASK23_MASK)) { 55 mask |= 0xff00; 56 } 57 58 if (env->v7m.ltpsize < 4 && 59 env->regs[14] <= (1 << (4 - env->v7m.ltpsize))) { 60 /* 61 * Tail predication active, and this is the last loop iteration. 62 * The element size is (1 << ltpsize), and we only want to process 63 * loopcount elements, so we want to retain the least significant 64 * (loopcount * esize) predicate bits and zero out bits above that. 65 */ 66 int masklen = env->regs[14] << env->v7m.ltpsize; 67 assert(masklen <= 16); 68 mask &= MAKE_64BIT_MASK(0, masklen); 69 } 70 71 if ((env->condexec_bits & 0xf) == 0) { 72 /* 73 * ECI bits indicate which beats are already executed; 74 * we handle this by effectively predicating them out. 75 */ 76 int eci = env->condexec_bits >> 4; 77 switch (eci) { 78 case ECI_NONE: 79 break; 80 case ECI_A0: 81 mask &= 0xfff0; 82 break; 83 case ECI_A0A1: 84 mask &= 0xff00; 85 break; 86 case ECI_A0A1A2: 87 case ECI_A0A1A2B0: 88 mask &= 0xf000; 89 break; 90 default: 91 g_assert_not_reached(); 92 } 93 } 94 95 return mask; 96 } 97 98 static void mve_advance_vpt(CPUARMState *env) 99 { 100 /* Advance the VPT and ECI state if necessary */ 101 uint32_t vpr = env->v7m.vpr; 102 unsigned mask01, mask23; 103 104 if ((env->condexec_bits & 0xf) == 0) { 105 env->condexec_bits = (env->condexec_bits == (ECI_A0A1A2B0 << 4)) ? 106 (ECI_A0 << 4) : (ECI_NONE << 4); 107 } 108 109 if (!(vpr & (R_V7M_VPR_MASK01_MASK | R_V7M_VPR_MASK23_MASK))) { 110 /* VPT not enabled, nothing to do */ 111 return; 112 } 113 114 mask01 = FIELD_EX32(vpr, V7M_VPR, MASK01); 115 mask23 = FIELD_EX32(vpr, V7M_VPR, MASK23); 116 if (mask01 > 8) { 117 /* high bit set, but not 0b1000: invert the relevant half of P0 */ 118 vpr ^= 0xff; 119 } 120 if (mask23 > 8) { 121 /* high bit set, but not 0b1000: invert the relevant half of P0 */ 122 vpr ^= 0xff00; 123 } 124 vpr = FIELD_DP32(vpr, V7M_VPR, MASK01, mask01 << 1); 125 vpr = FIELD_DP32(vpr, V7M_VPR, MASK23, mask23 << 1); 126 env->v7m.vpr = vpr; 127 } 128 129 130 #define DO_VLDR(OP, MSIZE, LDTYPE, ESIZE, TYPE) \ 131 void HELPER(mve_##OP)(CPUARMState *env, void *vd, uint32_t addr) \ 132 { \ 133 TYPE *d = vd; \ 134 uint16_t mask = mve_element_mask(env); \ 135 unsigned b, e; \ 136 /* \ 137 * R_SXTM allows the dest reg to become UNKNOWN for abandoned \ 138 * beats so we don't care if we update part of the dest and \ 139 * then take an exception. \ 140 */ \ 141 for (b = 0, e = 0; b < 16; b += ESIZE, e++) { \ 142 if (mask & (1 << b)) { \ 143 d[H##ESIZE(e)] = cpu_##LDTYPE##_data_ra(env, addr, GETPC()); \ 144 } \ 145 addr += MSIZE; \ 146 } \ 147 mve_advance_vpt(env); \ 148 } 149 150 #define DO_VSTR(OP, MSIZE, STTYPE, ESIZE, TYPE) \ 151 void HELPER(mve_##OP)(CPUARMState *env, void *vd, uint32_t addr) \ 152 { \ 153 TYPE *d = vd; \ 154 uint16_t mask = mve_element_mask(env); \ 155 unsigned b, e; \ 156 for (b = 0, e = 0; b < 16; b += ESIZE, e++) { \ 157 if (mask & (1 << b)) { \ 158 cpu_##STTYPE##_data_ra(env, addr, d[H##ESIZE(e)], GETPC()); \ 159 } \ 160 addr += MSIZE; \ 161 } \ 162 mve_advance_vpt(env); \ 163 } 164 165 DO_VLDR(vldrb, 1, ldub, 1, uint8_t) 166 DO_VLDR(vldrh, 2, lduw, 2, uint16_t) 167 DO_VLDR(vldrw, 4, ldl, 4, uint32_t) 168 169 DO_VSTR(vstrb, 1, stb, 1, uint8_t) 170 DO_VSTR(vstrh, 2, stw, 2, uint16_t) 171 DO_VSTR(vstrw, 4, stl, 4, uint32_t) 172 173 DO_VLDR(vldrb_sh, 1, ldsb, 2, int16_t) 174 DO_VLDR(vldrb_sw, 1, ldsb, 4, int32_t) 175 DO_VLDR(vldrb_uh, 1, ldub, 2, uint16_t) 176 DO_VLDR(vldrb_uw, 1, ldub, 4, uint32_t) 177 DO_VLDR(vldrh_sw, 2, ldsw, 4, int32_t) 178 DO_VLDR(vldrh_uw, 2, lduw, 4, uint32_t) 179 180 DO_VSTR(vstrb_h, 1, stb, 2, int16_t) 181 DO_VSTR(vstrb_w, 1, stb, 4, int32_t) 182 DO_VSTR(vstrh_w, 2, stw, 4, int32_t) 183 184 #undef DO_VLDR 185 #undef DO_VSTR 186 187 /* 188 * The mergemask(D, R, M) macro performs the operation "*D = R" but 189 * storing only the bytes which correspond to 1 bits in M, 190 * leaving other bytes in *D unchanged. We use _Generic 191 * to select the correct implementation based on the type of D. 192 */ 193 194 static void mergemask_ub(uint8_t *d, uint8_t r, uint16_t mask) 195 { 196 if (mask & 1) { 197 *d = r; 198 } 199 } 200 201 static void mergemask_sb(int8_t *d, int8_t r, uint16_t mask) 202 { 203 mergemask_ub((uint8_t *)d, r, mask); 204 } 205 206 static void mergemask_uh(uint16_t *d, uint16_t r, uint16_t mask) 207 { 208 uint16_t bmask = expand_pred_b_data[mask & 3]; 209 *d = (*d & ~bmask) | (r & bmask); 210 } 211 212 static void mergemask_sh(int16_t *d, int16_t r, uint16_t mask) 213 { 214 mergemask_uh((uint16_t *)d, r, mask); 215 } 216 217 static void mergemask_uw(uint32_t *d, uint32_t r, uint16_t mask) 218 { 219 uint32_t bmask = expand_pred_b_data[mask & 0xf]; 220 *d = (*d & ~bmask) | (r & bmask); 221 } 222 223 static void mergemask_sw(int32_t *d, int32_t r, uint16_t mask) 224 { 225 mergemask_uw((uint32_t *)d, r, mask); 226 } 227 228 static void mergemask_uq(uint64_t *d, uint64_t r, uint16_t mask) 229 { 230 uint64_t bmask = expand_pred_b_data[mask & 0xff]; 231 *d = (*d & ~bmask) | (r & bmask); 232 } 233 234 static void mergemask_sq(int64_t *d, int64_t r, uint16_t mask) 235 { 236 mergemask_uq((uint64_t *)d, r, mask); 237 } 238 239 #define mergemask(D, R, M) \ 240 _Generic(D, \ 241 uint8_t *: mergemask_ub, \ 242 int8_t *: mergemask_sb, \ 243 uint16_t *: mergemask_uh, \ 244 int16_t *: mergemask_sh, \ 245 uint32_t *: mergemask_uw, \ 246 int32_t *: mergemask_sw, \ 247 uint64_t *: mergemask_uq, \ 248 int64_t *: mergemask_sq)(D, R, M) 249 250 void HELPER(mve_vdup)(CPUARMState *env, void *vd, uint32_t val) 251 { 252 /* 253 * The generated code already replicated an 8 or 16 bit constant 254 * into the 32-bit value, so we only need to write the 32-bit 255 * value to all elements of the Qreg, allowing for predication. 256 */ 257 uint32_t *d = vd; 258 uint16_t mask = mve_element_mask(env); 259 unsigned e; 260 for (e = 0; e < 16 / 4; e++, mask >>= 4) { 261 mergemask(&d[H4(e)], val, mask); 262 } 263 mve_advance_vpt(env); 264 } 265 266 #define DO_1OP(OP, ESIZE, TYPE, FN) \ 267 void HELPER(mve_##OP)(CPUARMState *env, void *vd, void *vm) \ 268 { \ 269 TYPE *d = vd, *m = vm; \ 270 uint16_t mask = mve_element_mask(env); \ 271 unsigned e; \ 272 for (e = 0; e < 16 / ESIZE; e++, mask >>= ESIZE) { \ 273 mergemask(&d[H##ESIZE(e)], FN(m[H##ESIZE(e)]), mask); \ 274 } \ 275 mve_advance_vpt(env); \ 276 } 277 278 #define DO_CLS_B(N) (clrsb32(N) - 24) 279 #define DO_CLS_H(N) (clrsb32(N) - 16) 280 281 DO_1OP(vclsb, 1, int8_t, DO_CLS_B) 282 DO_1OP(vclsh, 2, int16_t, DO_CLS_H) 283 DO_1OP(vclsw, 4, int32_t, clrsb32) 284 285 #define DO_CLZ_B(N) (clz32(N) - 24) 286 #define DO_CLZ_H(N) (clz32(N) - 16) 287 288 DO_1OP(vclzb, 1, uint8_t, DO_CLZ_B) 289 DO_1OP(vclzh, 2, uint16_t, DO_CLZ_H) 290 DO_1OP(vclzw, 4, uint32_t, clz32) 291 292 DO_1OP(vrev16b, 2, uint16_t, bswap16) 293 DO_1OP(vrev32b, 4, uint32_t, bswap32) 294 DO_1OP(vrev32h, 4, uint32_t, hswap32) 295 DO_1OP(vrev64b, 8, uint64_t, bswap64) 296 DO_1OP(vrev64h, 8, uint64_t, hswap64) 297 DO_1OP(vrev64w, 8, uint64_t, wswap64) 298 299 #define DO_NOT(N) (~(N)) 300 301 DO_1OP(vmvn, 8, uint64_t, DO_NOT) 302 303 #define DO_ABS(N) ((N) < 0 ? -(N) : (N)) 304 #define DO_FABSH(N) ((N) & dup_const(MO_16, 0x7fff)) 305 #define DO_FABSS(N) ((N) & dup_const(MO_32, 0x7fffffff)) 306 307 DO_1OP(vabsb, 1, int8_t, DO_ABS) 308 DO_1OP(vabsh, 2, int16_t, DO_ABS) 309 DO_1OP(vabsw, 4, int32_t, DO_ABS) 310 311 /* We can do these 64 bits at a time */ 312 DO_1OP(vfabsh, 8, uint64_t, DO_FABSH) 313 DO_1OP(vfabss, 8, uint64_t, DO_FABSS) 314 315 #define DO_NEG(N) (-(N)) 316 #define DO_FNEGH(N) ((N) ^ dup_const(MO_16, 0x8000)) 317 #define DO_FNEGS(N) ((N) ^ dup_const(MO_32, 0x80000000)) 318 319 DO_1OP(vnegb, 1, int8_t, DO_NEG) 320 DO_1OP(vnegh, 2, int16_t, DO_NEG) 321 DO_1OP(vnegw, 4, int32_t, DO_NEG) 322 323 /* We can do these 64 bits at a time */ 324 DO_1OP(vfnegh, 8, uint64_t, DO_FNEGH) 325 DO_1OP(vfnegs, 8, uint64_t, DO_FNEGS) 326 327 #define DO_2OP(OP, ESIZE, TYPE, FN) \ 328 void HELPER(glue(mve_, OP))(CPUARMState *env, \ 329 void *vd, void *vn, void *vm) \ 330 { \ 331 TYPE *d = vd, *n = vn, *m = vm; \ 332 uint16_t mask = mve_element_mask(env); \ 333 unsigned e; \ 334 for (e = 0; e < 16 / ESIZE; e++, mask >>= ESIZE) { \ 335 mergemask(&d[H##ESIZE(e)], \ 336 FN(n[H##ESIZE(e)], m[H##ESIZE(e)]), mask); \ 337 } \ 338 mve_advance_vpt(env); \ 339 } 340 341 /* provide unsigned 2-op helpers for all sizes */ 342 #define DO_2OP_U(OP, FN) \ 343 DO_2OP(OP##b, 1, uint8_t, FN) \ 344 DO_2OP(OP##h, 2, uint16_t, FN) \ 345 DO_2OP(OP##w, 4, uint32_t, FN) 346 347 /* provide signed 2-op helpers for all sizes */ 348 #define DO_2OP_S(OP, FN) \ 349 DO_2OP(OP##b, 1, int8_t, FN) \ 350 DO_2OP(OP##h, 2, int16_t, FN) \ 351 DO_2OP(OP##w, 4, int32_t, FN) 352 353 /* 354 * "Long" operations where two half-sized inputs (taken from either the 355 * top or the bottom of the input vector) produce a double-width result. 356 * Here ESIZE, TYPE are for the input, and LESIZE, LTYPE for the output. 357 */ 358 #define DO_2OP_L(OP, TOP, ESIZE, TYPE, LESIZE, LTYPE, FN) \ 359 void HELPER(glue(mve_, OP))(CPUARMState *env, void *vd, void *vn, void *vm) \ 360 { \ 361 LTYPE *d = vd; \ 362 TYPE *n = vn, *m = vm; \ 363 uint16_t mask = mve_element_mask(env); \ 364 unsigned le; \ 365 for (le = 0; le < 16 / LESIZE; le++, mask >>= LESIZE) { \ 366 LTYPE r = FN((LTYPE)n[H##ESIZE(le * 2 + TOP)], \ 367 m[H##ESIZE(le * 2 + TOP)]); \ 368 mergemask(&d[H##LESIZE(le)], r, mask); \ 369 } \ 370 mve_advance_vpt(env); \ 371 } 372 373 #define DO_2OP_SAT(OP, ESIZE, TYPE, FN) \ 374 void HELPER(glue(mve_, OP))(CPUARMState *env, void *vd, void *vn, void *vm) \ 375 { \ 376 TYPE *d = vd, *n = vn, *m = vm; \ 377 uint16_t mask = mve_element_mask(env); \ 378 unsigned e; \ 379 bool qc = false; \ 380 for (e = 0; e < 16 / ESIZE; e++, mask >>= ESIZE) { \ 381 bool sat = false; \ 382 TYPE r = FN(n[H##ESIZE(e)], m[H##ESIZE(e)], &sat); \ 383 mergemask(&d[H##ESIZE(e)], r, mask); \ 384 qc |= sat & mask & 1; \ 385 } \ 386 if (qc) { \ 387 env->vfp.qc[0] = qc; \ 388 } \ 389 mve_advance_vpt(env); \ 390 } 391 392 /* provide unsigned 2-op helpers for all sizes */ 393 #define DO_2OP_SAT_U(OP, FN) \ 394 DO_2OP_SAT(OP##b, 1, uint8_t, FN) \ 395 DO_2OP_SAT(OP##h, 2, uint16_t, FN) \ 396 DO_2OP_SAT(OP##w, 4, uint32_t, FN) 397 398 /* provide signed 2-op helpers for all sizes */ 399 #define DO_2OP_SAT_S(OP, FN) \ 400 DO_2OP_SAT(OP##b, 1, int8_t, FN) \ 401 DO_2OP_SAT(OP##h, 2, int16_t, FN) \ 402 DO_2OP_SAT(OP##w, 4, int32_t, FN) 403 404 #define DO_AND(N, M) ((N) & (M)) 405 #define DO_BIC(N, M) ((N) & ~(M)) 406 #define DO_ORR(N, M) ((N) | (M)) 407 #define DO_ORN(N, M) ((N) | ~(M)) 408 #define DO_EOR(N, M) ((N) ^ (M)) 409 410 DO_2OP(vand, 8, uint64_t, DO_AND) 411 DO_2OP(vbic, 8, uint64_t, DO_BIC) 412 DO_2OP(vorr, 8, uint64_t, DO_ORR) 413 DO_2OP(vorn, 8, uint64_t, DO_ORN) 414 DO_2OP(veor, 8, uint64_t, DO_EOR) 415 416 #define DO_ADD(N, M) ((N) + (M)) 417 #define DO_SUB(N, M) ((N) - (M)) 418 #define DO_MUL(N, M) ((N) * (M)) 419 420 DO_2OP_U(vadd, DO_ADD) 421 DO_2OP_U(vsub, DO_SUB) 422 DO_2OP_U(vmul, DO_MUL) 423 424 DO_2OP_L(vmullbsb, 0, 1, int8_t, 2, int16_t, DO_MUL) 425 DO_2OP_L(vmullbsh, 0, 2, int16_t, 4, int32_t, DO_MUL) 426 DO_2OP_L(vmullbsw, 0, 4, int32_t, 8, int64_t, DO_MUL) 427 DO_2OP_L(vmullbub, 0, 1, uint8_t, 2, uint16_t, DO_MUL) 428 DO_2OP_L(vmullbuh, 0, 2, uint16_t, 4, uint32_t, DO_MUL) 429 DO_2OP_L(vmullbuw, 0, 4, uint32_t, 8, uint64_t, DO_MUL) 430 431 DO_2OP_L(vmulltsb, 1, 1, int8_t, 2, int16_t, DO_MUL) 432 DO_2OP_L(vmulltsh, 1, 2, int16_t, 4, int32_t, DO_MUL) 433 DO_2OP_L(vmulltsw, 1, 4, int32_t, 8, int64_t, DO_MUL) 434 DO_2OP_L(vmulltub, 1, 1, uint8_t, 2, uint16_t, DO_MUL) 435 DO_2OP_L(vmulltuh, 1, 2, uint16_t, 4, uint32_t, DO_MUL) 436 DO_2OP_L(vmulltuw, 1, 4, uint32_t, 8, uint64_t, DO_MUL) 437 438 /* 439 * Because the computation type is at least twice as large as required, 440 * these work for both signed and unsigned source types. 441 */ 442 static inline uint8_t do_mulh_b(int32_t n, int32_t m) 443 { 444 return (n * m) >> 8; 445 } 446 447 static inline uint16_t do_mulh_h(int32_t n, int32_t m) 448 { 449 return (n * m) >> 16; 450 } 451 452 static inline uint32_t do_mulh_w(int64_t n, int64_t m) 453 { 454 return (n * m) >> 32; 455 } 456 457 static inline uint8_t do_rmulh_b(int32_t n, int32_t m) 458 { 459 return (n * m + (1U << 7)) >> 8; 460 } 461 462 static inline uint16_t do_rmulh_h(int32_t n, int32_t m) 463 { 464 return (n * m + (1U << 15)) >> 16; 465 } 466 467 static inline uint32_t do_rmulh_w(int64_t n, int64_t m) 468 { 469 return (n * m + (1U << 31)) >> 32; 470 } 471 472 DO_2OP(vmulhsb, 1, int8_t, do_mulh_b) 473 DO_2OP(vmulhsh, 2, int16_t, do_mulh_h) 474 DO_2OP(vmulhsw, 4, int32_t, do_mulh_w) 475 DO_2OP(vmulhub, 1, uint8_t, do_mulh_b) 476 DO_2OP(vmulhuh, 2, uint16_t, do_mulh_h) 477 DO_2OP(vmulhuw, 4, uint32_t, do_mulh_w) 478 479 DO_2OP(vrmulhsb, 1, int8_t, do_rmulh_b) 480 DO_2OP(vrmulhsh, 2, int16_t, do_rmulh_h) 481 DO_2OP(vrmulhsw, 4, int32_t, do_rmulh_w) 482 DO_2OP(vrmulhub, 1, uint8_t, do_rmulh_b) 483 DO_2OP(vrmulhuh, 2, uint16_t, do_rmulh_h) 484 DO_2OP(vrmulhuw, 4, uint32_t, do_rmulh_w) 485 486 #define DO_MAX(N, M) ((N) >= (M) ? (N) : (M)) 487 #define DO_MIN(N, M) ((N) >= (M) ? (M) : (N)) 488 489 DO_2OP_S(vmaxs, DO_MAX) 490 DO_2OP_U(vmaxu, DO_MAX) 491 DO_2OP_S(vmins, DO_MIN) 492 DO_2OP_U(vminu, DO_MIN) 493 494 #define DO_ABD(N, M) ((N) >= (M) ? (N) - (M) : (M) - (N)) 495 496 DO_2OP_S(vabds, DO_ABD) 497 DO_2OP_U(vabdu, DO_ABD) 498 499 static inline uint32_t do_vhadd_u(uint32_t n, uint32_t m) 500 { 501 return ((uint64_t)n + m) >> 1; 502 } 503 504 static inline int32_t do_vhadd_s(int32_t n, int32_t m) 505 { 506 return ((int64_t)n + m) >> 1; 507 } 508 509 static inline uint32_t do_vhsub_u(uint32_t n, uint32_t m) 510 { 511 return ((uint64_t)n - m) >> 1; 512 } 513 514 static inline int32_t do_vhsub_s(int32_t n, int32_t m) 515 { 516 return ((int64_t)n - m) >> 1; 517 } 518 519 DO_2OP_S(vhadds, do_vhadd_s) 520 DO_2OP_U(vhaddu, do_vhadd_u) 521 DO_2OP_S(vhsubs, do_vhsub_s) 522 DO_2OP_U(vhsubu, do_vhsub_u) 523 524 #define DO_VSHLS(N, M) do_sqrshl_bhs(N, (int8_t)(M), sizeof(N) * 8, false, NULL) 525 #define DO_VSHLU(N, M) do_uqrshl_bhs(N, (int8_t)(M), sizeof(N) * 8, false, NULL) 526 #define DO_VRSHLS(N, M) do_sqrshl_bhs(N, (int8_t)(M), sizeof(N) * 8, true, NULL) 527 #define DO_VRSHLU(N, M) do_uqrshl_bhs(N, (int8_t)(M), sizeof(N) * 8, true, NULL) 528 529 DO_2OP_S(vshls, DO_VSHLS) 530 DO_2OP_U(vshlu, DO_VSHLU) 531 DO_2OP_S(vrshls, DO_VRSHLS) 532 DO_2OP_U(vrshlu, DO_VRSHLU) 533 534 #define DO_RHADD_S(N, M) (((int64_t)(N) + (M) + 1) >> 1) 535 #define DO_RHADD_U(N, M) (((uint64_t)(N) + (M) + 1) >> 1) 536 537 DO_2OP_S(vrhadds, DO_RHADD_S) 538 DO_2OP_U(vrhaddu, DO_RHADD_U) 539 540 static inline int32_t do_sat_bhw(int64_t val, int64_t min, int64_t max, bool *s) 541 { 542 if (val > max) { 543 *s = true; 544 return max; 545 } else if (val < min) { 546 *s = true; 547 return min; 548 } 549 return val; 550 } 551 552 #define DO_SQADD_B(n, m, s) do_sat_bhw((int64_t)n + m, INT8_MIN, INT8_MAX, s) 553 #define DO_SQADD_H(n, m, s) do_sat_bhw((int64_t)n + m, INT16_MIN, INT16_MAX, s) 554 #define DO_SQADD_W(n, m, s) do_sat_bhw((int64_t)n + m, INT32_MIN, INT32_MAX, s) 555 556 #define DO_UQADD_B(n, m, s) do_sat_bhw((int64_t)n + m, 0, UINT8_MAX, s) 557 #define DO_UQADD_H(n, m, s) do_sat_bhw((int64_t)n + m, 0, UINT16_MAX, s) 558 #define DO_UQADD_W(n, m, s) do_sat_bhw((int64_t)n + m, 0, UINT32_MAX, s) 559 560 #define DO_SQSUB_B(n, m, s) do_sat_bhw((int64_t)n - m, INT8_MIN, INT8_MAX, s) 561 #define DO_SQSUB_H(n, m, s) do_sat_bhw((int64_t)n - m, INT16_MIN, INT16_MAX, s) 562 #define DO_SQSUB_W(n, m, s) do_sat_bhw((int64_t)n - m, INT32_MIN, INT32_MAX, s) 563 564 #define DO_UQSUB_B(n, m, s) do_sat_bhw((int64_t)n - m, 0, UINT8_MAX, s) 565 #define DO_UQSUB_H(n, m, s) do_sat_bhw((int64_t)n - m, 0, UINT16_MAX, s) 566 #define DO_UQSUB_W(n, m, s) do_sat_bhw((int64_t)n - m, 0, UINT32_MAX, s) 567 568 /* 569 * For QDMULH and QRDMULH we simplify "double and shift by esize" into 570 * "shift by esize-1", adjusting the QRDMULH rounding constant to match. 571 */ 572 #define DO_QDMULH_B(n, m, s) do_sat_bhw(((int64_t)n * m) >> 7, \ 573 INT8_MIN, INT8_MAX, s) 574 #define DO_QDMULH_H(n, m, s) do_sat_bhw(((int64_t)n * m) >> 15, \ 575 INT16_MIN, INT16_MAX, s) 576 #define DO_QDMULH_W(n, m, s) do_sat_bhw(((int64_t)n * m) >> 31, \ 577 INT32_MIN, INT32_MAX, s) 578 579 #define DO_QRDMULH_B(n, m, s) do_sat_bhw(((int64_t)n * m + (1 << 6)) >> 7, \ 580 INT8_MIN, INT8_MAX, s) 581 #define DO_QRDMULH_H(n, m, s) do_sat_bhw(((int64_t)n * m + (1 << 14)) >> 15, \ 582 INT16_MIN, INT16_MAX, s) 583 #define DO_QRDMULH_W(n, m, s) do_sat_bhw(((int64_t)n * m + (1 << 30)) >> 31, \ 584 INT32_MIN, INT32_MAX, s) 585 586 DO_2OP_SAT(vqdmulhb, 1, int8_t, DO_QDMULH_B) 587 DO_2OP_SAT(vqdmulhh, 2, int16_t, DO_QDMULH_H) 588 DO_2OP_SAT(vqdmulhw, 4, int32_t, DO_QDMULH_W) 589 590 DO_2OP_SAT(vqrdmulhb, 1, int8_t, DO_QRDMULH_B) 591 DO_2OP_SAT(vqrdmulhh, 2, int16_t, DO_QRDMULH_H) 592 DO_2OP_SAT(vqrdmulhw, 4, int32_t, DO_QRDMULH_W) 593 594 DO_2OP_SAT(vqaddub, 1, uint8_t, DO_UQADD_B) 595 DO_2OP_SAT(vqadduh, 2, uint16_t, DO_UQADD_H) 596 DO_2OP_SAT(vqadduw, 4, uint32_t, DO_UQADD_W) 597 DO_2OP_SAT(vqaddsb, 1, int8_t, DO_SQADD_B) 598 DO_2OP_SAT(vqaddsh, 2, int16_t, DO_SQADD_H) 599 DO_2OP_SAT(vqaddsw, 4, int32_t, DO_SQADD_W) 600 601 DO_2OP_SAT(vqsubub, 1, uint8_t, DO_UQSUB_B) 602 DO_2OP_SAT(vqsubuh, 2, uint16_t, DO_UQSUB_H) 603 DO_2OP_SAT(vqsubuw, 4, uint32_t, DO_UQSUB_W) 604 DO_2OP_SAT(vqsubsb, 1, int8_t, DO_SQSUB_B) 605 DO_2OP_SAT(vqsubsh, 2, int16_t, DO_SQSUB_H) 606 DO_2OP_SAT(vqsubsw, 4, int32_t, DO_SQSUB_W) 607 608 /* 609 * This wrapper fixes up the impedance mismatch between do_sqrshl_bhs() 610 * and friends wanting a uint32_t* sat and our needing a bool*. 611 */ 612 #define WRAP_QRSHL_HELPER(FN, N, M, ROUND, satp) \ 613 ({ \ 614 uint32_t su32 = 0; \ 615 typeof(N) r = FN(N, (int8_t)(M), sizeof(N) * 8, ROUND, &su32); \ 616 if (su32) { \ 617 *satp = true; \ 618 } \ 619 r; \ 620 }) 621 622 #define DO_SQSHL_OP(N, M, satp) \ 623 WRAP_QRSHL_HELPER(do_sqrshl_bhs, N, M, false, satp) 624 #define DO_UQSHL_OP(N, M, satp) \ 625 WRAP_QRSHL_HELPER(do_uqrshl_bhs, N, M, false, satp) 626 #define DO_SQRSHL_OP(N, M, satp) \ 627 WRAP_QRSHL_HELPER(do_sqrshl_bhs, N, M, true, satp) 628 #define DO_UQRSHL_OP(N, M, satp) \ 629 WRAP_QRSHL_HELPER(do_uqrshl_bhs, N, M, true, satp) 630 631 DO_2OP_SAT_S(vqshls, DO_SQSHL_OP) 632 DO_2OP_SAT_U(vqshlu, DO_UQSHL_OP) 633 DO_2OP_SAT_S(vqrshls, DO_SQRSHL_OP) 634 DO_2OP_SAT_U(vqrshlu, DO_UQRSHL_OP) 635 636 /* 637 * Multiply add dual returning high half 638 * The 'FN' here takes four inputs A, B, C, D, a 0/1 indicator of 639 * whether to add the rounding constant, and the pointer to the 640 * saturation flag, and should do "(A * B + C * D) * 2 + rounding constant", 641 * saturate to twice the input size and return the high half; or 642 * (A * B - C * D) etc for VQDMLSDH. 643 */ 644 #define DO_VQDMLADH_OP(OP, ESIZE, TYPE, XCHG, ROUND, FN) \ 645 void HELPER(glue(mve_, OP))(CPUARMState *env, void *vd, void *vn, \ 646 void *vm) \ 647 { \ 648 TYPE *d = vd, *n = vn, *m = vm; \ 649 uint16_t mask = mve_element_mask(env); \ 650 unsigned e; \ 651 bool qc = false; \ 652 for (e = 0; e < 16 / ESIZE; e++, mask >>= ESIZE) { \ 653 bool sat = false; \ 654 if ((e & 1) == XCHG) { \ 655 TYPE r = FN(n[H##ESIZE(e)], \ 656 m[H##ESIZE(e - XCHG)], \ 657 n[H##ESIZE(e + (1 - 2 * XCHG))], \ 658 m[H##ESIZE(e + (1 - XCHG))], \ 659 ROUND, &sat); \ 660 mergemask(&d[H##ESIZE(e)], r, mask); \ 661 qc |= sat & mask & 1; \ 662 } \ 663 } \ 664 if (qc) { \ 665 env->vfp.qc[0] = qc; \ 666 } \ 667 mve_advance_vpt(env); \ 668 } 669 670 static int8_t do_vqdmladh_b(int8_t a, int8_t b, int8_t c, int8_t d, 671 int round, bool *sat) 672 { 673 int64_t r = ((int64_t)a * b + (int64_t)c * d) * 2 + (round << 7); 674 return do_sat_bhw(r, INT16_MIN, INT16_MAX, sat) >> 8; 675 } 676 677 static int16_t do_vqdmladh_h(int16_t a, int16_t b, int16_t c, int16_t d, 678 int round, bool *sat) 679 { 680 int64_t r = ((int64_t)a * b + (int64_t)c * d) * 2 + (round << 15); 681 return do_sat_bhw(r, INT32_MIN, INT32_MAX, sat) >> 16; 682 } 683 684 static int32_t do_vqdmladh_w(int32_t a, int32_t b, int32_t c, int32_t d, 685 int round, bool *sat) 686 { 687 int64_t m1 = (int64_t)a * b; 688 int64_t m2 = (int64_t)c * d; 689 int64_t r; 690 /* 691 * Architecturally we should do the entire add, double, round 692 * and then check for saturation. We do three saturating adds, 693 * but we need to be careful about the order. If the first 694 * m1 + m2 saturates then it's impossible for the *2+rc to 695 * bring it back into the non-saturated range. However, if 696 * m1 + m2 is negative then it's possible that doing the doubling 697 * would take the intermediate result below INT64_MAX and the 698 * addition of the rounding constant then brings it back in range. 699 * So we add half the rounding constant before doubling rather 700 * than adding the rounding constant after the doubling. 701 */ 702 if (sadd64_overflow(m1, m2, &r) || 703 sadd64_overflow(r, (round << 30), &r) || 704 sadd64_overflow(r, r, &r)) { 705 *sat = true; 706 return r < 0 ? INT32_MAX : INT32_MIN; 707 } 708 return r >> 32; 709 } 710 711 static int8_t do_vqdmlsdh_b(int8_t a, int8_t b, int8_t c, int8_t d, 712 int round, bool *sat) 713 { 714 int64_t r = ((int64_t)a * b - (int64_t)c * d) * 2 + (round << 7); 715 return do_sat_bhw(r, INT16_MIN, INT16_MAX, sat) >> 8; 716 } 717 718 static int16_t do_vqdmlsdh_h(int16_t a, int16_t b, int16_t c, int16_t d, 719 int round, bool *sat) 720 { 721 int64_t r = ((int64_t)a * b - (int64_t)c * d) * 2 + (round << 15); 722 return do_sat_bhw(r, INT32_MIN, INT32_MAX, sat) >> 16; 723 } 724 725 static int32_t do_vqdmlsdh_w(int32_t a, int32_t b, int32_t c, int32_t d, 726 int round, bool *sat) 727 { 728 int64_t m1 = (int64_t)a * b; 729 int64_t m2 = (int64_t)c * d; 730 int64_t r; 731 /* The same ordering issue as in do_vqdmladh_w applies here too */ 732 if (ssub64_overflow(m1, m2, &r) || 733 sadd64_overflow(r, (round << 30), &r) || 734 sadd64_overflow(r, r, &r)) { 735 *sat = true; 736 return r < 0 ? INT32_MAX : INT32_MIN; 737 } 738 return r >> 32; 739 } 740 741 DO_VQDMLADH_OP(vqdmladhb, 1, int8_t, 0, 0, do_vqdmladh_b) 742 DO_VQDMLADH_OP(vqdmladhh, 2, int16_t, 0, 0, do_vqdmladh_h) 743 DO_VQDMLADH_OP(vqdmladhw, 4, int32_t, 0, 0, do_vqdmladh_w) 744 DO_VQDMLADH_OP(vqdmladhxb, 1, int8_t, 1, 0, do_vqdmladh_b) 745 DO_VQDMLADH_OP(vqdmladhxh, 2, int16_t, 1, 0, do_vqdmladh_h) 746 DO_VQDMLADH_OP(vqdmladhxw, 4, int32_t, 1, 0, do_vqdmladh_w) 747 748 DO_VQDMLADH_OP(vqrdmladhb, 1, int8_t, 0, 1, do_vqdmladh_b) 749 DO_VQDMLADH_OP(vqrdmladhh, 2, int16_t, 0, 1, do_vqdmladh_h) 750 DO_VQDMLADH_OP(vqrdmladhw, 4, int32_t, 0, 1, do_vqdmladh_w) 751 DO_VQDMLADH_OP(vqrdmladhxb, 1, int8_t, 1, 1, do_vqdmladh_b) 752 DO_VQDMLADH_OP(vqrdmladhxh, 2, int16_t, 1, 1, do_vqdmladh_h) 753 DO_VQDMLADH_OP(vqrdmladhxw, 4, int32_t, 1, 1, do_vqdmladh_w) 754 755 DO_VQDMLADH_OP(vqdmlsdhb, 1, int8_t, 0, 0, do_vqdmlsdh_b) 756 DO_VQDMLADH_OP(vqdmlsdhh, 2, int16_t, 0, 0, do_vqdmlsdh_h) 757 DO_VQDMLADH_OP(vqdmlsdhw, 4, int32_t, 0, 0, do_vqdmlsdh_w) 758 DO_VQDMLADH_OP(vqdmlsdhxb, 1, int8_t, 1, 0, do_vqdmlsdh_b) 759 DO_VQDMLADH_OP(vqdmlsdhxh, 2, int16_t, 1, 0, do_vqdmlsdh_h) 760 DO_VQDMLADH_OP(vqdmlsdhxw, 4, int32_t, 1, 0, do_vqdmlsdh_w) 761 762 DO_VQDMLADH_OP(vqrdmlsdhb, 1, int8_t, 0, 1, do_vqdmlsdh_b) 763 DO_VQDMLADH_OP(vqrdmlsdhh, 2, int16_t, 0, 1, do_vqdmlsdh_h) 764 DO_VQDMLADH_OP(vqrdmlsdhw, 4, int32_t, 0, 1, do_vqdmlsdh_w) 765 DO_VQDMLADH_OP(vqrdmlsdhxb, 1, int8_t, 1, 1, do_vqdmlsdh_b) 766 DO_VQDMLADH_OP(vqrdmlsdhxh, 2, int16_t, 1, 1, do_vqdmlsdh_h) 767 DO_VQDMLADH_OP(vqrdmlsdhxw, 4, int32_t, 1, 1, do_vqdmlsdh_w) 768 769 #define DO_2OP_SCALAR(OP, ESIZE, TYPE, FN) \ 770 void HELPER(glue(mve_, OP))(CPUARMState *env, void *vd, void *vn, \ 771 uint32_t rm) \ 772 { \ 773 TYPE *d = vd, *n = vn; \ 774 TYPE m = rm; \ 775 uint16_t mask = mve_element_mask(env); \ 776 unsigned e; \ 777 for (e = 0; e < 16 / ESIZE; e++, mask >>= ESIZE) { \ 778 mergemask(&d[H##ESIZE(e)], FN(n[H##ESIZE(e)], m), mask); \ 779 } \ 780 mve_advance_vpt(env); \ 781 } 782 783 #define DO_2OP_SAT_SCALAR(OP, ESIZE, TYPE, FN) \ 784 void HELPER(glue(mve_, OP))(CPUARMState *env, void *vd, void *vn, \ 785 uint32_t rm) \ 786 { \ 787 TYPE *d = vd, *n = vn; \ 788 TYPE m = rm; \ 789 uint16_t mask = mve_element_mask(env); \ 790 unsigned e; \ 791 bool qc = false; \ 792 for (e = 0; e < 16 / ESIZE; e++, mask >>= ESIZE) { \ 793 bool sat = false; \ 794 mergemask(&d[H##ESIZE(e)], FN(n[H##ESIZE(e)], m, &sat), \ 795 mask); \ 796 qc |= sat & mask & 1; \ 797 } \ 798 if (qc) { \ 799 env->vfp.qc[0] = qc; \ 800 } \ 801 mve_advance_vpt(env); \ 802 } 803 804 /* provide unsigned 2-op scalar helpers for all sizes */ 805 #define DO_2OP_SCALAR_U(OP, FN) \ 806 DO_2OP_SCALAR(OP##b, 1, uint8_t, FN) \ 807 DO_2OP_SCALAR(OP##h, 2, uint16_t, FN) \ 808 DO_2OP_SCALAR(OP##w, 4, uint32_t, FN) 809 #define DO_2OP_SCALAR_S(OP, FN) \ 810 DO_2OP_SCALAR(OP##b, 1, int8_t, FN) \ 811 DO_2OP_SCALAR(OP##h, 2, int16_t, FN) \ 812 DO_2OP_SCALAR(OP##w, 4, int32_t, FN) 813 814 DO_2OP_SCALAR_U(vadd_scalar, DO_ADD) 815 DO_2OP_SCALAR_U(vsub_scalar, DO_SUB) 816 DO_2OP_SCALAR_U(vmul_scalar, DO_MUL) 817 DO_2OP_SCALAR_S(vhadds_scalar, do_vhadd_s) 818 DO_2OP_SCALAR_U(vhaddu_scalar, do_vhadd_u) 819 DO_2OP_SCALAR_S(vhsubs_scalar, do_vhsub_s) 820 DO_2OP_SCALAR_U(vhsubu_scalar, do_vhsub_u) 821 822 DO_2OP_SAT_SCALAR(vqaddu_scalarb, 1, uint8_t, DO_UQADD_B) 823 DO_2OP_SAT_SCALAR(vqaddu_scalarh, 2, uint16_t, DO_UQADD_H) 824 DO_2OP_SAT_SCALAR(vqaddu_scalarw, 4, uint32_t, DO_UQADD_W) 825 DO_2OP_SAT_SCALAR(vqadds_scalarb, 1, int8_t, DO_SQADD_B) 826 DO_2OP_SAT_SCALAR(vqadds_scalarh, 2, int16_t, DO_SQADD_H) 827 DO_2OP_SAT_SCALAR(vqadds_scalarw, 4, int32_t, DO_SQADD_W) 828 829 DO_2OP_SAT_SCALAR(vqsubu_scalarb, 1, uint8_t, DO_UQSUB_B) 830 DO_2OP_SAT_SCALAR(vqsubu_scalarh, 2, uint16_t, DO_UQSUB_H) 831 DO_2OP_SAT_SCALAR(vqsubu_scalarw, 4, uint32_t, DO_UQSUB_W) 832 DO_2OP_SAT_SCALAR(vqsubs_scalarb, 1, int8_t, DO_SQSUB_B) 833 DO_2OP_SAT_SCALAR(vqsubs_scalarh, 2, int16_t, DO_SQSUB_H) 834 DO_2OP_SAT_SCALAR(vqsubs_scalarw, 4, int32_t, DO_SQSUB_W) 835 836 DO_2OP_SAT_SCALAR(vqdmulh_scalarb, 1, int8_t, DO_QDMULH_B) 837 DO_2OP_SAT_SCALAR(vqdmulh_scalarh, 2, int16_t, DO_QDMULH_H) 838 DO_2OP_SAT_SCALAR(vqdmulh_scalarw, 4, int32_t, DO_QDMULH_W) 839 DO_2OP_SAT_SCALAR(vqrdmulh_scalarb, 1, int8_t, DO_QRDMULH_B) 840 DO_2OP_SAT_SCALAR(vqrdmulh_scalarh, 2, int16_t, DO_QRDMULH_H) 841 DO_2OP_SAT_SCALAR(vqrdmulh_scalarw, 4, int32_t, DO_QRDMULH_W) 842 843 /* 844 * Long saturating scalar ops. As with DO_2OP_L, TYPE and H are for the 845 * input (smaller) type and LESIZE, LTYPE, LH for the output (long) type. 846 * SATMASK specifies which bits of the predicate mask matter for determining 847 * whether to propagate a saturation indication into FPSCR.QC -- for 848 * the 16x16->32 case we must check only the bit corresponding to the T or B 849 * half that we used, but for the 32x32->64 case we propagate if the mask 850 * bit is set for either half. 851 */ 852 #define DO_2OP_SAT_SCALAR_L(OP, TOP, ESIZE, TYPE, LESIZE, LTYPE, FN, SATMASK) \ 853 void HELPER(glue(mve_, OP))(CPUARMState *env, void *vd, void *vn, \ 854 uint32_t rm) \ 855 { \ 856 LTYPE *d = vd; \ 857 TYPE *n = vn; \ 858 TYPE m = rm; \ 859 uint16_t mask = mve_element_mask(env); \ 860 unsigned le; \ 861 bool qc = false; \ 862 for (le = 0; le < 16 / LESIZE; le++, mask >>= LESIZE) { \ 863 bool sat = false; \ 864 LTYPE r = FN((LTYPE)n[H##ESIZE(le * 2 + TOP)], m, &sat); \ 865 mergemask(&d[H##LESIZE(le)], r, mask); \ 866 qc |= sat && (mask & SATMASK); \ 867 } \ 868 if (qc) { \ 869 env->vfp.qc[0] = qc; \ 870 } \ 871 mve_advance_vpt(env); \ 872 } 873 874 static inline int32_t do_qdmullh(int16_t n, int16_t m, bool *sat) 875 { 876 int64_t r = ((int64_t)n * m) * 2; 877 return do_sat_bhw(r, INT32_MIN, INT32_MAX, sat); 878 } 879 880 static inline int64_t do_qdmullw(int32_t n, int32_t m, bool *sat) 881 { 882 /* The multiply can't overflow, but the doubling might */ 883 int64_t r = (int64_t)n * m; 884 if (r > INT64_MAX / 2) { 885 *sat = true; 886 return INT64_MAX; 887 } else if (r < INT64_MIN / 2) { 888 *sat = true; 889 return INT64_MIN; 890 } else { 891 return r * 2; 892 } 893 } 894 895 #define SATMASK16B 1 896 #define SATMASK16T (1 << 2) 897 #define SATMASK32 ((1 << 4) | 1) 898 899 DO_2OP_SAT_SCALAR_L(vqdmullb_scalarh, 0, 2, int16_t, 4, int32_t, \ 900 do_qdmullh, SATMASK16B) 901 DO_2OP_SAT_SCALAR_L(vqdmullb_scalarw, 0, 4, int32_t, 8, int64_t, \ 902 do_qdmullw, SATMASK32) 903 DO_2OP_SAT_SCALAR_L(vqdmullt_scalarh, 1, 2, int16_t, 4, int32_t, \ 904 do_qdmullh, SATMASK16T) 905 DO_2OP_SAT_SCALAR_L(vqdmullt_scalarw, 1, 4, int32_t, 8, int64_t, \ 906 do_qdmullw, SATMASK32) 907 908 /* 909 * Long saturating ops 910 */ 911 #define DO_2OP_SAT_L(OP, TOP, ESIZE, TYPE, LESIZE, LTYPE, FN, SATMASK) \ 912 void HELPER(glue(mve_, OP))(CPUARMState *env, void *vd, void *vn, \ 913 void *vm) \ 914 { \ 915 LTYPE *d = vd; \ 916 TYPE *n = vn, *m = vm; \ 917 uint16_t mask = mve_element_mask(env); \ 918 unsigned le; \ 919 bool qc = false; \ 920 for (le = 0; le < 16 / LESIZE; le++, mask >>= LESIZE) { \ 921 bool sat = false; \ 922 LTYPE op1 = n[H##ESIZE(le * 2 + TOP)]; \ 923 LTYPE op2 = m[H##ESIZE(le * 2 + TOP)]; \ 924 mergemask(&d[H##LESIZE(le)], FN(op1, op2, &sat), mask); \ 925 qc |= sat && (mask & SATMASK); \ 926 } \ 927 if (qc) { \ 928 env->vfp.qc[0] = qc; \ 929 } \ 930 mve_advance_vpt(env); \ 931 } 932 933 DO_2OP_SAT_L(vqdmullbh, 0, 2, int16_t, 4, int32_t, do_qdmullh, SATMASK16B) 934 DO_2OP_SAT_L(vqdmullbw, 0, 4, int32_t, 8, int64_t, do_qdmullw, SATMASK32) 935 DO_2OP_SAT_L(vqdmullth, 1, 2, int16_t, 4, int32_t, do_qdmullh, SATMASK16T) 936 DO_2OP_SAT_L(vqdmulltw, 1, 4, int32_t, 8, int64_t, do_qdmullw, SATMASK32) 937 938 static inline uint32_t do_vbrsrb(uint32_t n, uint32_t m) 939 { 940 m &= 0xff; 941 if (m == 0) { 942 return 0; 943 } 944 n = revbit8(n); 945 if (m < 8) { 946 n >>= 8 - m; 947 } 948 return n; 949 } 950 951 static inline uint32_t do_vbrsrh(uint32_t n, uint32_t m) 952 { 953 m &= 0xff; 954 if (m == 0) { 955 return 0; 956 } 957 n = revbit16(n); 958 if (m < 16) { 959 n >>= 16 - m; 960 } 961 return n; 962 } 963 964 static inline uint32_t do_vbrsrw(uint32_t n, uint32_t m) 965 { 966 m &= 0xff; 967 if (m == 0) { 968 return 0; 969 } 970 n = revbit32(n); 971 if (m < 32) { 972 n >>= 32 - m; 973 } 974 return n; 975 } 976 977 DO_2OP_SCALAR(vbrsrb, 1, uint8_t, do_vbrsrb) 978 DO_2OP_SCALAR(vbrsrh, 2, uint16_t, do_vbrsrh) 979 DO_2OP_SCALAR(vbrsrw, 4, uint32_t, do_vbrsrw) 980 981 /* 982 * Multiply add long dual accumulate ops. 983 */ 984 #define DO_LDAV(OP, ESIZE, TYPE, XCHG, EVENACC, ODDACC) \ 985 uint64_t HELPER(glue(mve_, OP))(CPUARMState *env, void *vn, \ 986 void *vm, uint64_t a) \ 987 { \ 988 uint16_t mask = mve_element_mask(env); \ 989 unsigned e; \ 990 TYPE *n = vn, *m = vm; \ 991 for (e = 0; e < 16 / ESIZE; e++, mask >>= ESIZE) { \ 992 if (mask & 1) { \ 993 if (e & 1) { \ 994 a ODDACC \ 995 (int64_t)n[H##ESIZE(e - 1 * XCHG)] * m[H##ESIZE(e)]; \ 996 } else { \ 997 a EVENACC \ 998 (int64_t)n[H##ESIZE(e + 1 * XCHG)] * m[H##ESIZE(e)]; \ 999 } \ 1000 } \ 1001 } \ 1002 mve_advance_vpt(env); \ 1003 return a; \ 1004 } 1005 1006 DO_LDAV(vmlaldavsh, 2, int16_t, false, +=, +=) 1007 DO_LDAV(vmlaldavxsh, 2, int16_t, true, +=, +=) 1008 DO_LDAV(vmlaldavsw, 4, int32_t, false, +=, +=) 1009 DO_LDAV(vmlaldavxsw, 4, int32_t, true, +=, +=) 1010 1011 DO_LDAV(vmlaldavuh, 2, uint16_t, false, +=, +=) 1012 DO_LDAV(vmlaldavuw, 4, uint32_t, false, +=, +=) 1013 1014 DO_LDAV(vmlsldavsh, 2, int16_t, false, +=, -=) 1015 DO_LDAV(vmlsldavxsh, 2, int16_t, true, +=, -=) 1016 DO_LDAV(vmlsldavsw, 4, int32_t, false, +=, -=) 1017 DO_LDAV(vmlsldavxsw, 4, int32_t, true, +=, -=) 1018 1019 /* 1020 * Rounding multiply add long dual accumulate high: we must keep 1021 * a 72-bit internal accumulator value and return the top 64 bits. 1022 */ 1023 #define DO_LDAVH(OP, ESIZE, TYPE, XCHG, EVENACC, ODDACC, TO128) \ 1024 uint64_t HELPER(glue(mve_, OP))(CPUARMState *env, void *vn, \ 1025 void *vm, uint64_t a) \ 1026 { \ 1027 uint16_t mask = mve_element_mask(env); \ 1028 unsigned e; \ 1029 TYPE *n = vn, *m = vm; \ 1030 Int128 acc = int128_lshift(TO128(a), 8); \ 1031 for (e = 0; e < 16 / ESIZE; e++, mask >>= ESIZE) { \ 1032 if (mask & 1) { \ 1033 if (e & 1) { \ 1034 acc = ODDACC(acc, TO128(n[H##ESIZE(e - 1 * XCHG)] * \ 1035 m[H##ESIZE(e)])); \ 1036 } else { \ 1037 acc = EVENACC(acc, TO128(n[H##ESIZE(e + 1 * XCHG)] * \ 1038 m[H##ESIZE(e)])); \ 1039 } \ 1040 acc = int128_add(acc, int128_make64(1 << 7)); \ 1041 } \ 1042 } \ 1043 mve_advance_vpt(env); \ 1044 return int128_getlo(int128_rshift(acc, 8)); \ 1045 } 1046 1047 DO_LDAVH(vrmlaldavhsw, 4, int32_t, false, int128_add, int128_add, int128_makes64) 1048 DO_LDAVH(vrmlaldavhxsw, 4, int32_t, true, int128_add, int128_add, int128_makes64) 1049 1050 DO_LDAVH(vrmlaldavhuw, 4, uint32_t, false, int128_add, int128_add, int128_make64) 1051 1052 DO_LDAVH(vrmlsldavhsw, 4, int32_t, false, int128_add, int128_sub, int128_makes64) 1053 DO_LDAVH(vrmlsldavhxsw, 4, int32_t, true, int128_add, int128_sub, int128_makes64) 1054