src/float/ffloat/elem/cl_FF_scale.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
  15 namespace cln {
  16
  17 const cl_FF scale_float (const cl_FF& x, sintL delta)
  18 {
  19   // Methode:
  20   // x=0.0 -> x als Ergebnis
  21   // delta muß betragsmäßig <= FF_exp_high-FF_exp_low sein.
  22   // Neues FF mit um delta vergrößertem Exponenten bilden.
  23       // x entpacken:
  24       var cl_signean sign;
  25       var sintL exp;
  26       var uint32 mant;
  27       FF_decode(x, { return x; }, sign=,exp=,mant=);
  28       if (delta >= 0)
  29         // delta>=0
  30         { var uintL udelta = delta;
  31           if (udelta <= (uintL)(FF_exp_high-FF_exp_low))
  32             { exp = exp+udelta;
  33               return encode_FF(sign,exp,mant);
  34             }
  35             else
  36             { cl_error_floating_point_overflow(); }
  37         }
  38         else
  39         // delta<0
  40         { var uintL udelta = -delta;
  41           if (udelta <= (uintL)(FF_exp_high-FF_exp_low))
  42             { exp = exp-udelta;
  43               return encode_FF(sign,exp,mant);
  44             }
  45             else
  46             if (underflow_allowed())
  47               { cl_error_floating_point_underflow(); }
  48               else
  49               { return cl_FF_0; }
  50         }
  51 }
  52
  53 }  // namespace cln