remove unused code

Author: lperron

ortools/util/fp_utils.cc

@@ -25,7 +25,6 @@
 #include <vector>
 
 #include "absl/base/casts.h"
-#include "absl/base/internal/endian.h"
 #include "absl/log/check.h"
 #include "absl/types/span.h"
 #include "ortools/util/bitset.h"
@@ -230,33 +229,4 @@ int64_t ComputeGcdOfRoundedDoubles(absl::Span<const double> x,
   return gcd > 0 ? gcd : 1;
 }
 
-int fast_ilogb(double value) {
-  static_assert(CHAR_BIT == 8);
-  static_assert(sizeof(double) == 8);
-  // Get little-endian bit-representation of the floating point value.
-  const uint64_t bit_rep =
-      absl::little_endian::FromHost64(absl::bit_cast<uint64_t>(value));
-  return static_cast<int>((bit_rep >> 52) & 0x7FF) - 1023;
-}
-
-void fast_scalbn_inplace(double& mutable_value, int exponent) {
-  mutable_value = fast_scalbn(mutable_value, exponent);
-}
-
-double fast_scalbn(double value, int exponent) {
-  if (value == 0.0) return 0.0;
-  uint64_t bit_rep =
-      absl::little_endian::FromHost64(absl::bit_cast<uint64_t>(value));
-  // Binary representation is: (sign-bit)(11 exponent bits)(52 mantissa bits)
-  constexpr uint64_t kExponentMask(0x7FF0000000000000);
-  // This addition relies on the fact that signed numbers are written in
-  // two-s complement, and is correct as long as the sum does not
-  // overflow/underflow the result.
-  const uint64_t value_exponent =
-      (bit_rep + (static_cast<uint64_t>(exponent) << 52)) & kExponentMask;
-  bit_rep &= ~kExponentMask;
-  bit_rep |= value_exponent;
-  return absl::bit_cast<double>(absl::little_endian::ToHost64(bit_rep));
-}
-
 }  // namespace operations_research

ortools/util/fp_utils.h

@@ -253,19 +253,6 @@ inline FloatType Interpolate(FloatType x, FloatType y, FloatType alpha) {
   return alpha * x + (1 - alpha) * y;
 }
 
-// This is a fast implementation of the C99 function ilogb for normalized
-// doubles with the caveat that it returns -1023 for zero, and 1024 for infinity
-// an NaNs.
-int fast_ilogb(double value);
-
-// This is a fast implementation of the C99 function scalbn, with the caveat
-// that it works on normalized numbers and if the result underflows, overflows,
-// or is applied to a NaN or an +-infinity, the result is undefined behavior.
-// Note that the version of the function that takes a reference, modifies the
-// given value.
-double fast_scalbn(double value, int exponent);
-void fast_scalbn_inplace(double& mutable_value, int exponent);
-
 }  // namespace operations_research
 
 #endif  // OR_TOOLS_UTIL_FP_UTILS_H_