| ... | ... | @@ -2678,7 +2678,9 @@ pub const Const = struct { |
| 2678 | 2678 | |
| 2679 | 2679 | /// Calculate the base 2 logarithm, rounded down. |
| 2680 | 2680 | pub fn log2(a: Const) Limb { |
| 2681 | | return a.limbs.len * @bitSizeOf(Limb) - 1 - @clz(a.limbs[a.limbs.len - 1]); |
| 2681 | assert(a.positive); |
| 2682 | assert(!a.eqlZero()); |
| 2683 | return a.bitCountAbs() - 1; |
| 2682 | 2684 | } |
| 2683 | 2685 | |
| 2684 | 2686 | /// Calculate the base 10 logarithm, rounded down. |
| ... | ... | @@ -2695,13 +2697,9 @@ pub const Const = struct { |
| 2695 | 2697 | /// |
| 2696 | 2698 | /// `limbs_buffer` is used for temporary storage. The amount required is given by `calcLog10LimbsBufferLen`. |
| 2697 | 2699 | pub fn log10(a: Const, limbs_buffer: []Limb) Limb { |
| 2698 | | const max_digits_per_limb = std.math.log10(std.math.maxInt(Limb)); |
| 2699 | | const limb_base = comptime calc: { |
| 2700 | | var limb_base: comptime_int = 1; |
| 2701 | | for (0..max_digits_per_limb) |_| limb_base *= 10; |
| 2702 | | break :calc limb_base; |
| 2703 | | }; |
| 2704 | | const limb_base_as_bigint: Const = .{ .limbs = &.{limb_base}, .positive = true }; |
| 2700 | assert(a.positive); |
| 2701 | assert(!a.eqlZero()); |
| 2702 | const limb_base_as_bigint: Const = .{ .limbs = &.{constants.big_bases[10]}, .positive = true }; |
| 2705 | 2703 | |
| 2706 | 2704 | var q: Mutable = .{ |
| 2707 | 2705 | .limbs = limbs_buffer[0 .. a.limbs.len + 2], |
| ... | ... | @@ -2721,7 +2719,7 @@ pub const Const = struct { |
| 2721 | 2719 | var num_digits: Limb = 0; |
| 2722 | 2720 | while (q.len >= 2) { |
| 2723 | 2721 | q.divTrunc(&remainder, q.toConst(), limb_base_as_bigint, division_buf); |
| 2724 | | num_digits += max_digits_per_limb; |
| 2722 | num_digits += constants.digits_per_limb[10]; |
| 2725 | 2723 | } |
| 2726 | 2724 | var remaining_limb = q.limbs[0]; |
| 2727 | 2725 | while (remaining_limb != 0) { |