1 /* IEEE754 floating point arithmetic
5 * MIPS floating point support
6 * Copyright (C) 1994-2000 Algorithmics Ltd.
7 * Copyright (C) 2017 Imagination Technologies, Ltd.
8 * Author: Aleksandar Markovic <aleksandar.markovic@imgtec.com>
10 * This program is free software; you can distribute it and/or modify it
11 * under the terms of the GNU General Public License (Version 2) as
12 * published by the Free Software Foundation.
14 * This program is distributed in the hope it will be useful, but WITHOUT
15 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
16 * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19 * You should have received a copy of the GNU General Public License along
23 #include "ieee754sp.h"
25 union ieee754sp
ieee754sp_rint(union ieee754sp x
)
33 COMPXDP
; /* <-- DP needed for 64-bit mantissa tmp */
40 if (xc
== IEEE754_CLASS_SNAN
)
41 return ieee754sp_nanxcpt(x
);
43 if ((xc
== IEEE754_CLASS_QNAN
) ||
44 (xc
== IEEE754_CLASS_INF
) ||
45 (xc
== IEEE754_CLASS_ZERO
))
54 sticky
= residue
!= 0;
57 residue
= xm
<< (xe
+ 1);
58 residue
<<= 31 - SP_FBITS
;
59 round
= (residue
>> 31) != 0;
60 sticky
= (residue
<< 1) != 0;
64 odd
= (xm
& 0x1) != 0x0;
66 switch (ieee754_csr
.rm
) {
67 case FPU_CSR_RN
: /* toward nearest */
68 if (round
&& (sticky
|| odd
))
71 case FPU_CSR_RZ
: /* toward zero */
73 case FPU_CSR_RU
: /* toward +infinity */
74 if ((round
|| sticky
) && !xs
)
77 case FPU_CSR_RD
: /* toward -infinity */
78 if ((round
|| sticky
) && xs
)
84 ieee754_setcx(IEEE754_INEXACT
);
86 ret
= ieee754sp_flong(xm
);