<br><br><div class="gmail_quote">On Tue, Jul 24, 2012 at 6:07 PM, Rüdiger Kessel <span dir="ltr">&lt;<a href="mailto:ruediger.kessel@gmail.com" target="_blank">ruediger.kessel@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Rüdiger Kessel &lt;ruediger.kessel &lt;at&gt; <a href="http://gmail.com" target="_blank">gmail.com</a>&gt; writes:<br>
<br>
PS: With test_stable() one can tell implementations apart. If snrm2() is<br>
implemented completely in single precision then the result of norm(a) will be<br>
10000.0. If internally a double precision implementation is used but the result<br>
is rounded to single precision then the result of norm(a) is 10000.5. Both<br>
implementations are possible and probably present in different libs. Therefore<br>
the test should allow for both results. It might produce some printed comment<br>
about the implementation.<br></blockquote><div><br>Thanks Rudiger, that explanation was helpful. <a href="https://github.com/scipy/scipy/pull/286">https://github.com/scipy/scipy/pull/286</a> allows for both results in the test.<br>
<br>Ralf<br> <br></div></div>