t-ui_div.c
上传用户:qaz666999
上传日期:2022-08-06
资源大小:2570k
文件大小:3k
源码类别:

数学计算

开发平台:

Unix_Linux

  1. /* Test mpf_ui_div.
  2. Copyright 2004 Free Software Foundation, Inc.
  3. This file is part of the GNU MP Library.
  4. The GNU MP Library is free software; you can redistribute it and/or modify
  5. it under the terms of the GNU Lesser General Public License as published by
  6. the Free Software Foundation; either version 3 of the License, or (at your
  7. option) any later version.
  8. The GNU MP Library is distributed in the hope that it will be useful, but
  9. WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  10. or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  11. License for more details.
  12. You should have received a copy of the GNU Lesser General Public License
  13. along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  14. #include <stdio.h>
  15. #include <stdlib.h>
  16. #include "gmp.h"
  17. #include "gmp-impl.h"
  18. #include "tests.h"
  19. void
  20. check_one (const char *desc, mpf_ptr got, unsigned long u, mpf_srcptr v)
  21. {
  22.   mpf_t      uf;
  23.   mp_limb_t  ulimbs[2];
  24.   mp_size_t  usize;
  25.   ulimbs[0] = u & GMP_NUMB_MASK;
  26.   usize = (u != 0);
  27. #if BITS_PER_ULONG > GMP_NUMB_BITS
  28.   u >>= GMP_NUMB_BITS;
  29.   ulimbs[1] = u;
  30.   usize += (u != 0);
  31. #endif
  32.   PTR(uf) = ulimbs;
  33.   SIZ(uf) = usize;
  34.   EXP(uf) = usize;
  35.   if (! refmpf_validate_division ("mpf_ui_div", got, uf, v))
  36.     {
  37.       mp_trace_base = -16;
  38.       printf    ("  u 0x%lX  (%lu)n", u, u);
  39.       mpf_trace ("  v", v);
  40.       printf    ("  %sn", desc);
  41.       abort ();
  42.     }
  43. }
  44. void
  45. check_rand (void)
  46. {
  47.   unsigned long  min_prec = __GMPF_BITS_TO_PREC (1);
  48.   gmp_randstate_ptr  rands = RANDS;
  49.   unsigned long  prec, u;
  50.   mpf_t  got, v;
  51.   int    i;
  52.   mpf_init (got);
  53.   mpf_init (v);
  54.   for (i = 0; i < 200; i++)
  55.     {
  56.       /* got precision */
  57.       prec = min_prec + gmp_urandomm_ui (rands, 15L);
  58.       refmpf_set_prec_limbs (got, prec);
  59.       /* u */
  60.       prec = gmp_urandomm_ui (rands, BITS_PER_ULONG+1);
  61.       u = gmp_urandomb_ui (rands, prec);
  62.       /* v precision */
  63.       prec = min_prec + gmp_urandomm_ui (rands, 15L);
  64.       refmpf_set_prec_limbs (v, prec);
  65.       /* v, non-zero */
  66.       do {
  67.         mpf_random2 (v, PREC(v), (mp_exp_t) 20);
  68.       } while (SIZ(v) == 0);
  69.       /* v possibly negative */
  70.       if (gmp_urandomb_ui (rands, 1L))
  71.         mpf_neg (v, v);
  72.       if ((i % 2) == 0)
  73.         {
  74.           /* src != dst */
  75.           mpf_ui_div (got, u, v);
  76.           check_one ("separate", got, u, v);
  77.         }
  78.       else
  79.         {
  80.           /* src == dst */
  81.           prec = refmpf_set_overlap (got, v);
  82.           mpf_ui_div (got, u, got);
  83.           check_one ("overlap src==dst", got, u, v);
  84.           mpf_set_prec_raw (got, prec);
  85.         }
  86.     }
  87.   mpf_clear (got);
  88.   mpf_clear (v);
  89. }
  90. void
  91. check_various (void)
  92. {
  93.   mpf_t got, v;
  94.   mpf_init (got);
  95.   mpf_init (v);
  96.   /* 100/4 == 25 */
  97.   mpf_set_prec (got, 20L);
  98.   mpf_set_ui (v, 4L);
  99.   mpf_ui_div (got, 100L, v);
  100.   MPF_CHECK_FORMAT (got);
  101.   ASSERT_ALWAYS (mpf_cmp_ui (got, 25L) == 0);
  102.   {
  103.     /* 1/(2^n+1), a case where truncating the divisor would be wrong */
  104.     unsigned long  u = 1L;
  105.     mpf_set_prec (got, 500L);
  106.     mpf_set_prec (v, 900L);
  107.     mpf_set_ui (v, 1L);
  108.     mpf_mul_2exp (v, v, 800L);
  109.     mpf_add_ui (v, v, 1L);
  110.     mpf_ui_div (got, u, v);
  111.     check_one ("1/2^n+1, separate", got, u, v);
  112.   }
  113.   mpf_clear (got);
  114.   mpf_clear (v);
  115. }
  116. int
  117. main (void)
  118. {
  119.   tests_start ();
  120.   check_various ();
  121.   check_rand ();
  122.   tests_end ();
  123.   exit (0);
  124. }