src/real/input/cl_R_read.cc

   1 // read_real().
   2 // This file contains a slimmed down version of read_complex().
   3 // It does not pull in all the complex function code.
   4
   5 // General includes.
   6 #include "cl_sysdep.h"
   7
   8 // Specification.
   9 #include "cln/real_io.h"
  10
  11
  12 // Implementation.
  13
  14 #include <cstring>
  15 #include "cln/input.h"
  16 #include "cln/rational_io.h"
  17 #include "cln/integer_io.h"
  18 #include "cln/float_io.h"
  19 #include "cln/integer.h"
  20 #include "cl_I.h"
  21 #include "cl_F.h"
  22 #include "cln/abort.h"
  23
  24 #undef floor
  25 #include <cmath>
  26 #define floor cln_floor
  27
  28
  29 namespace cln {
  30
  31 // Step forward over all digits, to the end of string or to the next non-digit.
  32 static const char * skip_digits (const char * ptr, const char * string_limit, unsigned int base)
  33 {
  34         for ( ; ptr != string_limit; ptr++) {
  35                 var char ch = *ptr;
  36                 if ((ch >= '0') && (ch <= '9'))
  37                         if (ch < '0' + (int)base)
  38                                 continue;
  39                         else
  40                                 break;
  41                 else {
  42                         if (base <= 10)
  43                                 break;
  44                         if (((ch >= 'A') && (ch < 'A'-10+(int)base))
  45                             || ((ch >= 'a') && (ch < 'a'-10+(int)base))
  46                            )
  47                                 continue;
  48                         else
  49                                 break;
  50                 }
  51         }
  52         return ptr;
  53 }
  54
  55 #define at_end_of_parse(ptr)  \
  56   if (end_of_parse)                                                     \
  57     { *end_of_parse = (ptr); }                                          \
  58   else                                                                  \
  59     { if ((ptr) != string_limit) { read_number_junk((ptr),string,string_limit); } }
  60
  61 const cl_R read_real (const cl_read_flags& flags, const char * string, const char * string_limit, const char * * end_of_parse)
  62 {
  63         ASSERT((flags.syntax & ~(syntax_real|syntax_maybe_bad)) == 0);
  64         // If no string_limit is given, it defaults to the end of the string.
  65         if (!string_limit)
  66                 string_limit = string + ::strlen(string);
  67         if (flags.syntax & syntax_rational) {
  68                 // Check for rational number syntax.
  69                 var unsigned int rational_base = flags.rational_base;
  70                 var const char * ptr = string;
  71                 if (flags.lsyntax & lsyntax_commonlisp) {
  72                         if (ptr == string_limit) goto not_rational_syntax;
  73                         if (*ptr == '#') {
  74                                 // Check for #b, #o, #x, #nR syntax.
  75                                 ptr++;
  76                                 if (ptr == string_limit) goto not_rational_syntax;
  77                                 switch (*ptr) {
  78                                 case 'b': case 'B':
  79                                         rational_base = 2; break;
  80                                 case 'o': case 'O':
  81                                         rational_base = 8; break;
  82                                 case 'x': case 'X':
  83                                         rational_base = 16; break;
  84                                 default:
  85                                         var const char * base_end_ptr =
  86                                                 skip_digits(ptr,string_limit,10);
  87                                         if (base_end_ptr == ptr) goto not_rational_syntax;
  88                                         if (base_end_ptr == string_limit) goto not_rational_syntax;
  89                                         if (!((*base_end_ptr == 'r') || (*base_end_ptr == 'R')))
  90                                                 goto not_rational_syntax;
  91                                         var cl_I base = read_integer(10,0,ptr,0,base_end_ptr-ptr);
  92                                         if (!((base >= 2) && (base <= 36))) {
  93                                                 fprint(std::cerr, "Base must be an integer in the range from 2 to 36, not ");
  94                                                 fprint(std::cerr, base);
  95                                                 fprint(std::cerr, "\n");
  96                                                 cl_abort();
  97                                         }
  98                                         rational_base = FN_to_UV(base); ptr = base_end_ptr;
  99                                         break;
 100                                 }
 101                                 ptr++;
 102                         }
 103                 }
 104                 var const char * ptr_after_prefix = ptr;
 105                 var cl_signean sign = 0;
 106                 if (ptr == string_limit) goto not_rational_syntax;
 107                 switch (*ptr) {
 108                         case '-': sign = ~sign;
 109                         case '+': ptr++;
 110                         default: break;
 111                 }
 112                 var const char * ptr_after_sign = ptr;
 113                 if (flags.syntax & syntax_integer) {
 114                         // Check for integer syntax:  {'+'|'-'|} {digit}+ {'.'|}
 115                         // Allow final dot only in Common Lisp syntax if there was no #<base> prefix.
 116                         if ((flags.lsyntax & lsyntax_commonlisp) && (ptr_after_prefix == string)) {
 117                                 ptr = skip_digits(ptr_after_sign,string_limit,10);
 118                                 if (ptr != ptr_after_sign)
 119                                   if (ptr != string_limit)
 120                                     if (*ptr == '.') {
 121                                         ptr++;
 122                                         if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/'))) {
 123                                                 at_end_of_parse(ptr);
 124                                                 return read_integer(10,sign,ptr_after_sign,0,ptr-ptr_after_sign);
 125                                         }
 126                                 }
 127                         }
 128                         ptr = skip_digits(ptr_after_sign,string_limit,rational_base);
 129                         if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/'))) {
 130                                 at_end_of_parse(ptr);
 131                                 return read_integer(rational_base,sign,ptr_after_sign,0,ptr-ptr_after_sign);
 132                         }
 133                 }
 134                 if (flags.syntax & syntax_ratio) {
 135                         // Check for ratio syntax: {'+'|'-'|} {digit}+ '/' {digit}+
 136                         ptr = skip_digits(ptr_after_sign,string_limit,rational_base);
 137                         if (ptr != ptr_after_sign)
 138                           if (ptr != string_limit)
 139                             if (*ptr == '/') {
 140                                 var const char * ptr_at_slash = ptr;
 141                                 ptr = skip_digits(ptr_at_slash+1,string_limit,rational_base);
 142                                 if (ptr != ptr_at_slash+1)
 143                                   if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '_') || (*ptr == '/'))) {
 144                                         at_end_of_parse(ptr);
 145                                         return read_rational(rational_base,sign,ptr_after_sign,0,ptr_at_slash-ptr_after_sign,ptr-ptr_after_sign);
 146                                 }
 147                         }
 148                 }
 149         }
 150 not_rational_syntax:
 151         if (flags.syntax & syntax_float) {
 152                 // Check for floating-point number syntax:
 153                 // {'+'|'-'|} {digit}+ {'.' {digit}* | } expo {'+'|'-'|} {digit}+
 154                 // {'+'|'-'|} {digit}* '.' {digit}+ expo {'+'|'-'|} {digit}+
 155                 // {'+'|'-'|} {digit}* '.' {digit}+
 156                 var const char * ptr = string;
 157                 var const unsigned int float_base = 10;
 158                 var cl_signean sign = 0;
 159                 if (ptr == string_limit) goto not_float_syntax;
 160                 switch (*ptr) {
 161                         case '-': sign = ~sign;
 162                         case '+': ptr++;
 163                         default: break;
 164                 }
 165                 var const char * ptr_after_sign = ptr;
 166                 var const char * ptr_after_intpart = skip_digits(ptr_after_sign,string_limit,float_base);
 167                 var cl_boolean have_dot = cl_false;
 168                 var const char * ptr_before_fracpart = ptr_after_intpart;
 169                 var const char * ptr_after_fracpart = ptr_after_intpart;
 170                 ptr = ptr_after_intpart;
 171                 if (ptr != string_limit)
 172                   if (*ptr == '.') {
 173                         have_dot = cl_true;
 174                         ptr_before_fracpart = ptr+1;
 175                         ptr_after_fracpart = skip_digits(ptr_before_fracpart,string_limit,float_base);
 176                 }
 177                 ptr = ptr_after_fracpart;
 178                 var char exponent_marker;
 179                 var cl_boolean have_exponent;
 180                 var const char * ptr_in_exponent = ptr;
 181                 var const char * ptr_after_exponent = ptr;
 182                 if ((ptr == string_limit) || !(((*ptr >= '0') && (*ptr <= '9')) || ((*ptr >= 'A') && (*ptr <= 'Z') && (*ptr != 'I')) || ((*ptr >= 'a') && (*ptr <= 'z') && (*ptr != 'i')) || (*ptr == '.') || (*ptr == '/'))) {
 183                         // No exponent.
 184                         have_exponent = cl_false;
 185                         // Must have at least one fractional part digit.
 186                         if (ptr_after_fracpart == ptr_before_fracpart) goto not_float_syntax;
 187                         exponent_marker = 'E';
 188                 } else {
 189                         have_exponent = cl_true;
 190                         // Must have at least one digit.
 191                         if (ptr_after_sign == ptr_after_intpart)
 192                                 if (ptr_after_fracpart == ptr_before_fracpart)
 193                                         goto not_float_syntax;
 194                         exponent_marker = ((*ptr >= 'a') && (*ptr <= 'z') ? *ptr - 'a' + 'A' : *ptr);
 195                         switch (exponent_marker) {
 196                                 case 'E':
 197                                 case 'S': case 'F': case 'D': case 'L':
 198                                         break;
 199                                 default:
 200                                         goto not_float_syntax;
 201                         }
 202                 }
 203                 if (have_exponent) {
 204                         ptr++;
 205                         if (ptr == string_limit) goto not_float_syntax;
 206                         switch (*ptr) {
 207                                 case '-':
 208                                 case '+': ptr++;
 209                                 default: break;
 210                         }
 211                         ptr_in_exponent = ptr;
 212                         ptr_after_exponent = skip_digits(ptr_in_exponent,string_limit,10);
 213                         if (ptr_after_exponent == ptr_in_exponent) goto not_float_syntax;
 214                 }
 215                 ptr = ptr_after_exponent;
 216                 var const char * ptr_after_prec = ptr;
 217                 var float_format_t prec;
 218                 if ((ptr != string_limit) && (*ptr == '_')) {
 219                         ptr++;
 220                         ptr_after_prec = skip_digits(ptr,string_limit,10);
 221                         if (ptr_after_prec == ptr) goto not_float_syntax;
 222                         var cl_I prec1 = digits_to_I(ptr,ptr_after_prec-ptr,10);
 223                         var uintL prec2 = cl_I_to_UL(prec1);
 224                         prec = (float_base==10 ? float_format(prec2)
 225                                                : (float_format_t)((uintL)((1+prec2)*::log((double)float_base)*1.442695041)+1)
 226                                );
 227                 } else {
 228                         switch (exponent_marker) {
 229                                 case 'S': prec = float_format_sfloat; break;
 230                                 case 'F': prec = float_format_ffloat; break;
 231                                 case 'D': prec = float_format_dfloat; break;
 232                                 case 'L': prec = flags.float_flags.default_lfloat_format; break;
 233                                 case 'E': prec = flags.float_flags.default_float_format; break;
 234                                 default: NOTREACHED
 235                         }
 236                         if (flags.float_flags.mantissa_dependent_float_format) {
 237                                 // Count the number of significant digits.
 238                                 ptr = ptr_after_sign;
 239                                 while (ptr < ptr_after_fracpart && (*ptr == '0' || *ptr == '.')) ptr++;
 240                                 var uintL num_significant_digits =
 241                                   (ptr_after_fracpart - ptr) - (ptr_before_fracpart > ptr ? 1 : 0);
 242                                 var uintL prec2 = (num_significant_digits>=2 ? num_significant_digits-2 : 0);
 243                                 var float_format_t precx =
 244                                   (float_base==10 ? float_format(prec2)
 245                                                   : (float_format_t)((uintL)((1+prec2)*::log((double)float_base)*1.442695041)+1)
 246                                   );
 247                                 if ((uintL)precx > (uintL)prec)
 248                                         prec = precx;
 249                         }
 250                 }
 251                 floatformatcase(prec
 252                 ,       if (!(flags.syntax & syntax_sfloat)) goto not_float_syntax;
 253                 ,       if (!(flags.syntax & syntax_ffloat)) goto not_float_syntax;
 254                 ,       if (!(flags.syntax & syntax_dfloat)) goto not_float_syntax;
 255                 ,       unused len;
 256                         if (!(flags.syntax & syntax_lfloat)) goto not_float_syntax;
 257                 );
 258                 at_end_of_parse(ptr_after_prec);
 259                 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);
 260         }
 261 not_float_syntax:
 262         if (flags.syntax & syntax_maybe_bad) {
 263                 ASSERT(end_of_parse);
 264                 *end_of_parse = string;
 265                 return 0; // dummy return
 266         }
 267         read_number_bad_syntax(string,string_limit);
 268 }
 269
 270 }  // namespace cln