Path: utzoo!utgpu!news-server.csri.toronto.edu!cs.utexas.edu!sdd.hp.com!decwrl!bacchus.pa.dec.com!granite.pa.dec.com!mwm From: mwm@raven.pa.dec.com (Mike (Real Amigas have keyboard garages) Meyer) Newsgroups: comp.sys.amiga.tech Subject: Re: math in Lattice Message-ID: Date: 22 Aug 90 16:38:16 GMT References: <90233.133801KELLYDK@QUCDN.BITNET> <11700@wehi.dn.mu.oz> Sender: news@wrl.dec.com (News) Organization: Missionaria Phonibalonica Lines: 19 In-Reply-To: BAXTER_A@wehi.dn.mu.oz's message of 22 Aug 90 14:19:11 GMT In article <11700@wehi.dn.mu.oz> BAXTER_A@wehi.dn.mu.oz writes: Many of the examples from Numerical recipies don't work either. The explanation I got for it is that a lot of numerical work in C makes undue assumptions about the treatment of floats and doubles in functions. You could try #define double float to get a rough idea of what's going on. I gave up and hand coded. You can thank ANSI for this. K&R specified that all floating point arithmetic was done using doubles. ANSI changed that to allow the compiler to use the longest precision in the expression. Most compilers have a switch to force the old behavior, Lattice being among them - for some of it's math modes.