1 ! Copyright (C) 2012 John Benediktsson
2 ! See http://factorcode.org/license.txt for BSD license
4 USING: arrays kernel math math.extras math.ranges sequences
9 { { 1 -1/2 1/6 0 -1/30 0 1/42 0 -1/30 0 } }
10 [ 10 iota [ bernoulli ] map ] unit-test
12 { -1 } [ -1 7 jacobi ] unit-test
13 { 0 } [ 3 3 jacobi ] unit-test
14 { -1 } [ 127 703 jacobi ] unit-test
15 { 1 } [ -4 197 jacobi ] unit-test
17 { { 2 3 4 5 6 7 8 9 } } [ 10 [1,b] 3 moving-average ] unit-test
18 { { 1+1/2 2+1/2 3+1/2 4+1/2 5+1/2 6+1/2 7+1/2 8+1/2 9+1/2 } }
19 [ 10 [1,b] 2 moving-average ] unit-test
21 { { 1 1+1/2 2+1/4 3+1/8 4+1/16 5+1/32 } }
22 [ 6 [1,b] 1/2 exponential-moving-average ] unit-test
23 { { 1 3 3 5 5 7 7 9 9 11 } }
24 [ 10 [1,b] 2 exponential-moving-average ] unit-test
26 { { 2 5 5 4 3 } } [ { 1 2 5 6 1 4 3 } 3 moving-median ] unit-test
28 { { 2 1 1 2 1 1 } } [ { 1 1 2 3 5 8 13 } 2 [ odd? ] moving-count ] unit-test
30 { { } } [ { 0 0 } nonzero ] unit-test
31 { { 1 2 3 } } [ { 0 1 0 2 0 3 0 } nonzero ] unit-test
33 { { } } [ 0 bartlett ] unit-test
34 { { 1 } } [ 1 bartlett ] unit-test
35 { { 0 0 } } [ 2 bartlett ] unit-test
36 { { 0 1 0 } } [ 3 bartlett ] unit-test
37 { { 0 2/3 2/3 0 } } [ 4 bartlett ] unit-test
38 { { 0 1/2 1 1/2 0 } } [ 5 bartlett ] unit-test
39 { { 0 2/5 4/5 4/5 2/5 0 } } [ 6 bartlett ] unit-test
60 } concentration-coefficient
63 { 0 } [ { 1 } gini ] unit-test
64 { 0 } [ { 1 1 1 1 1 1 } gini ] unit-test
65 { 0 } [ { 10 10 10 10 } gini ] unit-test
66 { 0 } [ { } gini ] unit-test
68 { 0 } [ { 1 } concentration-coefficient ] unit-test
69 { 0 } [ { 1 1 1 1 1 1 } concentration-coefficient ] unit-test
70 { 0 } [ { 10 10 10 10 } concentration-coefficient ] unit-test
71 { 0 } [ { } concentration-coefficient ] unit-test
73 { 57/200 } [ { 80 60 10 20 30 } herfindahl ] unit-test
74 { 17/160 } [ { 80 60 10 20 30 } normalized-herfindahl ] unit-test
77 { -10 10 2 2.5 3 } { 1 2 3 4 5 } digitize]
81 { -10 10 2 2.5 3 } { 1 2 3 4 5 } digitize)
87 } [ { 1 2 1 1 3 4 } unique-indices ] unit-test
89 { { 1 8+4/5 16+3/5 24+2/5 32+1/5 } } [ 1 40 5 linspace[a,b) >array ] unit-test
90 { { 1 10+3/4 20+1/2 30+1/4 40 } } [ 1 40 5 linspace[a,b] >array ] unit-test
92 { f } [ { } majority ] unit-test
93 { 1 } [ { 1 } majority ] unit-test
94 { f } [ { 1 2 } majority ] unit-test
95 { 1 } [ { 1 1 2 } majority ] unit-test
96 { f } [ { 1 1 2 2 } majority ] unit-test
97 { 2 } [ { 1 1 2 2 2 } majority ] unit-test
98 { 3 } [ { 1 2 3 1 2 3 1 2 3 3 } majority ] unit-test
99 { CHAR: C } [ "AAACCBBCCCBCC" majority ] unit-test
101 { -5 } [ -4-3/5 round-to-even ] unit-test
102 { -4 } [ -4-1/2 round-to-even ] unit-test
103 { -4 } [ -4-2/5 round-to-even ] unit-test
104 { 5 } [ 4+3/5 round-to-even ] unit-test
105 { 4 } [ 4+1/2 round-to-even ] unit-test
106 { 4 } [ 4+2/5 round-to-even ] unit-test
108 { -5.0 } [ -4.6 round-to-even ] unit-test
109 { -4.0 } [ -4.5 round-to-even ] unit-test
110 { -4.0 } [ -4.4 round-to-even ] unit-test
111 { 5.0 } [ 4.6 round-to-even ] unit-test
112 { 4.0 } [ 4.5 round-to-even ] unit-test
113 { 4.0 } [ 4.4 round-to-even ] unit-test
115 { 0.0 } [ 0 2 round-to-decimal ] unit-test
116 { 1.0 } [ 1 2 round-to-decimal ] unit-test
117 { 1.23 } [ 1.2349 2 round-to-decimal ] unit-test
118 { 1.24 } [ 1.2350 2 round-to-decimal ] unit-test
119 { 1.24 } [ 1.2351 2 round-to-decimal ] unit-test
120 { -1.23 } [ -1.2349 2 round-to-decimal ] unit-test
121 { -1.24 } [ -1.2350 2 round-to-decimal ] unit-test
122 { -1.24 } [ -1.2351 2 round-to-decimal ] unit-test
125 0.0 0.0 10000.0 12000.0 12300.0 12350.0 12346.0 12345.7
126 12345.68 12345.679 12345.6789 12345.6789 12345.678901
127 12345.6789012 12345.67890123 12345.678901235
129 } [ 12345.67890123456 -6 9 [a,b] [ round-to-decimal ] with map ] unit-test
131 { 0 } [ 0 5 round-to-step ] unit-test
132 { 0 } [ 1 5 round-to-step ] unit-test
133 { 5 } [ 3 5 round-to-step ] unit-test
134 { 10 } [ 12 5 round-to-step ] unit-test
135 { 15 } [ 13 5 round-to-step ] unit-test
137 { 0 } [ 0 round-away-from-zero ] unit-test
138 { -1.0 } [ -0.1 round-away-from-zero ] unit-test
139 { 1.0 } [ 0.1 round-away-from-zero ] unit-test
140 { -2.0 } [ -1.9 round-away-from-zero ] unit-test
141 { 2.0 } [ 1.9 round-away-from-zero ] unit-test
143 { { 0 1 2 3 0 0 1 } } [ { 1 2 3 3 2 1 2 } [ <= ] monotonic-count ] unit-test
144 { 4 } [ { 1 2 3 1 2 3 4 5 } [ < ] max-monotonic-count ] unit-test
146 { 2470 } [ 20 iota sum-squares ] unit-test
147 { 2470 } [ 20 iota >array sum-squares ] unit-test
149 { 36100 } [ 20 iota sum-cubes ] unit-test
150 { 36100 } [ 20 iota >array sum-cubes ] unit-test
154 1 -1 -1 0 -1 1 -1 0 0 1 -1 0 -1 1 1
155 0 -1 0 -1 0 1 1 -1 0 0 1 0 0 -1 -1
158 30 [1,b] [ mobius ] map
161 { 1/5 } [ 3/5 1 kelly ] unit-test
162 { 0 } [ 1/2 1 kelly ] unit-test
163 { -1/5 } [ 2/5 1 kelly ] unit-test