[关于极限]为什么求π的程序会出现空数据？

Jackxwb

变量上已经用了long类型了，从数值上看也没有超出范围。
原程序地址：http://www.cnblogs.com/PegasusWang/archive/2013/03/17/2965189.html
关于printf命令用了中文社区 奈何col 的程序
地址：http://www.arduino.cn/thread-8366-1-1.html

移植版(本人copy后略微修改过)

1. int a = 10000, b, c = 1000,  e, g, ls, ls2, cs;
   
2. long d,  f[1001];
   
3. int serial_putc( char c, struct __file * )
   
4. {
   
5.   Serial.write( c );
   
6.   return c;
   
7. }
   
8. void printf_begin(void)
   
9. {
   
10.   fdevopen( &serial_putc, 0 );
    
11. }
    
12. void setup() {
    
13.   // put your setup code here, to run once:
    
14.   int i;
    
15.   Serial.begin(9600);
    
16.   printf_begin();
    
17.   Serial.println("pi =");
    
18.   for (i = 0; i < c; i++)
    
19.     f[i] = a / 5;
    
20.   while (c != 0)
    
21.   {
    
22.     d = 0;
    
23.     g = c * 2;
    
24.     b = c;
    
25.     while (1)
    
26.     {
    
27.       d = d + f[b] * a;
    
28.       g = g - 1;
    
29.       f[b] = d % g;
    
30.       d = d / g;
    
31.       g = g - 1;
    
32.       b = b - 1;
    
33. 
    
34.       if (b == 0) break;
    
35.       d = d * b;
    
36.     }
    
37.     c = c - 14;
    
38.     cs++;
    
39.     printf("(%.4d)%d", e + d / a, cs);//此次加括号仅为了好看。cs本为循环次数
    
40.     if (e + d / a == 0) break;
    
41.     e = d % a;
    
42.   }
    
43.   Serial.println("");
    
44.   Serial.println("run:" + ls);
    
45. }
    
46. 
    
47. void loop() {
    
48.   // put your main code here, to run repeatedly:
    
49. 
    
50. }
    

复制代码

Jackxwb

本人使用的是mega 2560，内存已使用了51%


如果编译时提示“没有足够的内存。”时，请把 f 的数组及 c 的数组改小一点。（c 的值等于 f 的数组数减一）

Super169

有檢查過 b 的值有可能 < 0 嗎?

在 設定

1. f[b] = d % g;

复制代码

之前, 加入以下的檢測看看:

1. if (b < 0 )Serial.println("Error: b < 0");

复制代码

janeasy

楼上正解，是由于c-14 造成的b小于0，c的数值应该是14的整数倍，就不会有此问题了。

Jackxwb

janeasy 发表于 2016-3-7 16:29
楼上正解，是由于c-14 造成的b小于0，c的数值应该是14的整数倍，就不会有此问题了。

谢谢！
.
一开始是想用来测试一下性能的，由于考虑到变量有效范围太窄，目前已经放弃了。
这么久没回复不好意思~

Jackxwb

Super169 发表于 2016-3-6 00:25
有檢查過 b 的值有可能 < 0 嗎?

在 設定之前, 加入以下的檢測看看:

谢谢！
虽然目前已经放弃这个项目了，但还是要感谢一下。
这么久没回复不好意思~