skill用蒙特卡罗法算圆周率（据说可以测性能）

chrysalis

来源于这里：
http://www.zhihu.com/question/27768721

提问的人表示，vba要570ms，而python要7秒。
我把vba的程序拿到我电脑上跑了下，大概410ms，我就估摸着用skill跑大概要多久，于是把这个算法移植到了skill上面。

代码：

1. 
   
2. 
   
3. axlCmdRegister("test3" 'test3)
   
4. 
   
5. defun(test3 ()
   
6.         prog((l_time x_rand1 x_rand2 x_radius x_loop x_count x_rand1Opp x_rand2Opp f_pi)
   
7.                 l_time = measureTime(
   
8.                         x_radius = 100
   
9.                         x_loop   = 2500000
   
10.                         x_count  = 0
    
11.                         for(i 1 x_loop
    
12.                                 x_rand1 = random(x_radius)
    
13.                                 x_rand2 = random(x_radius)
    
14.                                                        
    
15.                                 when((x_rand1 ** 2 + x_rand2 ** 2) <= (x_radius ** 2)
    
16.                                         x_count = x_count + 1
    
17.                                 );when
    
18.                                                        
    
19.                                 x_rand1Opp = x_radius - x_rand1
    
20.                                 x_rand2Opp = x_radius - x_rand2
    
21.                                 when((x_rand1Opp ** 2 + x_rand2Opp ** 2) <= (x_radius ** 2)
    
22.                                         x_count = x_count + 1
    
23.                                 );when
    
24.                         );for
    
25.                         f_pi = x_count / float(x_loop) * 2
    
26.                 )
    
27.                 printf("%f time:%f\n", f_pi nth(2 l_time))
    
28.         );prog
    
29. )
    
30. 
    
31. 
    
32. 
    
33. 
    
34. 
    
35. 
    
36.        

复制代码

实际跑了一下，大概要3秒左右，比我想象中快（我以为会和python差不多），看来skill优化还是挺好的。（或许是Lisp的优越性？）
-----

在自己的电脑上跑了下python代码，要4.24秒。这么一看skill也没比python快到哪儿去……