2005-02-08 22:02:44 -05:00
|
|
|
! Copyright (C) 2004, 2005 Slava Pestov.
|
|
|
|
! See http://factor.sf.net/license.txt for BSD license.
|
2004-12-24 02:52:02 -05:00
|
|
|
IN: math
|
2005-02-08 22:02:44 -05:00
|
|
|
USING: generic kernel kernel-internals math math-internals ;
|
2004-12-24 02:52:02 -05:00
|
|
|
|
|
|
|
GENERIC: numerator ( a/b -- a )
|
|
|
|
M: integer numerator ;
|
2004-12-31 02:17:45 -05:00
|
|
|
M: ratio numerator 0 slot %integer ;
|
2004-12-24 02:52:02 -05:00
|
|
|
|
|
|
|
GENERIC: denominator ( a/b -- b )
|
|
|
|
M: integer denominator drop 1 ;
|
2004-12-31 02:17:45 -05:00
|
|
|
M: ratio denominator 1 slot %integer ;
|
2004-12-24 02:52:02 -05:00
|
|
|
|
|
|
|
IN: math-internals
|
2004-07-16 02:26:21 -04:00
|
|
|
|
2004-12-18 23:18:32 -05:00
|
|
|
: 2>fraction ( a/b c/d -- a c b d )
|
|
|
|
[ swap numerator swap numerator ] 2keep
|
|
|
|
swap denominator swap denominator ; inline
|
2004-07-28 19:02:24 -04:00
|
|
|
|
2004-12-18 23:18:32 -05:00
|
|
|
M: ratio number= ( a/b c/d -- ? )
|
|
|
|
2>fraction number= [ number= ] [ 2drop f ] ifte ;
|
2004-07-16 02:26:21 -04:00
|
|
|
|
2004-12-18 23:18:32 -05:00
|
|
|
: scale ( a/b c/d -- a*d b*c )
|
2005-02-08 22:02:44 -05:00
|
|
|
2>fraction >r * swap r> * swap ;
|
2004-07-16 02:26:21 -04:00
|
|
|
|
2004-12-18 23:18:32 -05:00
|
|
|
: ratio+d ( a/b c/d -- b*d )
|
|
|
|
denominator swap denominator * ; inline
|
2004-07-16 02:26:21 -04:00
|
|
|
|
2004-12-18 23:18:32 -05:00
|
|
|
M: ratio < scale < ;
|
|
|
|
M: ratio <= scale <= ;
|
|
|
|
M: ratio > scale > ;
|
|
|
|
M: ratio >= scale >= ;
|
2004-07-16 02:26:21 -04:00
|
|
|
|
2004-12-18 23:18:32 -05:00
|
|
|
M: ratio + ( x y -- x+y ) 2dup scale + -rot ratio+d integer/ ;
|
|
|
|
M: ratio - ( x y -- x-y ) 2dup scale - -rot ratio+d integer/ ;
|
|
|
|
M: ratio * ( x y -- x*y ) 2>fraction * >r * r> integer/ ;
|
|
|
|
M: ratio / scale integer/ ;
|
|
|
|
M: ratio /i scale /i ;
|
|
|
|
M: ratio /f scale /f ;
|