1 USING: math math.functions kernel sequences io io.styles
2 prettyprint words hints ;
3 IN: benchmark.partial-sums
5 : summing ( n quot -- y )
6 [ >float ] swap [ + ] 3compose
7 0.0 -rot 1 -rot (each-integer) ; inline
9 : 2/3^k ( n -- y ) [ 2.0 3.0 / swap 1- ^ ] summing ;
13 : k^-0.5 ( n -- y ) [ -0.5 ^ ] summing ;
15 HINTS: k^-0.5 fixnum ;
17 : 1/k(k+1) ( n -- y ) [ dup 1+ * recip ] summing ;
19 HINTS: 1/k(k+1) fixnum ;
21 : cube ( x -- y ) dup dup * * ; inline
23 : flint-hills ( n -- y )
24 [ dup cube swap sin sq * recip ] summing ;
26 HINTS: flint-hills fixnum ;
28 : cookson-hills ( n -- y )
29 [ dup cube swap cos sq * recip ] summing ;
31 HINTS: cookson-hills fixnum ;
33 : harmonic ( n -- y ) [ recip ] summing ;
35 HINTS: harmonic fixnum ;
37 : riemann-zeta ( n -- y ) [ sq recip ] summing ;
39 HINTS: riemann-zeta fixnum ;
41 : -1^ 2 mod zero? 1 -1 ? ; inline
43 : alternating-harmonic ( n -- y ) [ dup -1^ swap / ] summing ;
45 HINTS: alternating-harmonic fixnum ;
47 : gregory ( n -- y ) [ dup -1^ swap 2 * 1- / ] summing ;
49 HINTS: gregory fixnum ;
52 { 2/3^k k^-0.5 1/k(k+1) flint-hills cookson-hills harmonic riemann-zeta alternating-harmonic gregory } ;
54 : partial-sums ( n -- )
55 standard-table-style [
57 [ tuck execute pprint-cell pprint-cell ] with-row
61 : partial-sums-main 2500000 partial-sums ;
63 MAIN: partial-sums-main