原来WINDOWS下面的super pi写的很差

wolfop

不知道是不是为了压CPU的原因。
LINUX下面的这个代码，计算1M PI，在2.5G的ATHLON64上面只要8－9s。牛啊。

1. 
   
2. /*
   
3. * he Gauss AGM algorithm using Schonhage variation. [1] Jorg Arndt,
   
4. * Christoph Haenel, "Pi Unleashed", pp. 93, Springer, 2001. (C) 2003
   
5. * Hironobu SUZUKI, licensed by GPL2
   
6. */
   
7. #include <stdio.h>
   
8. #include <gmp.h>
   
9. #include <time.h>
   
10. 
    
11. #define LOG_TEN_TWO  3.32192809488736234789
    
12. #define bprec(n) (int)(((n+10)*LOG_TEN_TWO)+2)
    
13. 
    
14. int
    
15. main(int ac, char *av[])
    
16. {
    
17.     long int        k,
    
18.                     loopmax;
    
19.     mpf_t           A,
    
20.                     B,
    
21.                     a,
    
22.                     b,
    
23.                     s,
    
24.                     s_1,
    
25.                     t,
    
26.                     t1,
    
27.                     t2,
    
28.                     t3,
    
29.                     c2;
    
30.     long int        prec,
    
31.                     dprec;
    
32.     dprec = 1000000L;           /* decimal precision */
    
33.     prec = bprec(dprec);        /* binary precision (plus alpha) */
    
34.     mpf_set_default_prec(prec);
    
35.     loopmax = 21;
    
36. 
    
37.     mpf_init(A);                /* big A */
    
38.     mpf_init(B);                /* big B */
    
39.     mpf_init(a);                /* a */
    
40.     mpf_init(b);                /* b */
    
41.     mpf_init(s);                /* s(n) */
    
42.     mpf_init(s_1);              /* s(n-1) */
    
43.     mpf_init(t);                /* temporary */
    
44.     mpf_init(t1);               /* temporary */
    
45.     mpf_init(t2);               /* temporary */
    
46.     mpf_init(t3);               /* temporary */
    
47.     mpf_init(c2);               /* 2 constant */
    
48.     mpf_set_ui(A, 1);
    
49.     mpf_set_ui(a, 1);
    
50.     mpf_set_ui(t1, 1);
    
51.     mpf_div_ui(B, t1, 2);
    
52.     mpf_div_ui(s, t1, 2);
    
53.     mpf_set_ui(t1, 10);
    
54.     mpf_set_ui(c2, 2);
    
55.     for (k = 1; k <= loopmax; k++) {
    
56.         mpf_add(t1, A, B);      /* (A+B) */
    
57.         mpf_div_ui(t, t1, 4);   /* t = (A+B)/4 */
    
58.         mpf_sqrt(b, B);         /* b = sqrt(B) */
    
59.         mpf_add(t1, a, b);      /* (a+b) */
    
60.         mpf_div_ui(a, t1, 2);   /* a = (a+b)/2 */
    
61.         mpf_mul(A, a, a);       /* A = a * a */
    
62.         mpf_sub(t1, A, t);      /* (A-t) */
    
63.         mpf_mul_ui(B, t1, 2);   /* B = (A - t) * 2 */
    
64.         mpf_sub(t1, B, A);      /* (B-A) */
    
65.         mpf_pow_ui(t2, c2, k);  /* 2^k */
    
66.         mpf_mul(t3, t1, t2);    /* (B-A) * 2^k */
    
67.         mpf_set(s_1, s);
    
68.         mpf_add(s, s_1, t3);    /* s = s + (B-A) * 2^k */
    
69.     }
    
70.     mpf_add(t1, a, b);
    
71.     mpf_pow_ui(t2, t1, 2);
    
72.     mpf_mul_ui(t1, s, 2);
    
73.     mpf_div(A, t2, t1);
    
74.     mpf_out_str(stdout, 10, dprec + 10, A);
    
75.     exit(0);
    
76. }
    

复制代码

bGFuZ3Vl

呵呵，早有人把 GMP 连接到 OpenSSL 的 engine 里去