2 // This file contains a slimmed down version of read_complex().
3 // It does not pull in all the complex function code.
9 #include "cln/real_io.h"
16 #include "cln/input.h"
17 #include "cln/rational_io.h"
18 #include "cln/integer_io.h"
19 #include "cln/float_io.h"
20 #include "cln/integer.h"
23 #include "cln/exception.h"
27 #define floor cln_floor
32 // Step forward over all digits, to the end of string or to the next non-digit.
33 static const char * skip_digits (const char * ptr, const char * string_limit, unsigned int base)
35 for ( ; ptr != string_limit; ptr++) {
37 if ((ch >= '0') && (ch <= '9'))
38 if (ch < '0' + (int)base)
45 if (((ch >= 'A') && (ch < 'A'-10+(int)base))
46 || ((ch >= 'a') && (ch < 'a'-10+(int)base))
56 #define at_end_of_parse(ptr) \
58 { *end_of_parse = (ptr); } \
60 { if ((ptr) != string_limit) { throw read_number_junk_exception((ptr),string,string_limit); } }
62 const cl_R read_real (const cl_read_flags& flags, const char * string, const char * string_limit, const char * * end_of_parse)
64 ASSERT((flags.syntax & ~(syntax_real|syntax_maybe_bad)) == 0);
65 // If no string_limit is given, it defaults to the end of the string.
67 string_limit = string + ::strlen(string);
68 if (flags.syntax & syntax_rational) {
69 // Check for rational number syntax.
70 var unsigned int rational_base = flags.rational_base;
71 var const char * ptr = string;
72 if (flags.lsyntax & lsyntax_commonlisp) {
73 if (ptr == string_limit) goto not_rational_syntax;
75 // Check for #b, #o, #x, #nR syntax.
77 if (ptr == string_limit) goto not_rational_syntax;
80 rational_base = 2; break;
82 rational_base = 8; break;
84 rational_base = 16; break;
86 var const char * base_end_ptr =
87 skip_digits(ptr,string_limit,10);
88 if (base_end_ptr == ptr) goto not_rational_syntax;
89 if (base_end_ptr == string_limit) goto not_rational_syntax;
90 if (!((*base_end_ptr == 'r') || (*base_end_ptr == 'R')))
91 goto not_rational_syntax;
92 var cl_I base = read_integer(10,0,ptr,0,base_end_ptr-ptr);
93 if (!((base >= 2) && (base <= 36))) {
94 std::ostringstream buf;
95 fprint(buf, "Base must be an integer in the range from 2 to 36, not ");
97 throw runtime_exception(buf.str());
99 rational_base = FN_to_UV(base); ptr = base_end_ptr;
105 var const char * ptr_after_prefix = ptr;
106 var cl_signean sign = 0;
107 if (ptr == string_limit) goto not_rational_syntax;
109 case '-': sign = ~sign;
113 var const char * ptr_after_sign = ptr;
114 if (flags.syntax & syntax_integer) {
115 // Check for integer syntax: {'+'|'-'|} {digit}+ {'.'|}
116 // Allow final dot only in Common Lisp syntax if there was no #<base> prefix.
117 if ((flags.lsyntax & lsyntax_commonlisp) && (ptr_after_prefix == string)) {
118 ptr = skip_digits(ptr_after_sign,string_limit,10);
119 if (ptr != ptr_after_sign)
120 if (ptr != string_limit)
123 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/'))) {
124 at_end_of_parse(ptr);
125 return read_integer(10,sign,ptr_after_sign,0,ptr-ptr_after_sign);
129 ptr = skip_digits(ptr_after_sign,string_limit,rational_base);
130 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/'))) {
131 at_end_of_parse(ptr);
132 return read_integer(rational_base,sign,ptr_after_sign,0,ptr-ptr_after_sign);
135 if (flags.syntax & syntax_ratio) {
136 // Check for ratio syntax: {'+'|'-'|} {digit}+ '/' {digit}+
137 ptr = skip_digits(ptr_after_sign,string_limit,rational_base);
138 if (ptr != ptr_after_sign)
139 if (ptr != string_limit)
141 var const char * ptr_at_slash = ptr;
142 ptr = skip_digits(ptr_at_slash+1,string_limit,rational_base);
143 if (ptr != ptr_at_slash+1)
144 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/'))) {
145 at_end_of_parse(ptr);
146 return read_rational(rational_base,sign,ptr_after_sign,0,ptr_at_slash-ptr_after_sign,ptr-ptr_after_sign);
152 if (flags.syntax & syntax_float) {
153 // Check for floating-point number syntax:
154 // {'+'|'-'|} {digit}+ {'.' {digit}* | } expo {'+'|'-'|} {digit}+
155 // {'+'|'-'|} {digit}* '.' {digit}+ expo {'+'|'-'|} {digit}+
156 // {'+'|'-'|} {digit}* '.' {digit}+
157 var const char * ptr = string;
158 var const unsigned int float_base = 10;
159 var cl_signean sign = 0;
160 if (ptr == string_limit) goto not_float_syntax;
162 case '-': sign = ~sign;
166 var const char * ptr_after_sign = ptr;
167 var const char * ptr_after_intpart = skip_digits(ptr_after_sign,string_limit,float_base);
168 var bool have_dot = false;
169 var const char * ptr_before_fracpart = ptr_after_intpart;
170 var const char * ptr_after_fracpart = ptr_after_intpart;
171 ptr = ptr_after_intpart;
172 if (ptr != string_limit)
175 ptr_before_fracpart = ptr+1;
176 ptr_after_fracpart = skip_digits(ptr_before_fracpart,string_limit,float_base);
178 ptr = ptr_after_fracpart;
179 var char exponent_marker;
180 var bool have_exponent;
181 var const char * ptr_in_exponent = ptr;
182 var const char * ptr_after_exponent = ptr;
183 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '/'))) {
185 have_exponent = false;
186 // Must have at least one fractional part digit.
187 if (ptr_after_fracpart == ptr_before_fracpart) goto not_float_syntax;
188 exponent_marker = 'E';
190 have_exponent = true;
191 // Must have at least one digit.
192 if (ptr_after_sign == ptr_after_intpart)
193 if (ptr_after_fracpart == ptr_before_fracpart)
194 goto not_float_syntax;
195 exponent_marker = ((*ptr >= 'a') && (*ptr <= 'z') ? *ptr - 'a' + 'A' : *ptr);
196 switch (exponent_marker) {
198 case 'S': case 'F': case 'D': case 'L':
201 goto not_float_syntax;
206 if (ptr == string_limit) goto not_float_syntax;
212 ptr_in_exponent = ptr;
213 ptr_after_exponent = skip_digits(ptr_in_exponent,string_limit,10);
214 if (ptr_after_exponent == ptr_in_exponent) goto not_float_syntax;
216 ptr = ptr_after_exponent;
217 var const char * ptr_after_prec = ptr;
218 var float_format_t prec;
219 if ((ptr != string_limit) && (*ptr == '_')) {
221 ptr_after_prec = skip_digits(ptr,string_limit,10);
222 if (ptr_after_prec == ptr) goto not_float_syntax;
223 var cl_I prec1 = digits_to_I(ptr,ptr_after_prec-ptr,10);
224 var uintC prec2 = cl_I_to_ulong(prec1);
225 prec = (float_base==10 ? float_format(prec2)
226 : (float_format_t)((uintC)((1+prec2)*::log((double)float_base)*1.442695041)+1)
229 switch (exponent_marker) {
230 case 'S': prec = float_format_sfloat; break;
231 case 'F': prec = float_format_ffloat; break;
232 case 'D': prec = float_format_dfloat; break;
233 case 'L': prec = flags.float_flags.default_lfloat_format; break;
234 case 'E': prec = flags.float_flags.default_float_format; break;
237 if (flags.float_flags.mantissa_dependent_float_format) {
238 // Count the number of significant digits.
239 ptr = ptr_after_sign;
240 while (ptr < ptr_after_fracpart && (*ptr == '0' || *ptr == '.')) ptr++;
241 var uintC num_significant_digits =
242 (ptr_after_fracpart - ptr) - (ptr_before_fracpart > ptr ? 1 : 0);
243 var uintC prec2 = (num_significant_digits>=2 ? num_significant_digits-2 : 0);
244 var float_format_t precx =
245 (float_base==10 ? float_format(prec2)
246 : (float_format_t)((uintC)((1+prec2)*::log((double)float_base)*1.442695041)+1)
248 if ((uintC)precx > (uintC)prec)
253 , if (!(flags.syntax & syntax_sfloat)) goto not_float_syntax;
254 , if (!(flags.syntax & syntax_ffloat)) goto not_float_syntax;
255 , if (!(flags.syntax & syntax_dfloat)) goto not_float_syntax;
257 if (!(flags.syntax & syntax_lfloat)) goto not_float_syntax;
259 at_end_of_parse(ptr_after_prec);
260 return 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);
263 if (flags.syntax & syntax_maybe_bad) {
264 ASSERT(end_of_parse);
265 *end_of_parse = string;
266 return 0; // dummy return
268 throw read_number_bad_syntax_exception(string,string_limit);