1 /* Software floating-point emulation: comparison.
2 Copyright (C) 1997-2022 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Lesser General Public
7 License as published by the Free Software Foundation; either
8 version 2.1 of the License, or (at your option) any later version.
9
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Lesser General Public License for more details.
14
15 You should have received a copy of the GNU Lesser General Public
16 License along with the GNU C Library. If not, see
17 <https://www.gnu.org/licenses/>. */
18
19 #include "local-soft-fp.h"
20
21 static long
internal_equality(long al,long ah,long bl,long bh,long neq)22 internal_equality (long al, long ah, long bl, long bh, long neq)
23 {
24 FP_DECL_EX;
25 FP_DECL_Q(A); FP_DECL_Q(B);
26 long r;
27
28 AXP_UNPACK_RAW_Q(A, a);
29 AXP_UNPACK_RAW_Q(B, b);
30
31 if ((A_e == _FP_EXPMAX_Q && !_FP_FRAC_ZEROP_2(A))
32 || (B_e == _FP_EXPMAX_Q && !_FP_FRAC_ZEROP_2(B)))
33 {
34 /* EQ and NE signal invalid operation only if either operand is SNaN. */
35 if (FP_ISSIGNAN_Q(A) || FP_ISSIGNAN_Q(B))
36 {
37 FP_SET_EXCEPTION(FP_EX_INVALID);
38 FP_HANDLE_EXCEPTIONS;
39 }
40 return -1;
41 }
42
43 r = (A_e == B_e
44 && _FP_FRAC_EQ_2 (A, B)
45 && (A_s == B_s || (!A_e && _FP_FRAC_ZEROP_2(A))));
46 r ^= neq;
47
48 return r;
49 }
50
51 long
_OtsEqlX(long al,long ah,long bl,long bh)52 _OtsEqlX (long al, long ah, long bl, long bh)
53 {
54 return internal_equality (al, ah, bl, bh, 0);
55 }
56
57 long
_OtsNeqX(long al,long ah,long bl,long bh)58 _OtsNeqX (long al, long ah, long bl, long bh)
59 {
60 return internal_equality (al, ah, bl, bh, 1);
61 }
62