SuperPi@home 高精度计算内核测试程序

JUST · 发表于 2004-5-22 00:00:00

试了一下，还不错：）
.bat应该是@%java% numbercruncher.program13_3.PiBorwein %digits%
原来的少了%

支持！

水鬼 · 发表于 2004-5-23 00:00:00

竟然要这么大的安装空间啊，晕倒。

hackerboy · 发表于 2004-5-25 00:00:00

以下是引用水鬼在2004-5-23 20:36:20的发言：
竟然要这么大的安装空间啊，晕倒。

正在精简……

wenmao · 发表于 2004-7-7 00:00:00

提示: 作者被禁止或删除内容自动屏蔽

count · 发表于 2004-7-7 00:00:00

楼主的程序时间复杂度是 O(n^2) 的。这可太不好了。

count · 发表于 2004-7-9 00:00:00

更精确的程序时间复杂度应是 O（(n log n)^2） (log 以2为底)
这说明了一个严重的问题：
hackerboy的程序中的乘法函数是用的“经典”方法。
就是类似于以下代码（C语言, invBas是进位制的倒数）

 void Mul(long*a, long Asize, long*b, long Bsize,
 long*c, long Csize);
{ long i, j, carry=0;
 double tmp;
 Csize=Asize+Bsize;
 for (i =0; i<Bsize; i++)
 {
 for (j = 0; i<Asize; j++)
 {
 tmp=a[j]*b[i]+carry
 carry=(long)(tmp*invBas)
 c[i+j]+=(long)(tmp-carry*Bas)
 }
 }
 if c[Asize+Bsize]=0 then Csize--;
 return 0;
}

有一句话说的好：当n大时，O(n^2)的汇编程序不如O(n)的Basic程序。

因此我建议楼主改用 FFT 乘法。这样其时间复杂度可降为 O(n(logn)^2)。
比如说，原来计算 5000 位要5秒,计算 10000 位要23秒。
 设改进后计算 5000 位也要5秒,则计算 10000 位只要11秒.
数据更大时，差距将极其明显。

hackerboy · 发表于 2004-10-12 00:00:00

感谢count的提示。此程序确实没用FFT算法，将在下一版中加入。

satellite · 发表于 2004-10-17 00:00:00

下载不了，为什么呢？

hackerboy · 发表于 2004-10-18 00:00:00

确实下载不了了。你需要的话我可以发给你。也可以看这里，我把他嵌如网页http://hackerboygn.91x.net/Pi.htm

satellite · 发表于 2004-10-20 00:00:00

能发给我吗？[email protected] 谢谢

hackerboy · 发表于 2004-10-20 00:00:00

以下是引用satellite在2004-10-20 14:14:12的发言：
能发给我吗？[email protected]
谢谢

请收邮件。没发对的话请到
http://www.equn.com/ofs/job.php?action=download&fid=1&tid=16&download=1_16_1084537888.rar%2CSuperPi%40home.rar%2C0%2C0%2Czip
http://www.equn.com/ofs/read.php?fid=1&tid=16&fpage=2
下载，不过要先注册！

[此贴子已经被作者于2004-10-20 18:48:52编辑过]

		自动登录	找回密码
密码			新注册用户

wenmao wenmao 当前离线积分 6505 UID 382 在线时间小时最后登录 1970-1-1 头像被屏蔽	发表于 2004-7-7 00:00:00 \| 显示全部楼层提示: 作者被禁止或删除内容自动屏蔽
	[建议]上传图片前先进来了解下相关经验
	回复使用道具举报