PageRenderTime 40ms CodeModel.GetById 18ms RepoModel.GetById 0ms app.codeStats 0ms

/cln-1.3.2/src/rational/division/cl_RA_floor2.cc

#
C++ | 39 lines | 21 code | 8 blank | 10 comment | 2 complexity | 0dcdc0bc0116601e4f8f64fc15391bbd MD5 | raw file
Possible License(s): GPL-2.0
  1. // floor2().
  2. // General includes.
  3. #include "base/cl_sysdep.h"
  4. // Specification.
  5. #include "cln/rational.h"
  6. // Implementation.
  7. #include "rational/cl_RA.h"
  8. #include "cln/integer.h"
  9. namespace cln {
  10. const cl_RA_div_t floor2 (const cl_RA& x)
  11. {
  12. // Methode:
  13. // x Integer -> (q,r) := (x,0)
  14. // x Ratio a/b ->
  15. // (floor a b) liefert q und r.
  16. // Liefere q und r/b (mit b>1 und ggT(r,b)=ggT(r+q*b,b)=ggT(a,b)=1).
  17. if (integerp(x)) {
  18. DeclareType(cl_I,x);
  19. // (q,r) := (x,0)
  20. return cl_RA_div_t(x,0);
  21. } else {
  22. DeclareType(cl_RT,x);
  23. var const cl_I& a = numerator(x);
  24. var const cl_I& b = denominator(x);
  25. var cl_I_div_t q_r = floor2(a,b);
  26. var cl_I& q = q_r.quotient;
  27. var cl_I& r = q_r.remainder;
  28. return cl_RA_div_t(q,I_I_to_RT(r,b));
  29. }
  30. }
  31. } // namespace cln