1  /* Software floating-point emulation: comparison.
       2     Copyright (C) 1997-2023 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
      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
      52  _OtsEqlX (long al, long ah, long bl, long bh)
      53  {
      54    return internal_equality (al, ah, bl, bh, 0);
      55  }
      56  
      57  long
      58  _OtsNeqX (long al, long ah, long bl, long bh)
      59  {
      60    return internal_equality (al, ah, bl, bh, 1);
      61  }