src/float/input/cl_F_readparsed.cc

   1 // read_float().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/float_io.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cln/integer.h"
  13 #include "cl_I.h"
  14 #include "cln/rational.h"
  15 #include "cl_RA.h"
  16 #include "cln/float.h"
  17 #include "cln/sfloat.h"
  18 #include "cln/ffloat.h"
  19 #include "cln/dfloat.h"
  20 #include "cln/lfloat.h"
  21 #include "cl_F.h"
  22 #include "cl_SF.h"
  23 #include "cl_FF.h"
  24 #include "cl_DF.h"
  25 #include "cl_LF.h"
  26
  27 namespace cln {
  28
  29 #if intCsize > intLsize
  30 const cl_F read_float (unsigned int base, float_format_t prec, cl_signean sign, const char * string, uintL index1, uintL index4, uintL index2, uintL index3)
  31 {
  32         // XXX: This signature is for binary compatibility with CLN-1.2.0 only.
  33         return read_float(base, prec, sign, string, uintC(index1), uintC(index4), uintC(index2), uintC(index3));
  34 }
  35 #endif
  36
  37 const cl_F read_float (unsigned int base, float_format_t prec, cl_signean sign, const char * string, uintC index1, uintC index4, uintC index2, uintC index3)
  38 {
  39         var cl_I exponent;
  40         {
  41                 var uintC exp_len = index2-index4; // Anzahl Stellen des Exponenten
  42                 if (exp_len > 0) {
  43                         var const char * ptr = &string[index4]; // zeigt auf den Exponentmarker
  44                         ptr++; exp_len--; // Exponentmarker überlesen
  45                         var cl_signean exp_sign = 0; // Exponenten-Vorzeichen
  46                         switch (*ptr) {
  47                                 case '-': exp_sign = ~exp_sign; // Vorzeichen := negativ
  48                                 case '+': ptr++; exp_len--; // Exponenten-Vorzeichen überlesen
  49                                 default: ;
  50                         }
  51                         exponent = digits_to_I(ptr,exp_len,(uintD)base); // Exponent in Integer umwandeln
  52                         if (!(exp_sign==0))
  53                                 exponent = -exponent; // incl. Vorzeichen
  54                 } else {
  55                         // kein Exponent da
  56                         exponent = 0;
  57                 }
  58         }
  59         // exponent = Exponent als Integer.
  60         var cl_RA base_power = expt(base,exponent-(index4-index3)); // zu multiplizierende Zehnerpotenz
  61         var cl_I mantisse = // Mantisse als Integer
  62           digits_to_I(&string[index1],index4-index1,(uintD)base);
  63         // Mantisse (Integer) und Zehnerpotenz (rational >0) unelegant zusammenmultiplizieren:
  64         var cl_RA ratvalue;
  65         if (integerp(base_power)) {
  66                 DeclareType(cl_I,base_power);
  67                 ratvalue = mantisse * base_power;
  68         } else {
  69                 // falls mantisse/=0, in exponent=1/10^i den Zähler durch mantisse
  70                 // ersetzen (liefert ungekürzten Bruch, Vorsicht!)
  71                 DeclareType(cl_RT,base_power);
  72                 if (zerop(mantisse))
  73                         ratvalue = 0;
  74                 else {
  75                         ASSERT(TheRatio(base_power)->refcount == 1);
  76                         TheRatio(base_power)->numerator = mantisse;
  77                         ratvalue = base_power;
  78                 }
  79         }
  80         // ratvalue = Mantisse * Zehnerpotenz, als ungekürzte rationale Zahl!
  81         floatformatcase(prec
  82         ,       // in Short-Float umwandeln
  83                 {
  84                         var cl_SF x = cl_RA_to_SF(ratvalue);
  85                         if (sign==0) { return x; } else { return -x; }
  86                 }
  87         ,       // in Single-Float umwandeln
  88                 {
  89                         var cl_FF x = cl_RA_to_FF(ratvalue);
  90                         if (sign==0) { return x; } else { return -x; }
  91                 }
  92         ,       // in Double-Float umwandeln
  93                 {
  94                         var cl_DF x = cl_RA_to_DF(ratvalue);
  95                         if (sign==0) { return x; } else { return -x; }
  96                 }
  97         ,       // in Long-Float umwandeln
  98                 {
  99                         var cl_LF x = cl_RA_to_LF(ratvalue,len);
 100                         if (sign==0) { return x; } else { return -x; }
 101                 }
 102         );
 103 }
 104
 105 }  // namespace cln