xref: /aosp_15_r20/external/skia/src/base/SkMathPriv.cpp (revision c8dee2aa9b3f27cf6c858bd81872bdeb2c07ed17)
1*c8dee2aaSAndroid Build Coastguard Worker /*
2*c8dee2aaSAndroid Build Coastguard Worker  * Copyright 2008 The Android Open Source Project
3*c8dee2aaSAndroid Build Coastguard Worker  *
4*c8dee2aaSAndroid Build Coastguard Worker  * Use of this source code is governed by a BSD-style license that can be
5*c8dee2aaSAndroid Build Coastguard Worker  * found in the LICENSE file.
6*c8dee2aaSAndroid Build Coastguard Worker  */
7*c8dee2aaSAndroid Build Coastguard Worker 
8*c8dee2aaSAndroid Build Coastguard Worker #include "src/base/SkMathPriv.h"
9*c8dee2aaSAndroid Build Coastguard Worker 
10*c8dee2aaSAndroid Build Coastguard Worker #include "include/private/base/SkAssert.h"
11*c8dee2aaSAndroid Build Coastguard Worker 
12*c8dee2aaSAndroid Build Coastguard Worker #include <cstdint>
13*c8dee2aaSAndroid Build Coastguard Worker 
14*c8dee2aaSAndroid Build Coastguard Worker ///////////////////////////////////////////////////////////////////////////////
15*c8dee2aaSAndroid Build Coastguard Worker 
16*c8dee2aaSAndroid Build Coastguard Worker /* www.worldserver.com/turk/computergraphics/FixedSqrt.pdf
17*c8dee2aaSAndroid Build Coastguard Worker */
SkSqrtBits(int32_t x,int count)18*c8dee2aaSAndroid Build Coastguard Worker int32_t SkSqrtBits(int32_t x, int count) {
19*c8dee2aaSAndroid Build Coastguard Worker     SkASSERT(x >= 0 && count > 0 && (unsigned)count <= 30);
20*c8dee2aaSAndroid Build Coastguard Worker 
21*c8dee2aaSAndroid Build Coastguard Worker     uint32_t    root = 0;
22*c8dee2aaSAndroid Build Coastguard Worker     uint32_t    remHi = 0;
23*c8dee2aaSAndroid Build Coastguard Worker     uint32_t    remLo = x;
24*c8dee2aaSAndroid Build Coastguard Worker 
25*c8dee2aaSAndroid Build Coastguard Worker     do {
26*c8dee2aaSAndroid Build Coastguard Worker         root <<= 1;
27*c8dee2aaSAndroid Build Coastguard Worker 
28*c8dee2aaSAndroid Build Coastguard Worker         remHi = (remHi<<2) | (remLo>>30);
29*c8dee2aaSAndroid Build Coastguard Worker         remLo <<= 2;
30*c8dee2aaSAndroid Build Coastguard Worker 
31*c8dee2aaSAndroid Build Coastguard Worker         uint32_t testDiv = (root << 1) + 1;
32*c8dee2aaSAndroid Build Coastguard Worker         if (remHi >= testDiv) {
33*c8dee2aaSAndroid Build Coastguard Worker             remHi -= testDiv;
34*c8dee2aaSAndroid Build Coastguard Worker             root++;
35*c8dee2aaSAndroid Build Coastguard Worker         }
36*c8dee2aaSAndroid Build Coastguard Worker     } while (--count >= 0);
37*c8dee2aaSAndroid Build Coastguard Worker 
38*c8dee2aaSAndroid Build Coastguard Worker     return root;
39*c8dee2aaSAndroid Build Coastguard Worker }
40*c8dee2aaSAndroid Build Coastguard Worker 
41*c8dee2aaSAndroid Build Coastguard Worker // Kernighan's method
SkPopCount_portable(uint32_t n)42*c8dee2aaSAndroid Build Coastguard Worker int SkPopCount_portable(uint32_t n) {
43*c8dee2aaSAndroid Build Coastguard Worker     int count = 0;
44*c8dee2aaSAndroid Build Coastguard Worker 
45*c8dee2aaSAndroid Build Coastguard Worker     while (n) {
46*c8dee2aaSAndroid Build Coastguard Worker         n &= (n - 1); // Remove the lowest bit in the integer.
47*c8dee2aaSAndroid Build Coastguard Worker         count++;
48*c8dee2aaSAndroid Build Coastguard Worker     }
49*c8dee2aaSAndroid Build Coastguard Worker     return count;
50*c8dee2aaSAndroid Build Coastguard Worker }
51*c8dee2aaSAndroid Build Coastguard Worker 
52*c8dee2aaSAndroid Build Coastguard Worker // Here we strip off the unwanted bits and then return the number of trailing zero bits
SkNthSet(uint32_t target,int n)53*c8dee2aaSAndroid Build Coastguard Worker int SkNthSet(uint32_t target, int n) {
54*c8dee2aaSAndroid Build Coastguard Worker     SkASSERT(n < SkPopCount(target));
55*c8dee2aaSAndroid Build Coastguard Worker 
56*c8dee2aaSAndroid Build Coastguard Worker     for (int i = 0; i < n; ++i) {
57*c8dee2aaSAndroid Build Coastguard Worker         target &= (target - 1); // Remove the lowest bit in the integer.
58*c8dee2aaSAndroid Build Coastguard Worker     }
59*c8dee2aaSAndroid Build Coastguard Worker 
60*c8dee2aaSAndroid Build Coastguard Worker     return SkCTZ(target);
61*c8dee2aaSAndroid Build Coastguard Worker }
62