Hi!

Does anyone have some high level source for FSIN and FCOS words? Of
course, FP words allowed ;-) BTW, FLN, F** or others are appreciated
as well.

Hans Bezemer

"The Beez'" bigfoot.com> writes:

: fsinn ( x -- {sin x}/x )
fsqr 1 >f fswap 1 >f 0 >r begin ( a x^2 t ; k )
1 r+! fover f* r 2* >f f/ r 2* 1+ >f f/ fnegate ( a x^2 t' ; k' )
frot fover fover f+ ( x^2 t' a a' ; k' )
ftuck f= >r f-rot r> until ( a' x^2 t' ; k' )
fdrop fdrop r> drop
;
: fsin fdup fsinn f* ;

: r+! rp@ +! ;
: >f s>d d>f ;
: f-rot frot frot ;

Obviously, it won't work for large arguments, you need another
way to calculate Taylor serie, but you want to take argument modulo pi
anyway.

cos(x) = sin (pi - x)

In article <4c1d96fb-b352-4cb1-a731-13f27c6926f0@d45g2000hsc.googlegroups.com>,
The Beez' bigfoot.com> wrote:

This is excerpted of the transputer sources. It will need some adaptation,
like removing the T in front of some words.
This is not rocket science: a standard Chebychev approximation is used
to arrive at most 2 bits from the exact value. (More on those on my
site below.) The binary representation of floats is IEEE.

The series have 8 terms. This is suboptimal for single precision.
_________________________________
TCREATE- SIN
#8 T,
$4A TC, $BB TC, $52 TC, $3C TC,
$2C TC, $26 TC, $06 TC,...

The Beez' wrote:

Consider "Computer Approximations" by John F Hart 1968 in the SIAM series.

"The Beez'" bigfoot.com> wrote in message
news:4c1d96fb-b352-4cb1-a731-13f27c6926f0@d45g2000hsc.googlegroups.com...

No.

_FIRST_ link Yahoo pulled up has source for them (Quartus):
http://www.quartus.net/files/PalmOS/Forth125/Snippets125/

You could try generating tables from algorithms here:
http://www.df.lth.se/~john_e/gems/gem0024.html
http://www.musicdsp.org/showone.php?id=9

Rod Pemberton

"Bill" wrote in message news:EPOdnab1x7Qqw47VnZ2dnUVZ_qbinZ2d@comcast.com...

This text is often quoted but unfortunately long out of print. Public
libraries and second-hand copies may be the only alternative.

It's ironic that algorithms and data useful to the public is ultimately
subject to considerations such as whether a publisher feels it
commercially viable to reprint.

On Apr 27, 10:52 am, "Ed" invalid.com> wrote:

In article news-01.bur.connect.com.au>,
Ed invalid.com> wrote: