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 | } |
