xref: /src/external/lgpl3/gmp/dist/mini-gmp/tests/t-mpq_addsub.c

/*

Copyright 2012, 2013, 2018 Free Software Foundation, Inc.

This file is part of the GNU MP Library test suite.

The GNU MP Library test suite is free software; you can redistribute it
and/or modify it under the terms of the GNU 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 MP Library test suite 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 General
Public License for more details.

You should have received a copy of the GNU General Public License along with
the GNU MP Library test suite.  If not, see https://www.gnu.org/licenses/.  */

#include <assert.h>
#include <stdlib.h>
#include <stdio.h>

#include "testutils.h"
#include "../mini-mpq.h"

#define MAXBITS 300
#define COUNT 10000

static void
_mpq_set_zz (mpq_t q, mpz_t n, mpz_t d)
{
  if (mpz_fits_ulong_p (d) && mpz_fits_slong_p (n))
    {
      mpq_set_si (q, mpz_get_si (n), mpz_get_ui (d));
    }
  else if (mpz_fits_ulong_p (d) && mpz_fits_ulong_p (n))
    {
      mpq_set_ui (q, mpz_get_ui (n), mpz_get_ui (d));
    }
  else
    {
      mpq_set_num (q, n);
      mpq_set_den (q, d);
    }
  mpq_canonicalize (q);
}

void
testcmpui ()
{
  unsigned d1, d2, n1, n2;
  mpq_t q1, q2;

  mpq_init (q1);
  mpq_init (q2);

  for (d1 = 1; d1 < 6; d1 += 2)
    for (n1 = 1; n1 < 6; n1 *= 2)
      {
	mpq_set_ui (q1, n1, d1);
	for (d2 = 1; d2 < 6; d2 += 2)
	  for (n2 = 1; n2 < 6; n2 *= 2)
	    {
	      int fres = mpq_cmp_ui (q1, n2, d2);
	      int ref = (d1*n2 < d2*n1) - (d1*n2 > d2*n1);

	      mpq_set_ui (q2, n2, d2);

	      if ((!ref) != mpq_equal (q1, q2))
		{
		  fprintf (stderr, "mpz_equal failed: %i / %i = %i / %i ? %i\n", n1, d1, n2, d2, ref);
		  abort ();
		}

	      if (ref != fres)
		{
		  fprintf (stderr, "mpz_cmp_ui failed: %i / %i = %i / %i ? %i != %i\n", n1, d1, n2, d2, ref, fres);
		  abort ();
		}
	    }
      }

  mpq_clear (q1);
  mpq_clear (q2);
}

void
testmain (int argc, char **argv)
{
  unsigned i;
  mpz_t a, b, q, r, c;
  mpq_t rr, ii, ff;
  int tst;

  testcmpui ();
  mpz_init (a);
  mpz_init (b);
  mpz_init (r);
  mpz_init (q);
  mpz_init (c);
  mpq_init (rr);
  mpq_init (ff);
  mpq_init (ii);

  for (i = 0; i < COUNT; i++)
    {
      mini_random_op4 (OP_TDIV, MAXBITS, a, b, q, r);

      _mpq_set_zz (rr, a, b);
      _mpq_set_zz (ff, r, b);

      mpq_set_z (ii, q);

      mpz_set_q (c, rr);
      if (mpz_cmp (c, q))
	{
	  fprintf (stderr, "mpz_set_q failed:\n");
	  dump ("a", a);
	  dump ("b", b);
	  dump ("c", c);
	  dump ("q", q);
	  abort ();
	}

      if ((mpz_sgn (r) != 0) ^ (mpz_cmp_ui (mpq_denref (rr), 1) != 0))
	{
	  fprintf (stderr, "mpq_canonicalize failed:\n");
	  dump ("a", a);
	  dump ("b", b);
	  dump ("r", r);
	  dump ("D", mpq_denref (rr));
	  abort ();
	}

      if (i & 1)
	{
	  if (mpz_fits_slong_p (q))
	    tst = mpq_cmp_si (rr, mpz_get_si (q), 1);
	  else if (mpz_fits_ulong_p (q))
	    tst = mpq_cmp_ui (rr, mpz_get_ui (q), 1);
	  else
	    tst = mpq_cmp_z (rr, q);
	  if (mpz_sgn (b) < 0)
	    tst = - tst;
	  if ((tst != mpz_sgn (r)) && ((tst < 0 && mpz_sgn (r) >= 0) || (tst > 0 && mpz_sgn (r) <= 0)))
	    {
	      fprintf (stderr, "mpq_cmp ii failed: %i %i\n", tst, mpz_sgn (r));
	      dump ("a", a);
	      dump ("b", b);
	      dump ("r", r);
	      dump ("q", q);
	      abort ();
	    }
	}
      else
	{
	  if (mpz_fits_ulong_p (b) && mpz_fits_slong_p (r))
	    tst = mpq_cmp_si (rr, mpz_get_si (r), mpz_get_ui (b));
	  else if (mpz_fits_ulong_p (b) && mpz_fits_ulong_p (r))
	    tst = mpq_cmp_ui (rr, mpz_get_ui (r), mpz_get_ui (b));
	  else
	    tst = mpq_cmp (rr, ff);
	  if ((tst != mpz_sgn (q)) && ((tst < 0 && mpz_sgn (q) >= 0) || (tst > 0 && mpz_sgn (q) <= 0)))
	    {
	      fprintf (stderr, "mpq_cmp ff failed: %i %i\n", tst, mpz_sgn (q));
	      dump ("a", a);
	      dump ("b", b);
	      dump ("r", r);
	      dump ("q", q);
	      abort ();
	    }
	}

      if (i & 1)
	{
	  mpq_sub (rr, rr, ff);
	}
      else
	{
	  mpq_neg (ff, ff);
	  mpq_add (rr, ff, rr);
	}

      if (!mpq_equal (ii, rr))
	{
	  fprintf (stderr, "mpq_%s failed:\n", (i & 1) ? "sub" : "add");
	  dump ("a", a);
	  dump ("b", b);
	  dump ("r", r);
	  dump ("q", q);
	  abort ();
	}
    }

  mpz_clear (a);
  mpz_clear (b);
  mpz_clear (r);
  mpz_clear (q);
  mpz_clear (c);
  mpq_clear (rr);
  mpq_clear (ff);
  mpq_clear (ii);
}