src/float/dfloat/elem/cl_DF_scale.cc

   1 // scale_float().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_dfloat.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_DF.h"
  13 #include "cl_F.h"
  14
  15 const cl_DF scale_float (const cl_DF& x, sintL delta)
  16 {
  17   // Methode:
  18   // x=0.0 -> x als Ergebnis
  19   // delta muß betragsmäßig <= DF_exp_high-DF_exp_low sein.
  20   // Neues DF mit um delta vergrößertem Exponenten bilden.
  21       // x entpacken:
  22       var cl_signean sign;
  23       var sintL exp;
  24 #if (cl_word_size==64)
  25       var uint64 mant;
  26       DF_decode(x, { return x; }, sign=,exp=,mant=);
  27 #else
  28       var uint32 manthi;
  29       var uint32 mantlo;
  30       DF_decode2(x, { return x; }, sign=,exp=,manthi=,mantlo=);
  31 #endif
  32       if (delta >= 0)
  33         // delta>=0
  34         { var uintL udelta = delta;
  35           if (udelta <= (uintL)(DF_exp_high-DF_exp_low))
  36             { exp = exp+udelta;
  37 #if (cl_word_size==64)
  38               return encode_DF(sign,exp,mant);
  39 #else
  40               return encode_DF(sign,exp,manthi,mantlo);
  41 #endif
  42             }
  43             else
  44             { cl_error_floating_point_overflow(); }
  45         }
  46         else
  47         // delta<0
  48         { var uintL udelta = -delta;
  49           if (udelta <= (uintL)(DF_exp_high-DF_exp_low))
  50             { exp = exp-udelta;
  51 #if (cl_word_size==64)
  52               return encode_DF(sign,exp,mant);
  53 #else
  54               return encode_DF(sign,exp,manthi,mantlo);
  55 #endif
  56             }
  57             else
  58             if (underflow_allowed())
  59               { cl_error_floating_point_underflow(); }
  60               else
  61               { return cl_DF_0; }
  62         }
  63 }