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

数学计算

开发平台:

Unix_Linux

  1. /* Test mpf_get_str and mpf_set_str.
  2. Copyright 1996, 2000, 2001, 2008 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 <string.h> /* for strlen */
  17. #include "gmp.h"
  18. #include "gmp-impl.h"
  19. #include "tests.h"
  20. #ifndef SIZE
  21. #define SIZE 10
  22. #endif
  23. #ifndef EXPO
  24. #define EXPO 200
  25. #endif
  26. int
  27. main (int argc, char **argv)
  28. {
  29.   mpf_t x, y;
  30.   int reps = 20000;
  31.   int i;
  32.   mp_size_t bprec = 100;
  33.   mpf_t d, rerr, max_rerr, limit_rerr;
  34.   char *str;
  35.   mp_exp_t bexp;
  36.   long size, exp;
  37.   int base;
  38.   char buf[SIZE * GMP_LIMB_BITS + 5];
  39.   tests_start ();
  40.   if (argc > 1)
  41.     {
  42.       reps = strtol (argv[1], 0, 0);
  43.       if (argc > 2)
  44. bprec = strtol (argv[2], 0, 0);
  45.     }
  46.   mpf_set_default_prec (bprec);
  47.   mpf_init_set_ui (limit_rerr, 1);
  48.   mpf_div_2exp (limit_rerr, limit_rerr, bprec);
  49. #if VERBOSE
  50.   mpf_dump (limit_rerr);
  51. #endif
  52.   mpf_init (rerr);
  53.   mpf_init_set_ui (max_rerr, 0);
  54.   mpf_init (x);
  55.   mpf_init (y);
  56.   mpf_init (d);
  57.   /* First test some specific values.  */
  58.   mpf_set_str (y, "1.23456e1000", 0);
  59.   mpf_set_str (x, "1.23456e1000", 10);
  60.   if (mpf_cmp (x, y) != 0)
  61.     abort ();
  62.   mpf_set_str (x, "1.23456e+1000", 0);
  63.   if (mpf_cmp (x, y) != 0)
  64.     abort ();
  65.   mpf_set_str (x, "1.23456e+1000", 10);
  66.   if (mpf_cmp (x, y) != 0)
  67.     abort ();
  68.   /* Now test random values.  */
  69.   for (i = 0; i < reps; i++)
  70.     {
  71.       if (i == 0)
  72.         {
  73.           /* exercise the special case in get_str for for x==0 */
  74.           mpf_set_ui (x, 0L);
  75.           base = 10;
  76.         }
  77.       else
  78.         {
  79.           size = urandom () % (2 * SIZE) - SIZE;
  80.           exp = urandom () % EXPO;
  81.           mpf_random2 (x, size, exp);
  82.           base = urandom () % 61 + 2;
  83.         }
  84.       str = mpf_get_str (0, &bexp, base, 0, x);
  85.       if (str[0] == '-')
  86. sprintf (buf, "-0.%s@%ld", str + 1, bexp);
  87.       else
  88. sprintf (buf, "0.%s@%ld", str, bexp);
  89.       mpf_set_str_or_abort (y, buf, -base);
  90.       (*__gmp_free_func) (str, strlen (str) + 1);
  91.       mpf_reldiff (rerr, x, y);
  92.       if (mpf_cmp (rerr, max_rerr) > 0)
  93. {
  94.   mpf_set (max_rerr, rerr);
  95. #if VERBOSE
  96.   mpf_dump (max_rerr);
  97. #endif
  98.   if (mpf_cmp (rerr, limit_rerr) > 0)
  99.     {
  100.       printf ("ERROR after %d testsn", i);
  101.       printf ("base = %dn", base);
  102.       printf ("   x = "); mpf_dump (x);
  103.       printf ("   y = "); mpf_dump (y);
  104.       abort ();
  105.     }
  106. }
  107.     }
  108.   mpf_clear (limit_rerr);
  109.   mpf_clear (rerr);
  110.   mpf_clear (max_rerr);
  111.   mpf_clear (x);
  112.   mpf_clear (y);
  113.   mpf_clear (d);
  114.   tests_end ();
  115.   exit (0);
  116. }