src/float/transcendental/cl_F_roundpi.cc

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