src/float/ffloat/elem/cl_FF_scale_I.cc

   1 // scale_float().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/ffloat.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_FF.h"
  13 #include "cl_F.h"
  14 #include "cl_I.h"
  15
  16 namespace cln {
  17
  18 const cl_FF scale_float (const cl_FF& x, const cl_I& delta)
  19 {
  20   // Methode:
  21   // x=0.0 -> x als Ergebnis
  22   // delta muß ein Fixnum betragsmäßig <= FF_exp_high-FF_exp_low sein.
  23   // Neues FF mit um delta vergrößertem Exponenten bilden.
  24       // x entpacken:
  25       var cl_signean sign;
  26       var sintL exp;
  27       var uint32 mant;
  28       FF_decode(x, { return x; }, sign=,exp=,mant=);
  29       if (!minusp(delta))
  30         // delta>=0
  31         { var uintL udelta;
  32           if (fixnump(delta)
  33               && ((udelta = FN_to_L(delta)) <= (uintL)(FF_exp_high-FF_exp_low))
  34              )
  35             { exp = exp+udelta;
  36               return encode_FF(sign,exp,mant);
  37             }
  38             else
  39             { cl_error_floating_point_overflow(); }
  40         }
  41         else
  42         // delta<0
  43         { var uintL udelta;
  44           if (fixnump(delta)
  45               && ((udelta = -FN_to_L(delta)) <= (uintL)(FF_exp_high-FF_exp_low))
  46               && ((cl_value_len+1<intLsize) || !(udelta==0))
  47              )
  48             { exp = exp-udelta;
  49               return encode_FF(sign,exp,mant);
  50             }
  51             else
  52             if (underflow_allowed())
  53               { cl_error_floating_point_underflow(); }
  54               else
  55               { return cl_FF_0; }
  56         }
  57 }
  58
  59 }  // namespace cln