如何从升压int256_t中获取以 2 为底的对数

Question

所以我想从int256_t取一个对数。

我发现了这个，但将其修改为大小int256_t不起作用。它会给出不正确的结果：https://stackoverflow.com/a/11376759

提升中是否有支持多精度的对数函数？

Answer 1

如果没有进一步的限制，我会写：

住在科里鲁

int256_t x("12345678901234567890");
std::cout << "log2(" << x << ") = " << log2(cpp_bin_float_100(x)) << "n";

指纹

log2(12345678901234567890) = 63.4206

简单

无论如何，如果您要对结果进行四舍五入，则不需要所有这些精度：

住在科里鲁

std::cout << "log2(" << x << ") = " << log2(x.convert_to<double>()) << "n";

临时

一个非常粗糙的形式可能看起来像：

住在科里鲁

template <typename T>
static inline int adhoc_log2(T i) {
    int n = 0;
    while(i /= 2) ++n;
    return n;
}

指纹

adhoc(12345678901234567890) = 63

最后，事实上，您可以回到位级别，并以@SamVarshavchik建议的方式应用位摆弄技巧

Answer 2

template<class T, std::size_t...Is>
std::array<uint8_t, bytes> decomponse( std::index_sequence<Is...>, T in ) {
  return {{ (uint8_t)( in >> (8*( (sizeof...(Is)) - Is ) )) }};
}
template<std::size_t bytes, class T>
std::array<int8_t, bytes> decomponse( T in ) {
  return decompose( in, std::make_index_sequence<bytes>{} );
}
template<std::size_t bytes, class T>
double big_log2( T in ) {
  auto bytes_array = decompose<bytes>(in);
  std::size_t zeros = 0;
  for (auto byte:bytes_array) {
    if (byte) break;
    ++zeros;
  }
  int32_t tmp = 0;
  for (auto i = zeros; (i < bytes) && i < (zeros+4); ++i) {
    tmp |= bytes_array[i] << (8*(3 - (i-zeros)));
  }
  return log2(tmp) + (bytes-zeros-4)*8;
}

其中log2(int32_t)生成 32 位值的 log2 并返回 double 。

根据字节序，反转bytes_array和修改算法可能会使其更快（理想情况下，它应该编译为内存或完全优化）。

这试图将 4 个最高有效性字节推入一个int32_t中，取其对数，然后添加正确的量以将其移动到所讨论的 N 字节整数的大小。