tests/mpq/t-aors.c

   1 /* Test mpq_add and mpq_sub.
   2
   3 Copyright 2001 Free Software Foundation, Inc.
   4
   5 This file is part of the GNU MP Library.
   6
   7 The GNU MP 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 3 of the License, or (at your
  10 option) any later version.
  11
  12 The GNU MP 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 GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  19
  20 #include "config.h"
  21
  22 #include <stdio.h>
  23 #include <stdlib.h>
  24 #include <string.h>
  25
  26 #include "gmp.h"
  27 #include "gmp-impl.h"
  28 #include "tests.h"
  29
  30
  31 void
  32 check_all (mpq_ptr x, mpq_ptr y, mpq_ptr want_add, mpq_ptr want_sub)
  33 {
  34   mpq_t  got;
  35   int    neg_x, neg_y, swap;
  36
  37   mpq_init (got);
  38
  39   MPQ_CHECK_FORMAT (want_add);
  40   MPQ_CHECK_FORMAT (want_sub);
  41   MPQ_CHECK_FORMAT (x);
  42   MPQ_CHECK_FORMAT (y);
  43
  44   for (swap = 0; swap <= 1; swap++)
  45     {
  46       for (neg_x = 0; neg_x <= 1; neg_x++)
  47         {
  48           for (neg_y = 0; neg_y <= 1; neg_y++)
  49             {
  50               mpq_add (got, x, y);
  51               MPQ_CHECK_FORMAT (got);
  52               if (! mpq_equal (got, want_add))
  53                 {
  54                   printf ("mpq_add wrong\n");
  55                   mpq_trace ("  x   ", x);
  56                   mpq_trace ("  y   ", y);
  57                   mpq_trace ("  got ", got);
  58                   mpq_trace ("  want", want_add);
  59                   abort ();
  60                 }
  61
  62               mpq_sub (got, x, y);
  63               MPQ_CHECK_FORMAT (got);
  64               if (! mpq_equal (got, want_sub))
  65                 {
  66                   printf ("mpq_sub wrong\n");
  67                   mpq_trace ("  x   ", x);
  68                   mpq_trace ("  y   ", y);
  69                   mpq_trace ("  got ", got);
  70                   mpq_trace ("  want", want_sub);
  71                   abort ();
  72                 }
  73
  74
  75               mpq_neg (y, y);
  76               mpq_swap (want_add, want_sub);
  77             }
  78
  79           mpq_neg (x, x);
  80           mpq_swap (want_add, want_sub);
  81           mpq_neg (want_add, want_add);
  82           mpq_neg (want_sub, want_sub);
  83         }
  84
  85       mpq_swap (x, y);
  86       mpq_neg (want_sub, want_sub);
  87     }
  88
  89   mpq_clear (got);
  90 }
  91
  92
  93 void
  94 check_data (void)
  95 {
  96   static const struct {
  97     const char  *x;
  98     const char  *y;
  99     const char  *want_add;
 100     const char  *want_sub;
 101
 102   } data[] = {
 103
 104     { "0", "0", "0", "0" },
 105     { "1", "0", "1", "1" },
 106     { "1", "1", "2", "0" },
 107
 108     { "1/2", "1/2", "1", "0" },
 109     { "5/6", "14/15", "53/30", "-1/10" },
 110   };
 111
 112   mpq_t  x, y, want_add, want_sub;
 113   int i;
 114
 115   mpq_init (x);
 116   mpq_init (y);
 117   mpq_init (want_add);
 118   mpq_init (want_sub);
 119
 120   for (i = 0; i < numberof (data); i++)
 121     {
 122       mpq_set_str_or_abort (x, data[i].x, 0);
 123       mpq_set_str_or_abort (y, data[i].y, 0);
 124       mpq_set_str_or_abort (want_add, data[i].want_add, 0);
 125       mpq_set_str_or_abort (want_sub, data[i].want_sub, 0);
 126
 127       check_all (x, y, want_add, want_sub);
 128     }
 129
 130   mpq_clear (x);
 131   mpq_clear (y);
 132   mpq_clear (want_add);
 133   mpq_clear (want_sub);
 134 }
 135
 136
 137 void
 138 check_rand (void)
 139 {
 140   mpq_t  x, y, want_add, want_sub;
 141   int i;
 142   gmp_randstate_ptr  rands = RANDS;
 143
 144   mpq_init (x);
 145   mpq_init (y);
 146   mpq_init (want_add);
 147   mpq_init (want_sub);
 148
 149   for (i = 0; i < 500; i++)
 150     {
 151       mpz_errandomb (mpq_numref(x), rands, 512L);
 152       mpz_errandomb_nonzero (mpq_denref(x), rands, 512L);
 153       mpq_canonicalize (x);
 154
 155       mpz_errandomb (mpq_numref(y), rands, 512L);
 156       mpz_errandomb_nonzero (mpq_denref(y), rands, 512L);
 157       mpq_canonicalize (y);
 158
 159       refmpq_add (want_add, x, y);
 160       refmpq_sub (want_sub, x, y);
 161
 162       check_all (x, y, want_add, want_sub);
 163     }
 164
 165   mpq_clear (x);
 166   mpq_clear (y);
 167   mpq_clear (want_add);
 168   mpq_clear (want_sub);
 169 }
 170
 171
 172 int
 173 main (void)
 174 {
 175   tests_start ();
 176
 177   check_data ();
 178   check_rand ();
 179
 180   tests_end ();
 181
 182   exit (0);
 183 }