Adding new ffi calls into the Hpc rts subsystem
[ghc-hetmet.git] / rts / gmp / mpz / inp_str.c
1 /* mpz_inp_str(dest_integer, stream, base) -- Input a number in base
2    BASE from stdio stream STREAM and store the result in DEST_INTEGER.
3
4 Copyright (C) 1991, 1993, 1994, 1996, 1998, 2000 Free Software Foundation, Inc.
5
6 This file is part of the GNU MP Library.
7
8 The GNU MP Library is free software; you can redistribute it and/or modify
9 it under the terms of the GNU Lesser General Public License as published by
10 the Free Software Foundation; either version 2.1 of the License, or (at your
11 option) any later version.
12
13 The GNU MP Library is distributed in the hope that it will be useful, but
14 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15 or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
16 License for more details.
17
18 You should have received a copy of the GNU Lesser General Public License
19 along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
20 the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
21 MA 02111-1307, USA. */
22
23 #include <stdio.h>
24 #include <ctype.h>
25 #include "gmp.h"
26 #include "gmp-impl.h"
27
28 static int
29 #if __STDC__
30 digit_value_in_base (int c, int base)
31 #else
32 digit_value_in_base (c, base)
33      int c;
34      int base;
35 #endif
36 {
37   int digit;
38
39   if (isdigit (c))
40     digit = c - '0';
41   else if (islower (c))
42     digit = c - 'a' + 10;
43   else if (isupper (c))
44     digit = c - 'A' + 10;
45   else
46     return -1;
47
48   if (digit < base)
49     return digit;
50   return -1;
51 }
52
53 size_t
54 #if __STDC__
55 mpz_inp_str (mpz_ptr x, FILE *stream, int base)
56 #else
57 mpz_inp_str (x, stream, base)
58      mpz_ptr x;
59      FILE *stream;
60      int base;
61 #endif
62 {
63   char *str;
64   size_t alloc_size, str_size;
65   int c;
66   int negative;
67   mp_size_t xsize;
68   size_t nread;
69
70   if (stream == 0)
71     stream = stdin;
72
73   nread = 0;
74
75   /* Skip whitespace.  */
76   do
77     {
78       c = getc (stream);
79       nread++;
80     }
81   while (isspace (c));
82
83   negative = 0;
84   if (c == '-')
85     {
86       negative = 1;
87       c = getc (stream);
88       nread++;
89     }
90
91   if (digit_value_in_base (c, base == 0 ? 10 : base) < 0)
92     return 0;                   /* error if no digits */
93
94   /* If BASE is 0, try to find out the base by looking at the initial
95      characters.  */
96   if (base == 0)
97     {
98       base = 10;
99       if (c == '0')
100         {
101           base = 8;
102           c = getc (stream);
103           nread++;
104           if (c == 'x' || c == 'X')
105             {
106               base = 16;
107               c = getc (stream);
108               nread++;
109             }
110           else if (c == 'b' || c == 'B')
111             {
112               base = 2;
113               c = getc (stream);
114               nread++;
115             }
116         }
117     }
118
119   /* Skip leading zeros.  */
120   while (c == '0')
121     {
122       c = getc (stream);
123       nread++;
124     }
125
126   alloc_size = 100;
127   str = (char *) (*_mp_allocate_func) (alloc_size);
128   str_size = 0;
129
130   for (;;)
131     {
132       int dig;
133       if (str_size >= alloc_size)
134         {
135           size_t old_alloc_size = alloc_size;
136           alloc_size = alloc_size * 3 / 2;
137           str = (char *) (*_mp_reallocate_func) (str, old_alloc_size, alloc_size);
138         }
139       dig = digit_value_in_base (c, base);
140       if (dig < 0)
141         break;
142       str[str_size++] = dig;
143       c = getc (stream);
144     }
145
146   ungetc (c, stream);
147
148   /* Make sure the string is not empty, mpn_set_str would fail.  */
149   if (str_size == 0)
150     {
151       x->_mp_size = 0;
152       (*_mp_free_func) (str, alloc_size);
153       return nread;
154     }
155
156   xsize = (((mp_size_t) (str_size / __mp_bases[base].chars_per_bit_exactly))
157            / BITS_PER_MP_LIMB + 2);
158   if (x->_mp_alloc < xsize)
159     _mpz_realloc (x, xsize);
160
161   /* Convert the byte array in base BASE to our bignum format.  */
162   xsize = mpn_set_str (x->_mp_d, (unsigned char *) str, str_size, base);
163   x->_mp_size = negative ? -xsize : xsize;
164
165   (*_mp_free_func) (str, alloc_size);
166   return str_size + nread;
167 }