Category

Development Platform
Home > SourceCode/Document > Mathimatics-Numerical algorithms > Algorithm > View Code
t-divis_2exp.c
Package: gmp-5.0.1.tar.gz　[view]
Upload User: qaz666999
Upload Date: 2022-08-06
Package Size: 2570k
Code Size: 3k
Category:

Algorithm

Development Platform:

Unix_Linux

t-divis_2exp.c：Code Content
  1. /* test mpz_divisible_2exp_p */
  3. /*
  4. Copyright 2001 Free Software Foundation, Inc.
  6. This file is part of the GNU MP Library.
  8. The GNU MP Library is free software; you can redistribute it and/or modify
  9. it under the terms of the GNU Lesser General Public License as published by
  10. the Free Software Foundation; either version 3 of the License, or (at your
  11. option) any later version.
  13. The GNU MP Library is distributed in the hope that it will be useful, but
  14. WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  15. or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  16. License for more details.
  18. You should have received a copy of the GNU Lesser General Public License
  19. along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  21. #include <stdio.h>
  22. #include <stdlib.h>
  24. #include "gmp.h"
  25. #include "gmp-impl.h"
  26. #include "tests.h"
  29. void
  30. check_one (mpz_srcptr a, unsigned long d, int want)
  31. {
  32.   int   got;
  34.   got = (mpz_divisible_2exp_p (a, d) != 0);
  35.   if (want != got)
  36.     {
  37.       printf ("mpz_divisible_2exp_p wrongn");
  38.       printf ("   expected %d got %dn", want, got);
  39.       mpz_trace ("   a", a);
  40.       printf    ("   d=%lun", d);
  41.       mp_trace_base = -16;
  42.       mpz_trace ("   a", a);
  43.       printf    ("   d=0x%lXn", d);
  44.       abort ();
  45.     }
  46. }
  48. void
  49. check_data (void)
  50. {
  51.   static const struct {
  52.     const char    *a;
  53.     unsigned long d;
  54.     int           want;
  56.   } data[] = {
  58.     { "0", 0, 1 },
  59.     { "0", 1, 1 },
  60.     { "0", 2, 1 },
  61.     { "0", 3, 1 },
  63.     { "1", 0, 1 },
  64.     { "1", 1, 0 },
  65.     { "1", 2, 0 },
  66.     { "1", 3, 0 },
  67.     { "1", 10000, 0 },
  69.     { "4", 0, 1 },
  70.     { "4", 1, 1 },
  71.     { "4", 2, 1 },
  72.     { "4", 3, 0 },
  73.     { "4", 4, 0 },
  74.     { "4", 10000, 0 },
  76.     { "0x80000000", 31, 1 },
  77.     { "0x80000000", 32, 0 },
  78.     { "0x80000000", 64, 0 },
  80.     { "0x100000000", 32, 1 },
  81.     { "0x100000000", 33, 0 },
  82.     { "0x100000000", 64, 0 },
  84.     { "0x8000000000000000", 63, 1 },
  85.     { "0x8000000000000000", 64, 0 },
  86.     { "0x8000000000000000", 128, 0 },
  88.     { "0x10000000000000000", 64, 1 },
  89.     { "0x10000000000000000", 65, 0 },
  90.     { "0x10000000000000000", 128, 0 },
  91.     { "0x10000000000000000", 256, 0 },
  93.     { "0x10000000000000000100000000", 32, 1 },
  94.     { "0x10000000000000000100000000", 33, 0 },
  95.     { "0x10000000000000000100000000", 64, 0 },
  97.     { "0x1000000000000000010000000000000000", 64, 1 },
  98.     { "0x1000000000000000010000000000000000", 65, 0 },
  99.     { "0x1000000000000000010000000000000000", 128, 0 },
  100.     { "0x1000000000000000010000000000000000", 256, 0 },
  101.     { "0x1000000000000000010000000000000000", 1024, 0 },
  103.   };
  105.   mpz_t   a, d;
  106.   int     i;
  108.   mpz_init (a);
  109.   mpz_init (d);
  111.   for (i = 0; i < numberof (data); i++)
  112.     {
  113.       mpz_set_str_or_abort (a, data[i].a, 0);
  114.       check_one (a, data[i].d, data[i].want);
  116.       mpz_neg (a, a);
  117.       check_one (a, data[i].d, data[i].want);
  118.     }
  120.   mpz_clear (a);
  121.   mpz_clear (d);
  122. }
  124. int
  125. main (int argc, char *argv[])
  126. {
  127.   tests_start ();
  129.   check_data ();
  131.   tests_end ();
  132.   exit (0);
  133. }