1 1.1 mrg /* mpf_cmp -- Compare two floats. 2 1.1 mrg 3 1.1.1.3 mrg Copyright 1993, 1994, 1996, 2001, 2015 Free Software Foundation, Inc. 4 1.1 mrg 5 1.1 mrg This file is part of the GNU MP Library. 6 1.1 mrg 7 1.1 mrg The GNU MP Library is free software; you can redistribute it and/or modify 8 1.1.1.3 mrg it under the terms of either: 9 1.1.1.3 mrg 10 1.1.1.3 mrg * the GNU Lesser General Public License as published by the Free 11 1.1.1.3 mrg Software Foundation; either version 3 of the License, or (at your 12 1.1.1.3 mrg option) any later version. 13 1.1.1.3 mrg 14 1.1.1.3 mrg or 15 1.1.1.3 mrg 16 1.1.1.3 mrg * the GNU General Public License as published by the Free Software 17 1.1.1.3 mrg Foundation; either version 2 of the License, or (at your option) any 18 1.1.1.3 mrg later version. 19 1.1.1.3 mrg 20 1.1.1.3 mrg or both in parallel, as here. 21 1.1 mrg 22 1.1 mrg The GNU MP Library is distributed in the hope that it will be useful, but 23 1.1 mrg WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 1.1.1.3 mrg or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 25 1.1.1.3 mrg for more details. 26 1.1 mrg 27 1.1.1.3 mrg You should have received copies of the GNU General Public License and the 28 1.1.1.3 mrg GNU Lesser General Public License along with the GNU MP Library. If not, 29 1.1.1.3 mrg see https://www.gnu.org/licenses/. */ 30 1.1 mrg 31 1.1 mrg #include "gmp-impl.h" 32 1.1 mrg 33 1.1 mrg int 34 1.1.1.2 mrg mpf_cmp (mpf_srcptr u, mpf_srcptr v) __GMP_NOTHROW 35 1.1 mrg { 36 1.1 mrg mp_srcptr up, vp; 37 1.1 mrg mp_size_t usize, vsize; 38 1.1 mrg mp_exp_t uexp, vexp; 39 1.1 mrg int cmp; 40 1.1 mrg int usign; 41 1.1 mrg 42 1.1.1.3 mrg usize = SIZ(u); 43 1.1.1.3 mrg vsize = SIZ(v); 44 1.1.1.3 mrg usign = usize >= 0 ? 1 : -1; 45 1.1 mrg 46 1.1 mrg /* 1. Are the signs different? */ 47 1.1 mrg if ((usize ^ vsize) >= 0) 48 1.1 mrg { 49 1.1 mrg /* U and V are both non-negative or both negative. */ 50 1.1 mrg if (usize == 0) 51 1.1 mrg /* vsize >= 0 */ 52 1.1 mrg return -(vsize != 0); 53 1.1 mrg if (vsize == 0) 54 1.1 mrg /* usize >= 0 */ 55 1.1 mrg return usize != 0; 56 1.1 mrg /* Fall out. */ 57 1.1 mrg } 58 1.1 mrg else 59 1.1 mrg { 60 1.1 mrg /* Either U or V is negative, but not both. */ 61 1.1.1.3 mrg return usign; 62 1.1 mrg } 63 1.1 mrg 64 1.1 mrg /* U and V have the same sign and are both non-zero. */ 65 1.1 mrg 66 1.1.1.3 mrg uexp = EXP(u); 67 1.1.1.3 mrg vexp = EXP(v); 68 1.1 mrg 69 1.1 mrg /* 2. Are the exponents different? */ 70 1.1 mrg if (uexp > vexp) 71 1.1 mrg return usign; 72 1.1 mrg if (uexp < vexp) 73 1.1 mrg return -usign; 74 1.1 mrg 75 1.1 mrg usize = ABS (usize); 76 1.1 mrg vsize = ABS (vsize); 77 1.1 mrg 78 1.1.1.3 mrg up = PTR (u); 79 1.1.1.3 mrg vp = PTR (v); 80 1.1 mrg 81 1.1 mrg #define STRICT_MPF_NORMALIZATION 0 82 1.1 mrg #if ! STRICT_MPF_NORMALIZATION 83 1.1 mrg /* Ignore zeroes at the low end of U and V. */ 84 1.1.1.3 mrg do { 85 1.1.1.3 mrg mp_limb_t tl; 86 1.1.1.3 mrg tl = up[0]; 87 1.1.1.3 mrg MPN_STRIP_LOW_ZEROS_NOT_ZERO (up, usize, tl); 88 1.1.1.3 mrg tl = vp[0]; 89 1.1.1.3 mrg MPN_STRIP_LOW_ZEROS_NOT_ZERO (vp, vsize, tl); 90 1.1.1.3 mrg } while (0); 91 1.1 mrg #endif 92 1.1 mrg 93 1.1 mrg if (usize > vsize) 94 1.1 mrg { 95 1.1 mrg cmp = mpn_cmp (up + usize - vsize, vp, vsize); 96 1.1.1.3 mrg /* if (cmp == 0) */ 97 1.1.1.3 mrg /* return usign; */ 98 1.1.1.3 mrg ++cmp; 99 1.1 mrg } 100 1.1 mrg else if (vsize > usize) 101 1.1 mrg { 102 1.1 mrg cmp = mpn_cmp (up, vp + vsize - usize, usize); 103 1.1.1.3 mrg /* if (cmp == 0) */ 104 1.1.1.3 mrg /* return -usign; */ 105 1.1 mrg } 106 1.1 mrg else 107 1.1 mrg { 108 1.1 mrg cmp = mpn_cmp (up, vp, usize); 109 1.1 mrg if (cmp == 0) 110 1.1 mrg return 0; 111 1.1 mrg } 112 1.1 mrg return cmp > 0 ? usign : -usign; 113 1.1 mrg } 114
Indexes created Sun Apr 26 00:22:38 UTC 2026