hc
2025-02-14 bbb9540dc49f70f6b703d1c8d1b85fa5f602d86e
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
// SPDX-License-Identifier: GPL-2.0-only
/* IEEE754 floating point arithmetic
 * single precision
 */
/*
 * MIPS floating point support
 * Copyright (C) 1994-2000 Algorithmics Ltd.
 */
 
#include "ieee754sp.h"
 
int ieee754sp_cmp(union ieee754sp x, union ieee754sp y, int cmp, int sig)
{
   int vx;
   int vy;
 
   COMPXSP;
   COMPYSP;
 
   EXPLODEXSP;
   EXPLODEYSP;
   FLUSHXSP;
   FLUSHYSP;
   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 ^ SP_SIGN_BIT;
       if (vy < 0)
           vy = -vy ^ SP_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;
   }
}