src/rational/elem/cl_RA_square.cc

   1 // square().
   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 square (const cl_RA& r)
  16 {
  17 // Methode:
  18 // r Integer -> klar.
  19 // r = a/b -> Ergebnis a^2/b^2
  20         if (integerp(r)) {
  21                 DeclareType(cl_I,r);
  22                 return square(r);
  23         } else {
  24                 DeclareType(cl_RT,r);
  25                 var const cl_I& a = numerator(r);
  26                 var const cl_I& b = denominator(r);
  27                 // Immer noch b^2>1 und ggT(a^2,b^2) = ggT(a,b)^2 = 1
  28                 return I_I_to_RT(square(a),square(b));
  29         }
  30 }