Based on a proposal by @gjlayde on https://www.mikrocontroller.net/topic/85256#6133370, speed of fp64_log could be improved by another 40%. Instead of approximating log(x) in the range [1;2[, it is faster to approximate it in [1/sqrt(2); sqrt(2)] with same accuracy.