Hex Artifact Content
Not logged in

Artifact d0581dee355d64b79198ddf81ee0c1bf9c5e6d90:


0000: 2f 2a 20 54 65 73 74 20 6d 70 7a 5f 67 65 74 5f  /* Test mpz_get_
0010: 64 2e 0a 0a 43 6f 70 79 72 69 67 68 74 20 32 30  d...Copyright 20
0020: 30 32 20 46 72 65 65 20 53 6f 66 74 77 61 72 65  02 Free Software
0030: 20 46 6f 75 6e 64 61 74 69 6f 6e 2c 20 49 6e 63   Foundation, Inc
0040: 2e 0a 0a 54 68 69 73 20 66 69 6c 65 20 69 73 20  ...This file is 
0050: 70 61 72 74 20 6f 66 20 74 68 65 20 47 4e 55 20  part of the GNU 
0060: 4d 50 20 4c 69 62 72 61 72 79 2e 0a 0a 54 68 65  MP Library...The
0070: 20 47 4e 55 20 4d 50 20 4c 69 62 72 61 72 79 20   GNU MP Library 
0080: 69 73 20 66 72 65 65 20 73 6f 66 74 77 61 72 65  is free software
0090: 3b 20 79 6f 75 20 63 61 6e 20 72 65 64 69 73 74  ; you can redist
00a0: 72 69 62 75 74 65 20 69 74 20 61 6e 64 2f 6f 72  ribute it and/or
00b0: 20 6d 6f 64 69 66 79 0a 69 74 20 75 6e 64 65 72   modify.it under
00c0: 20 74 68 65 20 74 65 72 6d 73 20 6f 66 20 74 68   the terms of th
00d0: 65 20 47 4e 55 20 4c 65 73 73 65 72 20 47 65 6e  e GNU Lesser Gen
00e0: 65 72 61 6c 20 50 75 62 6c 69 63 20 4c 69 63 65  eral Public Lice
00f0: 6e 73 65 20 61 73 20 70 75 62 6c 69 73 68 65 64  nse as published
0100: 20 62 79 0a 74 68 65 20 46 72 65 65 20 53 6f 66   by.the Free Sof
0110: 74 77 61 72 65 20 46 6f 75 6e 64 61 74 69 6f 6e  tware Foundation
0120: 3b 20 65 69 74 68 65 72 20 76 65 72 73 69 6f 6e  ; either version
0130: 20 32 2e 31 20 6f 66 20 74 68 65 20 4c 69 63 65   2.1 of the Lice
0140: 6e 73 65 2c 20 6f 72 20 28 61 74 20 79 6f 75 72  nse, or (at your
0150: 0a 6f 70 74 69 6f 6e 29 20 61 6e 79 20 6c 61 74  .option) any lat
0160: 65 72 20 76 65 72 73 69 6f 6e 2e 0a 0a 54 68 65  er version...The
0170: 20 47 4e 55 20 4d 50 20 4c 69 62 72 61 72 79 20   GNU MP Library 
0180: 69 73 20 64 69 73 74 72 69 62 75 74 65 64 20 69  is distributed i
0190: 6e 20 74 68 65 20 68 6f 70 65 20 74 68 61 74 20  n the hope that 
01a0: 69 74 20 77 69 6c 6c 20 62 65 20 75 73 65 66 75  it will be usefu
01b0: 6c 2c 20 62 75 74 0a 57 49 54 48 4f 55 54 20 41  l, but.WITHOUT A
01c0: 4e 59 20 57 41 52 52 41 4e 54 59 3b 20 77 69 74  NY WARRANTY; wit
01d0: 68 6f 75 74 20 65 76 65 6e 20 74 68 65 20 69 6d  hout even the im
01e0: 70 6c 69 65 64 20 77 61 72 72 61 6e 74 79 20 6f  plied warranty o
01f0: 66 20 4d 45 52 43 48 41 4e 54 41 42 49 4c 49 54  f MERCHANTABILIT
0200: 59 0a 6f 72 20 46 49 54 4e 45 53 53 20 46 4f 52  Y.or FITNESS FOR
0210: 20 41 20 50 41 52 54 49 43 55 4c 41 52 20 50 55   A PARTICULAR PU
0220: 52 50 4f 53 45 2e 20 20 53 65 65 20 74 68 65 20  RPOSE.  See the 
0230: 47 4e 55 20 4c 65 73 73 65 72 20 47 65 6e 65 72  GNU Lesser Gener
0240: 61 6c 20 50 75 62 6c 69 63 0a 4c 69 63 65 6e 73  al Public.Licens
0250: 65 20 66 6f 72 20 6d 6f 72 65 20 64 65 74 61 69  e for more detai
0260: 6c 73 2e 0a 0a 59 6f 75 20 73 68 6f 75 6c 64 20  ls...You should 
0270: 68 61 76 65 20 72 65 63 65 69 76 65 64 20 61 20  have received a 
0280: 63 6f 70 79 20 6f 66 20 74 68 65 20 47 4e 55 20  copy of the GNU 
0290: 4c 65 73 73 65 72 20 47 65 6e 65 72 61 6c 20 50  Lesser General P
02a0: 75 62 6c 69 63 20 4c 69 63 65 6e 73 65 0a 61 6c  ublic License.al
02b0: 6f 6e 67 20 77 69 74 68 20 74 68 65 20 47 4e 55  ong with the GNU
02c0: 20 4d 50 20 4c 69 62 72 61 72 79 3b 20 73 65 65   MP Library; see
02d0: 20 74 68 65 20 66 69 6c 65 20 43 4f 50 59 49 4e   the file COPYIN
02e0: 47 2e 4c 49 42 2e 20 20 49 66 20 6e 6f 74 2c 20  G.LIB.  If not, 
02f0: 77 72 69 74 65 20 74 6f 0a 74 68 65 20 46 72 65  write to.the Fre
0300: 65 20 53 6f 66 74 77 61 72 65 20 46 6f 75 6e 64  e Software Found
0310: 61 74 69 6f 6e 2c 20 49 6e 63 2e 2c 20 35 39 20  ation, Inc., 59 
0320: 54 65 6d 70 6c 65 20 50 6c 61 63 65 20 2d 20 53  Temple Place - S
0330: 75 69 74 65 20 33 33 30 2c 20 42 6f 73 74 6f 6e  uite 330, Boston
0340: 2c 0a 4d 41 20 30 32 31 31 31 2d 31 33 30 37 2c  ,.MA 02111-1307,
0350: 20 55 53 41 2e 20 2a 2f 0a 0a 23 69 6e 63 6c 75   USA. */..#inclu
0360: 64 65 20 3c 73 74 64 69 6f 2e 68 3e 0a 23 69 6e  de <stdio.h>.#in
0370: 63 6c 75 64 65 20 3c 73 74 64 6c 69 62 2e 68 3e  clude <stdlib.h>
0380: 0a 23 69 6e 63 6c 75 64 65 20 22 67 6d 70 2e 68  .#include "gmp.h
0390: 22 0a 23 69 6e 63 6c 75 64 65 20 22 67 6d 70 2d  ".#include "gmp-
03a0: 69 6d 70 6c 2e 68 22 0a 23 69 6e 63 6c 75 64 65  impl.h".#include
03b0: 20 22 74 65 73 74 73 2e 68 22 0a 0a 0a 76 6f 69   "tests.h"...voi
03c0: 64 0a 63 68 65 63 6b 5f 6f 6e 65 62 69 74 20 28  d.check_onebit (
03d0: 76 6f 69 64 29 0a 7b 0a 20 20 69 6e 74 20 20 20  void).{.  int   
03e0: 20 20 69 3b 0a 20 20 6d 70 7a 5f 74 20 20 20 7a    i;.  mpz_t   z
03f0: 3b 0a 20 20 64 6f 75 62 6c 65 20 20 67 6f 74 2c  ;.  double  got,
0400: 20 77 61 6e 74 3b 0a 20 20 2f 2a 20 46 49 58 4d   want;.  /* FIXM
0410: 45 3a 20 49 74 27 64 20 62 65 20 62 65 74 74 65  E: It'd be bette
0420: 72 20 74 6f 20 62 61 73 65 20 74 68 69 73 20 6f  r to base this o
0430: 6e 20 74 68 65 20 66 6c 6f 61 74 20 66 6f 72 6d  n the float form
0440: 61 74 2e 20 2a 2f 0a 23 69 66 64 65 66 20 5f 5f  at. */.#ifdef __
0450: 76 61 78 0a 20 20 69 6e 74 20 20 20 20 20 6c 69  vax.  int     li
0460: 6d 69 74 20 3d 20 31 32 37 3b 20 20 2f 2a 20 76  mit = 127;  /* v
0470: 61 78 20 66 70 20 6e 75 6d 62 65 72 73 20 68 61  ax fp numbers ha
0480: 76 65 20 6c 69 6d 69 74 65 64 20 72 61 6e 67 65  ve limited range
0490: 20 2a 2f 0a 23 65 6c 73 65 0a 20 20 69 6e 74 20   */.#else.  int 
04a0: 20 20 20 20 6c 69 6d 69 74 20 3d 20 35 31 32 3b      limit = 512;
04b0: 0a 23 65 6e 64 69 66 0a 0a 20 20 6d 70 7a 5f 69  .#endif..  mpz_i
04c0: 6e 69 74 20 28 7a 29 3b 0a 0a 20 20 6d 70 7a 5f  nit (z);..  mpz_
04d0: 73 65 74 5f 75 69 20 28 7a 2c 20 31 4c 29 3b 0a  set_ui (z, 1L);.
04e0: 20 20 77 61 6e 74 20 3d 20 31 2e 30 3b 0a 0a 20    want = 1.0;.. 
04f0: 20 66 6f 72 20 28 69 20 3d 20 30 3b 20 69 20 3c   for (i = 0; i <
0500: 20 6c 69 6d 69 74 3b 20 69 2b 2b 29 0a 20 20 20   limit; i++).   
0510: 20 7b 0a 20 20 20 20 20 20 67 6f 74 20 3d 20 6d   {.      got = m
0520: 70 7a 5f 67 65 74 5f 64 20 28 7a 29 3b 0a 0a 20  pz_get_d (z);.. 
0530: 20 20 20 20 20 69 66 20 28 67 6f 74 20 21 3d 20       if (got != 
0540: 77 61 6e 74 29 0a 20 20 20 20 20 20 20 20 7b 0a  want).        {.
0550: 20 20 20 20 20 20 20 20 20 20 70 72 69 6e 74 66            printf
0560: 20 20 20 20 28 22 6d 70 7a 5f 67 65 74 5f 64 20      ("mpz_get_d 
0570: 77 72 6f 6e 67 20 6f 6e 20 32 2a 2a 25 64 5c 6e  wrong on 2**%d\n
0580: 22 2c 20 69 29 3b 0a 20 20 20 20 20 20 20 20 20  ", i);.         
0590: 20 6d 70 7a 5f 74 72 61 63 65 20 28 22 20 20 20   mpz_trace ("   
05a0: 7a 20 20 20 20 22 2c 20 7a 29 3b 0a 20 20 20 20  z    ", z);.    
05b0: 20 20 20 20 20 20 70 72 69 6e 74 66 20 20 20 20        printf    
05c0: 28 22 20 20 20 77 61 6e 74 20 20 25 2e 32 30 67  ("   want  %.20g
05d0: 5c 6e 22 2c 20 77 61 6e 74 29 3b 0a 20 20 20 20  \n", want);.    
05e0: 20 20 20 20 20 20 70 72 69 6e 74 66 20 20 20 20        printf    
05f0: 28 22 20 20 20 67 6f 74 20 20 20 25 2e 32 30 67  ("   got   %.20g
0600: 5c 6e 22 2c 20 67 6f 74 29 3b 0a 20 20 20 20 20  \n", got);.     
0610: 20 20 20 20 20 61 62 6f 72 74 28 29 3b 0a 20 20       abort();.  
0620: 20 20 20 20 20 20 7d 0a 0a 20 20 20 20 20 20 6d        }..      m
0630: 70 7a 5f 6d 75 6c 5f 32 65 78 70 20 28 7a 2c 20  pz_mul_2exp (z, 
0640: 7a 2c 20 31 4c 29 3b 0a 20 20 20 20 20 20 77 61  z, 1L);.      wa
0650: 6e 74 20 2a 3d 20 32 2e 30 3b 0a 20 20 20 20 7d  nt *= 2.0;.    }
0660: 0a 20 20 6d 70 7a 5f 63 6c 65 61 72 20 28 7a 29  .  mpz_clear (z)
0670: 3b 0a 7d 0a 0a 0a 69 6e 74 0a 6d 61 69 6e 20 28  ;.}...int.main (
0680: 76 6f 69 64 29 0a 7b 0a 20 20 74 65 73 74 73 5f  void).{.  tests_
0690: 73 74 61 72 74 20 28 29 3b 0a 0a 20 20 63 68 65  start ();..  che
06a0: 63 6b 5f 6f 6e 65 62 69 74 20 28 29 3b 0a 0a 20  ck_onebit ();.. 
06b0: 20 74 65 73 74 73 5f 65 6e 64 20 28 29 3b 0a 20   tests_end ();. 
06c0: 20 65 78 69 74 20 28 30 29 3b 0a 7d 0a            exit (0);.}.