src/float/misc/cl_F_rational.cc

   1 // rational().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/float.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cln/integer.h"
  13 #include "cl_RA.h"
  14
  15 namespace cln {
  16
  17 const cl_RA rational (const cl_F& x)
  18 {
  19   // Methode:
  20   // Der mathematische Wert eines Float ist, wenn INTEGER-DECODE-FLOAT die
  21   // drei Zahlen m,e,s (Mantisse, Exponent, Vorzeichen) liefert,
  22   // = s * 2^e * m.
  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);
  31         if (!minusp(e))
  32                 return ash(n,e);
  33         else {
  34 #if 0
  35                 return I_posI_div_RA(n, ash(1,-e));
  36 #else // spart ggT
  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);
  42                 if (k >= _e)
  43                         // Kürze mit 2^(-e).
  44                         return ash(n,e);
  45                 else
  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));
  49 #endif
  50         }
  51 }
  52
  53 }  // namespace cln