extra/benchmark/partial-sums/partial-sums.factor

   1 USING: math math.functions kernel sequences io io.styles
   2 prettyprint words hints ;
   3 IN: benchmark.partial-sums
   4
   5 : summing ( n quot -- y )
   6     [ + ] compose 0.0 -rot 1 -rot (each-integer) ; inline
   7
   8 : 2/3^k ( n -- y ) [ 2.0 3.0 / swap 1- ^ ] summing ;
   9
  10 HINTS: 2/3^k fixnum ;
  11
  12 : k^-0.5 ( n -- y ) [ -0.5 ^ ] summing ;
  13
  14 HINTS: k^-0.5 fixnum ;
  15
  16 : 1/k(k+1) ( n -- y ) [ dup 1+ * recip ] summing ;
  17
  18 HINTS: 1/k(k+1) fixnum ;
  19
  20 : cube ( x -- y ) dup dup * * ; inline
  21
  22 : flint-hills ( n -- y )
  23     [ dup cube swap sin sq * recip ] summing ;
  24
  25 HINTS: flint-hills fixnum ;
  26
  27 : cookson-hills ( n -- y )
  28     [ dup cube swap cos sq * recip ] summing ;
  29
  30 HINTS: cookson-hills fixnum ;
  31
  32 : harmonic ( n -- y ) [ recip ] summing ;
  33
  34 HINTS: harmonic fixnum ;
  35
  36 : riemann-zeta ( n -- y ) [ sq recip ] summing ;
  37
  38 HINTS: riemann-zeta fixnum ;
  39
  40 : -1^ 2 mod zero? 1 -1 ? ; inline
  41
  42 : alternating-harmonic ( n -- y ) [ dup -1^ swap / ] summing ;
  43
  44 HINTS: alternating-harmonic fixnum ;
  45
  46 : gregory ( n -- y ) [ dup -1^ swap 2 * 1- / ] summing ;
  47
  48 HINTS: gregory fixnum ;
  49
  50 : functions
  51     { 2/3^k k^-0.5 1/k(k+1) flint-hills cookson-hills harmonic riemann-zeta alternating-harmonic gregory } ;
  52
  53 : partial-sums ( n -- )
  54     standard-table-style [
  55         functions [
  56             [ tuck execute pprint-cell pprint-cell ] with-row
  57         ] with each
  58     ] tabular-output ;
  59
  60 : partial-sums-main 2500000 partial-sums ;
  61
  62 MAIN: partial-sums-main