version 1.39, 2002/12/05 19:25:23
|
version 1.51, 2007/07/06 12:54:56
|
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,2004,2005,2006 Free Software Foundation, Inc. |
|
|
\ This file is part of Gforth. |
\ This file is part of Gforth. |
|
|
Line 51
|
Line 51
|
\G enough space to align it. |
\G enough space to align it. |
here dup dfaligned swap ?DO bl c, LOOP ; |
here dup dfaligned swap ?DO bl c, LOOP ; |
|
|
1 sfloats (Field) sfloat+ , ( sf-addr1 -- sf-addr2 ) \ float-ext s-float-plus |
(Field) sfloat+ ( sf-addr1 -- sf-addr2 ) \ float-ext s-float-plus |
\G @code{1 sfloats +}. |
\G @code{1 sfloats +}. |
|
1 sfloats , |
|
|
1 dfloats (Field) dfloat+ , ( df-addr1 -- df-addr2 ) \ float-ext d-float-plus |
(Field) dfloat+ ( df-addr1 -- df-addr2 ) \ float-ext d-float-plus |
\G @code{1 dfloats +}. |
\G @code{1 dfloats +}. |
|
1 dfloats , |
|
|
: f, ( f -- ) \ gforth |
: f, ( f -- ) \ gforth |
\G Reserve data space for one floating-point number and store |
\G Reserve data space for one floating-point number and store |
\G @i{f} in the space. |
\G @i{f} in the space. |
Line 127 DOES> ( -- r )
|
Line 129 DOES> ( -- r )
|
scratch over c@ emit '. emit 1 /string type |
scratch over c@ emit '. emit 1 /string type |
'E emit . ; |
'E emit . ; |
|
|
require debugs.fs |
|
|
|
: sfnumber ( c-addr u -- r true | false ) |
: sfnumber ( c-addr u -- r true | false ) |
2dup [CHAR] e scan ( c-addr u c-addr2 u2 ) |
2dup [CHAR] e scan ( c-addr u c-addr2 u2 ) |
dup 0= |
dup 0= |
Line 142 require debugs.fs
|
Line 142 require debugs.fs
|
2drop false |
2drop false |
THEN ; |
THEN ; |
|
|
:noname ( c-addr u -- ) |
[ifundef] compiler-notfound1 |
|
defer compiler-notfound1 |
|
' no.extensions IS compiler-notfound1 |
|
|
|
:noname compiler-notfound1 execute ; is compiler-notfound |
|
|
|
defer interpreter-notfound1 |
|
' no.extensions IS interpreter-notfound1 |
|
|
|
:noname interpreter-notfound1 execute ; is interpreter-notfound |
|
[then] |
|
|
|
:noname ( c-addr u -- ... xt ) |
2dup sfnumber |
2dup sfnumber |
IF |
IF |
2drop POSTPONE FLiteral |
2drop [comp'] FLiteral |
ELSE |
ELSE |
defers compiler-notfound |
defers compiler-notfound1 |
ENDIF ; |
ENDIF ; |
IS compiler-notfound |
IS compiler-notfound1 |
|
|
:noname ( c-addr u -- r ) |
:noname ( c-addr u -- ... xt ) |
2dup sfnumber |
2dup sfnumber |
IF |
IF |
2drop |
2drop ['] noop |
ELSE |
ELSE |
defers interpreter-notfound |
defers interpreter-notfound1 |
ENDIF ; |
ENDIF ; |
IS interpreter-notfound |
IS interpreter-notfound1 |
|
|
: fvariable ( "name" -- ) \ float f-variable |
: fvariable ( "name" -- ) \ float f-variable |
Create 0.0E0 f, ; |
Create 0.0E0 f, ; |
Line 211 set-current
|
Line 223 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 222 set-current
|
Line 234 set-current
|
fnegate f~rel |
fnegate f~rel |
THEN ; |
THEN ; |
|
|
1e 53 0 [do] f2* [loop] fconstant fround-offset \ 2^53 |
\ proposals from Krishna Myeni in <cjsp2d$47l$1@ngspool-d02.news.aol.com> |
|
\ not sure if they are a good idea |
|
|
: fround ( r1 -- r2 ) \ float f-round |
: FTRUNC ( r1 -- r2 ) |
\G Round to nearest integral value. Break ties with round-to-even. |
\ round towards 0 |
\ assumes IEEE DP FP in round-to-even mode. |
\ !! should be implemented properly |
\ for an explanation of this code read |
F>D D>F ; |
\ <2002Oct26.113823@a0.complang.tuwien.ac.at> ff. |
|
fdup f0> if |
: FMOD ( r1 r2 -- r ) |
fround-offset f- fround-offset f+ |
\ remainder of r1/r2 |
else fdup f0< if \ leave 0e and -0e as is |
FOVER FOVER F/ ftrunc F* F- ; |
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 ; |
|