1 ! Copyright (C) 2009 Doug Coleman.
2 ! See http://factorcode.org/license.txt for BSD license.
3 USING: accessors combinators.short-circuit kernel lexer math
4 math.functions math.parser parser sequences splitting
8 TUPLE: decimal { mantissa read-only } { exponent read-only } ;
10 : <decimal> ( mantissa exponent -- decimal ) decimal boa ;
12 : >decimal< ( decimal -- mantissa exponent )
13 [ mantissa>> ] [ exponent>> ] bi ; inline
15 : string>decimal ( string -- decimal )
17 [ [ CHAR: 0 = ] trim-head [ "0" ] when-empty ]
18 [ [ CHAR: 0 = ] trim-tail [ "" ] when-empty ] bi*
19 [ append string>number ] [ nip length neg ] 2bi <decimal> ;
21 : parse-decimal ( -- decimal ) scan-token string>decimal ;
23 SYNTAX: D: parse-decimal suffix! ;
25 : decimal>ratio ( decimal -- ratio ) >decimal< 10^ * ;
26 : decimal>float ( decimal -- ratio ) decimal>ratio >float ;
28 : scale-mantissas ( D1 D2 -- m1 m2 exp )
29 [ [ mantissa>> ] bi@ ]
35 [ 10^ [ * ] curry dip f ] if
39 : scale-decimals ( D1 D2 -- D1' D2' )
40 scale-mantissas [ <decimal> ] curry bi@ ;
42 ERROR: decimal-types-expected d1 d2 ;
44 : guard-decimals ( obj1 obj2 -- D1 D2 )
45 2dup [ decimal? ] both?
46 [ decimal-types-expected ] unless ;
50 [ [ decimal? ] both? ]
54 [ [ mantissa>> ] same? ]
55 [ [ exponent>> ] same? ]
61 guard-decimals scale-decimals
62 [ mantissa>> ] bi@ < ;
65 [ mantissa>> abs ] [ exponent>> ] bi <decimal> ;
68 guard-decimals scale-mantissas [ + ] dip <decimal> ;
71 guard-decimals scale-mantissas [ - ] dip <decimal> ;
74 guard-decimals [ >decimal< ] bi@ swapd + [ * ] dip <decimal> ;
76 :: D/ ( D1 D2 a -- D3 )
77 D1 D2 guard-decimals 2drop
78 D1 >decimal< :> ( m1 e1 )
79 D2 >decimal< :> ( m2 e2 )
87 2dup before? [ 2drop +lt+ ] [ equal? +eq+ +gt+ ? ] if ; inline