forked from Minki/linux
f6cc69f1f3
Based on 2 normalized pattern(s): it and or modify it under the terms of the gnu general public license version 2 as published by the free software foundation this program is distributed in the hope it will be useful but without any warranty without even the implied warranty of merchantability or fitness for a particular purpose see the gnu general public license for more details you should have received a copy of the gnu general public license along with this program this program is free software you can redistribute it and or modify it under the terms and conditions of the gnu general public license version 2 as published by the free software foundation this program is distributed in the hope it will be useful but without any warranty without even the implied warranty of merchantability or fitness for a particular purpose see the gnu general public license for more details you should have received a copy of the gnu general public license along with this program extracted by the scancode license scanner the SPDX license identifier GPL-2.0-only has been chosen to replace the boilerplate/reference in 11 file(s). Signed-off-by: Thomas Gleixner <tglx@linutronix.de> Reviewed-by: Alexios Zavras <alexios.zavras@intel.com> Reviewed-by: Allison Randal <allison@lohutok.net> Cc: linux-spdx@vger.kernel.org Link: https://lkml.kernel.org/r/20190530000434.249870634@linutronix.de Signed-off-by: Greg Kroah-Hartman <gregkh@linuxfoundation.org>
79 lines
1.4 KiB
C
79 lines
1.4 KiB
C
// SPDX-License-Identifier: GPL-2.0-only
|
|
/* IEEE754 floating point arithmetic
|
|
* double precision: common utilities
|
|
*/
|
|
/*
|
|
* MIPS floating point support
|
|
* Copyright (C) 1994-2000 Algorithmics Ltd.
|
|
* Copyright (C) 2017 Imagination Technologies, Ltd.
|
|
* Author: Aleksandar Markovic <aleksandar.markovic@imgtec.com>
|
|
*/
|
|
|
|
#include "ieee754dp.h"
|
|
|
|
union ieee754dp ieee754dp_rint(union ieee754dp x)
|
|
{
|
|
union ieee754dp ret;
|
|
u64 residue;
|
|
int sticky;
|
|
int round;
|
|
int odd;
|
|
|
|
COMPXDP;
|
|
|
|
ieee754_clearcx();
|
|
|
|
EXPLODEXDP;
|
|
FLUSHXDP;
|
|
|
|
if (xc == IEEE754_CLASS_SNAN)
|
|
return ieee754dp_nanxcpt(x);
|
|
|
|
if ((xc == IEEE754_CLASS_QNAN) ||
|
|
(xc == IEEE754_CLASS_INF) ||
|
|
(xc == IEEE754_CLASS_ZERO))
|
|
return x;
|
|
|
|
if (xe >= DP_FBITS)
|
|
return x;
|
|
|
|
if (xe < -1) {
|
|
residue = xm;
|
|
round = 0;
|
|
sticky = residue != 0;
|
|
xm = 0;
|
|
} else {
|
|
residue = xm << (64 - DP_FBITS + xe);
|
|
round = (residue >> 63) != 0;
|
|
sticky = (residue << 1) != 0;
|
|
xm >>= DP_FBITS - xe;
|
|
}
|
|
|
|
odd = (xm & 0x1) != 0x0;
|
|
|
|
switch (ieee754_csr.rm) {
|
|
case FPU_CSR_RN: /* toward nearest */
|
|
if (round && (sticky || odd))
|
|
xm++;
|
|
break;
|
|
case FPU_CSR_RZ: /* toward zero */
|
|
break;
|
|
case FPU_CSR_RU: /* toward +infinity */
|
|
if ((round || sticky) && !xs)
|
|
xm++;
|
|
break;
|
|
case FPU_CSR_RD: /* toward -infinity */
|
|
if ((round || sticky) && xs)
|
|
xm++;
|
|
break;
|
|
}
|
|
|
|
if (round || sticky)
|
|
ieee754_setcx(IEEE754_INEXACT);
|
|
|
|
ret = ieee754dp_flong(xm);
|
|
DPSIGN(ret) = xs;
|
|
|
|
return ret;
|
|
}
|