^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 1) // SPDX-License-Identifier: GPL-2.0-only
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 2) /* IEEE754 floating point arithmetic
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 3) * double precision: common utilities
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 4) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 5) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 6) * MIPS floating point support
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 7) * Copyright (C) 1994-2000 Algorithmics Ltd.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 8) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 9)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 10) #include "ieee754dp.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 11)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 12) int ieee754dp_cmp(union ieee754dp x, union ieee754dp y, int cmp, int sig)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 13) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 14) s64 vx;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 15) s64 vy;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 16)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 17) COMPXDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 18) COMPYDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 19)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 20) EXPLODEXDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 21) EXPLODEYDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 22) FLUSHXDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 23) FLUSHYDP;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 24) ieee754_clearcx(); /* Even clear inexact flag here */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 25)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 26) if (ieee754_class_nan(xc) || ieee754_class_nan(yc)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 27) if (sig ||
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 28) xc == IEEE754_CLASS_SNAN || yc == IEEE754_CLASS_SNAN)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 29) ieee754_setcx(IEEE754_INVALID_OPERATION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30) return (cmp & IEEE754_CUN) != 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 31) } else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 32) vx = x.bits;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 33) vy = y.bits;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 34)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 35) if (vx < 0)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 36) vx = -vx ^ DP_SIGN_BIT;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37) if (vy < 0)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) vy = -vy ^ DP_SIGN_BIT;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) if (vx < vy)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) return (cmp & IEEE754_CLT) != 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) else if (vx == vy)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) return (cmp & IEEE754_CEQ) != 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44) else
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) return (cmp & IEEE754_CGT) != 0;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) }