blob: 5e8a62f075eb9c7a7039d70665ca6a1bb784628c [file] [log] [blame]
Edward O'Callaghandabf71f2009-08-05 19:06:50 +00001/* ===-- popcountdi2.c - Implement __popcountdi2 ----------------------------===
2 *
3 * The LLVM Compiler Infrastructure
4 *
Howard Hinnant5b791f62010-11-16 22:13:33 +00005 * This file is dual licensed under the MIT and the University of Illinois Open
6 * Source Licenses. See LICENSE.TXT for details.
Edward O'Callaghandabf71f2009-08-05 19:06:50 +00007 *
8 * ===----------------------------------------------------------------------===
9 *
10 * This file implements __popcountdi2 for the compiler_rt library.
11 *
12 * ===----------------------------------------------------------------------===
13 */
Daniel Dunbarfd089992009-06-26 16:47:03 +000014
15#include "int_lib.h"
16
Edward O'Callaghandabf71f2009-08-05 19:06:50 +000017/* Returns: count of 1 bits */
Daniel Dunbarfd089992009-06-26 16:47:03 +000018
Anton Korobeynikove63da932011-04-19 17:52:09 +000019COMPILER_RT_ABI si_int
Daniel Dunbarfd089992009-06-26 16:47:03 +000020__popcountdi2(di_int a)
21{
22 du_int x2 = (du_int)a;
23 x2 = x2 - ((x2 >> 1) & 0x5555555555555555uLL);
Edward O'Callaghandabf71f2009-08-05 19:06:50 +000024 /* Every 2 bits holds the sum of every pair of bits (32) */
Daniel Dunbarfd089992009-06-26 16:47:03 +000025 x2 = ((x2 >> 2) & 0x3333333333333333uLL) + (x2 & 0x3333333333333333uLL);
Edward O'Callaghandabf71f2009-08-05 19:06:50 +000026 /* Every 4 bits holds the sum of every 4-set of bits (3 significant bits) (16) */
Daniel Dunbarfd089992009-06-26 16:47:03 +000027 x2 = (x2 + (x2 >> 4)) & 0x0F0F0F0F0F0F0F0FuLL;
Edward O'Callaghandabf71f2009-08-05 19:06:50 +000028 /* Every 8 bits holds the sum of every 8-set of bits (4 significant bits) (8) */
Daniel Dunbarfd089992009-06-26 16:47:03 +000029 su_int x = (su_int)(x2 + (x2 >> 32));
Edward O'Callaghandabf71f2009-08-05 19:06:50 +000030 /* The lower 32 bits hold four 16 bit sums (5 significant bits). */
31 /* Upper 32 bits are garbage */
Daniel Dunbarfd089992009-06-26 16:47:03 +000032 x = x + (x >> 16);
Edward O'Callaghandabf71f2009-08-05 19:06:50 +000033 /* The lower 16 bits hold two 32 bit sums (6 significant bits). */
34 /* Upper 16 bits are garbage */
35 return (x + (x >> 8)) & 0x0000007F; /* (7 significant bits) */
Daniel Dunbarfd089992009-06-26 16:47:03 +000036}