xref: /openbsd-src/gnu/llvm/compiler-rt/lib/builtins/popcountti2.c (revision 1f9cb04fc6f537ca6cf5a53c28927340cba218a2)
13cab2bb3Spatrick //===-- popcountti2.c - Implement __popcountti2
23cab2bb3Spatrick //----------------------------===//
33cab2bb3Spatrick //
43cab2bb3Spatrick // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
53cab2bb3Spatrick // See https://llvm.org/LICENSE.txt for license information.
63cab2bb3Spatrick // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
73cab2bb3Spatrick //
83cab2bb3Spatrick //===----------------------------------------------------------------------===//
93cab2bb3Spatrick //
103cab2bb3Spatrick // This file implements __popcountti2 for the compiler_rt library.
113cab2bb3Spatrick //
123cab2bb3Spatrick //===----------------------------------------------------------------------===//
133cab2bb3Spatrick 
143cab2bb3Spatrick #include "int_lib.h"
153cab2bb3Spatrick 
163cab2bb3Spatrick #ifdef CRT_HAS_128BIT
173cab2bb3Spatrick 
183cab2bb3Spatrick // Returns: count of 1 bits
193cab2bb3Spatrick 
__popcountti2(ti_int a)20*1f9cb04fSpatrick COMPILER_RT_ABI int __popcountti2(ti_int a) {
213cab2bb3Spatrick   tu_int x3 = (tu_int)a;
223cab2bb3Spatrick   x3 = x3 - ((x3 >> 1) &
233cab2bb3Spatrick              (((tu_int)0x5555555555555555uLL << 64) | 0x5555555555555555uLL));
243cab2bb3Spatrick   // Every 2 bits holds the sum of every pair of bits (64)
253cab2bb3Spatrick   x3 = ((x3 >> 2) &
263cab2bb3Spatrick         (((tu_int)0x3333333333333333uLL << 64) | 0x3333333333333333uLL)) +
273cab2bb3Spatrick        (x3 & (((tu_int)0x3333333333333333uLL << 64) | 0x3333333333333333uLL));
283cab2bb3Spatrick   // Every 4 bits holds the sum of every 4-set of bits (3 significant bits) (32)
293cab2bb3Spatrick   x3 = (x3 + (x3 >> 4)) &
303cab2bb3Spatrick        (((tu_int)0x0F0F0F0F0F0F0F0FuLL << 64) | 0x0F0F0F0F0F0F0F0FuLL);
313cab2bb3Spatrick   // Every 8 bits holds the sum of every 8-set of bits (4 significant bits) (16)
323cab2bb3Spatrick   du_int x2 = (du_int)(x3 + (x3 >> 64));
333cab2bb3Spatrick   // Every 8 bits holds the sum of every 8-set of bits (5 significant bits) (8)
343cab2bb3Spatrick   su_int x = (su_int)(x2 + (x2 >> 32));
353cab2bb3Spatrick   // Every 8 bits holds the sum of every 8-set of bits (6 significant bits) (4)
363cab2bb3Spatrick   x = x + (x >> 16);
373cab2bb3Spatrick   // Every 8 bits holds the sum of every 8-set of bits (7 significant bits) (2)
383cab2bb3Spatrick   //
393cab2bb3Spatrick   // Upper 16 bits are garbage
403cab2bb3Spatrick   return (x + (x >> 8)) & 0xFF; // (8 significant bits)
413cab2bb3Spatrick }
423cab2bb3Spatrick 
433cab2bb3Spatrick #endif // CRT_HAS_128BIT
44