Forcal运行速度测试例子

[复制链接]

字体大小: 正常放大

forcal

45 主题	3 听众	282 积分

升级 91%

TA的每日心情

	难过 2012-8-27 18:22

签到天数: 1 天

[LV.1]初来乍到

电梯直达

1^#

发表于 2010-9-27 19:28 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

Forcal的运行效率与Lua相当，是目前最快的脚本之一，但Forcal的扩充性能更好。
在大量函数调用时，Forcal更具速度优势。
以下是Forcal的运行速度测试例子，大家可转换为其他脚本代码（例如matlab）比较一下。

例1：Fibonacci递归程序（n取40）

SetRealStackMax(1000);
# L: @+ y9 @% g- o! R
F(n)= which{
$ b9 g/ f& R. S S; R
n == 0,' a- M+ R! ?$ d/ T! u
return(0),& v; ^5 q3 J! @7 W
n == 1,
: J1 ^/ S\" ?0 m( s0 W+ k6 H
return(1),
: B- ~& ?! G- d! q; G& u3 U
return [F(n - 1) + F(n - 2)]. Y1 B6 K8 P- j5 x- T$ X
};3 |: Y4 W3 C\" o0 S
main(:t,n)=7 `. k: m3 c/ z& l
{
N3 O- a% @0 u
t=sys::clock(),
9 h5 n' C5 p, W1 P% G' y
n=F(40),
, V5 j, Q$ B t; j% k
t=sys::clock()-t,
# a\" u, f0 L- ]9 m7 z4 ~
printff{"\r\nfibonacci={1,i}, fibonacci_Time={2,i}毫秒={3,i}秒\r\n",n,t,t/1000}' O' m. @) [6 q5 v/ ~
};
( q; o, g2 M. B8 k

复制代码

结果：

fibonacci=102334155, fibonacci_Time=41859毫秒=41秒

例2：八皇后问题

// 在运行不同的程序时，Forcal的速度，从接近C++到只有C++速度的几十分之一。 M2 J4 ~2 `. Y1 o
// Forcal的建议是：对运行时间较长的程序，如确有必要，设计成二级函数由Forcal调用，从而获得接近C++速度的性能。
% S9 ~4 k$ J! x' ]- f
// Forcal与C++是无缝链接的。故C++能实现的功能，借助二级函数，Forcal完全可以实现。
7 }4 \6 S- _+ n2 _
// 但没有Forcal支持的C++程序，将无法获得高效率地实时编译计算字符串表达式的功能。6 X$ B$ x6 y c; x- R( Z
// 据测定，以下八皇后问题，Forcal的运行速度约为C++的1/10。
/ o) f* O: ^7 p4 |0 o
// 八皇后问题是一个古老而著名的问题，是回溯算法的典型例题。
0 x: o+ T, c3 j- t0 c$ L' n% g( o4 }
// 该问题是19世纪著名的数学家高斯1850年提出：在8×8格的国际象棋盘上摆放8个皇后，使其不能互相攻击，即任意两个皇后都不能处于同一行、同一列或同一斜线上，问有多少种摆法。$ I6 e. B1 v; M* O. x1 c0 V* ]
// 高斯认为有76种方案。1854年在柏林的象棋杂志上不同的作者发表了40种不同的解，后来有人用图论的方法解出92种结果。$ Z3 z( `3 U6 p/ O% x
// 以下算法是从网上搜来的，该算法没有最终给出排列组合，仅仅给出有多少种组合，但是算法确实十分奥妙。
2 v: ~: o) c4 U) O! {2 U; q
//Forcal源程序' Z( z5 n j/ H0 b4 e: U R
i:(::sum,upperlim)= sum=0,upperlim=1,SetIntStackMax(1000);- T, s\" K6 g3 K! j, c/ X' _4 J
i:test(row, ld, rd : pos,p : sum,upperlim)=
) d$ t& Y% S D( z2 v7 G
{
& O2 ?9 [. s! G. i\" f% O1 g
which
! E\" u8 A, H) Q- o
{ row != upperlim,\" l( H1 V. `$ j5 f\" z3 o' u
{ pos = and{upperlim , not[row.or(ld).or(rd)]},5 t9 c8 m5 c/ I7 s2 x( o# h+ r
while{ pos,
9 X; \/ l5 j5 M' Z0 I* z, I% `
p = and(pos,-pos),
* Q7 j/ R1 k9 Y+ k\" ?1 d
pos = pos -p,0 f- `. k2 Y0 P9 H
test(row+p, shl(ld+p,1), shr(rd+p,1))4 |6 J3 V8 \! y
}
9 q3 G& p5 V- ~- n6 V
},
! x: P# O5 v: Y# l
sum++0 p3 j$ x0 l\" u7 a# B( }- ]
}: F- N* C, m# ~6 i5 H\" U/ I
};
# K2 @* x7 f- M* _2 t, E ^
i:main(:tm,n:sum,upperlim)=
\" h/ M. h+ r0 T/ I0 t$ i
{& K# e1 i; L\" j& _6 ^7 J8 C
n=15,
+ Y; l/ x, |# {& w9 N F) Y
tm=sys::clock(),% D Q. e* P0 {* V* Z7 H, a/ ^4 G' \
printff("Queens:{1,i}, ",n),\" a8 i) ?& l6 k2 ?
upperlim=shl(upperlim,n)-1,
; x9 v5 `0 _* J: e
test(0,0,0),
( \ g$ m: {' F* I( `# G$ g5 S
printff("sum:{1,i}, {2,i}毫秒.\r\n",sum,sys::clock()-tm)4 Y\" q( i, Q; D/ A2 A
};) F7 Y' N; v2 D) {

复制代码

Forcal运行结果：

Queens:15, sum:2279184, 59547毫秒.

例3：一段有趣的程序：此例Forcal速度约为C/C++的1/4左右

mandelbrot(x, y : cr,ci,zi,zr,i,temp,zr2,zi2) =
+ u/ Q( J1 Q+ ` L, D: r
{
8 E9 e- |\" G2 L3 Q
cr = y - 0.5,
) o- _. F2 m8 Q6 O i) P$ l' }- m
ci = x,
' E; F9 N8 I. ^8 }$ q
zi = 0.0,9 O5 i Y- o) d0 X: S/ N
zr = 0.0,
) [6 {* p. g' O- g
i = 0,2 D0 W C: ]& W+ k3 y
(1).while {6 ~5 V! x2 E, i6 `
i ++,. a/ D7 s% H2 a' ]
temp = zr * zi,
9 u, ^( d3 ` ]+ ^
zr2 = zr * zr,
0 t$ l, _& D' Y# H
zi2 = zi * zi,
& N5 j; U+ ?# {- F* ^2 ?9 q, K
zr = zr2 - zi2 + cr,0 k/ T# j4 s\" B+ x, l
zi = temp + temp + ci,$ t\" M. H6 a! P8 L& D
if [zi2 + zr2 > 16, return (i)],
/ Q! q' T& L, j8 d5 V! E
if [i > 1000, return (0)]
L8 _8 Y% m0 {3 \6 s' |+ H
}0 E8 V7 I; o$ m\" n! u( _% B. ~
: n3 v2 Z: l3 P
};
( N\" J( y# R' y% x# D
main (:i,x,y,old,now) = {8 N4 h! M7 D' V
old=sys::clock(),- J$ J9 i3 }! j
y = -39,; q& {: A7 j8 \- i, ]/ f& K8 R+ O
while{ y < 39,8 C$ i+ o1 p) r\" ]9 n7 S. o: {
printff("\r\n"),
5 C/ n! ~+ |8 l: Y7 R
x = -39,# @3 i\" d) n- W. z
while{ x < 39,
\" H; O$ V8 ?7 M! @; c. x& Y: s
i = mandelbrot(x/40.0, y/40.0),! P& b$ G& {# G
which{ i==0,
0 _\" g8 ], i2 S# r+ \( H: c! F0 X
printff("*"),7 Z: V% c0 F5 o8 o
printff(" ")
- |. B1 L1 p6 R. b* I: Q; C
},) Y1 }& t1 Q( H9 P* o9 r5 l
x++
7 G\" R/ ^6 i4 p0 k1 U
},8 W( A8 |0 j) E$ {
y++$ O4 ~9 f1 x* `* u9 E
},$ r. q/ {0 e0 L% r\" e/ q. l
now=sys::clock(),
' V s; p; l8 _* W1 F* r
printff("\r\nForcal Elapsed {1,r}\r\n",(now-old)/1000)
+ X2 v& |! \! R( R3 N
};+ V4 o( ^# \& m' d5 P2 h

复制代码

运行结果：

*
0 R8 R u% O3 K4 J( K
*# e+ [& u! v' G- x J b; D
* W0 Y% N2 i\" \4 }
*/ |! U/ R- @2 p& f: p
*+ |7 x3 y m6 S
***& ]3 H7 h8 ~. l# _ C/ r0 j
*****
+ h: o. w H2 n0 K4 y( r
*****2 }) J6 Z' C; \) Z+ d5 @
***
# O& B1 S9 Q8 P; W0 a, F
*
' |' r. Q# Q: t5 u
*********
! ^( h$ Q) K; c
*************- B\" T$ J5 S- t9 m' B
***************1 B- i9 n- c( J' w# V. Z/ g
*********************
' s\" c- M/ t1 T8 N3 Q
*********************
# h' o% g% C: C1 g8 N% w/ ^
*******************
5 m# ]0 r9 y9 c\" o/ r
*******************0 y\" ^& ]/ l7 M/ q
*******************$ J: V) ] J5 D D
*******************
) M% }$ D0 J1 a4 s6 X8 h' {
***********************, g4 v' v4 u; |8 D9 N x4 K4 Q
*******************0 t8 X1 {! O0 X# C# }3 ~: R; E1 l
*******************- k$ ]* |* }- G
*********************
1 K& P3 N0 X\" d\" C/ K$ ^
*******************
; `* N\" M% \- }' q. Q
*******************
. u+ y' {* ~: g8 p1 r
*****************
% D; y' z% Y6 ^) t+ t6 Q6 R
***************, ?) ~1 H: }* K& \6 K* `
*************
; Z\" v- x p7 [$ u$ ]
*********
) L9 P8 Y7 C- D( w0 v8 s* u2 z
*/ U: z( J) }* y% C6 F\" y1 n+ v
***************$ {! W4 `6 f! D; @, N9 E Z9 p
***********************
( x5 j9 {, h* N8 A0 l1 u
* ************************* * a+ N- i8 L; H) S% H! U
*****************************
& B( F1 |7 F6 w9 B5 Q4 ]
* ******************************* *; v7 p3 W, U2 s\" R4 q- q
*********************************/ p4 |0 g8 \7 [1 K, {$ H3 \4 l
***********************************
' U2 G% |- E% t! X* m9 g
***************************************+ k: r. |: e4 J3 `
*** ***************************************** ***9 e3 f: V# L. O
*************************************************
. z! _% F$ G' {3 G' u
***********************************************
/ {3 j. e8 T9 ]3 c- F/ W5 a
*********************************************
. k( }9 Z& @8 {4 Y8 Y
*********************************************3 S% r) G2 n! d0 V\" K
***********************************************- S# e, q\" {9 E m
***********************************************3 b; q8 Q- ?. M( R/ x5 z
***************************************************5 M1 x3 s' o8 o% ]
*************************************************6 X6 ~0 q! G8 m: W, ?
*************************************************- v) g) G S' K\" \# b
***************************************************
M- x! K/ n3 H8 \. s) a
***************************************************
, {5 H( V( {& G, r; f+ \( j/ o
* *************************************************** *
; F9 n% G/ k0 C$ c
***** *************************************************** *****9 V. [, @0 W7 O, k6 S
****** *************************************************** ******
- L+ E* }1 k& m, k\" S( T\" ?/ @
******* *************************************************** *******: K5 j! W R& D* [, F
***********************************************************************
1 I& U2 _) H( U9 O2 u
********* *************************************************** *********
) V3 y% W% K6 \* D
****** *************************************************** ******+ t* h; K+ X& C# w! ^* I. c& V4 J; v
***** *************************************************** *****5 Z7 Q1 C( r$ U& Q' Z\" P
***************************************************+ l, e' r5 t& i0 I, W* N# C
***************************************************
* Q' \' w+ f: h+ Y2 c: [) E/ c
***************************************************
& C, b3 E/ }, a; I Z+ d
***************************************************- u7 p4 S; m, X8 o5 r! _ g9 G
*************************************************
, q$ H9 D& _# |' A, }% `
*************************************************
+ A7 N3 }0 c/ X- }5 h2 d/ V2 _; J
***************************************************# S, B+ Z. F+ Q2 v5 f\" o0 | x
***********************************************7 c' s6 n* {# c z5 e7 p% T2 a
***********************************************& ^( n* `\" X8 u
*******************************************# j' ]* E V6 [9 L) q6 P1 U
*****************************************8 L9 T) M0 {\" ?\" s( j* i
*********************************************( I# K. D% u% s4 O* ^\" x% N
**** ****************** ****************** ****1 {% Z5 ?- T3 _ C; P) p
*** **************** **************** ***0 f+ o/ l7 _+ t1 u
* ************** ************** *, e# t3 V6 Q% t) T' K# E
*********** ***********
' T\" j X- o& D( x* @
** ***** ***** **
; Z4 s$ J( ^/ v* k# K9 y# ]# q( b
* * * *# x/ J7 X9 ?! q: `7 k7 N