X-Git-Url: https://ginac.de/CLN/cln.git//cln.git?a=blobdiff_plain;f=src%2Ffloat%2Ftranscendental%2Fcl_LF_cossin_aux.cc;h=86e630d300d75161ef7e75664adb7363bbec8d2a;hb=8b3d91dec77438c0fe679b10869ab29e6cdeba58;hp=943158be56257ef4e86c750fb3cd0cfbbc46537a;hpb=749cf017b2954ea256ce006b2d8b0b815c2ff131;p=cln.git diff --git a/src/float/transcendental/cl_LF_cossin_aux.cc b/src/float/transcendental/cl_LF_cossin_aux.cc index 943158b..86e630d 100644 --- a/src/float/transcendental/cl_LF_cossin_aux.cc +++ b/src/float/transcendental/cl_LF_cossin_aux.cc @@ -14,7 +14,7 @@ #include "cl_LF.h" #include "cln/integer.h" #include "cl_alloca.h" -#include "cln/abort.h" +#include "cln/exception.h" #undef floor #include @@ -30,7 +30,7 @@ const cl_LF_cos_sin_t cl_cossin_aux (const cl_I& p, uintE lq, uintC len) { { Mutable(cl_I,p); var uintE lp = integer_length(p); // now |p| < 2^lp. - if (!(lp <= lq)) cl_abort(); + if (!(lp <= lq)) throw runtime_exception(); lp = lq - lp; // now |p/2^lq| < 2^-lp. // Minimize lq (saves computation time). { @@ -112,7 +112,7 @@ const cl_LF_cos_sin_t cl_cossin_aux (const cl_I& p, uintE lq, uintC len) #else // TRIVIAL_SPEEDUP var cl_LF cossum = sqrt(cl_I_to_LF(1,actuallen) - square(sinsum)); #endif - return cl_LF_cos_sin_t(shorten(cossum,len),shorten(sinsum,len)); // verkürzen und fertig + return cl_LF_cos_sin_t(shorten(cossum,len),shorten(sinsum,len)); // verkürzen und fertig }} // Bit complexity (N = len, and if p has length O(log N) and ql = O(log N)): // O(log(N)*M(N)).