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 MAYBE_INLINE
  18 const cl_R sqrt (const cl_RA& x)
  19 {
  20         var cl_RA w;
  21         if (sqrtp(x,&w)) // auf Quadrat testen
  22                 return w; // war Quadrat, w ist die Wurzel
  23         else
  24                 // x in Float umwandeln, dann die Wurzel ziehen:
  25                 return sqrt(cl_float(x));
  26 }
  27
  28 }  // namespace cln