7 #include "cln/complex_io.h"
13 #include "cln/input.h"
14 #include "cln/real_io.h"
15 #include "cln/float_io.h"
16 #include "cln/rational_io.h"
17 #include "cln/integer_io.h"
18 #include "cln/integer.h"
22 #include "cln/abort.h"
26 #define floor cln_floor
27 // Ugh, some compilers #define stderr, confusing cln::stderr
35 // Step forward over all digits, to the end of string or to the next non-digit.
36 static const char * skip_digits (const char * ptr, const char * string_limit, unsigned int base)
38 for ( ; ptr != string_limit; ptr++) {
40 if ((ch >= '0') && (ch <= '9'))
41 if (ch < '0' + (int)base)
48 if (((ch >= 'A') && (ch < 'A'-10+(int)base))
49 || ((ch >= 'a') && (ch < 'a'-10+(int)base))
59 // Finish reading the "+yi" part of "x+yi" when "x" has already been read.
60 static const cl_N read_complex_number_rest (const cl_read_flags& flags, const char * string_rest, const char * string, const char * string_limit, const char * * end_of_parse, const cl_R& x);
62 #define at_end_of_parse(ptr) \
64 { *end_of_parse = (ptr); } \
66 { if ((ptr) != string_limit) { read_number_junk((ptr),string,string_limit); } }
68 const cl_N read_complex (const cl_read_flags& flags, const char * string, const char * string_limit, const char * * end_of_parse)
70 // If no string_limit is given, it defaults to the end of the string.
72 string_limit = string + ::strlen(string);
73 if (flags.syntax & syntax_rational) {
74 // Check for rational number syntax.
75 var unsigned int rational_base = flags.rational_base;
76 var const char * ptr = string;
77 if (flags.lsyntax & lsyntax_commonlisp) {
78 if (ptr == string_limit) goto not_rational_syntax;
80 // Check for #b, #o, #x, #nR syntax.
82 if (ptr == string_limit) goto not_rational_syntax;
85 rational_base = 2; break;
87 rational_base = 8; break;
89 rational_base = 16; break;
91 var const char * base_end_ptr =
92 skip_digits(ptr,string_limit,10);
93 if (base_end_ptr == ptr) goto not_rational_syntax;
94 if (base_end_ptr == string_limit) goto not_rational_syntax;
95 if (!((*base_end_ptr == 'r') || (*base_end_ptr == 'R')))
96 goto not_rational_syntax;
97 var cl_I base = read_integer(10,0,ptr,0,base_end_ptr-ptr);
98 if (!((base >= 2) && (base <= 36))) {
99 fprint(stderr, "Base must be an integer in the range from 2 to 36, not ");
100 fprint(stderr, base);
101 fprint(stderr, "\n");
104 rational_base = FN_to_UL(base); ptr = base_end_ptr;
110 var const char * ptr_after_prefix = ptr;
111 var cl_signean sign = 0;
112 if (ptr == string_limit) goto not_rational_syntax;
114 case '-': sign = ~sign;
118 var const char * ptr_after_sign = ptr;
119 if (flags.syntax & syntax_integer) {
120 // Check for integer syntax: {'+'|'-'|} {digit}+ {'.'|}
121 // Allow final dot only in Common Lisp syntax if there was no #<base> prefix.
122 if ((flags.lsyntax & lsyntax_commonlisp) && (ptr_after_prefix == string)) {
123 ptr = skip_digits(ptr_after_sign,string_limit,10);
124 if (ptr != ptr_after_sign)
125 if (ptr != string_limit)
128 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/')))
129 return read_complex_number_rest(flags,ptr,string,string_limit,end_of_parse,
130 read_integer(10,sign,ptr_after_sign,0,ptr-ptr_after_sign));
133 ptr = skip_digits(ptr_after_sign,string_limit,rational_base);
134 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/')))
135 return read_complex_number_rest(flags,ptr,string,string_limit,end_of_parse,
136 read_integer(rational_base,sign,ptr_after_sign,0,ptr-ptr_after_sign));
138 if (flags.syntax & syntax_ratio) {
139 // Check for ratio syntax: {'+'|'-'|} {digit}+ '/' {digit}+
140 ptr = skip_digits(ptr_after_sign,string_limit,rational_base);
141 if (ptr != ptr_after_sign)
142 if (ptr != string_limit)
144 var const char * ptr_at_slash = ptr;
145 ptr = skip_digits(ptr_at_slash+1,string_limit,rational_base);
146 if (ptr != ptr_at_slash+1)
147 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/')))
148 return read_complex_number_rest(flags,ptr,string,string_limit,end_of_parse,
149 read_rational(rational_base,sign,ptr_after_sign,0,ptr_at_slash-ptr_after_sign,ptr-ptr_after_sign));
154 if (flags.syntax & syntax_float) {
155 // Check for floating-point number syntax:
156 // {'+'|'-'|} {digit}+ {'.' {digit}* | } expo {'+'|'-'|} {digit}+
157 // {'+'|'-'|} {digit}* '.' {digit}+ expo {'+'|'-'|} {digit}+
158 // {'+'|'-'|} {digit}* '.' {digit}+
159 var const char * ptr = string;
160 var const unsigned int float_base = 10;
161 var cl_signean sign = 0;
162 if (ptr == string_limit) goto not_float_syntax;
164 case '-': sign = ~sign;
168 var const char * ptr_after_sign = ptr;
169 var const char * ptr_after_intpart = skip_digits(ptr_after_sign,string_limit,float_base);
170 var cl_boolean have_dot = cl_false;
171 var const char * ptr_before_fracpart = ptr_after_intpart;
172 var const char * ptr_after_fracpart = ptr_after_intpart;
173 ptr = ptr_after_intpart;
174 if (ptr != string_limit)
177 ptr_before_fracpart = ptr+1;
178 ptr_after_fracpart = skip_digits(ptr_before_fracpart,string_limit,float_base);
180 ptr = ptr_after_fracpart;
181 var char exponent_marker;
182 var cl_boolean have_exponent;
183 var const char * ptr_in_exponent = ptr;
184 var const char * ptr_after_exponent = ptr;
185 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '/'))) {
187 have_exponent = cl_false;
188 // Must have at least one fractional part digit.
189 if (ptr_after_fracpart == ptr_before_fracpart) goto not_float_syntax;
190 exponent_marker = 'E';
192 have_exponent = cl_true;
193 // Must have at least one digit.
194 if (ptr_after_sign == ptr_after_intpart)
195 if (ptr_after_fracpart == ptr_before_fracpart)
196 goto not_float_syntax;
197 exponent_marker = ((*ptr >= 'a') && (*ptr <= 'z') ? *ptr - 'a' + 'A' : *ptr);
198 switch (exponent_marker) {
200 case 'S': case 'F': case 'D': case 'L':
203 goto not_float_syntax;
208 if (ptr == string_limit) goto not_float_syntax;
214 ptr_in_exponent = ptr;
215 ptr_after_exponent = skip_digits(ptr_in_exponent,string_limit,10);
216 if (ptr_after_exponent == ptr_in_exponent) goto not_float_syntax;
218 ptr = ptr_after_exponent;
219 var const char * ptr_after_prec = ptr;
220 var float_format_t prec;
221 if ((ptr != string_limit) && (*ptr == '_')) {
223 ptr_after_prec = skip_digits(ptr,string_limit,10);
224 if (ptr_after_prec == ptr) goto not_float_syntax;
225 var cl_I prec1 = digits_to_I(ptr,ptr_after_prec-ptr,10);
226 var uintL prec2 = cl_I_to_UL(prec1);
227 prec = (float_base==10 ? float_format(prec2)
228 : (float_format_t)((uintL)((1+prec2)*::log((double)float_base)*1.442695041)+1)
231 switch (exponent_marker) {
232 case 'S': prec = float_format_sfloat; break;
233 case 'F': prec = float_format_ffloat; break;
234 case 'D': prec = float_format_dfloat; break;
235 case 'L': prec = flags.float_flags.default_lfloat_format; break;
236 case 'E': prec = flags.float_flags.default_float_format; break;
239 if (flags.float_flags.mantissa_dependent_float_format) {
240 // Count the number of significant digits.
241 ptr = ptr_after_sign;
242 while (ptr < ptr_after_fracpart && (*ptr == '0' || *ptr == '.')) ptr++;
243 var uintL num_significant_digits =
244 (ptr_after_fracpart - ptr) - (ptr_before_fracpart > ptr ? 1 : 0);
245 var uintL prec2 = (num_significant_digits>=2 ? num_significant_digits-2 : 0);
246 var float_format_t precx =
247 (float_base==10 ? float_format(prec2)
248 : (float_format_t)((uintL)((1+prec2)*::log((double)float_base)*1.442695041)+1)
250 if ((uintL)precx > (uintL)prec)
255 , if (!(flags.syntax & syntax_sfloat)) goto not_float_syntax;
256 , if (!(flags.syntax & syntax_ffloat)) goto not_float_syntax;
257 , if (!(flags.syntax & syntax_dfloat)) goto not_float_syntax;
259 if (!(flags.syntax & syntax_lfloat)) goto not_float_syntax;
261 return read_complex_number_rest(flags,ptr_after_prec,string,string_limit,end_of_parse,
262 read_float(float_base,prec,sign,ptr_after_sign,0,ptr_after_fracpart-ptr_after_sign,ptr_after_exponent-ptr_after_sign,ptr_before_fracpart-ptr_after_sign));
265 if ((flags.syntax & syntax_complex) && (flags.lsyntax & lsyntax_commonlisp)) {
266 // Check for complex number syntax:
267 // '#' {'C'|'c'} '(' realpart {' '}+ imagpart ')'
268 var const char * ptr = string;
269 if (ptr == string_limit) goto not_complex_syntax;
270 if (!(*ptr == '#')) goto not_complex_syntax;
272 if (ptr == string_limit) goto not_complex_syntax;
273 if (!((*ptr == 'C') || (*ptr == 'c'))) goto not_complex_syntax;
275 // Modified flags for parsing the realpart and imagpart:
276 var cl_read_flags flags_for_parts = flags;
277 flags_for_parts.syntax = (cl_read_syntax_t)((flags_for_parts.syntax & ~syntax_complex) | syntax_maybe_bad);
278 var const char * end_of_part;
279 if (ptr == string_limit) goto not_complex_syntax;
280 if (!(*ptr == '(')) goto not_complex_syntax;
282 var cl_R realpart = read_real(flags_for_parts,ptr,string_limit,&end_of_part);
283 if (end_of_part == ptr) goto not_complex_syntax;
285 if (ptr == string_limit) goto not_complex_syntax;
286 if (!(*ptr == ' ')) goto not_complex_syntax;
288 while ((ptr != string_limit) && (*ptr == ' ')) { ptr++; }
289 var cl_R imagpart = read_real(flags_for_parts,ptr,string_limit,&end_of_part);
290 if (end_of_part == ptr) goto not_complex_syntax;
292 if (ptr == string_limit) goto not_complex_syntax;
293 if (!(*ptr == ')')) goto not_complex_syntax;
295 at_end_of_parse(ptr);
296 return complex(realpart,imagpart);
300 if (flags.syntax & syntax_maybe_bad) {
301 ASSERT(end_of_parse);
302 *end_of_parse = string;
303 return 0; // dummy return
305 read_number_bad_syntax(string,string_limit);
308 static const cl_N read_complex_number_rest (const cl_read_flags& flags, const char * string_rest, const char * string, const char * string_limit, const char * * end_of_parse, const cl_R& x)
311 if ((flags.syntax & syntax_complex) && (flags.lsyntax & lsyntax_algebraic)) {
312 // Finish reading the "+yi" part of "x+yi".
313 // We allow "y" to begin with a '-'.
314 // We also allow the '+' to be replaced by '-', but in this case
315 // "y" may not begin with a '-'.
316 // We also allow the syntax "xi" (implicit realpart = 0).
317 var const char * ptr = string_rest;
318 if (ptr == string_limit) goto not_complex_syntax;
319 if ((*ptr == 'i') || (*ptr == 'I')) {
321 at_end_of_parse(ptr);
327 default: goto not_complex_syntax;
329 // Modified flags for parsing the imagpart:
330 var cl_read_flags flags_for_part = flags;
331 flags_for_part.syntax = (cl_read_syntax_t)((flags_for_part.syntax & ~syntax_complex) | syntax_maybe_bad);
332 var const char * end_of_part;
333 var const cl_R& realpart = x;
334 var cl_R imagpart = read_real(flags_for_part,ptr,string_limit,&end_of_part);
335 if (end_of_part == ptr) goto not_complex_syntax;
337 if (ptr == string_limit) goto not_complex_syntax;
338 if (!((*ptr == 'i') || (*ptr == 'I'))) goto not_complex_syntax;
340 at_end_of_parse(ptr);
341 return complex(realpart,imagpart);
344 at_end_of_parse(string_rest);