Home | History | Annotate | Line # | Download | only in mpq
t-aors.c revision 1.1 
      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 }
    184