I'll do it in a few days if noone else does, but right now I'm out of town for Thanksgiving and don't have access to my AMD computer.<br><br><div class="gmail_quote">On Thu, Nov 25, 2010 at 10:06 AM, Don Clugston <span dir="ltr"><<a href="mailto:dclugston@googlemail.com">dclugston@googlemail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">With regard to the AMD vs Intel bug:<br>
<br>
The previous results narrowed it down significantly, but I'm going to<br>
need more info to fix this.<br>
Please, somebody with AMD add this to line to<br>
std.math.internal.gammafunction line  624, add an import<br>
core.stdc.stdio; in the file,<br>
<br>
            if( x == 0.0L ) {<br>
                di = 0.5;<br>
<div class="im">                x = x0  +  di * (x1 - x0);<br>
</div>                if( x == 0.0 )<br>
                    goto under;<br>
<div class="im">            }<br>
            y = betaIncomplete( a, b, x );<br>
</div>+            printf("%d %La %La %La %La %La %La %La\n", i, x, y, di,<br>
yl, yh, a, b);<br>
            yp = (x1 - x0)/(x1 + x0);<br>
            if( fabs(yp) < dithresh )<br>
                goto newt;<br>
            yp = (y-y0)/y0;<br>
            if( fabs(yp) < dithresh )<br>
                goto newt;<br>
<br>
compile Phobos, and then run this program:<br>
<div class="im"><br>
import std.mathspecial;<br>
import std.stdio;<br>
<br>
void main()<br>
{<br>
</div><div class="im">writefln("%a", betaIncompleteInverse(0x1.ff1275ae5b939bcap-41, 4.6713e18,<br>
 0.0813601));<br>
}<br>
</div><div><div></div><div class="h5">_______________________________________________<br>
phobos mailing list<br>
<a href="mailto:phobos@puremagic.com">phobos@puremagic.com</a><br>
<a href="http://lists.puremagic.com/mailman/listinfo/phobos" target="_blank">http://lists.puremagic.com/mailman/listinfo/phobos</a><br>
</div></div></blockquote></div><br>