src/real/division/cl_R_trunc22.cc

   1 // truncate2().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_real.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_R.h"
  13 #include "cl_rational.h"
  14 #include "cl_R_div_t.h"
  15
  16 const cl_R_div_t truncate2 (const cl_R& x, const cl_R& y)
  17 {
  18 // Methode:
  19 // Beides rationale Zahlen -> truncate2(x,y).
  20 // Sonst: truncate2(x/y) -> (q,r). Liefere q und x-y*q=y*r.
  21         if (rationalp(x))
  22                 if (rationalp(y)) {
  23                         DeclareType(cl_RA,x);
  24                         DeclareType(cl_RA,y);
  25                         return truncate2(x,y);
  26                 }
  27         var cl_R_div_t q_r = truncate2(x/y);
  28         var cl_I& q = q_r.quotient;
  29         var cl_R& r = q_r.remainder;
  30         return cl_R_div_t(q,y*r);
  31 }