src/integer/bitwise/cl_I_logcount.cc

   1 // logcount().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_integer.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_I.h"
  13 #include "cl_DS.h"
  14 #include "cl_D.h"
  15 #include "cl_low.h"
  16
  17 uintL logcount (const cl_I& x)
  18 {
  19         if (fixnump(x))
  20           { var uint32 x32 = FN_to_L(x); // x als 32-Bit-Zahl
  21             if (FN_L_minusp(x,(sint32)x32)) { x32 = ~ x32; } // falls <0, komplementieren
  22             logcount_32(); // Bits von x32 zählen
  23             return x32;
  24           }
  25           else
  26           { var const uintD* MSDptr;
  27             var uintC len;
  28             BN_to_NDS_nocopy(x, MSDptr=,len=,); // DS zu x bilden, len>0.
  29             var uintL bitcount = 0; // Bitzähler
  30             var const uintD* ptr = MSDptr; // läuft durch die Digits durch
  31             var uintD sign = sign_of_sintD(mspref(ptr,0)); // Vorzeichen
  32             dotimespC(len,len,
  33               { bitcount += (uintL)logcountD(msprefnext(ptr) ^ sign); });
  34             // 0 <= bitcount < intDsize*2^intCsize.
  35             return bitcount;
  36           }
  37 }