是否有更好的方法来提取使用Pari/GP的实数的数字?

Question

这是我当前的代码，但它很丑，我担心非常大或很小的数字可能出现的边缘情况。有更好的方法吗?

real_to_int(n)={
if(n==floor(n),return(floor(n)));   \ If "n" is a whole number we're done
my(v=Vec(strprintf("%g",n)));       \ Convert "n" to a zero-padded character vector
my(d=sum(i=1,#v,i*(v[i]==".")));    \ Find the decimal point
my(t=eval(concat(v[^d])));          \ Delete the decimal point and reconvert to a number
my(z=valuation(t,10));              \ Count trailing zeroes
t/=10^z;                            \ Get rid of trailing zeroes
return(t)
}

Answer 1

您可以将输入的实数分成整数和小数部分，而无需查找点。

real_to_int(n) = {
my(intpart=digits(floor(n)));
my(fracpartrev=fromdigits(eval(Vecrev(Str(n))[1..-(2+#intpart)])));
fromdigits(concat(intpart, Vecrev(digits(fracpartrev))))
};
real_to_int(123456789.123456789009876543210000)
> 12345678912345678900987654321

注意，digits和fromdigits的组合为您消除了数字列表中所有的前导零。

Answer 2

这个问题没有很好地定义，因为从实数(内部存储为二进制)到十进制字符串的转换可能需要四舍五入，并且如何完成这取决于许多因素，例如format默认值或当前的bitprecision。

可以得到t_REAL的内部二进制表示为m * 2^e，其中m和e都是整数。

install(mantissa2nr, GL);
real_to_int(n) =
{
e = exponent(n) + 1 - bitprecision(n);
[mantissa2nr(n, 0), e];
}
? [m, e] = real_to_int(Pi)
%1 = [267257146016241686964920093290467695825, -126]
? m * 1. * 2^e
%2 = 3.1415926535897932384626433832795028842

使用[m, e]，我们获得了数字的精确(有理数)内部表示，并且两者都是定义良好的，即独立于所有设置。m是请求的十进制值的二进制等价物。