src/float/lfloat/elem/cl_LF_scale.cc

   1 // scale_float().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_lfloat.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_LF.h"
  13 #include "cl_LF_impl.h"
  14 #include "cl_F.h"
  15
  16 const cl_LF scale_float (const cl_LF& x, sintL delta)
  17 {
  18   // Methode:
  19   // delta=0 -> x als Ergebnis
  20   // x=0.0 -> x als Ergebnis
  21   // delta muß ein Integer betragsmäßig <= LF_exp_high-LF_exp_low sein.
  22   // Neues LF mit um delta vergrößertem Exponenten bilden.
  23       if (delta == 0) { return x; } // delta=0 -> x als Ergebnis
  24       var uintL uexp = TheLfloat(x)->expo;
  25       if (uexp==0) { return x; }
  26       var uintL udelta = delta;
  27       if (delta >= 0) {
  28         // udelta = delta >=0
  29         if (   ((uexp = uexp+udelta) < udelta) // Exponent-Überlauf?
  30             || (uexp > LF_exp_high) // oder Exponent zu groß?
  31            )
  32           { cl_error_floating_point_overflow(); }
  33       } else {
  34         // delta <0, udelta = 2^32+delta
  35         if (   ((uexp = uexp+udelta) >= udelta) // oder Exponent-Unterlauf?
  36             || (uexp < LF_exp_low) // oder Exponent zu klein?
  37            )
  38           { cl_error_floating_point_underflow(); }
  39       }
  40       var uintC len = TheLfloat(x)->len;
  41       return encode_LFu(TheLfloat(x)->sign,uexp,arrayMSDptr(TheLfloat(x)->data,len),len);
  42 }