diff options
Diffstat (limited to 'net/ipv4/tcp_cubic.c')
-rw-r--r-- | net/ipv4/tcp_cubic.c | 16 |
1 files changed, 5 insertions, 11 deletions
diff --git a/net/ipv4/tcp_cubic.c b/net/ipv4/tcp_cubic.c index 6f08adbda54e..0e6cdfeb207a 100644 --- a/net/ipv4/tcp_cubic.c +++ b/net/ipv4/tcp_cubic.c | |||
@@ -96,23 +96,17 @@ static void bictcp_init(struct sock *sk) | |||
96 | */ | 96 | */ |
97 | static u32 cubic_root(u64 a) | 97 | static u32 cubic_root(u64 a) |
98 | { | 98 | { |
99 | u32 x, x1; | 99 | u32 x; |
100 | 100 | ||
101 | /* Initial estimate is based on: | 101 | /* Initial estimate is based on: |
102 | * cbrt(x) = exp(log(x) / 3) | 102 | * cbrt(x) = exp(log(x) / 3) |
103 | */ | 103 | */ |
104 | x = 1u << (fls64(a)/3); | 104 | x = 1u << (fls64(a)/3); |
105 | 105 | ||
106 | /* | 106 | /* converges to 32 bits in 3 iterations */ |
107 | * Iteration based on: | 107 | x = (2 * x + (u32)div64_64(a, (u64)x*(u64)x)) / 3; |
108 | * 2 | 108 | x = (2 * x + (u32)div64_64(a, (u64)x*(u64)x)) / 3; |
109 | * x = ( 2 * x + a / x ) / 3 | 109 | x = (2 * x + (u32)div64_64(a, (u64)x*(u64)x)) / 3; |
110 | * k+1 k k | ||
111 | */ | ||
112 | do { | ||
113 | x1 = x; | ||
114 | x = (2 * x + (uint32_t) div64_64(a, x*x)) / 3; | ||
115 | } while (abs(x1 - x) > 1); | ||
116 | 110 | ||
117 | return x; | 111 | return x; |
118 | } | 112 | } |