src/libcore/num/flt2dec/estimator.rs

   1 // Copyright 2015 The Rust Project Developers. See the COPYRIGHT
   2 // file at the top-level directory of this distribution and at
   3 // http://rust-lang.org/COPYRIGHT.
   4 //
   5 // Licensed under the Apache License, Version 2.0 <LICENSE-APACHE or
   6 // http://www.apache.org/licenses/LICENSE-2.0> or the MIT license
   7 // <LICENSE-MIT or http://opensource.org/licenses/MIT>, at your
   8 // option. This file may not be copied, modified, or distributed
   9 // except according to those terms.
  10
  11 //! The exponent estimator.
  12
  13 /// Finds `k_0` such that `10^(k_0-1) < mant * 2^exp <= 10^(k_0+1)`.
  14 ///
  15 /// This is used to approximate `k = ceil(log_10 (mant * 2^exp))`;
  16 /// the true `k` is either `k_0` or `k_0+1`.
  17 #[doc(hidden)]
  18 pub fn estimate_scaling_factor(mant: u64, exp: i16) -> i16 {
  19     // 2^(nbits-1) < mant <= 2^nbits if mant > 0
  20     let nbits = 64 - (mant - 1).leading_zeros() as i64;
  21     // 1292913986 = floor(2^32 * log_10 2)
  22     // therefore this always underestimates (or is exact), but not much.
  23     (((nbits + exp as i64) * 1292913986) >> 32) as i16
  24 }