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 "cln/sfloat.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_SF.h"
  13 #include "cl_F.h"
  14
  15 namespace cln {
  16
  17 const cl_SF scale_float (const cl_SF& x, sintC delta)
  18 {
  19   // Methode:
  20   // x=0.0 -> x als Ergebnis
  21   // delta muß betragsmäßig <= SF_exp_high-SF_exp_low sein.
  22   // Neues SF mit um delta vergrößertem Exponenten bilden.
  23       // x entpacken:
  24       var cl_signean sign;
  25       var sintL exp;
  26       var uint32 mant;
  27       SF_decode(x, { return x; }, sign=,exp=,mant=);
  28       if (delta >= 0)
  29         // delta>=0
  30         { var uintC udelta = delta;
  31           if (udelta <= (uintL)(SF_exp_high-SF_exp_low))
  32             { exp = exp+udelta;
  33               return encode_SF(sign,exp,mant);
  34             }
  35             else
  36             { throw floating_point_overflow_exception(); }
  37         }
  38         else
  39         // delta<0
  40         { var uintC udelta = -delta;
  41           if (udelta <= (uintL)(SF_exp_high-SF_exp_low))
  42             { exp = exp-udelta;
  43               return encode_SF(sign,exp,mant);
  44             }
  45             else
  46             if (underflow_allowed())
  47               { throw floating_point_underflow_exception(); }
  48               else
  49               { return SF_0; }
  50         }
  51 }
  52
  53 }  // namespace cln