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

数学计算

开发平台:

Unix_Linux

  1. /* Test sub_ddmmss.
  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 "longlong.h"
  19. #include "tests.h"
  20. void
  21. check_data (void)
  22. {
  23. #define M  MP_LIMB_T_MAX
  24.   static const struct {
  25.     mp_limb_t  want_dh,want_dl, mh,ml, sh,sl;
  26.   } data[] = {
  27.     { 0,0,  0,0,  0,0 },
  28.     { 0,0,  0,1,  0,1 },
  29.     { 0,0,  1,2,  1,2 },
  30.     { 0,1,  0,2,  0,1 },
  31.     { 0,M,  1,0,  0,1 },
  32.     { M,M,  0,0,  0,1 },
  33.     { M,M,  0,M-1,  0,M },
  34.     { 0,0,  0,M-1,  0,M-1 },
  35.     { 0,1,  0,M-1,  0,M-2 },
  36.   };
  37.   int  i;
  38.   mp_limb_t  got_dh, got_dl;
  39.   for (i = 0; i < numberof (data); i++)
  40.     {
  41.       sub_ddmmss (got_dh,got_dl, data[i].mh,data[i].ml, data[i].sh,data[i].sl);
  42.       if (got_dh != data[i].want_dh || got_dl != data[i].want_dl)
  43.         {
  44.           printf ("check_data wrong at data[%d]n", i);
  45.           mp_limb_trace ("  mh", data[i].mh);
  46.           mp_limb_trace ("  ml", data[i].ml);
  47.           mp_limb_trace ("  sh", data[i].sh);
  48.           mp_limb_trace ("  sl", data[i].sl);
  49.           mp_limb_trace ("  want dh", data[i].want_dh);
  50.           mp_limb_trace ("  want dl", data[i].want_dl);
  51.           mp_limb_trace ("  got dh ", got_dh);
  52.           mp_limb_trace ("  got dl ", got_dl);
  53.           abort ();
  54.         }
  55.     }
  56. }
  57. void
  58. check_random (void)
  59. {
  60.   mp_limb_t  want_dh,want_dl, got_dh,got_dl, mh,ml, sh,sl;
  61.   int  i;
  62.   for (i = 0; i < 20; i++)
  63.     {
  64.       mh = urandom ();
  65.       ml = urandom ();
  66.       sh = urandom ();
  67.       sl = urandom ();
  68.       refmpn_sub_ddmmss (&want_dh,&want_dl, mh,ml, sh,sl);
  69.       sub_ddmmss (got_dh,got_dl, mh,ml, sh,sl);
  70.       if (got_dh != want_dh || got_dl != want_dl)
  71.         {
  72.           printf ("check_data wrong at data[%d]n", i);
  73.           mp_limb_trace ("  mh", mh);
  74.           mp_limb_trace ("  ml", ml);
  75.           mp_limb_trace ("  sh", sh);
  76.           mp_limb_trace ("  sl", sl);
  77.           mp_limb_trace ("  want dh", want_dh);
  78.           mp_limb_trace ("  want dl", want_dl);
  79.           mp_limb_trace ("  got dh ", got_dh);
  80.           mp_limb_trace ("  got dl ", got_dl);
  81.           abort ();
  82.         }
  83.     }
  84. }
  85. int
  86. main (void)
  87. {
  88.   tests_start ();
  89.   mp_trace_base = -16;
  90.   check_data ();
  91.   check_random ();
  92.   tests_end ();
  93.   exit (0);
  94. }