2005-04-02 02:39:33 -05:00
|
|
|
! Copyright (C) 2004, 2005 Slava Pestov.
|
2006-01-12 00:34:56 -05:00
|
|
|
! See http://factorcode.org/license.txt for BSD license.
|
2005-08-31 21:06:13 -04:00
|
|
|
IN: math
|
|
|
|
USING: errors generic kernel math-internals namespaces sequences
|
|
|
|
strings ;
|
2004-07-22 19:48:50 -04:00
|
|
|
|
2005-09-17 22:25:18 -04:00
|
|
|
: not-a-number "Not a number" throw ;
|
2004-07-22 19:48:50 -04:00
|
|
|
|
2005-12-24 23:01:49 -05:00
|
|
|
DEFER: base>
|
|
|
|
|
|
|
|
: string>ratio ( "a/b" radix -- a/b )
|
|
|
|
>r "/" split1 r> tuck base> >r base> r> / ;
|
|
|
|
|
2004-12-18 23:18:32 -05:00
|
|
|
GENERIC: digit> ( ch -- n )
|
|
|
|
M: digit digit> CHAR: 0 - ;
|
|
|
|
M: letter digit> CHAR: a - 10 + ;
|
|
|
|
M: LETTER digit> CHAR: A - 10 + ;
|
|
|
|
M: object digit> not-a-number ;
|
2004-07-22 19:48:50 -04:00
|
|
|
|
2004-12-16 18:36:26 -05:00
|
|
|
: digit+ ( num digit base -- num )
|
2005-09-24 15:21:17 -04:00
|
|
|
2dup < [ rot * + ] [ not-a-number ] if ;
|
2004-07-22 19:48:50 -04:00
|
|
|
|
2005-12-24 23:01:49 -05:00
|
|
|
: (string>integer) ( base str -- num )
|
2005-04-29 02:37:12 -04:00
|
|
|
dup empty? [
|
2004-08-06 18:40:44 -04:00
|
|
|
not-a-number
|
|
|
|
] [
|
2005-06-25 20:39:53 -04:00
|
|
|
0 [ digit> pick digit+ ] reduce nip
|
2005-09-24 15:21:17 -04:00
|
|
|
] if ;
|
2004-07-22 19:48:50 -04:00
|
|
|
|
2005-12-24 23:01:49 -05:00
|
|
|
: string>integer ( string -- n )
|
|
|
|
swap "-" ?head >r (string>integer) r> [ neg ] when ;
|
2004-12-16 18:36:26 -05:00
|
|
|
|
2005-12-24 23:01:49 -05:00
|
|
|
: base> ( string radix -- n )
|
2005-10-29 23:25:38 -04:00
|
|
|
{
|
2005-12-24 23:01:49 -05:00
|
|
|
{ [ CHAR: / pick member? ] [ string>ratio ] }
|
|
|
|
{ [ CHAR: . pick member? ] [ drop string>float ] }
|
|
|
|
{ [ t ] [ string>integer ] }
|
2005-10-29 23:25:38 -04:00
|
|
|
} cond ;
|
2004-08-04 22:43:58 -04:00
|
|
|
|
2005-12-24 23:01:49 -05:00
|
|
|
: string>number ( string -- num ) 10 base> ;
|
|
|
|
: bin> ( string -- num ) 2 base> ;
|
|
|
|
: oct> ( string -- num ) 8 base> ;
|
|
|
|
: hex> ( string -- num ) 16 base> ;
|
2005-08-21 20:50:14 -04:00
|
|
|
|
|
|
|
: >digit ( n -- ch )
|
2005-09-24 15:21:17 -04:00
|
|
|
dup 10 < [ CHAR: 0 + ] [ 10 - CHAR: a + ] if ;
|
2005-08-21 20:50:14 -04:00
|
|
|
|
|
|
|
: integer, ( num radix -- )
|
2005-08-31 21:06:13 -04:00
|
|
|
dup >r /mod >digit , dup 0 >
|
2005-09-24 15:21:17 -04:00
|
|
|
[ r> integer, ] [ r> 2drop ] if ;
|
2005-08-21 20:50:14 -04:00
|
|
|
|
2005-12-24 23:01:49 -05:00
|
|
|
G: >base ( num radix -- string ) [ over ] standard-combination ;
|
|
|
|
|
|
|
|
M: integer >base ( num radix -- string )
|
2005-08-21 20:50:14 -04:00
|
|
|
[
|
|
|
|
over 0 < [
|
|
|
|
swap neg swap integer, CHAR: - ,
|
|
|
|
] [
|
|
|
|
integer,
|
2005-09-24 15:21:17 -04:00
|
|
|
] if
|
2005-08-25 15:27:38 -04:00
|
|
|
] "" make reverse ;
|
2005-08-21 20:50:14 -04:00
|
|
|
|
2005-12-24 23:01:49 -05:00
|
|
|
M: ratio >base ( num radix -- string )
|
|
|
|
[
|
|
|
|
over numerator over >base %
|
|
|
|
CHAR: / ,
|
|
|
|
swap denominator swap >base %
|
|
|
|
] "" make ;
|
2005-08-21 20:50:14 -04:00
|
|
|
|
2006-01-26 23:01:14 -05:00
|
|
|
: fix-float
|
2005-12-24 23:01:49 -05:00
|
|
|
CHAR: . over member? [ ".0" append ] unless ;
|
|
|
|
|
2006-01-26 23:01:14 -05:00
|
|
|
M: float >base ( num radix -- string )
|
|
|
|
drop {
|
|
|
|
{ [ dup 1.0/0.0 = ] [ drop "1.0/0.0" ] }
|
|
|
|
{ [ dup -1.0/0.0 = ] [ drop "-1.0/0.0" ] }
|
2006-01-28 15:49:31 -05:00
|
|
|
{ [ dup fp-nan? ] [ drop "0.0/0.0" ] }
|
2006-01-26 23:01:14 -05:00
|
|
|
{ [ t ] [ float>string fix-float ] }
|
|
|
|
} cond ;
|
|
|
|
|
2005-12-24 23:01:49 -05:00
|
|
|
: number>string ( num -- string ) 10 >base ;
|
|
|
|
: >bin ( num -- string ) 2 >base ;
|
|
|
|
: >oct ( num -- string ) 8 >base ;
|
|
|
|
: >hex ( num -- string ) 16 >base ;
|