7 #include "cln/integer.h"
17 uintL power2p (const cl_I& x) // x > 0
19 // Methode 1: Wenn ord2(x) = integer_length(x)-1.
20 // Methode 2: Wenn logand(x,x-1) = 0.
21 // Methode 3: Wenn das erste Digit /=0 eine Zweierpotenz ist und alle weiteren
24 { var uintL x_ = FN_to_UL(x);
25 if (!((x_ & (x_-1)) == 0)) return 0; // keine Zweierpotenz
26 integerlength32(x_,return); // Zweierpotenz: n = integer_length(x)
29 { var const uintD* MSDptr;
31 var const uintD* LSDptr;
32 BN_to_NDS_nocopy(x, MSDptr=,len=,LSDptr=); // normalisierte DS zu x bilden.
33 var uintD msd = mspref(MSDptr,0);
34 if (msd==0) { msshrink(MSDptr); msd = mspref(MSDptr,0); len--; }
35 // len = Anzahl der Digits ab MSDptr, len>0, msd = erstes Digit (/=0)
36 if (!((msd & (msd-1)) == 0)) return 0; // erstes Digit muß Zweierpotenz sein
37 if (DS_test_loop(MSDptr mspop 1,len-1,LSDptr)) return 0; // danach alles Nullen
39 integerlengthD(msd, msdlen=);
40 return intDsize*(uintL)(len-1) + msdlen; // integer_length(x) als Ergebnis