2 // This file contains a slimmed down version of read_real().
3 // It does not pull in all the generic real number code.
9 #include "cln/float_io.h"
15 #include "cln/input.h"
16 #include "cln/integer.h"
22 #define floor cln_floor
26 // Step forward over all digits, to the end of string or to the next non-digit.
27 static const char * skip_digits (const char * ptr, const char * string_limit, unsigned int base)
29 for ( ; ptr != string_limit; ptr++) {
31 if ((ch >= '0') && (ch <= '9'))
32 if (ch < '0' + (int)base)
39 if (((ch >= 'A') && (ch < 'A'-10+(int)base))
40 || ((ch >= 'a') && (ch < 'a'-10+(int)base))
50 #define at_end_of_parse(ptr) \
52 { *end_of_parse = (ptr); } \
54 { if ((ptr) != string_limit) { read_number_junk((ptr),string,string_limit); } }
56 const cl_F read_float (const cl_read_flags& flags, const char * string, const char * string_limit, const char * * end_of_parse)
58 ASSERT((flags.syntax & ~(syntax_float|syntax_maybe_bad)) == 0);
59 // If no string_limit is given, it defaults to the end of the string.
61 string_limit = string + ::strlen(string);
62 if (flags.syntax & syntax_float) {
63 // Check for floating-point number syntax:
64 // {'+'|'-'|} {digit}+ {'.' {digit}* | } expo {'+'|'-'|} {digit}+
65 // {'+'|'-'|} {digit}* '.' {digit}+ expo {'+'|'-'|} {digit}+
66 // {'+'|'-'|} {digit}* '.' {digit}+
67 var const char * ptr = string;
68 var const unsigned int float_base = 10;
69 var cl_signean sign = 0;
70 if (ptr == string_limit) goto not_float_syntax;
72 case '-': sign = ~sign;
76 var const char * ptr_after_sign = ptr;
77 var const char * ptr_after_intpart = skip_digits(ptr_after_sign,string_limit,float_base);
78 var cl_boolean have_dot = cl_false;
79 var const char * ptr_before_fracpart = ptr_after_intpart;
80 var const char * ptr_after_fracpart = ptr_after_intpart;
81 ptr = ptr_after_intpart;
82 if (ptr != string_limit)
85 ptr_before_fracpart = ptr+1;
86 ptr_after_fracpart = skip_digits(ptr_before_fracpart,string_limit,float_base);
88 ptr = ptr_after_fracpart;
89 var char exponent_marker;
90 var cl_boolean have_exponent;
91 var const char * ptr_in_exponent = ptr;
92 var const char * ptr_after_exponent = ptr;
93 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '/'))) {
95 have_exponent = cl_false;
96 // Must have at least one fractional part digit.
97 if (ptr_after_fracpart == ptr_before_fracpart) goto not_float_syntax;
98 exponent_marker = 'E';
100 have_exponent = cl_true;
101 // Must have at least one digit.
102 if (ptr_after_sign == ptr_after_intpart)
103 if (ptr_after_fracpart == ptr_before_fracpart)
104 goto not_float_syntax;
105 exponent_marker = ((*ptr >= 'a') && (*ptr <= 'z') ? *ptr - 'a' + 'A' : *ptr);
106 switch (exponent_marker) {
108 case 'S': case 'F': case 'D': case 'L':
111 goto not_float_syntax;
116 if (ptr == string_limit) goto not_float_syntax;
122 ptr_in_exponent = ptr;
123 ptr_after_exponent = skip_digits(ptr_in_exponent,string_limit,10);
124 if (ptr_after_exponent == ptr_in_exponent) goto not_float_syntax;
126 ptr = ptr_after_exponent;
127 var const char * ptr_after_prec = ptr;
128 var float_format_t prec;
129 if ((ptr != string_limit) && (*ptr == '_')) {
131 ptr_after_prec = skip_digits(ptr,string_limit,10);
132 if (ptr_after_prec == ptr) goto not_float_syntax;
133 var cl_I prec1 = digits_to_I(ptr,ptr_after_prec-ptr,10);
134 var uintL prec2 = cl_I_to_UL(prec1);
135 prec = (float_base==10 ? float_format(prec2)
136 : (float_format_t)((uintL)((1+prec2)*::log((double)float_base)*1.442695041)+1)
139 switch (exponent_marker) {
140 case 'S': prec = float_format_sfloat; break;
141 case 'F': prec = float_format_ffloat; break;
142 case 'D': prec = float_format_dfloat; break;
143 case 'L': prec = flags.float_flags.default_lfloat_format; break;
144 case 'E': prec = flags.float_flags.default_float_format; break;
147 if (flags.float_flags.mantissa_dependent_float_format) {
148 // Count the number of significant digits.
149 ptr = ptr_after_sign;
150 while (ptr < ptr_after_fracpart && (*ptr == '0' || *ptr == '.')) ptr++;
151 var uintL num_significant_digits =
152 (ptr_after_fracpart - ptr) - (ptr_before_fracpart > ptr ? 1 : 0);
153 var uintL prec2 = (num_significant_digits>=2 ? num_significant_digits-2 : 0);
154 var float_format_t precx =
155 (float_base==10 ? float_format(prec2)
156 : (float_format_t)((uintL)((1+prec2)*::log((double)float_base)*1.442695041)+1)
158 if ((uintL)precx > (uintL)prec)
163 , if (!(flags.syntax & syntax_sfloat)) goto not_float_syntax;
164 , if (!(flags.syntax & syntax_ffloat)) goto not_float_syntax;
165 , if (!(flags.syntax & syntax_dfloat)) goto not_float_syntax;
167 if (!(flags.syntax & syntax_lfloat)) goto not_float_syntax;
169 at_end_of_parse(ptr_after_prec);
170 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);
174 if (flags.syntax & syntax_maybe_bad) {
175 ASSERT(end_of_parse);
176 *end_of_parse = string;
177 return cl_F(); // dummy return
179 read_number_bad_syntax(string,string_limit);