[BACK]Return to popcount.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / mpn / generic

Diff for /OpenXM_contrib/gmp/mpn/generic/Attic/popcount.c between version 1.1.1.2 and 1.1.1.3

version 1.1.1.2, 2000/09/09 14:12:26 version 1.1.1.3, 2003/08/25 16:06:20
Line 1 
Line 1 
 /* popcount.c  /* mpn_popcount -- mpn bit population count.
   
 Copyright (C) 1994, 1996, 2000 Free Software Foundation, Inc.  Copyright 1994, 1996, 2000, 2001 Free Software Foundation, Inc.
   
 This file is part of the GNU MP Library.  This file is part of the GNU MP Library.
   
Line 22  MA 02111-1307, USA. */
Line 22  MA 02111-1307, USA. */
 #include "gmp.h"  #include "gmp.h"
 #include "gmp-impl.h"  #include "gmp-impl.h"
   
 #if defined __GNUC__  
 /* No processor claiming to be SPARC v9 compliant seem to  
    implement the POPC instruction.  Disable pattern for now.  */  
 #if 0 && defined __sparc_v9__ && BITS_PER_MP_LIMB == 64  
 #define popc_limb(a) \  
   ({                                                                    \  
     DItype __res;                                                       \  
     asm ("popc %1,%0" : "=r" (__res) : "rI" (a));                       \  
     __res;                                                              \  
   })  
 #endif  
 #endif  
   
 #ifndef popc_limb  unsigned long
   mpn_popcount (mp_srcptr p, mp_size_t size)
 /* Cool population count of a mp_limb_t.  
    You have to figure out how this works, I won't tell you!  */  
   
 static inline unsigned int  
 #if __STDC__  
 popc_limb (mp_limb_t x)  
 #else  
 popc_limb (x)  
      mp_limb_t x;  
 #endif  
 {  {
 #if BITS_PER_MP_LIMB == 64    unsigned long  result = 0;
   /* We have to go into some trouble to define these constants.    unsigned long  n;
      (For mp_limb_t being `long long'.)  */  
   mp_limb_t cnst;  
   cnst = 0xaaaaaaaaL | ((mp_limb_t) 0xaaaaaaaaL << BITS_PER_MP_LIMB/2);  
   x -= (x & cnst) >> 1;  
   cnst = 0x33333333L | ((mp_limb_t) 0x33333333L << BITS_PER_MP_LIMB/2);  
   x = ((x & ~cnst) >> 2) + (x & cnst);  
   cnst = 0x0f0f0f0fL | ((mp_limb_t) 0x0f0f0f0fL << BITS_PER_MP_LIMB/2);  
   x = ((x >> 4) + x) & cnst;  
   x = ((x >> 8) + x);  
   x = ((x >> 16) + x);  
   x = ((x >> 32) + x) & 0xff;  
 #endif  
 #if BITS_PER_MP_LIMB == 32  
   x -= (x & 0xaaaaaaaa) >> 1;  
   x = ((x >> 2) & 0x33333333L) + (x & 0x33333333L);  
   x = ((x >> 4) + x) & 0x0f0f0f0fL;  
   x = ((x >> 8) + x);  
   x = ((x >> 16) + x) & 0xff;  
 #endif  
   return x;  
 }  
 #endif  
   
 unsigned long int    ASSERT (size >= 1);
 #if __STDC__  
 mpn_popcount (register mp_srcptr p, register mp_size_t size)  
 #else  
 mpn_popcount (p, size)  
      register mp_srcptr p;  
      register mp_size_t size;  
 #endif  
 {  
   unsigned long int popcnt;  
   mp_size_t i;  
   
   popcnt = 0;    do {
   for (i = 0; i < size; i++)      popc_limb (n, *p++);
     popcnt += popc_limb (p[i]);      result += n;
     } while (--size != 0);
   
   return popcnt;    return result;
 }  }
Legend:
Removed from v.1.1.1.2  
changed lines
  Added in v.1.1.1.3
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>