Blame tests/mpq/t-aors.c

Packit 5c3484
/* Test mpq_add and mpq_sub.
Packit 5c3484
Packit 5c3484
Copyright 2001 Free Software Foundation, Inc.
Packit 5c3484
Packit 5c3484
This file is part of the GNU MP Library test suite.
Packit 5c3484
Packit 5c3484
The GNU MP Library test suite is free software; you can redistribute it
Packit 5c3484
and/or modify it under the terms of the GNU General Public License as
Packit 5c3484
published by the Free Software Foundation; either version 3 of the License,
Packit 5c3484
or (at your option) any later version.
Packit 5c3484
Packit 5c3484
The GNU MP Library test suite is distributed in the hope that it will be
Packit 5c3484
useful, but WITHOUT ANY WARRANTY; without even the implied warranty of
Packit 5c3484
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General
Packit 5c3484
Public License for more details.
Packit 5c3484
Packit 5c3484
You should have received a copy of the GNU General Public License along with
Packit 5c3484
the GNU MP Library test suite.  If not, see https://www.gnu.org/licenses/.  */
Packit 5c3484
Packit 5c3484
#include "config.h"
Packit 5c3484
Packit 5c3484
#include <stdio.h>
Packit 5c3484
#include <stdlib.h>
Packit 5c3484
#include <string.h>
Packit 5c3484
Packit 5c3484
#include "gmp.h"
Packit 5c3484
#include "gmp-impl.h"
Packit 5c3484
#include "tests.h"
Packit 5c3484
Packit 5c3484
Packit 5c3484
void
Packit 5c3484
check_all (mpq_ptr x, mpq_ptr y, mpq_ptr want_add, mpq_ptr want_sub)
Packit 5c3484
{
Packit 5c3484
  mpq_t  got;
Packit 5c3484
  int    neg_x, neg_y, swap;
Packit 5c3484
Packit 5c3484
  mpq_init (got);
Packit 5c3484
Packit 5c3484
  MPQ_CHECK_FORMAT (want_add);
Packit 5c3484
  MPQ_CHECK_FORMAT (want_sub);
Packit 5c3484
  MPQ_CHECK_FORMAT (x);
Packit 5c3484
  MPQ_CHECK_FORMAT (y);
Packit 5c3484
Packit 5c3484
  for (swap = 0; swap <= 1; swap++)
Packit 5c3484
    {
Packit 5c3484
      for (neg_x = 0; neg_x <= 1; neg_x++)
Packit 5c3484
        {
Packit 5c3484
          for (neg_y = 0; neg_y <= 1; neg_y++)
Packit 5c3484
            {
Packit 5c3484
              mpq_add (got, x, y);
Packit 5c3484
              MPQ_CHECK_FORMAT (got);
Packit 5c3484
              if (! mpq_equal (got, want_add))
Packit 5c3484
                {
Packit 5c3484
                  printf ("mpq_add wrong\n");
Packit 5c3484
                  mpq_trace ("  x   ", x);
Packit 5c3484
                  mpq_trace ("  y   ", y);
Packit 5c3484
                  mpq_trace ("  got ", got);
Packit 5c3484
                  mpq_trace ("  want", want_add);
Packit 5c3484
                  abort ();
Packit 5c3484
                }
Packit 5c3484
Packit 5c3484
              mpq_sub (got, x, y);
Packit 5c3484
              MPQ_CHECK_FORMAT (got);
Packit 5c3484
              if (! mpq_equal (got, want_sub))
Packit 5c3484
                {
Packit 5c3484
                  printf ("mpq_sub wrong\n");
Packit 5c3484
                  mpq_trace ("  x   ", x);
Packit 5c3484
                  mpq_trace ("  y   ", y);
Packit 5c3484
                  mpq_trace ("  got ", got);
Packit 5c3484
                  mpq_trace ("  want", want_sub);
Packit 5c3484
                  abort ();
Packit 5c3484
                }
Packit 5c3484
Packit 5c3484
Packit 5c3484
              mpq_neg (y, y);
Packit 5c3484
              mpq_swap (want_add, want_sub);
Packit 5c3484
            }
Packit 5c3484
Packit 5c3484
          mpq_neg (x, x);
Packit 5c3484
          mpq_swap (want_add, want_sub);
Packit 5c3484
          mpq_neg (want_add, want_add);
Packit 5c3484
          mpq_neg (want_sub, want_sub);
Packit 5c3484
        }
Packit 5c3484
Packit 5c3484
      mpq_swap (x, y);
Packit 5c3484
      mpq_neg (want_sub, want_sub);
Packit 5c3484
    }
Packit 5c3484
Packit 5c3484
  mpq_clear (got);
Packit 5c3484
}
Packit 5c3484
Packit 5c3484
Packit 5c3484
void
Packit 5c3484
check_data (void)
Packit 5c3484
{
Packit 5c3484
  static const struct {
Packit 5c3484
    const char  *x;
Packit 5c3484
    const char  *y;
Packit 5c3484
    const char  *want_add;
Packit 5c3484
    const char  *want_sub;
Packit 5c3484
Packit 5c3484
  } data[] = {
Packit 5c3484
Packit 5c3484
    { "0", "0", "0", "0" },
Packit 5c3484
    { "1", "0", "1", "1" },
Packit 5c3484
    { "1", "1", "2", "0" },
Packit 5c3484
Packit 5c3484
    { "1/2", "1/2", "1", "0" },
Packit 5c3484
    { "5/6", "14/15", "53/30", "-1/10" },
Packit 5c3484
  };
Packit 5c3484
Packit 5c3484
  mpq_t  x, y, want_add, want_sub;
Packit 5c3484
  int i;
Packit 5c3484
Packit 5c3484
  mpq_init (x);
Packit 5c3484
  mpq_init (y);
Packit 5c3484
  mpq_init (want_add);
Packit 5c3484
  mpq_init (want_sub);
Packit 5c3484
Packit 5c3484
  for (i = 0; i < numberof (data); i++)
Packit 5c3484
    {
Packit 5c3484
      mpq_set_str_or_abort (x, data[i].x, 0);
Packit 5c3484
      mpq_set_str_or_abort (y, data[i].y, 0);
Packit 5c3484
      mpq_set_str_or_abort (want_add, data[i].want_add, 0);
Packit 5c3484
      mpq_set_str_or_abort (want_sub, data[i].want_sub, 0);
Packit 5c3484
Packit 5c3484
      check_all (x, y, want_add, want_sub);
Packit 5c3484
    }
Packit 5c3484
Packit 5c3484
  mpq_clear (x);
Packit 5c3484
  mpq_clear (y);
Packit 5c3484
  mpq_clear (want_add);
Packit 5c3484
  mpq_clear (want_sub);
Packit 5c3484
}
Packit 5c3484
Packit 5c3484
Packit 5c3484
void
Packit 5c3484
check_rand (void)
Packit 5c3484
{
Packit 5c3484
  mpq_t  x, y, want_add, want_sub;
Packit 5c3484
  int i;
Packit 5c3484
  gmp_randstate_ptr  rands = RANDS;
Packit 5c3484
Packit 5c3484
  mpq_init (x);
Packit 5c3484
  mpq_init (y);
Packit 5c3484
  mpq_init (want_add);
Packit 5c3484
  mpq_init (want_sub);
Packit 5c3484
Packit 5c3484
  for (i = 0; i < 500; i++)
Packit 5c3484
    {
Packit 5c3484
      mpz_errandomb (mpq_numref(x), rands, 512L);
Packit 5c3484
      mpz_errandomb_nonzero (mpq_denref(x), rands, 512L);
Packit 5c3484
      mpq_canonicalize (x);
Packit 5c3484
Packit 5c3484
      mpz_errandomb (mpq_numref(y), rands, 512L);
Packit 5c3484
      mpz_errandomb_nonzero (mpq_denref(y), rands, 512L);
Packit 5c3484
      mpq_canonicalize (y);
Packit 5c3484
Packit 5c3484
      refmpq_add (want_add, x, y);
Packit 5c3484
      refmpq_sub (want_sub, x, y);
Packit 5c3484
Packit 5c3484
      check_all (x, y, want_add, want_sub);
Packit 5c3484
    }
Packit 5c3484
Packit 5c3484
  mpq_clear (x);
Packit 5c3484
  mpq_clear (y);
Packit 5c3484
  mpq_clear (want_add);
Packit 5c3484
  mpq_clear (want_sub);
Packit 5c3484
}
Packit 5c3484
Packit 5c3484
Packit 5c3484
int
Packit 5c3484
main (void)
Packit 5c3484
{
Packit 5c3484
  tests_start ();
Packit 5c3484
Packit 5c3484
  check_data ();
Packit 5c3484
  check_rand ();
Packit 5c3484
Packit 5c3484
  tests_end ();
Packit 5c3484
Packit 5c3484
  exit (0);
Packit 5c3484
}