<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On Dec 29, 2012, at 10:04 PM, David Rush wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr">On 28 December 2012 16:53, Bradley Lucier <span dir="ltr"><<a href="mailto:lucier@math.purdue.edu" target="_blank">lucier@math.purdue.edu</a>></span> wrote:<br><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Inexact numbers are the same as C doubles (64-bit).  Extended precision floating-point numbers (80-bit numbers) are used in extremis on machines that do not support 64-bit floating-point numbers natively.  (Don't ask.)<br>
</div></blockquote><div><br></div><div style="">Too late. Besides i have a penchant for quirky/bizarre architectures. Are we talking about x86 here?</div></div></div></div></div></blockquote><br></div><div>x86 without the SSE extension.  In which case the precision of the mantissa of operations is set to 53 bits, but the exponent range of operations is still the exponent range of extended reals.  But inexact numbers are still stored as doubles.</div><div><br></div><div>I don't remember all the details; most of the time, results are as expected, but not always---for subnormal numbers double rounding can sometimes occur, for example.</div><div><br></div><div>Brad</div></body></html>