media: vimc: constify structures stored in fields of v4l2_subdev_ops structure
[linux-2.6-block.git] / lib / int_sqrt.c
CommitLineData
b2441318 1// SPDX-License-Identifier: GPL-2.0
30493cc9
DB
2/*
3 * Copyright (C) 2013 Davidlohr Bueso <davidlohr.bueso@hp.com>
4 *
5 * Based on the shift-and-subtract algorithm for computing integer
6 * square root from Guy L. Steele.
7 */
1da177e4
LT
8
9#include <linux/kernel.h>
8bc3bcc9 10#include <linux/export.h>
f8ae107e 11#include <linux/bitops.h>
1da177e4
LT
12
13/**
e813a614 14 * int_sqrt - computes the integer square root
1da177e4
LT
15 * @x: integer of which to calculate the sqrt
16 *
e813a614 17 * Computes: floor(sqrt(x))
1da177e4
LT
18 */
19unsigned long int_sqrt(unsigned long x)
20{
30493cc9 21 unsigned long b, m, y = 0;
1da177e4 22
30493cc9
DB
23 if (x <= 1)
24 return x;
1da177e4 25
f8ae107e 26 m = 1UL << (__fls(x) & ~1UL);
30493cc9
DB
27 while (m != 0) {
28 b = y + m;
29 y >>= 1;
1da177e4 30
30493cc9
DB
31 if (x >= b) {
32 x -= b;
33 y += m;
1da177e4 34 }
30493cc9 35 m >>= 2;
1da177e4 36 }
30493cc9
DB
37
38 return y;
1da177e4
LT
39}
40EXPORT_SYMBOL(int_sqrt);
47a36163
CM
41
42#if BITS_PER_LONG < 64
43/**
44 * int_sqrt64 - strongly typed int_sqrt function when minimum 64 bit input
45 * is expected.
46 * @x: 64bit integer of which to calculate the sqrt
47 */
48u32 int_sqrt64(u64 x)
49{
50 u64 b, m, y = 0;
51
52 if (x <= ULONG_MAX)
53 return int_sqrt((unsigned long) x);
54
55 m = 1ULL << (fls64(x) & ~1ULL);
56 while (m != 0) {
57 b = y + m;
58 y >>= 1;
59
60 if (x >= b) {
61 x -= b;
62 y += m;
63 }
64 m >>= 2;
65 }
66
67 return y;
68}
69EXPORT_SYMBOL(int_sqrt64);
70#endif