@@ -86,41 +86,14 @@ void compute_outputs(const uint32_t first_row, const uint32_t num_rows) {
8686
8787 const FLOAT_TYPE d = FLOAT_TYPE(data_a[ib0 + i].d);
8888
89- uint32_t ql0_u32 = uint32_t(blkcache[ix].blk[ql_offset / 2]) | (uint32_t(blkcache[ix].blk[ql_offset / 2 + 1]) << 16);
90- uint32_t ql32_u32 = uint32_t(blkcache[ix].blk[ql_offset / 2 + 16]) | (uint32_t(blkcache[ix].blk[ql_offset / 2 + 17]) << 16);
91-
92- uint32_t ql0_u32_lo4 = ql0_u32 & 0x0F0F0F0F;
93- uint32_t ql0_u32_hi4 = (ql0_u32 >> 4) & 0x0F0F0F0F;
94- uint32_t ql32_u32_lo4 = ql32_u32 & 0x0F0F0F0F;
95- uint32_t ql32_u32_hi4 = (ql32_u32 >> 4) & 0x0F0F0F0F;
96-
97- uint32_t qh_u32 = uint32_t(blkcache[ix].blk[64 + qh_offset / 2]) | (uint32_t(blkcache[ix].blk[64 + qh_offset / 2 + 1]) << 16);
98- uint32_t qh0_u32 = (qh_u32 & 0x03030303) << 4;
99- uint32_t qh2_u32 = (qh_u32 & 0x0C0C0C0C) << 2;
100- uint32_t qh4_u32 = (qh_u32 & 0x30303030);
101- uint32_t qh6_u32 = (qh_u32 & 0xC0C0C0C0) >> 2;
102-
103- uint32_t q0_u32 = ql0_u32_lo4 | qh0_u32;
104- uint32_t q1_u32 = ql32_u32_lo4 | qh2_u32;
105- uint32_t q2_u32 = ql0_u32_hi4 | qh4_u32;
106- uint32_t q3_u32 = ql32_u32_hi4 | qh6_u32;
107-
108- uvec4 q0 = uvec4(unpack8(q0_u32));
109- uvec4 q1 = uvec4(unpack8(q1_u32));
110- uvec4 q2 = uvec4(unpack8(q2_u32));
111- uvec4 q3 = uvec4(unpack8(q3_u32));
112-
113- FLOAT_TYPE sum[4] = {0, 0, 0, 0};
114- [[unroll]] for (uint l = 0; l < 4; ++l) {
115- sum[0] = fma(FLOAT_TYPE(bycache[ix][128*v_im + l0 + l]), FLOAT_TYPE(int8_t(q0[l]) - 32), sum[0]);
116- sum[1] = fma(FLOAT_TYPE(bycache[ix][128*v_im + l0 + 32 + l]), FLOAT_TYPE(int8_t(q1[l]) - 32), sum[1]);
117- sum[2] = fma(FLOAT_TYPE(bycache[ix][128*v_im + l0 + 64 + l]), FLOAT_TYPE(int8_t(q2[l]) - 32), sum[2]);
118- sum[3] = fma(FLOAT_TYPE(bycache[ix][128*v_im + l0 + 96 + l]), FLOAT_TYPE(int8_t(q3[l]) - 32), sum[3]);
119- }
120-
121- [[unroll]] for (uint l = 0; l < 4; ++l)
122- sum[l] *= sccache[ix][s_offset + l*2];
123- temp[n] += (sum[0] + sum[1] + sum[2] + sum[3]) * d;
89+ FLOAT_TYPE dq[16];
90+ FLOAT_TYPE sum = 0;
91+ [[unroll]] for (uint l = 0; l < 16; ++l) {
92+ dq[l] = bitfieldExtract(blkcache[ix].blk[l/2], 1, 4) | (bitfieldExtract(blkcache[ix].blk[64 + l], 2, 2) << 4);
93+ sum = fma(FLOAT_TYPE(bycache[ix][16*itid + l]), dq[l], sum);
94+ }
95+
96+ temp[n] += sum * sccache[ix][itid] * d;
12497 }
12598 }
12699
0 commit comments