tests/comparisons.c

   1 /* Comparison functions.
   2
   3 Copyright (C) 2008, 2009 Philippe Th\'eveny
   4
   5 This file is part of the MPC Library.
   6
   7 The MPC Library is free software; you can redistribute it and/or modify
   8 it under the terms of the GNU Lesser General Public License as published by
   9 the Free Software Foundation; either version 2.1 of the License, or (at your
  10 option) any later version.
  11
  12 The MPC Library is distributed in the hope that it will be useful, but
  13 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  14 or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  15 License for more details.
  16
  17 You should have received a copy of the GNU Lesser General Public License
  18 along with the MPC Library; see the file COPYING.LIB.  If not, write to
  19 the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
  20 MA 02111-1307, USA. */
  21
  22 #include "mpc-tests.h"
  23
  24 /* comparisons, see description in mpc-tests.h */
  25 int
  26 same_mpfr_value (mpfr_ptr got, mpfr_ptr ref, int known_sign)
  27 {
  28    /* The sign of zeroes and infinities is checked only when
  29       known_sign is true.                                    */
  30    if (mpfr_nan_p (got))
  31       return mpfr_nan_p (ref);
  32    if (mpfr_inf_p (got))
  33       return mpfr_inf_p (ref) &&
  34             (!known_sign || mpfr_signbit (got) == mpfr_signbit (ref));
  35    if (mpfr_zero_p (got))
  36       return mpfr_zero_p (ref) &&
  37             (!known_sign || mpfr_signbit (got) == mpfr_signbit (ref));
  38    return mpfr_cmp (got, ref) == 0;
  39 }
  40
  41 int
  42 same_mpc_value (mpc_ptr got, mpc_ptr ref, known_signs_t known_signs)
  43 {
  44    return    same_mpfr_value (MPC_RE (got), MPC_RE (ref), known_signs.re)
  45           && same_mpfr_value (MPC_IM (got), MPC_IM (ref), known_signs.im);
  46 }