12 #include "cln/integer.h"
17 const cl_RA rational (const cl_F& x)
20 // Der mathematische Wert eines Float ist, wenn INTEGER-DECODE-FLOAT die
21 // drei Zahlen m,e,s (Mantisse, Exponent, Vorzeichen) liefert,
23 // n:=m. Falls s<0, setze n:=-m.
24 // Falls e>=0, ist (ash n e) das Ergebnis,
25 // sonst ist die rationale Zahl (/ n (ash 1 (- e))) das Ergebnis.
26 var cl_idecoded_float x_decoded = integer_decode_float(x);
27 var cl_I& m = x_decoded.mantissa;
28 var cl_I& e = x_decoded.exponent;
29 var cl_I& s = x_decoded.sign;
30 var cl_I n = (!minusp(s) ? m : -m);
35 return I_posI_div_RA(n, ash(1,-e));
37 // n /= 0, -e > 0. Kürze mit ggT(n,2^(-e)) = 2^min(ord2(n),-e).
38 // 0 < -e <= LF_exp_mid-LF_exp_low + intDsize*len < 2^32,
39 var cl_I minus_e = -e;
40 var uintL _e = cl_I_to_UL(minus_e); // daher kein Überlauf
41 var uintC k = ord2(n);
46 // Kürze mit 2^k, 0 <= k < -e.
47 return I_I_to_RT(ash(n,-(sintC)k),
48 ash(1,minus_e-(cl_I)(unsigned long)k));