1*4882a593Smuzhiyun // SPDX-License-Identifier: GPL-2.0-only 2*4882a593Smuzhiyun /* IEEE754 floating point arithmetic 3*4882a593Smuzhiyun * single precision 4*4882a593Smuzhiyun */ 5*4882a593Smuzhiyun /* 6*4882a593Smuzhiyun * MIPS floating point support 7*4882a593Smuzhiyun * Copyright (C) 1994-2000 Algorithmics Ltd. 8*4882a593Smuzhiyun */ 9*4882a593Smuzhiyun 10*4882a593Smuzhiyun #include "ieee754sp.h" 11*4882a593Smuzhiyun ieee754sp_cmp(union ieee754sp x,union ieee754sp y,int cmp,int sig)12*4882a593Smuzhiyunint ieee754sp_cmp(union ieee754sp x, union ieee754sp y, int cmp, int sig) 13*4882a593Smuzhiyun { 14*4882a593Smuzhiyun int vx; 15*4882a593Smuzhiyun int vy; 16*4882a593Smuzhiyun 17*4882a593Smuzhiyun COMPXSP; 18*4882a593Smuzhiyun COMPYSP; 19*4882a593Smuzhiyun 20*4882a593Smuzhiyun EXPLODEXSP; 21*4882a593Smuzhiyun EXPLODEYSP; 22*4882a593Smuzhiyun FLUSHXSP; 23*4882a593Smuzhiyun FLUSHYSP; 24*4882a593Smuzhiyun ieee754_clearcx(); /* Even clear inexact flag here */ 25*4882a593Smuzhiyun 26*4882a593Smuzhiyun if (ieee754_class_nan(xc) || ieee754_class_nan(yc)) { 27*4882a593Smuzhiyun if (sig || 28*4882a593Smuzhiyun xc == IEEE754_CLASS_SNAN || yc == IEEE754_CLASS_SNAN) 29*4882a593Smuzhiyun ieee754_setcx(IEEE754_INVALID_OPERATION); 30*4882a593Smuzhiyun return (cmp & IEEE754_CUN) != 0; 31*4882a593Smuzhiyun } else { 32*4882a593Smuzhiyun vx = x.bits; 33*4882a593Smuzhiyun vy = y.bits; 34*4882a593Smuzhiyun 35*4882a593Smuzhiyun if (vx < 0) 36*4882a593Smuzhiyun vx = -vx ^ SP_SIGN_BIT; 37*4882a593Smuzhiyun if (vy < 0) 38*4882a593Smuzhiyun vy = -vy ^ SP_SIGN_BIT; 39*4882a593Smuzhiyun 40*4882a593Smuzhiyun if (vx < vy) 41*4882a593Smuzhiyun return (cmp & IEEE754_CLT) != 0; 42*4882a593Smuzhiyun else if (vx == vy) 43*4882a593Smuzhiyun return (cmp & IEEE754_CEQ) != 0; 44*4882a593Smuzhiyun else 45*4882a593Smuzhiyun return (cmp & IEEE754_CGT) != 0; 46*4882a593Smuzhiyun } 47*4882a593Smuzhiyun } 48