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