basis/math/functions/functions-tests.factor

   1 USING: kernel math math.constants math.functions math.libm
   2 math.order math.ranges math.private sequences tools.test ;
   3
   4 IN: math.functions.tests
   5
   6 [ t ] [ 4 4 .00000001 ~ ] unit-test
   7 [ t ] [ 4.0000001 4.0000001 .000001 ~ ] unit-test
   8 [ f ] [ -4.0000001 4.0000001 .00001 ~ ] unit-test
   9 [ t ] [ -.0000000000001 0 .0000000001 ~ ] unit-test
  10 [ t ] [ 100 101 -.9 ~ ] unit-test
  11 [ f ] [ 100 120 -.09 ~ ] unit-test
  12 [ t ] [ 0 0 -.9 ~ ] unit-test
  13 [ f ] [ 0 10 -.9 ~ ] unit-test
  14
  15 ! Lets get the argument order correct, eh?
  16 [ 0.0 ] [ 0.0 1.0 fatan2 ] unit-test
  17 [ 0.25 ] [ 2.0 -2.0 fpow ] unit-test
  18
  19 [ 4.0 ] [ 16 sqrt ] unit-test
  20 [ 2.0 ] [ 4.0 0.5 ^ ] unit-test
  21 [ C{ 0.0 4.0 } ] [ -16 sqrt ] unit-test
  22
  23 [ 4 ] [ 2 2 ^ ] unit-test
  24 [ 1/4 ] [ 2 -2 ^ ] unit-test
  25 [ t ] [ 2 0.5 ^ 2 ^ 2 2.00001 between? ] unit-test
  26 [ t ] [ e pi i* ^ real-part -1.0 = ] unit-test
  27 [ t ] [ e pi i* ^ imaginary-part -0.00001 0.00001 between? ] unit-test
  28
  29 [ 1/0. ] [ 2.0 1024 ^ ] unit-test
  30 [ HEX: 1.0p-1024 ] [ 2.0 -1024 ^ ] unit-test
  31
  32 [ t ] [ 0 0 ^ fp-nan? ] unit-test
  33 [ 0.0 ] [ 0.0 1.0 ^ ] unit-test
  34 [ 1/0. ] [ 0 -2 ^ ] unit-test
  35 [ t ] [ 0 0.0 ^ fp-nan? ] unit-test
  36 [ t ] [ 0.0 0.0 ^ fp-nan? ] unit-test
  37 [ 1/0. ] [ 0 -2.0 ^ ] unit-test
  38 [ 0 ] [ 0 3.0 ^ ] unit-test
  39 [ 0 ] [ 0 3 ^ ] unit-test
  40
  41 : factorial ( n -- n! ) [ 1 ] [ [1,b] 1 [ * ] reduce ] if-zero ;
  42
  43 [ 0.0 0 ] [ 0 frexp ] unit-test
  44 [ 0.5 1 ] [ 1 frexp ] unit-test
  45 [ -0.5 1 ] [ -1 frexp ] unit-test
  46 [ 0.5 2 ] [ 2 frexp ] unit-test
  47 [ -0.5 2 ] [ -2 frexp ] unit-test
  48 [ 0.64 -6 ] [ 0.01 frexp ] unit-test
  49 [ -0.64 -6 ] [ -0.01 frexp ] unit-test
  50 [ 0.75 0 ] [ 0.75 frexp ] unit-test
  51 [ -0.75 0 ] [ -0.75 frexp ] unit-test
  52 [ 1/0. 0 ] [ 1/0. frexp ] unit-test
  53 [ -1/0. 0 ] [ -1/0. frexp ] unit-test
  54 [ 0.6588418960767314 8530 t ] [ 1000 factorial [ frexp ] [ bignum? ] bi ] unit-test
  55 [ -0.6588418960767314 8530 t ] [ 1000 factorial neg [ frexp ] [ bignum? ] bi ] unit-test
  56
  57 [ 0.0 ] [ 1 log ] unit-test
  58 [ 0.0 ] [ 1.0 log ] unit-test
  59 [ 1.0 ] [ e log ] unit-test
  60 [ 5912.128178488163 t ] [ 1000 factorial [ log ] [ bignum? ] bi ] unit-test
  61 [ C{ 5912.128178488163 3.141592653589793 } t ] [ 1000 factorial neg [ log ] [ bignum? ] bi ] unit-test
  62
  63 [ 0.0 ] [ 1.0 log10 ] unit-test
  64 [ 1.0 ] [ 10.0 log10 ] unit-test
  65 [ 2.0 ] [ 100.0 log10 ] unit-test
  66 [ 3.0 ] [ 1000.0 log10 ] unit-test
  67 [ 4.0 ] [ 10000.0 log10 ] unit-test
  68 [ 2567.604644222133 t ] [ 1000 factorial [ log10 ] [ bignum? ] bi ] unit-test
  69
  70 [ t ] [ 1 exp e 1.e-10 ~ ] unit-test
  71 [ f ] [ 1 exp 0/0. 1.e-10 ~ ] unit-test
  72 [ f ] [ 0/0. 1 exp 1.e-10 ~ ] unit-test
  73 [ t ] [ 1.0 exp e 1.e-10 ~ ] unit-test
  74 [ t ] [ -1 exp e * 1.0 1.e-10 ~ ] unit-test
  75 [ f ] [ 1/0. 1/0. 1.e-10 ~ ] unit-test
  76 [ f ] [ 1/0. -1/0. 1.e-10 ~ ] unit-test
  77 [ f ] [ 1/0. 0/0. 1.e-10 ~ ] unit-test
  78 [ f ] [ 0/0. -1/0. 1.e-10 ~ ] unit-test
  79
  80 [ 1.0 ] [ 0 cosh ] unit-test
  81 [ 1.0 ] [ 0.0 cosh ] unit-test
  82 [ 0.0 ] [ 1 acosh ] unit-test
  83 [ 0.0 ] [ 1.0 acosh ] unit-test
  84
  85 [ 1.0 ] [ 0 cos ] unit-test
  86 [ 1.0 ] [ 0.0 cos ] unit-test
  87 [ 0.0 ] [ 1 acos ] unit-test
  88 [ 0.0 ] [ 1.0 acos ] unit-test
  89
  90 [ 0.0 ] [ 0 sinh ] unit-test
  91 [ 0.0 ] [ 0.0 sinh ] unit-test
  92 [ 0.0 ] [ 0 asinh ] unit-test
  93 [ 0.0 ] [ 0.0 asinh ] unit-test
  94
  95 [ 0.0 ] [ 0 sin ] unit-test
  96 [ 0.0 ] [ 0.0 sin ] unit-test
  97 [ 0.0 ] [ 0 asin ] unit-test
  98 [ 0.0 ] [ 0.0 asin ] unit-test
  99
 100 [ 0.0 ] [ 0 tan ] unit-test
 101 [ t ] [ pi 2 / tan 1.e10 > ] unit-test
 102
 103 [ t ] [ 10 atan real? ] unit-test
 104 [ t ] [ 10.0 atan real? ] unit-test
 105 [ f ] [ 10 atanh real? ] unit-test
 106 [ f ] [ 10.0 atanh real? ] unit-test
 107
 108 [ t ] [ 10 asin sin 10 1.e-10 ~ ] unit-test
 109 [ t ] [ -1 sqrt neg dup acos cos 1.e-10 ~ ] unit-test
 110 [ t ] [ -100 atan tan -100 1.e-10 ~ ] unit-test
 111 [ t ] [ 10 asinh sinh 10 1.e-10 ~ ] unit-test
 112 [ t ] [ 10 atanh tanh 10 1.e-10 ~ ] unit-test
 113 [ t ] [ 0.5 atanh tanh 0.5 1.e-10 ~ ] unit-test
 114
 115 [ 100 ] [ 100 100 gcd nip ] unit-test
 116 [ 100 ] [ 1000 100 gcd nip ] unit-test
 117 [ 100 ] [ 100 1000 gcd nip ] unit-test
 118 [ 4 ] [ 132 64 gcd nip ] unit-test
 119 [ 4 ] [ -132 64 gcd nip ] unit-test
 120 [ 4 ] [ -132 -64 gcd nip ] unit-test
 121 [ 4 ] [ 132 -64 gcd nip ] unit-test
 122 [ 4 ] [ -132 -64 gcd nip ] unit-test
 123
 124 [ 100 ] [ 100 >bignum 100 >bignum gcd nip ] unit-test
 125 [ 100 ] [ 1000 >bignum 100 >bignum gcd nip ] unit-test
 126 [ 100 ] [ 100 >bignum 1000 >bignum gcd nip ] unit-test
 127 [ 4 ] [ 132 >bignum 64 >bignum gcd nip ] unit-test
 128 [ 4 ] [ -132 >bignum 64 >bignum gcd nip ] unit-test
 129 [ 4 ] [ -132 >bignum -64 >bignum gcd nip ] unit-test
 130 [ 4 ] [ 132 >bignum -64 >bignum gcd nip ] unit-test
 131 [ 4 ] [ -132 >bignum -64 >bignum gcd nip ] unit-test
 132
 133 [ 6 ] [
 134     1326264299060955293181542400000006
 135     1591517158873146351817850880000000
 136     gcd nip
 137 ] unit-test
 138
 139 [ 11 ] [
 140     13262642990609552931815424
 141     159151715887314635181785
 142     gcd nip
 143 ] unit-test
 144
 145 [ 3 ] [
 146     13262642990609552931
 147     1591517158873146351
 148     gcd nip
 149 ] unit-test
 150
 151 [ 26525285981219 ] [
 152     132626429906095
 153     159151715887314
 154     gcd nip
 155 ] unit-test
 156
 157
 158 : verify-gcd ( a b -- ? )
 159     2dup gcd
 160     [ rot * swap rem ] dip = ;
 161
 162 [ t ] [ 123 124 verify-gcd ] unit-test
 163 [ t ] [ 50 120 verify-gcd ] unit-test
 164
 165 [ t ] [ 0 42 divisor? ] unit-test
 166 [ t ] [ 42 7 divisor? ] unit-test
 167 [ t ] [ 42 -7 divisor? ] unit-test
 168 [ t ] [ 42 42 divisor? ] unit-test
 169 [ f ] [ 42 16 divisor? ] unit-test
 170
 171 [ 3 ] [ 5 7 mod-inv ] unit-test
 172 [ 78572682077 ] [ 234829342 342389423843 mod-inv ] unit-test
 173
 174 [ 2 10 mod-inv ] must-fail
 175
 176 [ t ] [ 0 0 ^ fp-nan? ] unit-test
 177 [ 1 ] [ 10 0 ^ ] unit-test
 178 [ 1/8 ] [ 1/2 3 ^ ] unit-test
 179 [ 1/8 ] [ 2 -3 ^ ] unit-test
 180 [ t ] [ 1 100 shift 2 100 ^ = ] unit-test
 181
 182 [ 1 ] [ 7/8 ceiling ] unit-test
 183 [ 2 ] [ 3/2 ceiling ] unit-test
 184 [ 0 ] [ -7/8 ceiling ] unit-test
 185 [ -1 ] [ -3/2 ceiling ] unit-test
 186
 187 [ 4.0 ] [ 4.5 truncate ] unit-test
 188 [ 4.0 ] [ 4.5 floor ] unit-test
 189 [ 5.0 ] [ 4.5 ceiling ] unit-test
 190
 191 [ -4.0 ] [ -4.5 truncate ] unit-test
 192 [ -5.0 ] [ -4.5 floor ] unit-test
 193 [ -4.0 ] [ -4.5 ceiling ] unit-test
 194
 195 [ -4.0 ] [ -4.0 truncate ] unit-test
 196 [ -4.0 ] [ -4.0 floor ] unit-test
 197 [ -4.0 ] [ -4.0 ceiling ] unit-test
 198
 199 [ -5.0 ] [ -4.5 round ] unit-test
 200 [ -4.0 ] [ -4.4 round ] unit-test
 201 [ 5.0 ] [ 4.5 round ] unit-test
 202 [ 4.0 ] [ 4.4 round ] unit-test
 203
 204 [ 6 59967 ] [ 3837888 factor-2s ] unit-test
 205 [ 6 -59967 ] [ -3837888 factor-2s ] unit-test
 206
 207 [ 1 ] [
 208     183009416410801897
 209     1067811677921310779
 210     2135623355842621559
 211     ^mod
 212 ] unit-test
 213
 214 [ 1 ] [
 215     183009416410801897
 216     1067811677921310779
 217     2135623355842621559
 218     [ >bignum ] tri@ ^mod
 219 ] unit-test
 220
 221 [ 1.0  ] [ 1.0 2.5 0.0 lerp ] unit-test
 222 [ 2.5  ] [ 1.0 2.5 1.0 lerp ] unit-test
 223 [ 1.75 ] [ 1.0 2.5 0.5 lerp ] unit-test
 224
 225 [ C{ 1 2 } ] [ C{ 1 2 } 1 ^ ] unit-test