/* * VIS op helpers * * Copyright (c) 2003-2005 Fabrice Bellard * * This library is free software; you can redistribute it and/or * modify it under the terms of the GNU Lesser General Public * License as published by the Free Software Foundation; either * version 2 of the License, or (at your option) any later version. * * This library is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * Lesser General Public License for more details. * * You should have received a copy of the GNU Lesser General Public * License along with this library; if not, see <http://www.gnu.org/licenses/>. */ #include "cpu.h" #include "helper.h" /* This function uses non-native bit order */ #define GET_FIELD(X, FROM, TO) \ ((X) >> (63 - (TO)) & ((1ULL << ((TO) - (FROM) + 1)) - 1)) /* This function uses the order in the manuals, i.e. bit 0 is 2^0 */ #define GET_FIELD_SP(X, FROM, TO) \ GET_FIELD(X, 63 - (TO), 63 - (FROM)) target_ulong helper_array8(target_ulong pixel_addr, target_ulong cubesize) { return (GET_FIELD_SP(pixel_addr, 60, 63) << (17 + 2 * cubesize)) | (GET_FIELD_SP(pixel_addr, 39, 39 + cubesize - 1) << (17 + cubesize)) | (GET_FIELD_SP(pixel_addr, 17 + cubesize - 1, 17) << 17) | (GET_FIELD_SP(pixel_addr, 56, 59) << 13) | (GET_FIELD_SP(pixel_addr, 35, 38) << 9) | (GET_FIELD_SP(pixel_addr, 13, 16) << 5) | (((pixel_addr >> 55) & 1) << 4) | (GET_FIELD_SP(pixel_addr, 33, 34) << 2) | GET_FIELD_SP(pixel_addr, 11, 12); } #ifdef HOST_WORDS_BIGENDIAN #define VIS_B64(n) b[7 - (n)] #define VIS_W64(n) w[3 - (n)] #define VIS_SW64(n) sw[3 - (n)] #define VIS_L64(n) l[1 - (n)] #define VIS_B32(n) b[3 - (n)] #define VIS_W32(n) w[1 - (n)] #else #define VIS_B64(n) b[n] #define VIS_W64(n) w[n] #define VIS_SW64(n) sw[n] #define VIS_L64(n) l[n] #define VIS_B32(n) b[n] #define VIS_W32(n) w[n] #endif typedef union { uint8_t b[8]; uint16_t w[4]; int16_t sw[4]; uint32_t l[2]; uint64_t ll; float64 d; } VIS64; typedef union { uint8_t b[4]; uint16_t w[2]; uint32_t l; float32 f; } VIS32; uint64_t helper_fpmerge(uint64_t src1, uint64_t src2) { VIS64 s, d; s.ll = src1; d.ll = src2; /* Reverse calculation order to handle overlap */ d.VIS_B64(7) = s.VIS_B64(3); d.VIS_B64(6) = d.VIS_B64(3); d.VIS_B64(5) = s.VIS_B64(2); d.VIS_B64(4) = d.VIS_B64(2); d.VIS_B64(3) = s.VIS_B64(1); d.VIS_B64(2) = d.VIS_B64(1); d.VIS_B64(1) = s.VIS_B64(0); /* d.VIS_B64(0) = d.VIS_B64(0); */ return d.ll; } uint64_t helper_fmul8x16(uint64_t src1, uint64_t src2) { VIS64 s, d; uint32_t tmp; s.ll = src1; d.ll = src2; #define PMUL(r) \ tmp = (int32_t)d.VIS_SW64(r) * (int32_t)s.VIS_B64(r); \ if ((tmp & 0xff) > 0x7f) { \ tmp += 0x100; \ } \ d.VIS_W64(r) = tmp >> 8; PMUL(0); PMUL(1); PMUL(2); PMUL(3); #undef PMUL return d.ll; } uint64_t helper_fmul8x16al(uint64_t src1, uint64_t src2) { VIS64 s, d; uint32_t tmp; s.ll = src1; d.ll = src2; #define PMUL(r) \ tmp = (int32_t)d.VIS_SW64(1) * (int32_t)s.VIS_B64(r); \ if ((tmp & 0xff) > 0x7f) { \ tmp += 0x100; \ } \ d.VIS_W64(r) = tmp >> 8; PMUL(0); PMUL(1); PMUL(2); PMUL(3); #undef PMUL return d.ll; } uint64_t helper_fmul8x16au(uint64_t src1, uint64_t src2) { VIS64 s, d; uint32_t tmp; s.ll = src1; d.ll = src2; #define PMUL(r) \ tmp = (int32_t)d.VIS_SW64(0) * (int32_t)s.VIS_B64(r); \ if ((tmp & 0xff) > 0x7f) { \ tmp += 0x100; \ } \ d.VIS_W64(r) = tmp >> 8; PMUL(0); PMUL(1); PMUL(2); PMUL(3); #undef PMUL return d.ll; } uint64_t helper_fmul8sux16(uint64_t src1, uint64_t src2) { VIS64 s, d; uint32_t tmp; s.ll = src1; d.ll = src2; #define PMUL(r) \ tmp = (int32_t)d.VIS_SW64(r) * ((int32_t)s.VIS_SW64(r) >> 8); \ if ((tmp & 0xff) > 0x7f) { \ tmp += 0x100; \ } \ d.VIS_W64(r) = tmp >> 8; PMUL(0); PMUL(1); PMUL(2); PMUL(3); #undef PMUL return d.ll; } uint64_t helper_fmul8ulx16(uint64_t src1, uint64_t src2) { VIS64 s, d; uint32_t tmp; s.ll = src1; d.ll = src2; #define PMUL(r) \ tmp = (int32_t)d.VIS_SW64(r) * ((uint32_t)s.VIS_B64(r * 2)); \ if ((tmp & 0xff) > 0x7f) { \ tmp += 0x100; \ } \ d.VIS_W64(r) = tmp >> 8; PMUL(0); PMUL(1); PMUL(2); PMUL(3); #undef PMUL return d.ll; } uint64_t helper_fmuld8sux16(uint64_t src1, uint64_t src2) { VIS64 s, d; uint32_t tmp; s.ll = src1; d.ll = src2; #define PMUL(r) \ tmp = (int32_t)d.VIS_SW64(r) * ((int32_t)s.VIS_SW64(r) >> 8); \ if ((tmp & 0xff) > 0x7f) { \ tmp += 0x100; \ } \ d.VIS_L64(r) = tmp; /* Reverse calculation order to handle overlap */ PMUL(1); PMUL(0); #undef PMUL return d.ll; } uint64_t helper_fmuld8ulx16(uint64_t src1, uint64_t src2) { VIS64 s, d; uint32_t tmp; s.ll = src1; d.ll = src2; #define PMUL(r) \ tmp = (int32_t)d.VIS_SW64(r) * ((uint32_t)s.VIS_B64(r * 2)); \ if ((tmp & 0xff) > 0x7f) { \ tmp += 0x100; \ } \ d.VIS_L64(r) = tmp; /* Reverse calculation order to handle overlap */ PMUL(1); PMUL(0); #undef PMUL return d.ll; } uint64_t helper_fexpand(uint64_t src1, uint64_t src2) { VIS32 s; VIS64 d; s.l = (uint32_t)src1; d.ll = src2; d.VIS_W64(0) = s.VIS_B32(0) << 4; d.VIS_W64(1) = s.VIS_B32(1) << 4; d.VIS_W64(2) = s.VIS_B32(2) << 4; d.VIS_W64(3) = s.VIS_B32(3) << 4; return d.ll; } #define VIS_HELPER(name, F) \ uint64_t name##16(uint64_t src1, uint64_t src2) \ { \ VIS64 s, d; \ \ s.ll = src1; \ d.ll = src2; \ \ d.VIS_W64(0) = F(d.VIS_W64(0), s.VIS_W64(0)); \ d.VIS_W64(1) = F(d.VIS_W64(1), s.VIS_W64(1)); \ d.VIS_W64(2) = F(d.VIS_W64(2), s.VIS_W64(2)); \ d.VIS_W64(3) = F(d.VIS_W64(3), s.VIS_W64(3)); \ \ return d.ll; \ } \ \ uint32_t name##16s(uint32_t src1, uint32_t src2) \ { \ VIS32 s, d; \ \ s.l = src1; \ d.l = src2; \ \ d.VIS_W32(0) = F(d.VIS_W32(0), s.VIS_W32(0)); \ d.VIS_W32(1) = F(d.VIS_W32(1), s.VIS_W32(1)); \ \ return d.l; \ } \ \ uint64_t name##32(uint64_t src1, uint64_t src2) \ { \ VIS64 s, d; \ \ s.ll = src1; \ d.ll = src2; \ \ d.VIS_L64(0) = F(d.VIS_L64(0), s.VIS_L64(0)); \ d.VIS_L64(1) = F(d.VIS_L64(1), s.VIS_L64(1)); \ \ return d.ll; \ } \ \ uint32_t name##32s(uint32_t src1, uint32_t src2) \ { \ VIS32 s, d; \ \ s.l = src1; \ d.l = src2; \ \ d.l = F(d.l, s.l); \ \ return d.l; \ } #define FADD(a, b) ((a) + (b)) #define FSUB(a, b) ((a) - (b)) VIS_HELPER(helper_fpadd, FADD) VIS_HELPER(helper_fpsub, FSUB) #define VIS_CMPHELPER(name, F) \ uint64_t name##16(uint64_t src1, uint64_t src2) \ { \ VIS64 s, d; \ \ s.ll = src1; \ d.ll = src2; \ \ d.VIS_W64(0) = F(s.VIS_W64(0), d.VIS_W64(0)) ? 1 : 0; \ d.VIS_W64(0) |= F(s.VIS_W64(1), d.VIS_W64(1)) ? 2 : 0; \ d.VIS_W64(0) |= F(s.VIS_W64(2), d.VIS_W64(2)) ? 4 : 0; \ d.VIS_W64(0) |= F(s.VIS_W64(3), d.VIS_W64(3)) ? 8 : 0; \ d.VIS_W64(1) = d.VIS_W64(2) = d.VIS_W64(3) = 0; \ \ return d.ll; \ } \ \ uint64_t name##32(uint64_t src1, uint64_t src2) \ { \ VIS64 s, d; \ \ s.ll = src1; \ d.ll = src2; \ \ d.VIS_L64(0) = F(s.VIS_L64(0), d.VIS_L64(0)) ? 1 : 0; \ d.VIS_L64(0) |= F(s.VIS_L64(1), d.VIS_L64(1)) ? 2 : 0; \ d.VIS_L64(1) = 0; \ \ return d.ll; \ } #define FCMPGT(a, b) ((a) > (b)) #define FCMPEQ(a, b) ((a) == (b)) #define FCMPLE(a, b) ((a) <= (b)) #define FCMPNE(a, b) ((a) != (b)) VIS_CMPHELPER(helper_fcmpgt, FCMPGT) VIS_CMPHELPER(helper_fcmpeq, FCMPEQ) VIS_CMPHELPER(helper_fcmple, FCMPLE) VIS_CMPHELPER(helper_fcmpne, FCMPNE) uint64_t helper_pdist(uint64_t sum, uint64_t src1, uint64_t src2) { int i; for (i = 0; i < 8; i++) { int s1, s2; s1 = (src1 >> (56 - (i * 8))) & 0xff; s2 = (src2 >> (56 - (i * 8))) & 0xff; /* Absolute value of difference. */ s1 -= s2; if (s1 < 0) { s1 = -s1; } sum += s1; } return sum; } uint32_t helper_fpack16(uint64_t gsr, uint64_t rs2) { int scale = (gsr >> 3) & 0xf; uint32_t ret = 0; int byte; for (byte = 0; byte < 4; byte++) { uint32_t val; int16_t src = rs2 >> (byte * 16); int32_t scaled = src << scale; int32_t from_fixed = scaled >> 7; val = (from_fixed < 0 ? 0 : from_fixed > 255 ? 255 : from_fixed); ret |= val << (8 * byte); } return ret; } uint64_t helper_fpack32(uint64_t gsr, uint64_t rs1, uint64_t rs2) { int scale = (gsr >> 3) & 0x1f; uint64_t ret = 0; int word; ret = (rs1 << 8) & ~(0x000000ff000000ffULL); for (word = 0; word < 2; word++) { uint64_t val; int32_t src = rs2 >> (word * 32); int64_t scaled = (int64_t)src << scale; int64_t from_fixed = scaled >> 23; val = (from_fixed < 0 ? 0 : (from_fixed > 255) ? 255 : from_fixed); ret |= val << (32 * word); } return ret; } uint32_t helper_fpackfix(uint64_t gsr, uint64_t rs2) { int scale = (gsr >> 3) & 0x1f; uint32_t ret = 0; int word; for (word = 0; word < 2; word++) { uint32_t val; int32_t src = rs2 >> (word * 32); int64_t scaled = src << scale; int64_t from_fixed = scaled >> 16; val = (from_fixed < -32768 ? -32768 : from_fixed > 32767 ? 32767 : from_fixed); ret |= (val & 0xffff) << (word * 16); } return ret; } uint64_t helper_bshuffle(uint64_t gsr, uint64_t src1, uint64_t src2) { union { uint64_t ll[2]; uint8_t b[16]; } s; VIS64 r; uint32_t i, mask, host; /* Set up S such that we can index across all of the bytes. */ #ifdef HOST_WORDS_BIGENDIAN s.ll[0] = src1; s.ll[1] = src2; host = 0; #else s.ll[1] = src1; s.ll[0] = src2; host = 15; #endif mask = gsr >> 32; for (i = 0; i < 8; ++i) { unsigned e = (mask >> (28 - i*4)) & 0xf; r.VIS_B64(i) = s.b[e ^ host]; } return r.ll; }