40# 没——问题
:) ,好久不见~~
mpfr里面有现成的
在线文档:
http://www.mpfr.org/mpfr-current/mpfr.html#Special-Functions
MPFR很好很强大
必装
GMP, MPFR, MPC,LiDIA必装
恩,编译很顺利~效果很强大~
学习学习
购买后怎么不让下载呢?还需要金币吗
log(x)=pi/2*(1/agm(1,10^-n)-1/agm(1,x*10^-n))
还是有点不明白。请教几点
10^-n的n是所需精度吗?用2^-n行么?
假如我需要的是固定精度的log(x),pi/2*(1/agm(1,10^-n)是不是个可以预先计算的常数?
1. 我认为如果有你欲精确到n bits,你能使用2^-n
2. 如果你总是需要精确到n位有效数字,你能够预先计算pi/2*(1/agm(1,10^-n)
突然又想到个问题:我欲精确到n bits,每次计算agm要多计算到几位,
才能保证开方累计的截断误差不至于影响最终结果。
和13#的结论不同,我实验得到
计算过程中至少需要保持 p=1.5 n 个bit小数