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

数学计算

开发平台:

Unix_Linux

  1. /* Exercise mpf_mul_ui.
  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, mpf_srcptr u, unsigned long v)
  21. {
  22.   mp_size_t  usize, usign;
  23.   mp_ptr     wp;
  24.   mpf_t      want;
  25.   MPF_CHECK_FORMAT (got);
  26.   /* this code not nailified yet */
  27.   ASSERT_ALWAYS (BITS_PER_ULONG <= GMP_NUMB_BITS);
  28.   usign = SIZ (u);
  29.   usize = ABS (usign);
  30.   wp = refmpn_malloc_limbs (usize + 1);
  31.   wp[usize] = mpn_mul_1 (wp, PTR(u), usize, (mp_limb_t) v);
  32.   PTR(want) = wp;
  33.   SIZ(want) = (usign >= 0 ? usize+1 : -(usize+1));
  34.   EXP(want) = EXP(u) + 1;
  35.   refmpf_normalize (want);
  36.   if (! refmpf_validate ("mpf_mul_ui", got, want))
  37.     {
  38.       mp_trace_base = -16;
  39.       printf    ("  %sn", desc);
  40.       mpf_trace ("  u", u);
  41.       printf    ("  v %ld  0x%lXn", v, v);
  42.       abort ();
  43.     }
  44.   free (wp);
  45. }
  46. void
  47. check_rand (void)
  48. {
  49.   unsigned long  min_prec = __GMPF_BITS_TO_PREC (1);
  50.   gmp_randstate_ptr  rands = RANDS;
  51.   mpf_t              got, u;
  52.   unsigned long      prec, v;
  53.   int                i;
  54.   /* The nails code in mpf_mul_ui currently isn't exact, so suppress these
  55.      tests for now.  */
  56.   if (BITS_PER_ULONG > GMP_NUMB_BITS)
  57.     return;
  58.   mpf_init (got);
  59.   mpf_init (u);
  60.   for (i = 0; i < 200; i++)
  61.     {
  62.       /* got precision */
  63.       prec = min_prec + gmp_urandomm_ui (rands, 15L);
  64.       refmpf_set_prec_limbs (got, prec);
  65.       /* u precision */
  66.       prec = min_prec + gmp_urandomm_ui (rands, 15L);
  67.       refmpf_set_prec_limbs (u, prec);
  68.       /* u, possibly negative */
  69.       mpf_random2 (u, PREC(u), (mp_exp_t) 20);
  70.       if (gmp_urandomb_ui (rands, 1L))
  71.         mpf_neg (u, u);
  72.       /* v, 0 to BITS_PER_ULONG bits (inclusive) */
  73.       prec = gmp_urandomm_ui (rands, BITS_PER_ULONG+1);
  74.       v = gmp_urandomb_ui (rands, prec);
  75.       if ((i % 2) == 0)
  76.         {
  77.           /* separate */
  78.           mpf_mul_ui (got, u, v);
  79.           check_one ("separate", got, u, v);
  80.         }
  81.       else
  82.         {
  83.           /* overlap */
  84.           prec = refmpf_set_overlap (got, u);
  85.           mpf_mul_ui (got, got, v);
  86.           check_one ("overlap src==dst", got, u, v);
  87.           mpf_set_prec_raw (got, prec);
  88.         }
  89.     }
  90.   mpf_clear (got);
  91.   mpf_clear (u);
  92. }
  93. void
  94. check_various (void)
  95. {
  96.   mpf_t  u, got, want;
  97.   char   *s;
  98.   mpf_init2 (u,    2*8*sizeof(long));
  99.   mpf_init2 (got,  2*8*sizeof(long));
  100.   mpf_init2 (want, 2*8*sizeof(long));
  101.   s = "0 * ULONG_MAX";
  102.   mpf_set_ui (u, 0L);
  103.   mpf_mul_ui (got, u, ULONG_MAX);
  104.   MPF_CHECK_FORMAT (got);
  105.   mpf_set_ui (want, 0L);
  106.   if (mpf_cmp (got, want) != 0)
  107.     {
  108.     error:
  109.       printf ("Wrong result from %sn", s);
  110.       mpf_trace ("u   ", u);
  111.       mpf_trace ("got ", got);
  112.       mpf_trace ("want", want);
  113.       abort ();
  114.     }
  115.   s = "1 * ULONG_MAX";
  116.   mpf_set_ui (u, 1L);
  117.   mpf_mul_ui (got, u, ULONG_MAX);
  118.   MPF_CHECK_FORMAT (got);
  119.   mpf_set_ui (want, ULONG_MAX);
  120.   if (mpf_cmp (got, want) != 0)
  121.     goto error;
  122.   mpf_clear (u);
  123.   mpf_clear (got);
  124.   mpf_clear (want);
  125. }
  126. int
  127. main (void)
  128. {
  129.   tests_start ();
  130.   check_various ();
  131.   check_rand ();
  132.   tests_end ();
  133.   exit (0);
  134. }