<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Apr 23, 2014 at 9:43 AM, Nathaniel Smith <span dir="ltr"><<a href="mailto:njs@pobox.com" target="_blank">njs@pobox.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="">On Wed, Apr 23, 2014 at 6:22 AM, Matthew Brett <<a href="mailto:matthew.brett@gmail.com">matthew.brett@gmail.com</a>> wrote:<br>

> Hi,<br>
><br>
> I'm exploring Mingw-w64 for numpy building, and I've found it gives a<br>
> slightly different answer for 'exp' than - say - gcc on OSX.<br>
><br>
> The difference is of the order of the eps value for the output number<br>
> (2 * eps for a result of ~2.0).<br>
><br>
> Is accuracy somewhere specified for C functions like exp? ¬†Or is<br>
> accuracy left as an implementation detail for the C library author?<br>
<br>
</div>C99 says (sec 5.2.4.2.2) that "The accuracy of the floating point<br>
operations ... and of the library functions in <math.h> and<br>
<complex.h> that return floating point results is implemenetation<br>
defined. The implementation may state that the accuracy is unknown."<br>
(This last sentence is basically saying that with regard to some<br>
higher up clauses that required all conforming implementations to<br>
document this stuff, saying "eh, who knows" counts as documenting it.<br>
Hooray for standards!)<br>
<br>
Presumably the accuracy in this case is a function of the C library<br>
anyway, not the compiler? Numpy has its own implementations for a<br>
bunch of the math functions, and it's been unclear in the past whether<br>
numpy or the libc implementations were better in any particular case.<br></blockquote><div><br></div><div>In the case of MS runtime, at least 9 (as shipped in VS 2008), our implementation is likely to be better (most of the code was taken from the sun math library when the license allowed it).<br>
<br>David <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span class="HOEnZb"><font color="#888888"><br>
-n<br>
<br>
--<br>
Nathaniel J. Smith<br>
Postdoctoral researcher - Informatics - University of Edinburgh<br>
<a href="http://vorpus.org" target="_blank">http://vorpus.org</a><br>
</font></span><div class="HOEnZb"><div class="h5">_______________________________________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@scipy.org">NumPy-Discussion@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://mail.scipy.org/mailman/listinfo/numpy-discussion</a><br>
</div></div></blockquote></div><br></div></div>