src/rational/division/cl_RA_ceil2.cc

   1 // ceiling2().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_rational.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_RA.h"
  13 #include "cl_integer.h"
  14
  15 const cl_RA_div_t ceiling2 (const cl_RA& x)
  16 {
  17 // Methode:
  18 // x Integer -> (q,r) := (x,0)
  19 // x Ratio a/b ->
  20 //   (ceiling a b) liefert q und r.
  21 //   Liefere q und r/b (mit b>1 und ggT(r,b)=ggT(r+q*b,b)=ggT(a,b)=1).
  22         if (integerp(x)) {
  23                 DeclareType(cl_I,x);
  24                 // (q,r) := (x,0)
  25                 return cl_RA_div_t(x,0);
  26         } else {
  27                 DeclareType(cl_RT,x);
  28                 var const cl_I& a = numerator(x);
  29                 var const cl_I& b = denominator(x);
  30                 var cl_I_div_t q_r = ceiling2(a,b);
  31                 var cl_I& q = q_r.quotient;
  32                 var cl_I& r = q_r.remainder;
  33                 return cl_RA_div_t(q,I_I_to_RT(r,b));
  34         }
  35 }