src/float/transcendental/cl_F_roundpi2.cc

   1 // cl_round_pi2().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_F_tran.h"
   8
   9
  10 // Implementation.
  11
  12 namespace cln {
  13
  14 const cl_F_div_t cl_round_pi2 (const cl_F& x)
  15 {
  16         if (float_exponent(x) < 0)
  17                 // Exponent <0 -> |x|<1/2 -> |x/(pi/2)| < 1/2, also Division unnötig
  18                 return cl_F_div_t(0,x); // Quotient 0, Rest x
  19         else
  20                 // x durch pi/2 (mit hinreichender Genauigkeit) dividieren
  21                 return round2(x,scale_float(pi(x),-1));
  22 }
  23
  24 }  // namespace cln