Category

Development Platform
Home > SourceCode/Document > Mathimatics-Numerical algorithms > Algorithm > View Code
t-dm2exp.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-dm2exp.c：Code Content
  1. /* Test mpf_div, mpf_div_2exp, mpf_mul_2exp.
  3. Copyright 1996, 2000, 2001 Free Software Foundation, Inc.
  5. This file is part of the GNU MP Library.
  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.
  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.
  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/.  */
  20. #include <stdio.h>
  21. #include <stdlib.h>
  23. #include "gmp.h"
  24. #include "gmp-impl.h"
  25. #include "tests.h"
  27. #ifndef SIZE
  28. #define SIZE 16
  29. #endif
  31. int
  32. main (int argc, char **argv)
  33. {
  34.   int reps = 100000;
  35.   int i;
  36.   mpf_t u, v, w1, w2, w3;
  37.   mp_size_t bprec = 100;
  38.   mpf_t rerr, limit_rerr;
  39.   mp_size_t un;
  40.   mp_exp_t ue;
  42.   tests_start ();
  44.   if (argc > 1)
  45.     {
  46.       reps = strtol (argv[1], 0, 0);
  47.       if (argc > 2)
  48. bprec = strtol (argv[2], 0, 0);
  49.     }
  51.   mpf_set_default_prec (bprec);
  53.   mpf_init (rerr);
  54.   mpf_init (limit_rerr);
  56.   mpf_init (u);
  57.   mpf_init (v);
  58.   mpf_init (w1);
  59.   mpf_init (w2);
  60.   mpf_init (w3);
  62.   for (i = 0; i < reps; i++)
  63.     {
  64.       unsigned long int res_prec;
  65.       unsigned long int pow2;
  67.       res_prec = urandom () % (bprec + 100);
  68.       mpf_set_prec (w1, res_prec);
  69.       mpf_set_prec (w2, res_prec);
  70.       mpf_set_prec (w3, res_prec);
  72.       mpf_set_ui (limit_rerr, 1);
  73.       mpf_div_2exp (limit_rerr, limit_rerr, res_prec);
  75.       pow2 = urandom () % 0x10000;
  76.       mpf_set_ui (v, 1);
  77.       mpf_mul_2exp (v, v, pow2);
  79.       un = urandom () % (2 * SIZE) - SIZE;
  80.       ue = urandom () % SIZE;
  81.       mpf_random2 (u, un, ue);
  83.       mpf_div_2exp (w1, u, pow2);
  84.       mpf_div (w2, u, v);
  85.       mpf_reldiff (rerr, w1, w2);
  86.       if (mpf_cmp (rerr, limit_rerr) > 0)
  87. {
  88.   printf ("ERROR in mpf_div or mpf_div_2exp after %d testsn", i);
  89.   printf ("   u = "); mpf_dump (u);
  90.   printf ("   v = "); mpf_dump (v);
  91.   printf ("  w1 = "); mpf_dump (w1);
  92.   printf ("  w2 = "); mpf_dump (w2);
  93.   abort ();
  94. }
  95.       mpf_mul_2exp (w3, w1, pow2);
  96.       mpf_reldiff (rerr, u, w3);
  97.       if (mpf_cmp (rerr, limit_rerr) > 0)
  98. {
  99.   printf ("ERROR in mpf_mul_2exp after %d testsn", i);
  100.   printf ("   u = "); mpf_dump (u);
  101.   printf ("   v = "); mpf_dump (v);
  102.   printf ("  w1 = "); mpf_dump (w1);
  103.   printf ("  w3 = "); mpf_dump (w3);
  104.   abort ();
  105. }
  106.     }
  108.   mpf_clear (rerr);
  109.   mpf_clear (limit_rerr);
  111.   mpf_clear (u);
  112.   mpf_clear (v);
  113.   mpf_clear (w1);
  114.   mpf_clear (w2);
  115.   mpf_clear (w3);
  117.   tests_end ();
  118.   exit (0);
  119. }