]> www.ginac.de Git - cln.git/blob - src/float/input/cl_F_readparsed.cc
a33e9682e4b3ad8d90a3caf2d6c192ea261d63c6
[cln.git] / src / float / input / cl_F_readparsed.cc
1 // read_float().
2
3 // General includes.
4 #include "cl_sysdep.h"
5
6 // Specification.
7 #include "cl_float_io.h"
8
9
10 // Implementation.
11
12 #include "cl_integer.h"
13 #include "cl_I.h"
14 #include "cl_rational.h"
15 #include "cl_RA.h"
16 #include "cl_float.h"
17 #include "cl_sfloat.h"
18 #include "cl_ffloat.h"
19 #include "cl_dfloat.h"
20 #include "cl_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 const cl_F read_float (unsigned int base, cl_float_format_t prec, cl_signean sign, const char * string, uintL index1, uintL index4, uintL index2, uintL index3)
28 {
29         var cl_I exponent;
30         {
31                 var uintL exp_len = index2-index4; // Anzahl Stellen des Exponenten
32                 if (exp_len > 0) {
33                         var const char * ptr = &string[index4]; // zeigt auf den Exponentmarker
34                         ptr++; exp_len--; // Exponentmarker überlesen
35                         var cl_signean exp_sign = 0; // Exponenten-Vorzeichen
36                         switch (*ptr) {
37                                 case '-': exp_sign = ~exp_sign; // Vorzeichen := negativ
38                                 case '+': ptr++; exp_len--; // Exponenten-Vorzeichen überlesen
39                                 default: ;
40                         }
41                         exponent = digits_to_I(ptr,exp_len,(uintD)base); // Exponent in Integer umwandeln
42                         if (!(exp_sign==0))
43                                 exponent = -exponent; // incl. Vorzeichen
44                 } else {
45                         // kein Exponent da
46                         exponent = 0;
47                 }
48         }
49         // exponent = Exponent als Integer.
50         var cl_RA base_power = expt(base,exponent-(index4-index3)); // zu multiplizierende Zehnerpotenz
51         var cl_I mantisse = // Mantisse als Integer
52           digits_to_I(&string[index1],index4-index1,(uintD)base);
53         // Mantisse (Integer) und Zehnerpotenz (rational >0) unelegant zusammenmultiplizieren:
54         var cl_RA ratvalue;
55         if (integerp(base_power)) {
56                 DeclareType(cl_I,base_power);
57                 ratvalue = mantisse * base_power;
58         } else {
59                 // falls mantisse/=0, in exponent=1/10^i den Zähler durch mantisse
60                 // ersetzen (liefert ungekürzten Bruch, Vorsicht!)
61                 DeclareType(cl_RT,base_power);
62                 if (zerop(mantisse))
63                         ratvalue = 0;
64                 else {
65                         ASSERT(TheRatio(base_power)->refcount == 1);
66                         TheRatio(base_power)->numerator = mantisse;
67                         ratvalue = base_power;
68                 }
69         }
70         // ratvalue = Mantisse * Zehnerpotenz, als ungekürzte rationale Zahl!
71         floatformatcase(prec
72         ,       // in Short-Float umwandeln
73                 {
74                         var cl_SF x = cl_RA_to_SF(ratvalue);
75                         if (sign==0) { return x; } else { return -x; }
76                 }
77         ,       // in Single-Float umwandeln
78                 {
79                         var cl_FF x = cl_RA_to_FF(ratvalue);
80                         if (sign==0) { return x; } else { return -x; }
81                 }
82         ,       // in Double-Float umwandeln
83                 {
84                         var cl_DF x = cl_RA_to_DF(ratvalue);
85                         if (sign==0) { return x; } else { return -x; }
86                 }
87         ,       // in Long-Float umwandeln
88                 {
89                         var cl_LF x = cl_RA_to_LF(ratvalue,len);
90                         if (sign==0) { return x; } else { return -x; }
91                 }
92         );
93 }