%Qf in gmp_printf (was: Re: Wrong division ?)
gmpdiscuss at formauri.es
Sat Apr 2 16:55:05 UTC 2016
Torbjörn Granlund wrote, On 2016-04-02 11:09:
> Wouldn't a prospective Qf in essence perform the same division that is
> done for mpf_set_q, and thus use the same time and memory?
I guess it depends on the implementation. A naïve algorithm to output it in decimal, simplified to always output a given number of decimals, is:
mpz_t div = num / den
mpz_t rem = num % den
do (decimals) times:
div = rem * 10 / den
rem = rem * 10 % den
An obvious way to improve it would be to use base 10^a for some adequate a, not necessarily 1. I can't tell if a method similar to the one used to output integers could be used here (I haven't checked that part of the code in many years, but I recall discussions of a recursive algorithm that I guess got implemented).
An advantage of this approach is that one doesn't need to estimate or calculate the minimum required precision in advance. To be honest, I haven't checked the process that would be followed when using mpf_set_q followed by gmp_printf("%Ff"), so I don't actually know whether there would be a gain. I presumed there would, but I might be wrong.
Note that this is not intended to address the OP's need for the result as a float; it's a different use case, namely the problem of outputting a fraction as a decimal number with the required number of decimals.
More information about the gmp-discuss