src/real/algebraic/cl_RA_sqrt.cc

   1 // sqrt().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/real.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cln/rational.h"
  13 #include "cln/float.h"
  14
  15 namespace cln {
  16
  17 CL_INLINE const cl_R CL_INLINE_DECL(sqrt) (const cl_RA& x)
  18 {
  19         var cl_RA w;
  20         if (sqrtp(x,&w)) // auf Quadrat testen
  21                 return w; // war Quadrat, w ist die Wurzel
  22         else
  23                 // x in Float umwandeln, dann die Wurzel ziehen:
  24                 return sqrt(cl_float(x));
  25 }
  26
  27 }  // namespace cln