src/float/sfloat/elem/cl_SF_scale.cc

   1 // scale_float().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_sfloat.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_SF.h"
  13 #include "cl_F.h"
  14
  15 const cl_SF scale_float (const cl_SF& x, sintL delta)
  16 {
  17   // Methode:
  18   // x=0.0 -> x als Ergebnis
  19   // delta muß betragsmäßig <= SF_exp_high-SF_exp_low sein.
  20   // Neues SF mit um delta vergrößertem Exponenten bilden.
  21       // x entpacken:
  22       var cl_signean sign;
  23       var sintL exp;
  24       var uint32 mant;
  25       SF_decode(x, { return x; }, sign=,exp=,mant=);
  26       if (delta >= 0)
  27         // delta>=0
  28         { var uintL udelta = delta;
  29           if (udelta <= (uintL)(SF_exp_high-SF_exp_low))
  30             { exp = exp+udelta;
  31               return encode_SF(sign,exp,mant);
  32             }
  33             else
  34             { cl_error_floating_point_overflow(); }
  35         }
  36         else
  37         // delta<0
  38         { var uintL udelta = -delta;
  39           if (udelta <= (uintL)(SF_exp_high-SF_exp_low))
  40             { exp = exp-udelta;
  41               return encode_SF(sign,exp,mant);
  42             }
  43             else
  44             if (underflow_allowed())
  45               { cl_error_floating_point_underflow(); }
  46               else
  47               { return SF_0; }
  48         }
  49 }