Hex Artifact Content
Not logged in

Artifact 103e94153710e31d89042aa1d4ecf31de2af2f0a:


0000: 2f 2a 20 6d 70 66 5f 69 6e 70 5f 73 74 72 28 64  /* mpf_inp_str(d
0010: 65 73 74 5f 66 6c 6f 61 74 2c 20 73 74 72 65 61  est_float, strea
0020: 6d 2c 20 62 61 73 65 29 20 2d 2d 20 49 6e 70 75  m, base) -- Inpu
0030: 74 20 61 20 6e 75 6d 62 65 72 20 69 6e 20 62 61  t a number in ba
0040: 73 65 0a 20 20 20 42 41 53 45 20 66 72 6f 6d 20  se.   BASE from 
0050: 73 74 64 69 6f 20 73 74 72 65 61 6d 20 53 54 52  stdio stream STR
0060: 45 41 4d 20 61 6e 64 20 73 74 6f 72 65 20 74 68  EAM and store th
0070: 65 20 72 65 73 75 6c 74 20 69 6e 20 44 45 53 54  e result in DEST
0080: 5f 46 4c 4f 41 54 2e 0a 0a 43 6f 70 79 72 69 67  _FLOAT...Copyrig
0090: 68 74 20 31 39 39 39 2c 20 32 30 30 31 20 46 72  ht 1999, 2001 Fr
00a0: 65 65 20 53 6f 66 74 77 61 72 65 20 46 6f 75 6e  ee Software Foun
00b0: 64 61 74 69 6f 6e 2c 20 49 6e 63 2e 0a 28 43 6f  dation, Inc..(Co
00c0: 70 69 65 64 20 66 72 6f 6d 20 47 4d 50 2c 20 66  pied from GMP, f
00d0: 69 6c 65 20 6d 70 66 2f 69 6e 70 5f 73 74 72 2e  ile mpf/inp_str.
00e0: 63 29 0a 0a 54 68 69 73 20 66 69 6c 65 20 69 73  c)..This file is
00f0: 20 70 61 72 74 20 6f 66 20 74 68 65 20 4d 50 46   part of the MPF
0100: 52 20 4c 69 62 72 61 72 79 2e 0a 0a 54 68 65 20  R Library...The 
0110: 4d 50 46 52 20 4c 69 62 72 61 72 79 20 69 73 20  MPFR Library is 
0120: 66 72 65 65 20 73 6f 66 74 77 61 72 65 3b 20 79  free software; y
0130: 6f 75 20 63 61 6e 20 72 65 64 69 73 74 72 69 62  ou can redistrib
0140: 75 74 65 20 69 74 20 61 6e 64 2f 6f 72 20 6d 6f  ute it and/or mo
0150: 64 69 66 79 0a 69 74 20 75 6e 64 65 72 20 74 68  dify.it under th
0160: 65 20 74 65 72 6d 73 20 6f 66 20 74 68 65 20 47  e terms of the G
0170: 4e 55 20 4c 65 73 73 65 72 20 47 65 6e 65 72 61  NU Lesser Genera
0180: 6c 20 50 75 62 6c 69 63 20 4c 69 63 65 6e 73 65  l Public License
0190: 20 61 73 20 70 75 62 6c 69 73 68 65 64 20 62 79   as published by
01a0: 0a 74 68 65 20 46 72 65 65 20 53 6f 66 74 77 61  .the Free Softwa
01b0: 72 65 20 46 6f 75 6e 64 61 74 69 6f 6e 3b 20 65  re Foundation; e
01c0: 69 74 68 65 72 20 76 65 72 73 69 6f 6e 20 32 2e  ither version 2.
01d0: 31 20 6f 66 20 74 68 65 20 4c 69 63 65 6e 73 65  1 of the License
01e0: 2c 20 6f 72 20 28 61 74 20 79 6f 75 72 0a 6f 70  , or (at your.op
01f0: 74 69 6f 6e 29 20 61 6e 79 20 6c 61 74 65 72 20  tion) any later 
0200: 76 65 72 73 69 6f 6e 2e 0a 0a 54 68 65 20 4d 50  version...The MP
0210: 46 52 20 4c 69 62 72 61 72 79 20 69 73 20 64 69  FR Library is di
0220: 73 74 72 69 62 75 74 65 64 20 69 6e 20 74 68 65  stributed in the
0230: 20 68 6f 70 65 20 74 68 61 74 20 69 74 20 77 69   hope that it wi
0240: 6c 6c 20 62 65 20 75 73 65 66 75 6c 2c 20 62 75  ll be useful, bu
0250: 74 0a 57 49 54 48 4f 55 54 20 41 4e 59 20 57 41  t.WITHOUT ANY WA
0260: 52 52 41 4e 54 59 3b 20 77 69 74 68 6f 75 74 20  RRANTY; without 
0270: 65 76 65 6e 20 74 68 65 20 69 6d 70 6c 69 65 64  even the implied
0280: 20 77 61 72 72 61 6e 74 79 20 6f 66 20 4d 45 52   warranty of MER
0290: 43 48 41 4e 54 41 42 49 4c 49 54 59 0a 6f 72 20  CHANTABILITY.or 
02a0: 46 49 54 4e 45 53 53 20 46 4f 52 20 41 20 50 41  FITNESS FOR A PA
02b0: 52 54 49 43 55 4c 41 52 20 50 55 52 50 4f 53 45  RTICULAR PURPOSE
02c0: 2e 20 20 53 65 65 20 74 68 65 20 47 4e 55 20 4c  .  See the GNU L
02d0: 65 73 73 65 72 20 47 65 6e 65 72 61 6c 20 50 75  esser General Pu
02e0: 62 6c 69 63 0a 4c 69 63 65 6e 73 65 20 66 6f 72  blic.License for
02f0: 20 6d 6f 72 65 20 64 65 74 61 69 6c 73 2e 0a 0a   more details...
0300: 59 6f 75 20 73 68 6f 75 6c 64 20 68 61 76 65 20  You should have 
0310: 72 65 63 65 69 76 65 64 20 61 20 63 6f 70 79 20  received a copy 
0320: 6f 66 20 74 68 65 20 47 4e 55 20 4c 65 73 73 65  of the GNU Lesse
0330: 72 20 47 65 6e 65 72 61 6c 20 50 75 62 6c 69 63  r General Public
0340: 20 4c 69 63 65 6e 73 65 0a 61 6c 6f 6e 67 20 77   License.along w
0350: 69 74 68 20 74 68 65 20 4d 50 46 52 20 4c 69 62  ith the MPFR Lib
0360: 72 61 72 79 3b 20 73 65 65 20 74 68 65 20 66 69  rary; see the fi
0370: 6c 65 20 43 4f 50 59 49 4e 47 2e 4c 49 42 2e 20  le COPYING.LIB. 
0380: 20 49 66 20 6e 6f 74 2c 20 77 72 69 74 65 20 74   If not, write t
0390: 6f 0a 74 68 65 20 46 72 65 65 20 53 6f 66 74 77  o.the Free Softw
03a0: 61 72 65 20 46 6f 75 6e 64 61 74 69 6f 6e 2c 20  are Foundation, 
03b0: 49 6e 63 2e 2c 20 35 39 20 54 65 6d 70 6c 65 20  Inc., 59 Temple 
03c0: 50 6c 61 63 65 20 2d 20 53 75 69 74 65 20 33 33  Place - Suite 33
03d0: 30 2c 20 42 6f 73 74 6f 6e 2c 0a 4d 41 20 30 32  0, Boston,.MA 02
03e0: 31 31 31 2d 31 33 30 37 2c 20 55 53 41 2e 20 2a  111-1307, USA. *
03f0: 2f 0a 0a 23 69 6e 63 6c 75 64 65 20 3c 73 74 64  /..#include <std
0400: 69 6f 2e 68 3e 0a 23 69 6e 63 6c 75 64 65 20 3c  io.h>.#include <
0410: 63 74 79 70 65 2e 68 3e 0a 23 69 6e 63 6c 75 64  ctype.h>.#includ
0420: 65 20 22 67 6d 70 2e 68 22 0a 23 69 6e 63 6c 75  e "gmp.h".#inclu
0430: 64 65 20 22 67 6d 70 2d 69 6d 70 6c 2e 68 22 0a  de "gmp-impl.h".
0440: 23 69 6e 63 6c 75 64 65 20 22 6d 70 66 72 2e 68  #include "mpfr.h
0450: 22 0a 23 69 6e 63 6c 75 64 65 20 22 6d 70 66 72  ".#include "mpfr
0460: 2d 69 6d 70 6c 2e 68 22 0a 0a 73 69 7a 65 5f 74  -impl.h"..size_t
0470: 0a 6d 70 66 72 5f 69 6e 70 5f 73 74 72 20 28 6d  .mpfr_inp_str (m
0480: 70 66 72 5f 70 74 72 20 72 6f 70 2c 20 46 49 4c  pfr_ptr rop, FIL
0490: 45 20 2a 73 74 72 65 61 6d 2c 20 69 6e 74 20 62  E *stream, int b
04a0: 61 73 65 2c 20 6d 70 5f 72 6e 64 5f 74 20 72 6e  ase, mp_rnd_t rn
04b0: 64 5f 6d 6f 64 65 29 0a 7b 0a 20 20 63 68 61 72  d_mode).{.  char
04c0: 20 2a 73 74 72 3b 0a 20 20 73 69 7a 65 5f 74 20   *str;.  size_t 
04d0: 61 6c 6c 6f 63 5f 73 69 7a 65 2c 20 73 74 72 5f  alloc_size, str_
04e0: 73 69 7a 65 3b 0a 20 20 69 6e 74 20 63 3b 0a 20  size;.  int c;. 
04f0: 20 73 69 7a 65 5f 74 20 72 65 74 76 61 6c 3b 0a   size_t retval;.
0500: 20 20 73 69 7a 65 5f 74 20 6e 72 65 61 64 3b 0a    size_t nread;.
0510: 0a 20 20 4d 50 46 52 5f 43 4c 45 41 52 5f 46 4c  .  MPFR_CLEAR_FL
0520: 41 47 53 28 72 6f 70 29 3b 20 0a 20 20 69 66 20  AGS(rop); .  if 
0530: 28 73 74 72 65 61 6d 20 3d 3d 20 30 29 0a 20 20  (stream == 0).  
0540: 20 20 73 74 72 65 61 6d 20 3d 20 73 74 64 69 6e    stream = stdin
0550: 3b 0a 0a 20 20 61 6c 6c 6f 63 5f 73 69 7a 65 20  ;..  alloc_size 
0560: 3d 20 31 30 30 3b 0a 20 20 73 74 72 20 3d 20 28  = 100;.  str = (
0570: 63 68 61 72 20 2a 29 20 28 2a 5f 5f 67 6d 70 5f  char *) (*__gmp_
0580: 61 6c 6c 6f 63 61 74 65 5f 66 75 6e 63 29 20 28  allocate_func) (
0590: 61 6c 6c 6f 63 5f 73 69 7a 65 29 3b 0a 20 20 73  alloc_size);.  s
05a0: 74 72 5f 73 69 7a 65 20 3d 20 30 3b 0a 20 20 6e  tr_size = 0;.  n
05b0: 72 65 61 64 20 3d 20 30 3b 0a 0a 20 20 2f 2a 20  read = 0;..  /* 
05c0: 53 6b 69 70 20 77 68 69 74 65 73 70 61 63 65 2e  Skip whitespace.
05d0: 20 20 2a 2f 0a 20 20 64 6f 0a 20 20 20 20 7b 0a    */.  do.    {.
05e0: 20 20 20 20 20 20 63 20 3d 20 67 65 74 63 20 28        c = getc (
05f0: 73 74 72 65 61 6d 29 3b 0a 20 20 20 20 20 20 6e  stream);.      n
0600: 72 65 61 64 2b 2b 3b 0a 20 20 20 20 7d 0a 20 20  read++;.    }.  
0610: 77 68 69 6c 65 20 28 69 73 73 70 61 63 65 20 28  while (isspace (
0620: 63 29 29 3b 0a 0a 20 20 66 6f 72 20 28 3b 3b 29  c));..  for (;;)
0630: 0a 20 20 20 20 7b 0a 20 20 20 20 20 20 69 66 20  .    {.      if 
0640: 28 73 74 72 5f 73 69 7a 65 20 3e 3d 20 61 6c 6c  (str_size >= all
0650: 6f 63 5f 73 69 7a 65 29 0a 09 7b 0a 09 20 20 73  oc_size)..{..  s
0660: 69 7a 65 5f 74 20 6f 6c 64 5f 61 6c 6c 6f 63 5f  ize_t old_alloc_
0670: 73 69 7a 65 20 3d 20 61 6c 6c 6f 63 5f 73 69 7a  size = alloc_siz
0680: 65 3b 0a 09 20 20 61 6c 6c 6f 63 5f 73 69 7a 65  e;..  alloc_size
0690: 20 3d 20 61 6c 6c 6f 63 5f 73 69 7a 65 20 2a 20   = alloc_size * 
06a0: 33 20 2f 20 32 3b 0a 09 20 20 73 74 72 20 3d 20  3 / 2;..  str = 
06b0: 28 63 68 61 72 20 2a 29 20 28 2a 5f 5f 67 6d 70  (char *) (*__gmp
06c0: 5f 72 65 61 6c 6c 6f 63 61 74 65 5f 66 75 6e 63  _reallocate_func
06d0: 29 20 28 73 74 72 2c 20 6f 6c 64 5f 61 6c 6c 6f  ) (str, old_allo
06e0: 63 5f 73 69 7a 65 2c 20 61 6c 6c 6f 63 5f 73 69  c_size, alloc_si
06f0: 7a 65 29 3b 0a 09 7d 0a 20 20 20 20 20 20 69 66  ze);..}.      if
0700: 20 28 63 20 3d 3d 20 45 4f 46 20 7c 7c 20 69 73   (c == EOF || is
0710: 73 70 61 63 65 20 28 63 29 29 0a 09 62 72 65 61  space (c))..brea
0720: 6b 3b 0a 20 20 20 20 20 20 73 74 72 5b 73 74 72  k;.      str[str
0730: 5f 73 69 7a 65 2b 2b 5d 20 3d 20 63 3b 0a 20 20  _size++] = c;.  
0740: 20 20 20 20 63 20 3d 20 67 65 74 63 20 28 73 74      c = getc (st
0750: 72 65 61 6d 29 3b 0a 20 20 20 20 7d 0a 20 20 75  ream);.    }.  u
0760: 6e 67 65 74 63 20 28 63 2c 20 73 74 72 65 61 6d  ngetc (c, stream
0770: 29 3b 0a 0a 20 20 69 66 20 28 73 74 72 5f 73 69  );..  if (str_si
0780: 7a 65 20 3e 3d 20 61 6c 6c 6f 63 5f 73 69 7a 65  ze >= alloc_size
0790: 29 0a 20 20 20 20 7b 0a 20 20 20 20 20 20 73 69  ).    {.      si
07a0: 7a 65 5f 74 20 6f 6c 64 5f 61 6c 6c 6f 63 5f 73  ze_t old_alloc_s
07b0: 69 7a 65 20 3d 20 61 6c 6c 6f 63 5f 73 69 7a 65  ize = alloc_size
07c0: 3b 0a 20 20 20 20 20 20 61 6c 6c 6f 63 5f 73 69  ;.      alloc_si
07d0: 7a 65 20 3d 20 61 6c 6c 6f 63 5f 73 69 7a 65 20  ze = alloc_size 
07e0: 2a 20 33 20 2f 20 32 3b 0a 20 20 20 20 20 20 73  * 3 / 2;.      s
07f0: 74 72 20 3d 20 28 63 68 61 72 20 2a 29 20 28 2a  tr = (char *) (*
0800: 5f 5f 67 6d 70 5f 72 65 61 6c 6c 6f 63 61 74 65  __gmp_reallocate
0810: 5f 66 75 6e 63 29 20 28 73 74 72 2c 20 6f 6c 64  _func) (str, old
0820: 5f 61 6c 6c 6f 63 5f 73 69 7a 65 2c 20 61 6c 6c  _alloc_size, all
0830: 6f 63 5f 73 69 7a 65 29 3b 0a 20 20 20 20 7d 0a  oc_size);.    }.
0840: 20 20 73 74 72 5b 73 74 72 5f 73 69 7a 65 5d 20    str[str_size] 
0850: 3d 20 30 3b 0a 0a 20 20 72 65 74 76 61 6c 20 3d  = 0;..  retval =
0860: 20 6d 70 66 72 5f 73 65 74 5f 73 74 72 20 28 72   mpfr_set_str (r
0870: 6f 70 2c 20 73 74 72 2c 20 62 61 73 65 2c 20 72  op, str, base, r
0880: 6e 64 5f 6d 6f 64 65 29 3b 0a 20 20 28 2a 5f 5f  nd_mode);.  (*__
0890: 67 6d 70 5f 66 72 65 65 5f 66 75 6e 63 29 20 28  gmp_free_func) (
08a0: 73 74 72 2c 20 61 6c 6c 6f 63 5f 73 69 7a 65 29  str, alloc_size)
08b0: 3b 0a 0a 20 20 69 66 20 28 72 65 74 76 61 6c 20  ;..  if (retval 
08c0: 3d 3d 20 2d 31 29 0a 20 20 20 20 72 65 74 75 72  == -1).    retur
08d0: 6e 20 30 3b 09 09 09 2f 2a 20 65 72 72 6f 72 20  n 0;.../* error 
08e0: 2a 2f 0a 0a 20 20 72 65 74 75 72 6e 20 73 74 72  */..  return str
08f0: 5f 73 69 7a 65 20 2b 20 6e 72 65 61 64 3b 0a 7d  _size + nread;.}
0900: 0a                                               .