1234567891011121314151617181920212223242526272829303132333435363738394041424344454647 |
- // 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.
- */
- #include "ieee754dp.h"
- int ieee754dp_cmp(union ieee754dp x, union ieee754dp y, int cmp, int sig)
- {
- s64 vx;
- s64 vy;
- COMPXDP;
- COMPYDP;
- EXPLODEXDP;
- EXPLODEYDP;
- FLUSHXDP;
- FLUSHYDP;
- ieee754_clearcx(); /* Even clear inexact flag here */
- if (ieee754_class_nan(xc) || ieee754_class_nan(yc)) {
- if (sig ||
- xc == IEEE754_CLASS_SNAN || yc == IEEE754_CLASS_SNAN)
- ieee754_setcx(IEEE754_INVALID_OPERATION);
- return (cmp & IEEE754_CUN) != 0;
- } else {
- vx = x.bits;
- vy = y.bits;
- if (vx < 0)
- vx = -vx ^ DP_SIGN_BIT;
- if (vy < 0)
- vy = -vy ^ DP_SIGN_BIT;
- if (vx < vy)
- return (cmp & IEEE754_CLT) != 0;
- else if (vx == vy)
- return (cmp & IEEE754_CEQ) != 0;
- else
- return (cmp & IEEE754_CGT) != 0;
- }
- }
|