/* Test file for mpfr_cmpabs and mpfr_cmpabs_ui.

Copyright 2004-2023 Free Software Foundation, Inc.
Contributed by the AriC and Caramba projects, INRIA.

This file is part of the GNU MPFR Library.

The GNU MPFR Library is free software; you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation; either version 3 of the License, or (at your
option) any later version.

The GNU MPFR Library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
License for more details.

You should have received a copy of the GNU Lesser General Public License
along with the GNU MPFR Library; see the file COPYING.LESSER.  If not, see
https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc.,
51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */

#include "mpfr-test.h"

#define PRINT_ERROR(s) do { printf ("Error: %s\n", s); exit (1); } while (0)

static int
cmpabs (mpfr_srcptr x, mpfr_srcptr y)
{
  unsigned int i;
  int r[4];
  mpfr_flags_t f1, f2, flags[2] = { 0, MPFR_FLAGS_ALL };
  mpfr_exp_t emin, emax;

  emin = mpfr_get_emin ();
  emax = mpfr_get_emax ();

  for (i = 0; i < 4; i++)
    {
      if (i & 2)
        {
          mpfr_exp_t ex = MPFR_IS_SINGULAR (x) ? emax : MPFR_GET_EXP (x);
          mpfr_exp_t ey = MPFR_IS_SINGULAR (y) ? emax : MPFR_GET_EXP (y);
          set_emin (ex < ey ? ex : ey);
          set_emax (ex < ey ? ey : ex);
        }

      __gmpfr_flags = f1 = flags[i % 2];
      r[i] = mpfr_cmpabs (x, y);
      f2 = __gmpfr_flags;
      if (MPFR_IS_NAN (x) || MPFR_IS_NAN (y))
        f1 |= MPFR_FLAGS_ERANGE;

      if (i & 2)
        {
          set_emin (emin);
          set_emax (emax);
        }

      if (f1 != f2)
        {
          printf ("Flags error in mpfr_cmpabs for i = %u\n  x = ", i);
          mpfr_dump (x);
          printf ("  y = ");
          mpfr_dump (y);
          printf ("Expected flags = ");
          flags_out (f1);
          printf ("Obtained flags = ");
          flags_out (f2);
          exit (1);
        }

      if (i > 0)
        MPFR_ASSERTN (r[i] == r[0]);
    }

  return r[0];
}

static void
test_cmpabs (void)
{
  mpfr_t xx, yy;

  mpfr_init2 (xx, 2);
  mpfr_init2 (yy, 2);

  mpfr_clear_erangeflag ();
  MPFR_SET_NAN (xx);
  MPFR_SET_NAN (yy);
  if (cmpabs (xx, yy) != 0)
    PRINT_ERROR ("mpfr_cmpabs (NAN,NAN) returns non-zero");
  if (!mpfr_erangeflag_p ())
    PRINT_ERROR ("mpfr_cmpabs (NAN,NAN) doesn't set erange flag");

  mpfr_set_si (yy, -1, MPFR_RNDN);
  if (cmpabs (xx, yy) != 0)
    PRINT_ERROR ("mpfr_cmpabs (NAN,-1) returns non-zero");
  if (cmpabs (yy, xx) != 0)
    PRINT_ERROR ("mpfr_cmpabs (-1,NAN) returns non-zero");

  mpfr_set_str_binary (xx, "0.10E0");
  mpfr_set_str_binary (yy, "-0.10E0");
  if (cmpabs (xx, yy) != 0)
    PRINT_ERROR ("mpfr_cmpabs (xx, yy) returns non-zero for prec=2");

  mpfr_set_prec (xx, 65);
  mpfr_set_prec (yy, 65);
  mpfr_set_str_binary (xx, "-0.10011010101000110101010000000011001001001110001011101011111011101E623");
  mpfr_set_str_binary (yy, "0.10011010101000110101010000000011001001001110001011101011111011100E623");
  if (cmpabs (xx, yy) <= 0)
    PRINT_ERROR ("Error (1) in mpfr_cmpabs");

  mpfr_set_str_binary (xx, "-0.10100010001110110111000010001000010011111101000100011101000011100");
  mpfr_set_str_binary (yy, "-0.10100010001110110111000010001000010011111101000100011101000011011");
  if (cmpabs (xx, yy) <= 0)
    PRINT_ERROR ("Error (2) in mpfr_cmpabs");

  mpfr_set_prec (xx, 160);
  mpfr_set_prec (yy, 160);
  mpfr_set_str_binary (xx, "0.1E1");
  mpfr_set_str_binary (yy, "-0.1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111100000110001110100");
  if (cmpabs (xx, yy) <= 0)
    PRINT_ERROR ("Error (3) in mpfr_cmpabs");

  mpfr_set_prec(xx, 53);
  mpfr_set_prec(yy, 200);
  mpfr_set_ui (xx, 1, (mpfr_rnd_t) 0);
  mpfr_set_ui (yy, 1, (mpfr_rnd_t) 0);
  if (cmpabs (xx, yy) != 0)
    PRINT_ERROR ("Error in mpfr_cmpabs: 1.0 != 1.0");

  mpfr_set_prec (yy, 31);
  mpfr_set_str (xx, "-1.0000000002", 10, (mpfr_rnd_t) 0);
  mpfr_set_ui (yy, 1, (mpfr_rnd_t) 0);
  if (cmpabs (xx, yy) <= 0)
    PRINT_ERROR ("Error in mpfr_cmpabs: not 1.0000000002 > 1.0");
  mpfr_set_prec(yy, 53);

  mpfr_set_ui(xx, 0, MPFR_RNDN);
  mpfr_set_str (yy, "-0.1", 10, MPFR_RNDN);
  if (cmpabs (xx, yy) >= 0)
    PRINT_ERROR ("Error in mpfr_cmpabs(0.0, 0.1)");

  mpfr_set_inf (xx, -1);
  mpfr_set_str (yy, "23489745.0329", 10, MPFR_RNDN);
  if (cmpabs (xx, yy) <= 0)
    PRINT_ERROR ("Error in mpfr_cmpabs(-Inf, 23489745.0329)");

  mpfr_set_inf (xx, 1);
  mpfr_set_inf (yy, -1);
  if (cmpabs (xx, yy) != 0)
    PRINT_ERROR ("Error in mpfr_cmpabs(Inf, -Inf)");

  mpfr_set_inf (yy, -1);
  mpfr_set_str (xx, "2346.09234", 10, MPFR_RNDN);
  if (cmpabs (xx, yy) >= 0)
    PRINT_ERROR ("Error in mpfr_cmpabs(-Inf, 2346.09234)");

  mpfr_set_prec (xx, 2);
  mpfr_set_prec (yy, 128);
  mpfr_set_str_binary (xx, "0.1E10");
  mpfr_set_str_binary (yy,
                       "0.100000000000000000000000000000000000000000000000"
                       "00000000000000000000000000000000000000000000001E10");
  if (cmpabs (xx, yy) >= 0)
    PRINT_ERROR ("Error in mpfr_cmpabs(10.235, 2346.09234)");
  mpfr_swap (xx, yy);
  if (cmpabs (xx, yy) <= 0)
    PRINT_ERROR ("Error in mpfr_cmpabs(2346.09234, 10.235)");
  mpfr_swap (xx, yy);

  mpfr_clear (xx);
  mpfr_clear (yy);
}

static int
cmpabs_ui (mpfr_srcptr x, unsigned long u)
{
  unsigned int i;
  int r[4];
  mpfr_flags_t f1, f2, flags[2] = { 0, MPFR_FLAGS_ALL };
  mpfr_exp_t emin, emax;

  emin = mpfr_get_emin ();
  emax = mpfr_get_emax ();

  for (i = 0; i < 4; i++)
    {
      if (i & 2)
        {
          mpfr_exp_t e = MPFR_IS_SINGULAR (x) ? emax : MPFR_GET_EXP (x);
          set_emin (e);
          set_emax (e);
        }

      __gmpfr_flags = f1 = flags[i % 2];
      r[i] = mpfr_cmpabs_ui (x, u);
      f2 = __gmpfr_flags;
      if (MPFR_IS_NAN (x))
        f1 |= MPFR_FLAGS_ERANGE;

      if (i & 2)
        {
          set_emin (emin);
          set_emax (emax);
        }

      if (f1 != f2)
        {
          printf ("Flags error in mpfr_cmpabs_ui for i = %u\n  x = ", i);
          mpfr_dump (x);
          printf ("  u = %lu\n", u);
          printf ("Expected flags = ");
          flags_out (f1);
          printf ("Obtained flags = ");
          flags_out (f2);
          exit (1);
        }

      if (i > 0)
        MPFR_ASSERTN (r[i] == r[0]);
    }

  return r[0];
}

static void
test_cmpabs_ui (void)
{
  mpfr_t x;

  mpfr_init2 (x, 64); /* should be enough so that ULONG_MAX fits */
  /* check NaN */
  mpfr_set_nan (x);
  if (cmpabs_ui (x, 17) != 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (NaN, 17)");
  /* check -Inf */
  mpfr_set_inf (x, -1);
  if (cmpabs_ui (x, 17) <= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (-Inf, 17)");
  /* check +Inf */
  mpfr_set_inf (x, +1);
  if (cmpabs_ui (x, 17) <= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (+Inf, 17)");
  /* check -1 */
  mpfr_set_si (x, -1, MPFR_RNDN);
  if (cmpabs_ui (x, 0) <= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (-1, 0)");
  if (cmpabs_ui (x, 17) >= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (-1, 17)");
  /* check -0 */
  mpfr_set_zero (x, -1);
  if (cmpabs_ui (x, 0) != 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (-0, 0)");
  if (cmpabs_ui (x, 17) >= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (-0, 17)");
  /* check +0 */
  mpfr_set_zero (x, +1);
  if (cmpabs_ui (x, 0) != 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (+0, 0)");
  if (cmpabs_ui (x, 17) >= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (+0, 17)");
  /* check 1 */
  mpfr_set_ui (x, 1, MPFR_RNDN);
  if (cmpabs_ui (x, 0) <= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (1, 0)");
  if (cmpabs_ui (x, 1) != 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (1, 1)");
  if (cmpabs_ui (x, 17) >= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (1, 17)");
  /* check ULONG_MAX */
  mpfr_set_ui (x, ULONG_MAX, MPFR_RNDN);
  if (cmpabs_ui (x, 0) <= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (ULONG_MAX, 0)");
  if (cmpabs_ui (x, 1) <= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (ULONG_MAX, 1)");
  if (cmpabs_ui (x, 17) <= 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (ULONG_MAX, 17)");
  if (cmpabs_ui (x, ULONG_MAX) != 0)
    PRINT_ERROR ("mpfr_cmpabs_ui (ULONG_MAX, ULONG_MAX)");
  mpfr_clear (x);
}

int
main (void)
{
  tests_start_mpfr ();

  test_cmpabs ();
  test_cmpabs_ui ();

  tests_end_mpfr ();
  return 0;
}
