| Index: net/quic/congestion_control/cube_root.cc
|
| diff --git a/net/quic/congestion_control/cube_root.cc b/net/quic/congestion_control/cube_root.cc
|
| deleted file mode 100644
|
| index 191cccc1f9124cfee9d081c08b16942a6428707c..0000000000000000000000000000000000000000
|
| --- a/net/quic/congestion_control/cube_root.cc
|
| +++ /dev/null
|
| @@ -1,86 +0,0 @@
|
| -// Copyright (c) 2013 The Chromium Authors. All rights reserved.
|
| -// Use of this source code is governed by a BSD-style license that can be
|
| -// found in the LICENSE file.
|
| -
|
| -#include "net/quic/congestion_control/cube_root.h"
|
| -
|
| -#include "base/logging.h"
|
| -
|
| -namespace {
|
| -
|
| -// Find last bit in a 64-bit word.
|
| -int FindMostSignificantBit(uint64 x) {
|
| - if (!x) {
|
| - return 0;
|
| - }
|
| - int r = 0;
|
| - if (x & 0xffffffff00000000ull) {
|
| - x >>= 32;
|
| - r += 32;
|
| - }
|
| - if (x & 0xffff0000u) {
|
| - x >>= 16;
|
| - r += 16;
|
| - }
|
| - if (x & 0xff00u) {
|
| - x >>= 8;
|
| - r += 8;
|
| - }
|
| - if (x & 0xf0u) {
|
| - x >>= 4;
|
| - r += 4;
|
| - }
|
| - if (x & 0xcu) {
|
| - x >>= 2;
|
| - r += 2;
|
| - }
|
| - if (x & 0x02u) {
|
| - x >>= 1;
|
| - r++;
|
| - }
|
| - if (x & 0x01u) {
|
| - r++;
|
| - }
|
| - return r;
|
| -}
|
| -
|
| -// 6 bits table [0..63]
|
| -const uint32 cube_root_table[] = {
|
| - 0, 54, 54, 54, 118, 118, 118, 118, 123, 129, 134, 138, 143, 147, 151,
|
| - 156, 157, 161, 164, 168, 170, 173, 176, 179, 181, 185, 187, 190, 192, 194,
|
| - 197, 199, 200, 202, 204, 206, 209, 211, 213, 215, 217, 219, 221, 222, 224,
|
| - 225, 227, 229, 231, 232, 234, 236, 237, 239, 240, 242, 244, 245, 246, 248,
|
| - 250, 251, 252, 254
|
| -};
|
| -} // namespace
|
| -
|
| -namespace net {
|
| -
|
| -// Calculate the cube root using a table lookup followed by one Newton-Raphson
|
| -// iteration.
|
| -uint32 CubeRoot::Root(uint64 a) {
|
| - uint32 msb = FindMostSignificantBit(a);
|
| - DCHECK_LE(msb, 64u);
|
| -
|
| - if (msb < 7) {
|
| - // MSB in our table.
|
| - return ((cube_root_table[a]) + 31) >> 6;
|
| - }
|
| - // MSB 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, ...
|
| - // cubic_shift 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, ...
|
| - uint32 cubic_shift = (msb - 4);
|
| - cubic_shift = ((cubic_shift * 342) >> 10); // Div by 3, biased high.
|
| -
|
| - // 4 to 6 bits accuracy depending on MSB.
|
| - uint64 root =
|
| - ((cube_root_table[a >> (cubic_shift * 3)] + 10) << cubic_shift) >> 6;
|
| -
|
| - // Make one Newton-Raphson iteration.
|
| - // Since x has an error (inaccuracy due to the use of fix point) we get a
|
| - // more accurate result by doing x * (x - 1) instead of x * x.
|
| - root = 2 * root + (a / (root * (root - 1)));
|
| - root = ((root * 341) >> 10); // Div by 3, biased low.
|
| - return static_cast<uint32>(root);
|
| -}
|
| -
|
| -} // namespace net
|
|
|