#### /cln-1.3.2/src/integer/algebraic/cl_I_sqrt.cc

#
C++ | 43 lines | 30 code | 8 blank | 5 comment | 1 complexity | 9c7fb750967bff60aba495aac7d13ee1 MD5 | raw file
``````
// isqrt().

// General includes.
#include "base/cl_sysdep.h"

// Specification.
#include "cln/integer.h"

// Implementation.

#include "cln/number.h"
#include "cln/io.h"
#include "cln/integer_io.h"
#include "cln/exception.h"
#include "integer/cl_I.h"
#include "base/digitseq/cl_DS.h"
#include <sstream>

namespace cln {

bool isqrt (const cl_I& x, cl_I* w)
{
if (minusp(x)) {
std::ostringstream buf;
fprint(buf, "isqrt: applied to negative number: ");
fprint(buf, x);
throw runtime_exception(buf.str());
}
CL_ALLOCA_STACK;
var const uintD* x_MSDptr;
var uintC x_len;
var const uintD* x_LSDptr;
I_to_NDS_nocopy(x, x_MSDptr=,x_len=,x_LSDptr=,true,); // Digit sequence >=0 zu x
var DS y;
var bool squarep;
UDS_sqrt(x_MSDptr,x_len,x_LSDptr, &y, squarep=); // Wurzel ziehen
*w = NUDS_to_I(y.MSDptr,y.len); // als Integer
return squarep;
}
// Bit complexity (x of length N): O(M(N)).

}  // namespace cln

``````