Hex Artifact Content
Not logged in

Artifact c3d915f0988a33450d0eb692e2fc958f6397268c:


0000: 2f 2a 20 54 65 73 74 20 66 69 6c 65 20 66 6f 72  /* Test file for
0010: 20 6d 70 66 72 5f 63 6f 6e 73 74 5f 65 75 6c 65   mpfr_const_eule
0020: 72 2e 0a 0a 43 6f 70 79 72 69 67 68 74 20 32 30  r...Copyright 20
0030: 30 31 20 46 72 65 65 20 53 6f 66 74 77 61 72 65  01 Free Software
0040: 20 46 6f 75 6e 64 61 74 69 6f 6e 2e 0a 0a 54 68   Foundation...Th
0050: 69 73 20 66 69 6c 65 20 69 73 20 70 61 72 74 20  is file is part 
0060: 6f 66 20 74 68 65 20 4d 50 46 52 20 4c 69 62 72  of the MPFR Libr
0070: 61 72 79 2e 0a 0a 54 68 65 20 4d 50 46 52 20 4c  ary...The MPFR L
0080: 69 62 72 61 72 79 20 69 73 20 66 72 65 65 20 73  ibrary is free s
0090: 6f 66 74 77 61 72 65 3b 20 79 6f 75 20 63 61 6e  oftware; you can
00a0: 20 72 65 64 69 73 74 72 69 62 75 74 65 20 69 74   redistribute it
00b0: 20 61 6e 64 2f 6f 72 20 6d 6f 64 69 66 79 0a 69   and/or modify.i
00c0: 74 20 75 6e 64 65 72 20 74 68 65 20 74 65 72 6d  t under the term
00d0: 73 20 6f 66 20 74 68 65 20 47 4e 55 20 4c 65 73  s of the GNU Les
00e0: 73 65 72 20 47 65 6e 65 72 61 6c 20 50 75 62 6c  ser General Publ
00f0: 69 63 20 4c 69 63 65 6e 73 65 20 61 73 20 70 75  ic License as pu
0100: 62 6c 69 73 68 65 64 20 62 79 0a 74 68 65 20 46  blished by.the F
0110: 72 65 65 20 53 6f 66 74 77 61 72 65 20 46 6f 75  ree Software Fou
0120: 6e 64 61 74 69 6f 6e 3b 20 65 69 74 68 65 72 20  ndation; either 
0130: 76 65 72 73 69 6f 6e 20 32 2e 31 20 6f 66 20 74  version 2.1 of t
0140: 68 65 20 4c 69 63 65 6e 73 65 2c 20 6f 72 20 28  he License, or (
0150: 61 74 20 79 6f 75 72 0a 6f 70 74 69 6f 6e 29 20  at your.option) 
0160: 61 6e 79 20 6c 61 74 65 72 20 76 65 72 73 69 6f  any later versio
0170: 6e 2e 0a 0a 54 68 65 20 4d 50 46 52 20 4c 69 62  n...The MPFR Lib
0180: 72 61 72 79 20 69 73 20 64 69 73 74 72 69 62 75  rary is distribu
0190: 74 65 64 20 69 6e 20 74 68 65 20 68 6f 70 65 20  ted in the hope 
01a0: 74 68 61 74 20 69 74 20 77 69 6c 6c 20 62 65 20  that it will be 
01b0: 75 73 65 66 75 6c 2c 20 62 75 74 0a 57 49 54 48  useful, but.WITH
01c0: 4f 55 54 20 41 4e 59 20 57 41 52 52 41 4e 54 59  OUT ANY WARRANTY
01d0: 3b 20 77 69 74 68 6f 75 74 20 65 76 65 6e 20 74  ; without even t
01e0: 68 65 20 69 6d 70 6c 69 65 64 20 77 61 72 72 61  he implied warra
01f0: 6e 74 79 20 6f 66 20 4d 45 52 43 48 41 4e 54 41  nty of MERCHANTA
0200: 42 49 4c 49 54 59 0a 6f 72 20 46 49 54 4e 45 53  BILITY.or FITNES
0210: 53 20 46 4f 52 20 41 20 50 41 52 54 49 43 55 4c  S FOR A PARTICUL
0220: 41 52 20 50 55 52 50 4f 53 45 2e 20 20 53 65 65  AR PURPOSE.  See
0230: 20 74 68 65 20 47 4e 55 20 4c 65 73 73 65 72 20   the GNU Lesser 
0240: 47 65 6e 65 72 61 6c 20 50 75 62 6c 69 63 0a 4c  General Public.L
0250: 69 63 65 6e 73 65 20 66 6f 72 20 6d 6f 72 65 20  icense for more 
0260: 64 65 74 61 69 6c 73 2e 0a 0a 59 6f 75 20 73 68  details...You sh
0270: 6f 75 6c 64 20 68 61 76 65 20 72 65 63 65 69 76  ould have receiv
0280: 65 64 20 61 20 63 6f 70 79 20 6f 66 20 74 68 65  ed a copy of the
0290: 20 47 4e 55 20 4c 65 73 73 65 72 20 47 65 6e 65   GNU Lesser Gene
02a0: 72 61 6c 20 50 75 62 6c 69 63 20 4c 69 63 65 6e  ral Public Licen
02b0: 73 65 0a 61 6c 6f 6e 67 20 77 69 74 68 20 74 68  se.along with th
02c0: 65 20 4d 50 46 52 20 4c 69 62 72 61 72 79 3b 20  e MPFR Library; 
02d0: 73 65 65 20 74 68 65 20 66 69 6c 65 20 43 4f 50  see the file COP
02e0: 59 49 4e 47 2e 4c 49 42 2e 20 20 49 66 20 6e 6f  YING.LIB.  If no
02f0: 74 2c 20 77 72 69 74 65 20 74 6f 0a 74 68 65 20  t, write to.the 
0300: 46 72 65 65 20 53 6f 66 74 77 61 72 65 20 46 6f  Free Software Fo
0310: 75 6e 64 61 74 69 6f 6e 2c 20 49 6e 63 2e 2c 20  undation, Inc., 
0320: 35 39 20 54 65 6d 70 6c 65 20 50 6c 61 63 65 20  59 Temple Place 
0330: 2d 20 53 75 69 74 65 20 33 33 30 2c 20 42 6f 73  - Suite 330, Bos
0340: 74 6f 6e 2c 0a 4d 41 20 30 32 31 31 31 2d 31 33  ton,.MA 02111-13
0350: 30 37 2c 20 55 53 41 2e 20 2a 2f 0a 0a 23 69 6e  07, USA. */..#in
0360: 63 6c 75 64 65 20 3c 73 74 64 69 6f 2e 68 3e 0a  clude <stdio.h>.
0370: 23 69 6e 63 6c 75 64 65 20 3c 73 74 64 6c 69 62  #include <stdlib
0380: 2e 68 3e 0a 23 69 6e 63 6c 75 64 65 20 22 67 6d  .h>.#include "gm
0390: 70 2e 68 22 0a 23 69 6e 63 6c 75 64 65 20 22 6d  p.h".#include "m
03a0: 70 66 72 2e 68 22 0a 23 69 6e 63 6c 75 64 65 20  pfr.h".#include 
03b0: 22 6d 70 66 72 2d 69 6d 70 6c 2e 68 22 0a 0a 69  "mpfr-impl.h"..i
03c0: 6e 74 0a 6d 61 69 6e 20 28 69 6e 74 20 61 72 67  nt.main (int arg
03d0: 63 2c 20 63 68 61 72 20 2a 61 72 67 76 5b 5d 29  c, char *argv[])
03e0: 0a 7b 0a 20 20 6d 70 66 72 5f 74 20 67 61 6d 6d  .{.  mpfr_t gamm
03f0: 61 2c 20 79 2c 20 7a 2c 20 74 3b 0a 20 20 75 6e  a, y, z, t;.  un
0400: 73 69 67 6e 65 64 20 69 6e 74 20 65 72 72 2c 20  signed int err, 
0410: 70 72 65 63 2c 20 79 70 72 65 63 2c 20 70 30 20  prec, yprec, p0 
0420: 3d 20 32 2c 20 70 31 20 3d 20 32 30 30 3b 0a 20  = 2, p1 = 200;. 
0430: 20 6d 70 5f 72 6e 64 5f 74 20 72 6e 64 3b 0a 0a   mp_rnd_t rnd;..
0440: 20 20 70 72 65 63 20 3d 20 28 61 72 67 63 20 3c    prec = (argc <
0450: 20 32 29 20 3f 20 35 33 20 3a 20 61 74 6f 69 28   2) ? 53 : atoi(
0460: 61 72 67 76 5b 31 5d 29 3b 0a 0a 20 20 69 66 20  argv[1]);..  if 
0470: 28 61 72 67 63 20 3e 20 31 29 0a 20 20 20 20 7b  (argc > 1).    {
0480: 0a 20 20 20 20 20 20 6d 70 66 72 5f 69 6e 69 74  .      mpfr_init
0490: 32 20 28 67 61 6d 6d 61 2c 20 70 72 65 63 29 3b  2 (gamma, prec);
04a0: 0a 20 20 20 20 20 20 6d 70 66 72 5f 63 6f 6e 73  .      mpfr_cons
04b0: 74 5f 65 75 6c 65 72 20 28 67 61 6d 6d 61 2c 20  t_euler (gamma, 
04c0: 47 4d 50 5f 52 4e 44 4e 29 3b 0a 20 20 20 20 20  GMP_RNDN);.     
04d0: 20 70 72 69 6e 74 66 28 22 67 61 6d 6d 61 3d 22   printf("gamma="
04e0: 29 3b 20 6d 70 66 72 5f 6f 75 74 5f 73 74 72 20  ); mpfr_out_str 
04f0: 28 73 74 64 6f 75 74 2c 20 31 30 2c 20 30 2c 20  (stdout, 10, 0, 
0500: 67 61 6d 6d 61 2c 20 47 4d 50 5f 52 4e 44 44 29  gamma, GMP_RNDD)
0510: 3b 0a 20 20 20 20 20 20 70 75 74 63 68 61 72 20  ;.      putchar 
0520: 28 27 5c 6e 27 29 3b 0a 20 20 20 20 20 20 6d 70  ('\n');.      mp
0530: 66 72 5f 63 6c 65 61 72 20 28 67 61 6d 6d 61 29  fr_clear (gamma)
0540: 3b 0a 20 20 20 20 20 20 72 65 74 75 72 6e 20 30  ;.      return 0
0550: 3b 0a 20 20 20 20 7d 0a 0a 20 20 6d 70 66 72 5f  ;.    }..  mpfr_
0560: 69 6e 69 74 20 28 79 29 3b 0a 20 20 6d 70 66 72  init (y);.  mpfr
0570: 5f 69 6e 69 74 20 28 7a 29 3b 0a 20 20 6d 70 66  _init (z);.  mpf
0580: 72 5f 69 6e 69 74 20 28 74 29 3b 0a 0a 20 20 66  r_init (t);..  f
0590: 6f 72 20 28 70 72 65 63 20 3d 20 70 30 3b 20 70  or (prec = p0; p
05a0: 72 65 63 20 3c 3d 20 70 31 3b 20 70 72 65 63 2b  rec <= p1; prec+
05b0: 2b 29 0a 20 20 20 20 7b 0a 20 20 20 20 20 20 6d  +).    {.      m
05c0: 70 66 72 5f 73 65 74 5f 70 72 65 63 20 28 7a 2c  pfr_set_prec (z,
05d0: 20 70 72 65 63 29 3b 0a 20 20 20 20 20 20 6d 70   prec);.      mp
05e0: 66 72 5f 73 65 74 5f 70 72 65 63 20 28 74 2c 20  fr_set_prec (t, 
05f0: 70 72 65 63 29 3b 0a 20 20 20 20 20 20 79 70 72  prec);.      ypr
0600: 65 63 20 3d 20 70 72 65 63 20 2b 20 31 30 3b 0a  ec = prec + 10;.
0610: 0a 20 20 20 20 20 20 66 6f 72 20 28 72 6e 64 3d  .      for (rnd=
0620: 30 3b 20 72 6e 64 3c 34 3b 20 72 6e 64 2b 2b 29  0; rnd<4; rnd++)
0630: 0a 09 7b 0a 09 20 20 6d 70 66 72 5f 73 65 74 5f  ..{..  mpfr_set_
0640: 70 72 65 63 20 28 79 2c 20 79 70 72 65 63 29 3b  prec (y, yprec);
0650: 0a 09 20 20 6d 70 66 72 5f 63 6f 6e 73 74 5f 65  ..  mpfr_const_e
0660: 75 6c 65 72 20 28 79 2c 20 72 6e 64 29 3b 0a 09  uler (y, rnd);..
0670: 20 20 65 72 72 20 3d 20 28 72 6e 64 20 3d 3d 20    err = (rnd == 
0680: 47 4d 50 5f 52 4e 44 4e 29 20 3f 20 79 70 72 65  GMP_RNDN) ? ypre
0690: 63 20 2b 20 31 20 3a 20 79 70 72 65 63 3b 0a 09  c + 1 : yprec;..
06a0: 20 20 69 66 20 28 6d 70 66 72 5f 63 61 6e 5f 72    if (mpfr_can_r
06b0: 6f 75 6e 64 20 28 79 2c 20 65 72 72 2c 20 72 6e  ound (y, err, rn
06c0: 64 2c 20 72 6e 64 2c 20 70 72 65 63 29 29 0a 09  d, rnd, prec))..
06d0: 20 20 20 20 7b 0a 09 20 20 20 20 20 20 6d 70 66      {..      mpf
06e0: 72 5f 73 65 74 20 28 74 2c 20 79 2c 20 72 6e 64  r_set (t, y, rnd
06f0: 29 3b 0a 09 20 20 20 20 20 20 6d 70 66 72 5f 63  );..      mpfr_c
0700: 6f 6e 73 74 5f 65 75 6c 65 72 20 28 7a 2c 20 72  onst_euler (z, r
0710: 6e 64 29 3b 0a 09 20 20 20 20 20 20 69 66 20 28  nd);..      if (
0720: 6d 70 66 72 5f 63 6d 70 20 28 74 2c 20 7a 29 29  mpfr_cmp (t, z))
0730: 0a 09 09 7b 0a 09 09 20 20 70 72 69 6e 74 66 20  ...{...  printf 
0740: 28 22 72 65 73 75 6c 74 73 20 64 69 66 66 65 72  ("results differ
0750: 20 66 6f 72 20 70 72 65 63 3d 25 75 20 72 6e 64   for prec=%u rnd
0760: 5f 6d 6f 64 65 3d 25 73 5c 6e 22 2c 20 70 72 65  _mode=%s\n", pre
0770: 63 2c 0a 09 09 09 20 20 6d 70 66 72 5f 70 72 69  c,....  mpfr_pri
0780: 6e 74 5f 72 6e 64 5f 6d 6f 64 65 20 28 72 6e 64  nt_rnd_mode (rnd
0790: 29 29 3b 0a 09 09 20 20 70 72 69 6e 74 66 20 28  ));...  printf (
07a0: 22 20 20 20 67 6f 74 20 20 20 20 20 20 22 29 3b  "   got      ");
07b0: 0a 09 09 20 20 6d 70 66 72 5f 6f 75 74 5f 73 74  ...  mpfr_out_st
07c0: 72 20 28 73 74 64 6f 75 74 2c 20 32 2c 20 70 72  r (stdout, 2, pr
07d0: 65 63 2c 20 7a 2c 20 47 4d 50 5f 52 4e 44 4e 29  ec, z, GMP_RNDN)
07e0: 3b 0a 09 09 20 20 70 75 74 63 68 61 72 20 28 27  ;...  putchar ('
07f0: 5c 6e 27 29 3b 0a 09 09 20 20 70 72 69 6e 74 66  \n');...  printf
0800: 20 28 22 20 20 20 65 78 70 65 63 74 65 64 20 22   ("   expected "
0810: 29 3b 0a 09 09 20 20 6d 70 66 72 5f 6f 75 74 5f  );...  mpfr_out_
0820: 73 74 72 20 28 73 74 64 6f 75 74 2c 20 32 2c 20  str (stdout, 2, 
0830: 70 72 65 63 2c 20 74 2c 20 47 4d 50 5f 52 4e 44  prec, t, GMP_RND
0840: 4e 29 3b 0a 09 09 20 20 70 75 74 63 68 61 72 20  N);...  putchar 
0850: 28 27 5c 6e 27 29 3b 0a 09 09 20 20 70 72 69 6e  ('\n');...  prin
0860: 74 66 20 28 22 20 20 20 61 70 70 72 6f 78 69 6d  tf ("   approxim
0870: 61 74 69 6f 6e 20 77 61 73 20 22 29 3b 0a 09 09  ation was ");...
0880: 20 20 6d 70 66 72 5f 70 72 69 6e 74 5f 62 69 6e    mpfr_print_bin
0890: 61 72 79 20 28 79 29 3b 0a 09 09 20 20 70 75 74  ary (y);...  put
08a0: 63 68 61 72 20 28 27 5c 6e 27 29 3b 0a 09 09 20  char ('\n');... 
08b0: 20 65 78 69 74 20 28 31 29 3b 0a 09 09 7d 0a 09   exit (1);...}..
08c0: 20 20 20 20 7d 0a 09 7d 0a 20 20 20 20 7d 0a 0a      }..}.    }..
08d0: 20 20 6d 70 66 72 5f 63 6c 65 61 72 20 28 79 29    mpfr_clear (y)
08e0: 3b 0a 20 20 6d 70 66 72 5f 63 6c 65 61 72 20 28  ;.  mpfr_clear (
08f0: 7a 29 3b 0a 20 20 6d 70 66 72 5f 63 6c 65 61 72  z);.  mpfr_clear
0900: 20 28 74 29 3b 0a 0a 20 20 72 65 74 75 72 6e 20   (t);..  return 
0910: 30 3b 0a 7d 0a                                   0;.}.