gforth/float.fs - diff

Return to float.fs CVS log

Up to [gforth] / gforth

version 1.53, 2007/12/31 18:40:24	version 1.57, 2010/03/22 17:08:06
Line 1	Line 1
\ High level floating point 14jan94py	\ High level floating point 14jan94py

\ Copyright (C) 1995,1997,2003,2004,2005,2006,2007 Free Software Foundation, Inc.	\ Copyright (C) 1995,1997,2003,2004,2005,2006,2007,2009 Free Software Foundation, Inc.

\ This file is part of Gforth.	\ This file is part of Gforth.

Line 233 set-current	Line 233 set-current
fnegate f~rel	fnegate f~rel
THEN ;	THEN ;

	-0e 8 0 [do] fp@ [i] + c@ $80 = [if] [i] constant fsign-offset [then] [loop]

	: fcopysign ( r1 r2 -- r3 ) \ gforth
	\G r3 takes its absolute value from r1 and its sign from r2
	\ !! implementation relies on IEEE DP format
	fp@ fsign-offset + dup c@ $80 and >r ( r1 r2 addr-r1sign )
	float+ dup c@ $7f and r> or swap c!
	fdrop ;

\ proposals from Krishna Myeni in <cjsp2d$47l$1@ngspool-d02.news.aol.com>	\ proposals from Krishna Myeni in <cjsp2d$47l$1@ngspool-d02.news.aol.com>
\ not sure if they are a good idea	\ not sure if they are a good idea

: FTRUNC ( r1 -- r2 )	: ftrunc ( r1 -- r2 ) \ X:ftrunc
\ round towards 0	\ round towards 0
\ !! should be implemented properly	fdup fabs floor fswap fcopysign ;
F>D D>F ;

: FMOD ( r1 r2 -- r )	: FMOD ( r1 r2 -- r )
\ remainder of r1/r2	\ remainder of r1/r2
FOVER FOVER F/ ftrunc F* F- ;	FOVER FOVER F/ ftrunc F* F- ;