wayne
发表于 2010-8-11 09:36:16
40# 没——问题
:) ,好久不见~~
mpfr里面有现成的
wayne
发表于 2010-8-11 09:39:00
在线文档:
http://www.mpfr.org/mpfr-current/mpfr.html#Special-Functions
无心人
发表于 2010-8-11 10:16:28
MPFR很好很强大
必装
GMP, MPFR, MPC,LiDIA必装
没——问题
发表于 2010-8-14 12:36:32
恩,编译很顺利~效果很强大~
forlorm
发表于 2011-1-23 12:45:22
学习学习
forlorm
发表于 2011-1-23 12:50:33
购买后怎么不让下载呢?还需要金币吗
zeroieme
发表于 2011-3-28 22:56:36
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)是不是个可以预先计算的常数?
liangbch
发表于 2011-3-29 16:18:46
1. 我认为如果有你欲精确到n bits,你能使用2^-n
2. 如果你总是需要精确到n位有效数字,你能够预先计算pi/2*(1/agm(1,10^-n)
zeroieme
发表于 2011-7-1 08:24:54
突然又想到个问题:我欲精确到n bits,每次计算agm要多计算到几位,
才能保证开方累计的截断误差不至于影响最终结果。
zeroieme
发表于 2011-7-8 17:39:59
和13#的结论不同,我实验得到
计算过程中至少需要保持 p=1.5 n 个bit小数