Back to index

glibc  2.9
ots_cmp.c
Go to the documentation of this file.
00001 /* Software floating-point emulation: comparison.
00002    Copyright (C) 1997,1999,2004 Free Software Foundation, Inc.
00003    This file is part of the GNU C Library.
00004    Contributed by Richard Henderson (rth@cygnus.com) and
00005                 Jakub Jelinek (jj@ultra.linux.cz).
00006 
00007    The GNU C Library is free software; you can redistribute it and/or
00008    modify it under the terms of the GNU Lesser General Public
00009    License as published by the Free Software Foundation; either
00010    version 2.1 of the License, or (at your option) any later version.
00011 
00012    The GNU C Library is distributed in the hope that it will be useful,
00013    but WITHOUT ANY WARRANTY; without even the implied warranty of
00014    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00015    Lesser General Public License for more details.
00016 
00017    You should have received a copy of the GNU Lesser General Public
00018    License along with the GNU C Library; if not, write to the Free
00019    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
00020    02111-1307 USA.  */
00021 
00022 #include "local-soft-fp.h"
00023 
00024 static long
00025 internal_equality (long al, long ah, long bl, long bh, long neq)
00026 {
00027   FP_DECL_EX;
00028   FP_DECL_Q(A); FP_DECL_Q(B);
00029   long r;
00030 
00031   FP_UNPACK_RAW_Q(A, a);
00032   FP_UNPACK_RAW_Q(B, b);
00033 
00034   if ((A_e == _FP_EXPMAX_Q && !_FP_FRAC_ZEROP_2(A))
00035        || (B_e == _FP_EXPMAX_Q && !_FP_FRAC_ZEROP_2(B)))
00036     {
00037       /* EQ and NE signal invalid operation only if either operand is SNaN.  */
00038       if (FP_ISSIGNAN_Q(A) || FP_ISSIGNAN_Q(B))
00039        {
00040          FP_SET_EXCEPTION(FP_EX_INVALID);
00041          FP_HANDLE_EXCEPTIONS;
00042        }
00043       return -1;
00044     }
00045 
00046   r = (A_e == B_e
00047        && _FP_FRAC_EQ_2 (A, B)
00048        && (A_s == B_s || (!A_e && _FP_FRAC_ZEROP_2(A))));
00049   r ^= neq;
00050 
00051   return r;
00052 }
00053 
00054 long
00055 _OtsEqlX (long al, long ah, long bl, long bh)
00056 {
00057   return internal_equality (al, ah, bl, bh, 0);
00058 }
00059 
00060 long
00061 _OtsNeqX (long al, long ah, long bl, long bh)
00062 {
00063   return internal_equality (al, ah, bl, bh, 1);
00064 }