There is something I don't understand with bigfloat precision. I have some
task where I need to compute exactly with a thousand significant digit. I
wrote a small benchmark, just by taking a power of two and its logarithm.
julia> with_bigfloat_precision(10_000) do
log(2,big(2)^100)
end
1.000000000000000033456773379254076470278870835616524546375 [long
string of random digits]
9219723375649773e+02 with 10000 bits of precision
I would have expected to obtain 1.000... with a far better precision (at
least 600 zeros after the point). What I am doing wrong here?