src/float/input/cl_F_read.cc

   1 // read_float().
   2 // This file contains a slimmed down version of read_real().
   3 // It does not pull in all the generic real number code.
   4
   5 // General includes.
   6 #include "cl_sysdep.h"
   7
   8 // Specification.
   9 #include "cln/float_io.h"
  10
  11
  12 // Implementation.
  13
  14 #include <string.h>
  15 #include "cln/input.h"
  16 #include "cln/integer.h"
  17 #include "cl_I.h"
  18 #include "cl_F.h"
  19
  20 #undef floor
  21 #include <cmath>
  22 #define floor cln_floor
  23
  24 namespace cln {
  25
  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)
  28 {
  29         for ( ; ptr != string_limit; ptr++) {
  30                 var char ch = *ptr;
  31                 if ((ch >= '0') && (ch <= '9'))
  32                         if (ch < '0' + (int)base)
  33                                 continue;
  34                         else
  35                                 break;
  36                 else {
  37                         if (base <= 10)
  38                                 break;
  39                         if (((ch >= 'A') && (ch < 'A'-10+(int)base))
  40                             || ((ch >= 'a') && (ch < 'a'-10+(int)base))
  41                            )
  42                                 continue;
  43                         else
  44                                 break;
  45                 }
  46         }
  47         return ptr;
  48 }
  49
  50 #define at_end_of_parse(ptr)  \
  51   if (end_of_parse)                                                     \
  52     { *end_of_parse = (ptr); }                                          \
  53   else                                                                  \
  54     { if ((ptr) != string_limit) { read_number_junk((ptr),string,string_limit); } }
  55
  56 const cl_F read_float (const cl_read_flags& flags, const char * string, const char * string_limit, const char * * end_of_parse)
  57 {
  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.
  60         if (!string_limit)
  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;
  71                 switch (*ptr) {
  72                         case '-': sign = ~sign;
  73                         case '+': ptr++;
  74                         default: break;
  75                 }
  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)
  83                   if (*ptr == '.') {
  84                         have_dot = cl_true;
  85                         ptr_before_fracpart = ptr+1;
  86                         ptr_after_fracpart = skip_digits(ptr_before_fracpart,string_limit,float_base);
  87                 }
  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 == '/'))) {
  94                         // No exponent.
  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';
  99                 } else {
 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) {
 107                                 case 'E':
 108                                 case 'S': case 'F': case 'D': case 'L':
 109                                         break;
 110                                 default:
 111                                         goto not_float_syntax;
 112                         }
 113                 }
 114                 if (have_exponent) {
 115                         ptr++;
 116                         if (ptr == string_limit) goto not_float_syntax;
 117                         switch (*ptr) {
 118                                 case '-':
 119                                 case '+': ptr++;
 120                                 default: break;
 121                         }
 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;
 125                 }
 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 == '_')) {
 130                         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)
 137                                );
 138                 } else {
 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;
 145                                 default: NOTREACHED
 146                         }
 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)
 157                                   );
 158                                 if ((uintL)precx > (uintL)prec)
 159                                         prec = precx;
 160                         }
 161                 }
 162                 floatformatcase(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;
 166                 ,       unused len;
 167                         if (!(flags.syntax & syntax_lfloat)) goto not_float_syntax;
 168                 );
 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);
 171         }
 172 not_float_syntax:
 173 bad_syntax:
 174         if (flags.syntax & syntax_maybe_bad) {
 175                 ASSERT(end_of_parse);
 176                 *end_of_parse = string;
 177                 return cl_F(); // dummy return
 178         }
 179         read_number_bad_syntax(string,string_limit);
 180 }
 181
 182 }  // namespace cln