7 #include "cl_rational.h"
13 #include "cl_integer.h"
15 cl_boolean sqrtp (const cl_RA& x, cl_RA* w)
18 // Bei Integers: klar.
19 // Bei Brüchen a/b : muß a=c^2 und b=d^2 sein. Dann ist die Wurzel = c/d
20 // (mit ggT(c,d)=1 und d>1).
23 return sqrtp(x,(cl_I*)w);
27 var const cl_I& b = denominator(x);
29 if (!sqrtp(b,&d)) // Nenner auf Quadratzahl testen
31 var const cl_I& a = numerator(x);
33 if (!sqrtp(a,&c)) // Zähler auf Quadratzahl testen
35 // beides Quadratzahlen -> Quotient der Wurzeln bilden
36 *w = I_I_to_RT(c,d); return cl_true;