blob: a0c7f99a6cb980495b5e3645306212c9cab903cb (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
// vsmul: Signed saturating and rounding fractional multiply
VRM xrm = P.VU.get_vround_mode();
uint64_t int_max = (uint64_t(1) << (P.VU.vsew - 1)) - 1;
uint64_t int_min = - (1 << (P.VU.vsew - 1));
uint64_t sign_mask = uint64_t(1) << (P.VU.vsew - 1);
VI_VV_ULOOP
({
uint64_t vs1_sign;
uint64_t vs2_sign;
uint64_t result_sign;
vs1_sign = vs1 & sign_mask;
vs2_sign = vs2 & sign_mask;
bool overflow = vs1 == vs2 && vs1 == int_min;
uint128_t result = (uint128_t)vs1 * (uint128_t)vs2;
result &= ((uint128_t)1llu << ((sew * 2) - 2)) - 1;
result_sign = (vs1_sign ^ vs2_sign) & sign_mask;
// rounding
INT_ROUNDING(result, xrm, sew - 1);
// unsigned shifting
result = result >> (sew - 1);
// saturation
if (overflow) {
result = int_max;
P.VU.vxsat = 1;
} else {
result |= result_sign;
}
vd = result;
})
|