version 1.39, 2002/12/05 19:25:23
|
version 1.43, 2004/05/08 17:14:30
|
Line 1
|
Line 1
|
\ High level floating point 14jan94py |
\ High level floating point 14jan94py |
|
|
\ Copyright (C) 1995,1997 Free Software Foundation, Inc. |
\ Copyright (C) 1995,1997,2003 Free Software Foundation, Inc. |
|
|
\ This file is part of Gforth. |
\ This file is part of Gforth. |
|
|
Line 211 set-current
|
Line 211 set-current
|
\G @code{f~abs}; r3=0: bitwise comparison; r3<0: @code{fnegate f~rel}. |
\G @code{f~abs}; r3=0: bitwise comparison; r3<0: @code{fnegate f~rel}. |
fdup f0= |
fdup f0= |
IF \ bitwise comparison |
IF \ bitwise comparison |
fp@ float+ 1 floats over float+ -text 0= |
fp@ float+ 1 floats over float+ over str= |
fdrop fdrop fdrop |
fdrop fdrop fdrop |
EXIT |
EXIT |
THEN |
THEN |
Line 221 set-current
|
Line 221 set-current
|
ELSE |
ELSE |
fnegate f~rel |
fnegate f~rel |
THEN ; |
THEN ; |
|
|
1e 53 0 [do] f2* [loop] fconstant fround-offset \ 2^53 |
|
|
|
: fround ( r1 -- r2 ) \ float f-round |
|
\G Round to nearest integral value. Break ties with round-to-even. |
|
\ assumes IEEE DP FP in round-to-even mode. |
|
\ for an explanation of this code read |
|
\ <2002Oct26.113823@a0.complang.tuwien.ac.at> ff. |
|
fdup f0> if |
|
fround-offset f- fround-offset f+ |
|
else fdup f0< if \ leave 0e and -0e as is |
|
fround-offset f+ fround-offset f- |
|
then then ; |
|
|
|
: f.s ( -- ) \ gforth f-dot-s |
|
\G Display the number of items on the floating-point stack, |
|
\G followed by a list of the items; TOS is the right-most item. |
|
." <" fdepth 0 .r ." > " fdepth 0 max maxdepth-.s @ min dup 0 |
|
?DO dup i - 1- floats fp@ + f@ f. LOOP drop ; |
|