[phobos] phobos commit, revision 2186

[David Simcha]

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.

On Thu, Nov 25, 2010 at 10:06 AM, Don Clugston <dclugston at googlemail.com>wrote:

> With regard to the AMD vs Intel bug:
>
> The previous results narrowed it down significantly, but I'm going to
> need more info to fix this.
> Please, somebody with AMD add this to line to
> std.math.internal.gammafunction line  624, add an import
> core.stdc.stdio; in the file,
>
>            if( x == 0.0L ) {
>                di = 0.5;
>                 x = x0  +  di * (x1 - x0);
>                 if( x == 0.0 )
>                    goto under;
>             }
>            y = betaIncomplete( a, b, x );
> +            printf("%d %La %La %La %La %La %La %La\n", i, x, y, di,
> yl, yh, a, b);
>            yp = (x1 - x0)/(x1 + x0);
>            if( fabs(yp) < dithresh )
>                goto newt;
>            yp = (y-y0)/y0;
>            if( fabs(yp) < dithresh )
>                goto newt;
>
> compile Phobos, and then run this program:
>
> import std.mathspecial;
> import std.stdio;
>
> void main()
> {
> writefln("%a", betaIncompleteInverse(0x1.ff1275ae5b939bcap-41, 4.6713e18,
>  0.0813601));
> }
> _______________________________________________
> phobos mailing list
> phobos at puremagic.com
> http://lists.puremagic.com/mailman/listinfo/phobos
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.puremagic.com/pipermail/phobos/attachments/20101125/7d1c02bf/attachment-0001.html>