steve
/
factor
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

fix almost=

darcs
erg 2006-09-28 06:05:43 +00:00
parent ce8ed04543
commit f0f291d8ff
2 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
contrib/math/test.factor
View File

@ -47,15 +47,15 @@ USING: kernel math test sequences math-contrib ;
[ 1 ] [ 2 0 nCk ] unit-test
[ 1 ] [ 2 0 nPk ] unit-test
[ t ] [ -9000000000000000000000000000000000000000000 gamma 1/0. = ] unit-test
[ t ] [ -1.5 gamma 2.36327 almost= ] unit-test
[ t ] [ -1.5 gamma 2.363271801207344 almost= ] unit-test
[ t ] [ -1 gamma 1/0. = ] unit-test
[ t ] [ -0.5 gamma -3.5449 almost= ] unit-test
[ t ] [ -0.5 gamma -3.544907701811 almost= ] unit-test
[ t ] [ 0 gamma 1/0. = ] unit-test
[ t ] [ .5 gamma 1.7724538 almost= ] unit-test
[ t ] [ 1 gamma 1 almost= ] unit-test
[ t ] [ 2 gamma 1 almost= ] unit-test
[ t ] [ 3 gamma 2 almost= ] unit-test
[ t ] [ 11 gamma 3628800 almost= ] unit-test
[ t ] [ 11 gamma 3628800.000015679 almost= ] unit-test
[ t ] [ 90000000000000000000000000000000000000000000 gamma 1/0. = ] unit-test
! some fun identities
[ t ] [ 2/3 gamma 2 pi * 3 sqrt 1/3 gamma * / almost= ] unit-test
@ -71,7 +71,7 @@ USING: kernel math test sequences math-contrib ;
[ t ] [ 1 gammaln 0 almost= ] unit-test
[ t ] [ 2 gammaln 0 almost= ] unit-test
[ t ] [ 3 gammaln 0.693147180 almost= ] unit-test
[ t ] [ 11 gammaln 15.1044 almost= ] unit-test
[ t ] [ 11 gammaln 15.1044125 almost= ] unit-test
[ t ] [ 9000000000000000000000000000000000000000000 gammaln 8.811521863477754e+44 almost= ] unit-test
[ 1 ] [ qi norm ] unit-test

2
contrib/math/utils.factor
View File

@ -64,7 +64,7 @@ USING: errors kernel sequences math sequences-internals namespaces arrays ;
#! find the absolute values of the min and max of a seq in one pass
minmax 2dup [ abs ] 2apply > [ swap ] when ;
SYMBOL: almost=-precision .0001 almost=-precision set-global
SYMBOL: almost=-precision .000001 almost=-precision set-global
: almost= ( a b -- bool )
- abs almost=-precision get < ;