src/rational/elem/cl_RA_square.cc

   1 // square().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/rational.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_RA.h"
  13 #include "cln/integer.h"
  14
  15 namespace cln {
  16
  17 const cl_RA square (const cl_RA& r)
  18 {
  19 // Methode:
  20 // r Integer -> klar.
  21 // r = a/b -> Ergebnis a^2/b^2
  22         if (integerp(r)) {
  23                 DeclareType(cl_I,r);
  24                 return square(r);
  25         } else {
  26                 DeclareType(cl_RT,r);
  27                 var const cl_I& a = numerator(r);
  28                 var const cl_I& b = denominator(r);
  29                 // Immer noch b^2>1 und ggT(a^2,b^2) = ggT(a,b)^2 = 1
  30                 return I_I_to_RT(square(a),square(b));
  31         }
  32 }
  33
  34 }  // namespace cln